Statistical Analysis of the Spatial Distribution of Multi-elements in an Island Arc Region: Complicating Factors and Transfer by Water Currents

The compositions and transfer processes affecting coastal sea sediments from the Seto Inland Sea and the Pacific Ocean are examined through the construction of comprehensive terrestrial and marine geochemical maps for western Japan. Two-way analysis of variance (ANOVA) suggests that the elemental concentrations of marine sediments vary with particle size, and that this has a greater effect than the regional provenance of the terrestrial material. Cluster analysis is employed to reveal similarities and differences in the geochemistry of coastal sea and stream sediments. This analysis suggests that the geochemical features of fine sands and silts in the marine environment reflect those of stream sediments in the adjacent terrestrial areas. However, gravels and coarse sands do not show this direct relationship, which is likely a result of mineral segregation by strong tidal currents and the denudation of old basement rocks. Finally, the transport processes for the fine-grained sediments are discussed, using the spatial distribution patterns of outliers for those elements enriched in silt and clay. Silty and clayey sediments are found to be transported and dispersed widely by a periodic current in the inner sea, and are selectively deposited at the boundary of different water masses in the outer sea.


Introduction
The Geological Survey of Japan, part of the National Institute of Advanced Industrial Science and Technology (AIST), provides nationwide geochemical maps of elements within stream and marine sediments [1].These geochemical maps have been utilized to explore mineral occurrences and determine the natural abundance of elements, and nationwide and cross-boundary geochemical maps have been developed for such purposes in many countries (e.g., [2][3][4][5][6]).In addition, Japanese geochemical maps have also been created specifically for the purposes of environmental assessment.Ohta and Imai [7] demonstrate an additional use of land and marine geochemical maps, examining particle transfer processes from the land to the sea, or within the marine environment.However, determining provenances with geochemical maps is typically challenging, due to the mixing and homogenization of marine sediments, especially fine sediments such as silt, during transport (e.g., [8]).
In this study, we sought to use geochemical maps of the provenance and transfer analyses of coastal sea sediments in the Chugoku and Shikoku regions, including the Seto Inland Sea and nearby Pacific Ocean as inner and outer seas, respectively (Figure 1a).This region is appropriate for clarifying the influence of terrestrial source materials on the adjacent marine environment, because most rock types found in Japan are distributed in the study area.The Seto Inland Sea is subject to significant tidal variation, while the Pacific Ocean off the Shikoku region is influenced by the Kuroshio Current and the Kuroshio Counter-Current.A number of previous studies of the Seto Inland Sea have analyzed the seafloor topography [9], periodic currents [10], marine organization [11], surface sediments [12], particle transport by marine currents [13,14], and contamination processes [15].Therefore, this region is well-understood and is considered suitable for investigations into water-current transport processes, especially those between an inner and outer sea.The present study is intended to objectively reveal the complex factors affecting the spatial distributions of elements in coastal sea sediments using various statistical analyses.
Water 2017, 9, 37 2 of 26 significant tidal variation, while the Pacific Ocean off the Shikoku region is influenced by the Kuroshio Current and the Kuroshio Counter-Current.A number of previous studies of the Seto Inland Sea have analyzed the seafloor topography [9], periodic currents [10], marine organization [11], surface sediments [12], particle transport by marine currents [13,14], and contamination processes [15].Therefore, this region is well-understood and is considered suitable for investigations into water-current transport processes, especially those between an inner and outer sea.The present study is intended to objectively reveal the complex factors affecting the spatial distributions of elements in coastal sea sediments using various statistical analyses.A1.The blue colored area is the Seto Inland Sea.1: Kuroshio Current; 1': Kuroshio Counter-Current; 2: Tsushima-Current; 3: Oyashio-Current; 4: Liman-Current.Bathymetric depth contours are delineated using a dataset provided by the Japan Oceanographic Data Center.

Riverine System
Figure 1b presents a generalized location map of the study area showing the 24 major rivers.The region is mainly mountainous with small tracts of flat land.The bed slope of the streams on Shikoku Island and the Kii Peninsula is very steep, so rainwater is immediately discharged through the rivers to the sea.Table 1 presents the potential sediment yield from each terrestrial region, which is calculated using the data of Akimoto et al. [16] (Appendix A).Several rivers in the Chugoku and Kyushu regions that flow directly into the Sea of Japan, and from which marine sediments were collected, are not discussed in this study.In this regard, the sediment yield data of these rivers are excluded from Table 1.The Yodo, Yoshino, and Shimanto Rivers are the three largest rivers, each with a high sediment yield (Yd, Ys, and Sm, respectively, in Figure 1b), and the Chugoku, Kinki, Shikoku, and Kyushu regions supply 15%, 5%, 5%, and 1% of the total sediment yield to the Seto Inland Sea, respectively.Therefore, it can be seen that most sediments in the Seto Inland Sea originate from the Chugoku region.The Yodo River flows through several large cities, including Kyoto and Osaka, and discharges into Osaka Bay.Together, the Yodo and Yamato Rivers supply 19% of the total sediment yield to Osaka Bay.Rivers on Shikoku Island are the primary source of sediments to the Pacific Ocean side with 22% of the total sediment yield supplied to Tosa Bay.The Kii Channel also receives a high sediment discharge of 23% from Shikoku Island and the Kii Peninsula.In contrast, little sediment is discharged from rivers into the Bungo Channel.

Geology and Terrestrial Metalliferous Deposits
Figure 2 presents a geological map of the study area, simplified from the Geological Map of Japan 1:1,000,000 [17].The geology varies considerably between the Chugoku, Kinki, Shikoku, and Kyusyu regions.Rhyolitic-dacitic volcanic rocks and granitic rocks of Cretaceous and Paleogene age are widely distributed in the Chugoku region.In addition, Permian accretionary complexes associated with large limestone blocks, and Triassic high-pressure Sangun metamorphic rocks outcrop in western and central Chugoku.On Shikoku Island, the zonal arrangement of rock units from north to south is as follows: (1) Cretaceous granite-granodiorite; (2) Cretaceous sedimentary rocks; (3) Jurassic-Cretaceous high-pressure Sambagawa metamorphic rocks, comprising quartz schist and greenschist, associated with Mikabu greenstones consisting of basaltic, pyroclastic and ultramafic rocks; and (4) sedimentary rocks of accretionary complexes dated mainly to the Cretaceous-Paleogene.These rock types also outcrop on southern Kyushu Island and on the Kii Peninsula, while andesitic volcanic rocks and debris-pyroclastic rocks of Neogene-Quaternary age extensively outcrop in the northeast of Kyushu Island.Sedimentary rocks of accretionary complexes dating mostly to the Jurassic-Cretaceous outcrop in the northern Kinki region.Unconsolidated Quaternary sediments are restricted in distribution, occurring mainly in the Kinki region.The Osaka Plain is formed of these Quaternary sediments and is the widest plain in the entire study area, with a population of approximately 10 million people.
Figure 2 also shows the locations of some major economic metalliferous deposits.The Besshi mine is the largest Copper (Cu) mine in Japan, and the Kaneuchi and Ohtani mines are large-scale Tungsten-Tin (W-Sn) mines.The Ikuno and Akenobe mines are the largest polymetallic mineralization mines, and the Ichinokawa mine has the highest levels of Antimony (Sb) in Japan.

Figure 2.
Geological map of the study area simplified from the Geological Map of Japan 1:1,000,000 [17].

Marine Topography, Hydrographic Condition, and Geology
The Seto Inland Sea is mostly less than 60 m deep, and passes through to the Pacific Ocean through the Bungo and Kii Channels, and to the Sea of Japan through the Kan-mon Strait.The Seto Inland Sea is divided into several regions, termed "Nada", for example, Suo-Nada, Iyo-Nada, Aki-Nada, Bingo-Nada, Hiuchi-Nada, Bisan-Seto-Nada, and Harima-Nada (Figure 1b).Periodic currents flow from the Pacific Ocean through the Bungo and Kii Channels, meeting at Hiuchi-Nada and concentrating fine sediments in the water mass there [10,18].The marine geology of the Seto Inland Sea comprises mainly glacial to Holocene sediments [12].
The Pacific side of Shikoku Island is characterized by a narrow continental shelf.Topographic highs include Ashizuri Spar, Muroto Spar, Tosa-Bae, Ashizuri Sea Knoll, and Muroto Sea Knoll, and consist of Miocene-Pliocene siltstones, which include benthic foraminifer fossils that are distributed in the deep sub-bottom of the basin [19] In contrast, Tosa Basin, Hyuga Basin, and Muroto Trough are deep-sea basins in which the water depth exceeds 1000 m.Tosa Basin is covered thickly by Quaternary turbidite deposits.The continental shelf and slope of Tosa Bay, the Kii Channel, and the Bungo Channel comprise delta deposits formed during Quaternary regression and transgression cycles [19].Within the water mass, the Kuroshio Current flows off the coast of Shikoku Island from the southwest to the northeast, while the Kuroshio Counter-Current flows anticlockwise in Tosa Bay and Hyuga Basin [20,21] (Figure 1b).

Samples, Sampling Methods, and Processing
Terrestrial and marine sampling locations are presented in Figure 3.A total of 554 stream samples were collected from the study area for a regional geochemical mapping project during 1999-2004 [22].In addition, for the current study, 22 stream sediments were collected in 2008 from small islands, Geological map of the study area simplified from the Geological Map of Japan 1:1,000,000 [17].

Marine Topography, Hydrographic Condition, and Geology
The Seto Inland Sea is mostly less than 60 m deep, and passes through to the Pacific Ocean through the Bungo and Kii Channels, and to the Sea of Japan through the Kan-mon Strait.The Seto Inland Sea is divided into several regions, termed "Nada", for example, Suo-Nada, Iyo-Nada, Aki-Nada, Bingo-Nada, Hiuchi-Nada, Bisan-Seto-Nada, and Harima-Nada (Figure 1b).Periodic currents flow from the Pacific Ocean through the Bungo and Kii Channels, meeting at Hiuchi-Nada and concentrating fine sediments in the water mass there [10,18].The marine geology of the Seto Inland Sea comprises mainly glacial to Holocene sediments [12].
The Pacific side of Shikoku Island is characterized by a narrow continental shelf.Topographic highs include Ashizuri Spar, Muroto Spar, Tosa-Bae, Ashizuri Sea Knoll, and Muroto Sea Knoll, and consist of Miocene-Pliocene siltstones, which include benthic foraminifer fossils that are distributed in the deep sub-bottom of the basin [19] In contrast, Tosa Basin, Hyuga Basin, and Muroto Trough are deep-sea basins in which the water depth exceeds 1000 m.Tosa Basin is covered thickly by Quaternary turbidite deposits.The continental shelf and slope of Tosa Bay, the Kii Channel, and the Bungo Channel comprise delta deposits formed during Quaternary regression and transgression cycles [19].Within the water mass, the Kuroshio Current flows off the coast of Shikoku Island from the southwest to the northeast, while the Kuroshio Counter-Current flows anticlockwise in Tosa Bay and Hyuga Basin [20,21] (Figure 1b).

Samples, Sampling Methods, and Processing
Terrestrial and marine sampling locations are presented in Figure 3.A total of 554 stream samples were collected from the study area for a regional geochemical mapping project during 1999-2004 [22].In addition, for the current study, 22 stream sediments were collected in 2008 from small islands, including Awaji Island in the Seto Inland Sea.Sediment sample were air-dried and sieved using a 180 µm (83-mesh) screen.Magnetic minerals were removed from the sieved samples using a magnet, to minimize the effect of magnetic mineral accumulation [22,23].
Water 2017, 9, 37 5 of 26 including Awaji Island in the Seto Inland Sea.Sediment sample were air-dried and sieved using a 180 μm (83-mesh) screen.Magnetic minerals were removed from the sieved samples using a magnet, to minimize the effect of magnetic mineral accumulation [22,23].A total of 366 marine samples were collected from the Pacific Ocean during cruises GH82-1 and GH83-2, in 1982 and 1983, respectively [24][25][26].In addition, 97 samples were collected from the Seto Inland Sea and around the Kii Peninsula in 2005.The total 463 samples were collected using a K-grab sampler, and the uppermost 3 cm of each sediment sample was separated, air-dried, ground with an agate mortar and pestle, and retained for chemical analysis.Some samples were composed mainly of rock fragments and gravels, in which case, the sandy infilling materials were collected.The particle sizes of 223 marine sediments were determined, based on the median particle diameter of the surface sediments, and were classified as: coarse sediments, comprising lithic fragments, gravels, and coarse sand; medium sand; fine sand; silt; and clay (Figure 3).The median particle diameter was not measured for the remaining 240 samples, so their classification is based solely on a visual inspection of texture.
Figure 3 shows that silty sediments occur in Suo-Nada, Hiuchi-Nada, Bingo-Nada, and Harima-Nada, as well as basins affected by small periodic current flows or eddies [10].These localities are also associated with high turbidity [12,27].Coarse sands were collected from Aki-Nada and Bisan-Seto-Nada, and in northeastern Iyo-Nada.Rocks, gravels, and coarse sands occur in the Bungo Channel.In contrast, the Kii Channel is dominated by silty sediments, despite the presence of high velocity periodic currents similar to those in the Bungo Channel.This is likely due to the greater terrestrial sediment supply to the Kii Channel (18%) compared with the Bungo Channel (1%) (Table 1).On the Pacific side, silty sediments are widely distributed obliquely across the depth contour in the western part of Tosa Bay.Fine sands are found: (1) near to river mouths; (2) at water depths of 100-500 m on the shelf and slope off Tosa Bay; and (3) around Ashizuri Knoll and Muroto Knoll where the water depth is 500-1600 m.Silty and clayey sediments are found in Tosa Basin, Hyuga Basin, and A total of 366 marine samples were collected from the Pacific Ocean during cruises GH82-1 and GH83-2, in 1982 and 1983, respectively [24][25][26].In addition, 97 samples were collected from the Seto Inland Sea and around the Kii Peninsula in 2005.The total 463 samples were collected using a K-grab sampler, and the uppermost 3 cm of each sediment sample was separated, air-dried, ground with an agate mortar and pestle, and retained for chemical analysis.Some samples were composed mainly of rock fragments and gravels, in which case, the sandy infilling materials were collected.The particle sizes of 223 marine sediments were determined, based on the median particle diameter of the surface sediments, and were classified as: coarse sediments, comprising lithic fragments, gravels, and coarse sand; medium sand; fine sand; silt; and clay (Figure 3).The median particle diameter was not measured for the remaining 240 samples, so their classification is based solely on a visual inspection of texture.
Figure 3 shows that silty sediments occur in Suo-Nada, Hiuchi-Nada, Bingo-Nada, and Harima-Nada, as well as basins affected by small periodic current flows or eddies [10].These localities are also associated with high turbidity [12,27].Coarse sands were collected from Aki-Nada and Bisan-Seto-Nada, and in northeastern Iyo-Nada.Rocks, gravels, and coarse sands occur in the Bungo Channel.In contrast, the Kii Channel is dominated by silty sediments, despite the presence of high velocity periodic currents similar to those in the Bungo Channel.This is likely due to the greater terrestrial sediment supply to the Kii Channel (18%) compared with the Bungo Channel (1%) (Table 1).On the Pacific side, silty sediments are widely distributed obliquely across the depth contour in the western part of Tosa Bay.Fine sands are found: (1) near to river mouths; (2) at water depths of 100-500 m on the shelf and slope off Tosa Bay; and (3) around Ashizuri Knoll and Muroto Knoll where the water depth is 500-1600 m.Silty and clayey sediments are found in Tosa Basin, Hyuga Basin, and Muroto Trough, where the water depth exceeds 1000 m, and rock, gravels, and coarse sands were sampled from the topographic highs of Ashizuri Spar, Muroto Spar, and Tosa-Bae.

Geochemical and Spatial Analyses
The geochemical analytical method is detailed by Imai [28].Each 0.2 g sample was digested using HF, HNO 3 , and HClO 4 solutions at 120 • C for 2 h.The degraded solution was evaporated to dryness at 200 • C, before the residue was dissolved in 100 mL of 0.35 M HNO 3 solution.Concentrations of Na 2 O, MgO, Al 2 O 3 , P 2 O 5 , K 2 O, CaO, MnO, total (T-) Fe 2 O 3 , V, Sr, and Ba were determined using ICP atomic emission spectrometry, while contents of Li, Be, Sc, Cr, Co, Ni, Cu, Zn, Ga, Rb, Nb, Y, Mo, Cd, Sn, Sb, Cs, lanthanides (Ln: La-Lu), Ta, Tl, Pb, Bi, Th, and U were measured using ICP mass spectrometry.Analyses of As in all samples and Hg in stream sediment samples were subcontracted to ALS Chemex in Vancouver, B.C. Levels of Hg in marine sediments were determined using an atomic absorption spectrometer that measured the quantity of Hg vapor generated by direct thermal decomposition.Table 2 presents a summary of the analytical results obtained from the marine and stream sediments.The Na 2 O component of marine sediments should be used only as a guide, since these marine sediments were not desalinated.The spatial distributions of elemental concentrations in both the terrestrial and marine environments were plotted using geographic information system software (ArcGIS 10.3; Environmental Systems Research Institute (ESRI) Japan Corporation, Tokyo, Japan) after Ohta et al. [29].Different classes of elemental concentrations are applied to the terrestrial and marine environments, because the chemical and mineralogical compositions of sediments differ according to particle size see the details in [7,8].This simply improves geovisualization of the geochemical maps.If the same classification is applied to the geochemical maps of both terrestrial and marine environments, the regional geochemical differences (as depicted by color variation) of most of the elements in the land or the sea will be obscured [7,8].For example, for CaO and Sr, which show higher concentrations in the marine sediments (Table 2), their geochemical differences in the terrestrial areas would be obscured in resulting maps using the same classification, as shown in the Legend.Percentile ranges are used for the selection of elemental concentration intervals in the geochemical maps: 0 ≤ x ≤ 5, 5 < x ≤ 10, 10 < x ≤ 25, 25 < x ≤ 50, 50 < x ≤ 75, 75 < x ≤ 90, 90 < x ≤ 95, and 95 < x ≤ 100%, where x represents the elemental concentration [30].This class selection is advantageous in that the same range of percentiles (e.g., 90%-95%) implies the same statistical weight, even at different numerical scales [7,8,30].Subsequent statistical analysis of geochemical data is performed using EXCEL TOUKEI 7.0 (ESUMI Co. Ltd., Tokyo, Japan).

Spatial Distribution of Terrestrial Elemental Concentrations
Figure 4 shows the terrestrial and marine geochemical maps for 12 elements.Mikoshiba et al. [23] and Ohta et al. [31,32] have examined the factors controlling spatial distribution patterns of terrestrial elemental concentrations.Be, Na 2 O, K 2 O, Rb, Y, Nb, Ba, Ln, Ta, Th, and U are abundant in the Chugoku region and in the northern part of Shikoku Island where predominantly granitic and felsic volcanic rocks outcrop.Stream sediments derived from granitic rocks are enriched in Li, Be, Na 2 O, K 2 O, Rb, Y, Nb, Ln, Ta, Th, and U. Felsic volcanic rocks also elevate the K 2 O, Rb, and Ba contents of stream sediments.Li and Cs are abundant in muddy sedimentary rocks and the mélange matrix of accretionary complexes in northern Kinki.In contrast, the sandstone-dominated sedimentary rocks of accretionary complexes in the southern part of Shikoku Island and the Kii Peninsula contain fewer elements.Mafic volcanic rocks associated with debris flows and pyroclastic rocks around Aso volcano, and metamorphic rocks, dominantly greenschist, on Shikoku Island and the Kii Peninsula, elevate MgO, Al 2 O 3 , P 2 O 5 , CaO, Sc, V, TiO 2 , MnO, T-Fe 2 O 3 , Cr, Co, Ni, and Cu concentrations in the stream sediments.Additionally, extreme enrichments of MgO, Cr, Co, and Ni are caused by the presence of mafic-ultramafic rocks associated with accretionary complexes.Finally, elevated concentrations of P 2 O 5 , Cu, Zn, Mo, Cd, Sn, Sb, Hg, Pb and Bi are found in Osaka Plain, which is polluted by anthropogenic activity [32].
Therefore, lithology is considered to be the main factor affecting the concentrations of elements in stream sediments.For reference, the median elemental concentrations of stream sediments according to their representative rock type, which outcrops over more than half of the drainage basin, are summarized in Appendix B. It should be noted that Cu, Zn, As, Mo, Cd, Sn, Sb, Hg, Pb, and Bi are extremely enriched near to known metalliferous deposits, however, the effects of these mineral deposits on the geochemical maps are restricted to limited areas.

Spatial Distribution of Elemental Concentrations in Coastal Sea Sediments
Within the Seto Inland Sea, high levels of K 2 O, Rb, Ba, and Tl are found in the coarse sands of Aki-Nada and Bisan-Seto-Nada, and in northeastern Iyo-Nada (Figure 4).Coarse sediments with rock fragments, gravels, and coarse sands in western Iyo-Nada and the Bungo Channel are enriched in MgO, P 2 O 5 , CaO, TiO 2 , MnO, Fe 2 O 3 , Sc, V, Co, Zn, and Sr.In contrast, the sandy sediments in eastern Iyo-Nada and the Bungo Channel are depleted in almost all of the measured elements, except CaO and Sr. Silt and clay in Hiuchi-Nada, Harima-Nada, and the Kii Channel are abundant in Li, Be, MgO, P 2 O 5 , 3d transition metals, Ga, Y, Nb, Sn, Cs, Ln, Ta, Hg, Pb, Bi, and Th.Suo-Nada is enriched in Li, Zn, Ga, Sn, Mo, Cd, Ta, Pb, Bi, and U. Concentrations of P 2 O 5 , Cr, Ni, Cu, Zn, Mo, Cd, Sn, Hg, Pb, and Bi are especially elevated in the silt of Osaka Bay.Finally, MnO and Tl are abundant in all sediments within the Seto Inland Sea, irrespective of particle size.
On the Pacific Ocean side, coarse sediments with rock fragments, gravels, and coarse sands distributed around the topographic highs of Ashizuri Spar, Muroto Spar, Tosa-Bae, and the Muroto Sea Knoll are abundant in MgO, CaO, Al 2 O 3 , P 2 O 5 , TiO 2 , MnO, T-Fe 2 O 3 , Sc, V, Co, As, Sr, Mo, Sb, and U. Sediments abundant in TiO 2 , Cr, Co, Ni, Nb, light REEs (lanthanides from La to Sm), and Ta, occur in Tosa Bay.The silty and clayey sediments of Tosa Basin, Hyuga Basin, and the Muroto trough are enriched in Li, Cu, Nb, Cd, Sb, Cs, Hg, and U. Fine sands collected from the Ashizuri and Muroto Sea Knolls are richer in P 2 O 5 , CaO, Sc, TiO 2 , V, MnO, Fe 2 O 3 , Co, Cu, As, Sr, Mo, Cd, Sb, Cs, heavy REEs (Y and lanthanides from Gd to Lu), and Bi than the fine sands distributed across the continental shelf.

Analysis of Variance (ANOVA) to Reveal the Effects of Particle Size and Regional Difference
Coastal sea sediments typically originate from their adjacent terrestrial materials.Therefore, the geochemical features of terrestrial lithology are fundamentally reflected in those of nearby marine sediments.This phenomenon is referred to herein as "regional difference".Additionally, elemental concentrations in marine sediments are known to vary with particle size, in a phenomenon referred to herein as the "particle size effect".This variation is caused simply by the dilution effects of quartz and biogenic calcareous materials, which are abundant in coarse grains but less abundant in fine grains.Thus, we apply a two-way analysis of variance (ANOVA) to determine the factor that most significantly controls the chemical compositions of marine sediments: regional differences or the particle size effect [8].The procedure is detailed by [8] (Appendix C).
The marine samples are grouped into those from the Seto Inland Sea region, including the Kii and Bungo Channels, and those from the Pacific Ocean region.They are further classified into: coarse sediments, including rock fragments, gravels, coarse sands, and medium sands; fine sands; and silt-clay.We assume that the coarse sediments, fine sands, and silt-clay are moved by water power over short, middle, and long distances, respectively.Median elemental concentrations of marine sediments are calculated for each region and particle size (Table 3).Table 3 also presents the median concentrations of stream sediments from two regions, for reference.In this regard, 96 stream sediment samples were collected from rivers that flow directly into the Sea of Japan, and so these samples are excluded from Table 3.  Table 4 presents the results of a two-way ANOVA: the variance ratios (F), probabilities (p), and effect size (η 2 ) due to regional difference (factor A), particle size (factor B), and the interaction effect (factor A × B).When the estimated probability is lower than 0.01, we conclude that the factor makes a significant difference to the chemical compositions, and in fact, each factor was statistically significant (p < 0.01) in most cases (Table 4).This is because a statistical test with a large amount of input data, such as n = 463 in this study, is highly sensitive to very small differences.Therefore, we calculate η 2 to form a plausible estimation of p value irrespective of sample number [33,34].The η 2 is an easy-to-understand metric that is calculated as the ratio of the sum of squares (SS) for each factor and the total SS.Magnitudes of η 2 < 0.01, 0.01 ≤ η 2 < 0.06, 0.06 ≤ η 2 < 0.14, and η 2 ≥ 0.14 can be considered to represent no effect, a small effect, an intermediate effect, and a large effect, respectively [33,34].In this study, we conclude that each factor with η 2 ≥ 0.06 has a significant effect on the elemental concentrations of the sediments.Furthermore, this indicator can be used to decide the most dominant factor, on the basis of the highest η 2 score [33].Notes: The factor A, B, and A × B indicate the regional difference effect, particle size effect and interaction effect, respectively.MF means major factor.The boldface in η 2 scores indicates that each factor with η 2 ≥ 0.06 has a significant effect on the elemental concentrations of the sediments.
Results of the ANOVA suggest that particle size is the dominant influencing factor for many elements in the sediments (Table 4).The regional difference effect is more significant than the particle size effect only for Al 2 O 3 , CaO, MnO, Sc, V, Mo, and Pb.This suggests that consideration of the particle size effect must be made in order to fully understand the influence of terrestrial materials on coastal sea sediments.However, it should be noted that the interaction effect is more influential than either of the two main factors alone for MgO, P 2 O 5 , Fe 2 O 3 , Sc, Co, Rb, Ba, and Eu, and is significant for 14 elements, including K 2 O and TiO 2 .This interaction effect refers to the effect of one factor on the other factor, such that the two main factors synergistically affect the data e.g., [35].To visualize the interaction effect, variations in the median values of Al 2 O 3 , MgO, K 2 O, Co, Cs, and Ba are presented in Figure 5.When the interaction effect is insignificant, the median values vary approximately in parallel between the Seto Inland Sea and the Pacific Ocean, irrespective of particle size (Al 2 O 3 and Cs in Figure 5).In the cases of large η 2 values for the interaction effect, such as for MgO and Co, the elemental concentrations in the Pacific Ocean sediments decrease with decreasing particle size, but increase in the sediments of Water 2017, 9, 37 13 of 26 the Seto Inland Sea with decreasing particle size.Variations in the K 2 O and Ba concentrations with particle size are exactly opposite to those of MgO and Co, and still indicate a significant contribution of the interaction effect.These results suggest that regional difference and particle size effects are not independent for many elements.Therefore, we need to evaluate the geochemical similarities and differences among groups subdivided by region and particle size.
for 14 elements, including K2O and TiO2.This interaction effect refers to the effect of one factor on the other factor, such that the two main factors synergistically affect the data e.g., [35].To visualize the interaction effect, variations in the median values of Al2O3, MgO, K2O, Co, Cs, and Ba are presented in Figure 5.When the interaction effect is insignificant, the median values vary approximately in parallel between the Seto Inland Sea and the Pacific Ocean, irrespective of particle size (Al2O3 and Cs in Figure 5).In the cases of large η 2 values for the interaction effect, such as for MgO and Co, the elemental concentrations in the Pacific Ocean sediments decrease with decreasing particle size, but increase in the sediments of the Seto Inland Sea with decreasing particle size.Variations in the K2O and Ba concentrations with particle size are exactly opposite to those of MgO and Co, and still indicate a significant contribution of the interaction effect.These results suggest that regional difference and particle size effects are not independent for many elements.Therefore, we need to evaluate the geochemical similarities and differences among groups subdivided by region and particle size.A, B, and C represent coarse sediment, fine sand, and silty-clayey sediments, respectively.η 2 (A), η 2 (B), and η 2 (A × B) are the size effects of regional differences, particle sizes, and the interaction effect, respectively, obtained from the two-way ANOVA (Table 4).

Comparison of Elemental Abundances in Coastal Sea and Stream Sediments between the Seto Inland Sea and the Pacific Ocean
The chemical and mineralogical compositions of marine sediments differ from those of stream sediments [7,8].Furthermore, the range in concentration completely differs for each element: major and minor elements in wt % oxide amounts, and trace elements in mg/kg or μg/kg amounts.Therefore, to eliminate the effect of scale and units, the data have to be normalized to elucidate differences amongst the elements or the samples.Thus, in this study, enrichment factors (EF) of the marine and stream sediments compared with the upper continental crust (UCC) [36] are employed to allow direct comparison (Figure 6).The EF is calculated as follows:  η 2 (A), η 2 (B), and η 2 (A × B) are the size effects of regional differences, particle sizes, and the interaction effect, respectively, obtained from the two-way ANOVA (Table 4).

Comparison of Elemental Abundances in Coastal Sea and Stream Sediments between the Seto Inland Sea and the Pacific Ocean
The chemical and mineralogical compositions of marine sediments differ from those of stream sediments [7,8].Furthermore, the range in concentration completely differs for each element: major and minor elements in wt % oxide amounts, and trace elements in mg/kg or µg/kg amounts.Therefore, to eliminate the effect of scale and units, the data have to be normalized to elucidate differences amongst the elements or the samples.Thus, in this study, enrichment factors (EF) of the marine and stream sediments compared with the upper continental crust (UCC) [36] are employed to allow direct comparison (Figure 6).The EF is calculated as follows: where [C] and [Al 2 O 3 ] are the concentrations of a given element and Al 2 O 3 , respectively.Al 2 O 3 is a major element group, and a dominant constituent of many minerals such as plagioclase.Therefore by using it to normalize the EFs, we can effectively remove the effects of quartz, calcareous materials, and organic matter [8].The systematically high EFs of Li, As, Sb, Hg, and Bi, and low EFs of Nb, Sn, and Ta for all samples, as well as CaO and Sr in stream sediments, are fundamental geochemical signatures of upper crust materials in an island arc setting [37] (Figure 6a,b).The high abundance ratios of CaO and Sr in marine sediments indicate the presence of biogenic calcareous materials, such as shell fragments and foraminifera tests.In the Seto Inland Sea, the UCC normalized patterns of fine sands are similar to those of stream sediments (Figure 6a).Abundance ratios of most elements in marine sediments decrease with increasing grain size.CaO, K 2 O, Rb, Sr, and Ba are especially dominant in coarse sediments, whereas silt and clay have higher abundance ratios for Li, Cr, Co, Ni, Zn, Mo, Cd, Cs, Hg, Pb, and Bi.In the Pacific Ocean, the UCC normalized patterns of fine sands and silt-clay samples are most similar to those of stream sediments (Figure 6b).The variations in EF of most elements with particle size of the marine sediments in the Pacific Ocean are likely smaller than those for samples in the Seto Inland Sea, except for MnO, As, Ce, Eu, and Tl.Coarse sediments have relatively high abundance ratios for MnO, Fe 2 O 3 , Co, and As, and low abundance ratios for Li, Cu, Nb, Sn, Cs, light REEs, Ta, Hg, Tl, Th, and U.The silty and clayey sediments are particularly enriched in Cr, Ni, Cu, Cd, Sb, Hg, and Bi.
major element group, and a dominant constituent of many minerals such as plagioclase.Therefore by using it to normalize the EFs, we can effectively remove the effects of quartz, calcareous materials, and organic matter [8].The systematically high EFs of Li, As, Sb, Hg, and Bi, and low EFs of Nb, Sn, and Ta for all samples, as well as CaO and Sr in stream sediments, are fundamental geochemical signatures of upper crust materials in an island arc setting [37] (Figure 6a,b).The high abundance ratios of CaO and Sr in marine sediments indicate the presence of biogenic calcareous materials, such as shell fragments and foraminifera tests.In the Seto Inland Sea, the UCC normalized patterns of fine sands are similar to those of stream sediments (Figure 6a).Abundance ratios of most elements in marine sediments decrease with increasing grain size.CaO, K2O, Rb, Sr, and Ba are especially dominant in coarse sediments, whereas silt and clay have higher abundance ratios for Li, Cr, Co, Ni, Zn, Mo, Cd, Cs, Hg, Pb, and Bi.In the Pacific Ocean, the UCC normalized patterns of fine sands and silt-clay samples are most similar to those of stream sediments (Figure 6b).The variations in EF of most elements with particle size of the marine sediments in the Pacific Ocean are likely smaller than those for samples in the Seto Inland Sea, except for MnO, As, Ce, Eu, and Tl.Coarse sediments have relatively high abundance ratios for MnO, Fe2O3, Co, and As, and low abundance ratios for Li, Cu, Nb, Sn, Cs, light REEs, Ta, Hg, Tl, Th, and U.The silty and clayey sediments are particularly enriched in Cr, Ni, Cu, Cd, Sb, Hg, and Bi.

Discrimination of Coastal Sea Sediments Using Cluster Analysis
Coarse sediments have UCC normalized patterns that notably differ from those of terrestrial materials (Figure 6a

Discrimination of Coastal Sea Sediments Using Cluster Analysis
Coarse sediments have UCC normalized patterns that notably differ from those of terrestrial materials (Figure 6a,b), which implies that they have different origins to the modern terrestrial materials or are supplied by a particular lithology.Thus, we refine the identification of probable source materials for marine sediments using cluster analysis.The classification of marine sediment samples is as presented in Table 3. Stream sediments grouped according to the seven rock types summarized in Table B1 are used for comparison.The distance between data points was calculated as a Euclidian distance, requiring the geochemical data from marine and stream sediments to be standardized or transformed appropriately.This is because the compositional data constitute "closed" data that sum to a constant [38,39].To mitigate biased analysis arising from data closure, there is a requirement for data transformation to yield data with normal distributions suitable for parametric statistical analysis.Additive log-ratio, the centered log-ratio, and the isometric log-ratio transformation of compositional data have been proposed [38,39].We used a log transformation of the EF data to generate logarithmic EFs for each sample, which were used as the input data for cluster analysis.The logarithmic EF corresponds to the additive log-ratio transformation of the compositional data [38,39]: where alr(C) stands for additive log-ratio transformation (log {[C]/[Al 2 O 3 ]}).We have confirmed that the EF is effective in eliminating the effect of scale and units, and in removing the dilution effect, although Al 2 O 3 concentration is empirically used for denominator in Equation ( 1).The logarithmic EF also incorporates data normality and is thus amenable to multivariate parametric tests such as cluster analysis.
Figure 7 presents dendrograms that show the distances between two points.The distances between clusters were obtained by Ward's method [40].We calculated two kinds of dendrograms using different combinations of elements.Figure 7a uses a dataset that excludes sea salt (Na 2 O), biogenic carbonate materials (CaO and Sr), and heavy metals related to mining and anthropogenic activities (As, Mo, Cd, Sn, Sb, Hg, Pb, and Bi), which do not reflect the geochemistry of the parent lithology.Figure 7b uses the EFs of immobile elements (Sc, TiO 2 , Cr, Nb, Ta, Y, Ln, and Th), which are not strongly influenced by weathering processes.The obtained dendrograms are mutually similar (Figure 7a,b).Coarse sediments from the Pacific Ocean plot in the same group as stream sediments derived from mafic volcanic and metamorphic rocks in Figure 7a.However, Figure 7b shows that all marine sediments in the Pacific Ocean cluster together with coarse sediments from Seto Inland Sea and accretionary complexes.Fine sand and silt-clay in the Pacific Ocean and coarse sands in the Seto Inland Sea are always grouped together with stream sediments derived from sedimentary rocks of accretionary complexes.These results for fine sands and silt-clay are reasonable, because such sedimentary rocks outcrop extensively across the southern part of Shikoku Island (Figure 2), which supplies the highest sediment yield to the Pacific Ocean (Table 1).In contrast, the fine sands and silt-clay in the Seto Inland Sea plot with those stream sediments derived from granitic and felsic volcanic rocks.This result is also expected, due to the extensive outcropping of these rock types across the Chugoku region, which supplies the highest sediment yield to the Seto Inland Sea (Table 1).materials or are supplied by a particular lithology.Thus, we refine the identification of probable source materials for marine sediments using cluster analysis.The classification of marine sediment samples is as presented in Table 3. Stream sediments grouped according to the seven rock types summarized in Table B1 are used for comparison.The distance between data points was calculated as a Euclidian distance, requiring the geochemical data from marine and stream sediments to be standardized or transformed appropriately.This is because the compositional data constitute "closed" data that sum to a constant [38,39].To mitigate biased analysis arising from data closure, there is a requirement for data transformation to yield data with normal distributions suitable for parametric statistical analysis.Additive log-ratio, the centered log-ratio, and the isometric log-ratio transformation of compositional data have been proposed [38,39].We used a log transformation of the EF data to generate logarithmic EFs for each sample, which were used as the input data for cluster analysis.The logarithmic EF corresponds to the additive log-ratio transformation of the compositional data [38,39]: where alr(C) stands for additive log-ratio transformation (log {[C]/[Al2O3]}).We have confirmed that the EF is effective in eliminating the effect of scale and units, and in removing the dilution effect, although Al2O3 concentration is empirically used for denominator in Equation ( 1).The logarithmic EF also incorporates data normality and is thus amenable to multivariate parametric tests such as cluster analysis.
Figure 7 presents dendrograms that show the distances between two points.The distances between clusters were obtained by Ward's method [40].We calculated two kinds of dendrograms using different combinations of elements.Figure 7a uses a dataset that excludes sea salt (Na2O), biogenic carbonate materials (CaO and Sr), and heavy metals related to mining and anthropogenic activities (As, Mo, Cd, Sn, Sb, Hg, Pb, and Bi), which do not reflect the geochemistry of the parent lithology.Figure 7b uses the EFs of immobile elements (Sc, TiO2, Cr, Nb, Ta, Y, Ln, and Th), which are not strongly influenced by weathering processes.The obtained dendrograms are mutually similar (Figure 7a,b).Coarse sediments from the Pacific Ocean plot in the same group as stream sediments derived from mafic volcanic and metamorphic rocks in Figure 7a.However, Figure 7b shows that all marine sediments in the Pacific Ocean cluster together with coarse sediments from Seto Inland Sea and accretionary complexes.Fine sand and silt-clay in the Pacific Ocean and coarse sands in the Seto Inland Sea are always grouped together with stream sediments derived from sedimentary rocks of accretionary complexes.These results for fine sands and silt-clay are reasonable, because such sedimentary rocks outcrop extensively across the southern part of Shikoku Island (Figure 2), which supplies the highest sediment yield to the Pacific Ocean (Table 1).In contrast, the fine sands and siltclay in the Seto Inland Sea plot with those stream sediments derived from granitic and felsic volcanic rocks.This result is also expected, due to the extensive outcropping of these rock types across the Chugoku region, which supplies the highest sediment yield to the Seto Inland Sea (Table 1).Sed: sediments and sedimentary rocks; Acc: sedimentary rocks of accretionary complexes; Mv: mafic volcanic rocks; Fv: felsic volcanic rocks; Gr: granitic rocks; Meta: metamorphic rocks.Cluster analysis suggests that the coarse sediments in the Seto Inland Sea are related to stream sediments sourced from sedimentary rocks of accretionary complexes.However, granitic and felsic volcanic rocks are geographically more suitable sources for the coarse sediments, in analogy with the fine sands.To reveal their geochemical features in more detail, Figure 8a shows the abundance patterns of elements in coarse sediments collected from the narrow sea channels (Aki-Nada, Bisan-Seto-Nada and around Awaji Island), Iyo-Nada, and Bungo Channel, normalized to the stream sediments in the Seto Inland Sea region.

Fractionation of Mineralogical Compositions in Coarse Sediments by a Strong Tidal Current
Cluster analysis suggests that the coarse sediments in the Seto Inland Sea are related to stream sediments sourced from sedimentary rocks of accretionary complexes.However, granitic and felsic volcanic rocks are geographically more suitable sources for the coarse sediments, in analogy with the fine sands.To reveal their geochemical features in more detail, Figure 8a shows the abundance patterns of elements in coarse sediments collected from the narrow sea channels (Aki-Nada, Bisan-Seto-Nada and around Awaji Island), Iyo-Nada, and Bungo Channel, normalized to the stream sediments in the Seto Inland Sea region.Coarse sands in the narrow sea channels are particularly rich in K2O, CaO, Rb, Sr, Ba, and Tl, and have a distinctive positive Eu anomaly, but are extremely depleted in many elements, such as Sc and TiO2 (Figure 8a).These features are indicative of the selective accumulation of biogenic calcareous materials, quartz, and alkali-feldspars.Quartz and alkali-feldspar are abundant in granitic and felsic volcanic rocks, and are highly resistant to weathering processes, so are preserved as coarse particles e.g., [41].In contrast, the coarse sediments in Iyo-Nada and the Bungo Channel are relatively abundant in MgO, P2O5, CaO, Sc, TiO2, V, Fe2O3, Co, and Zn in addition to the elements mentioned above (Figure 8a).These additional elements are found abundantly in mafic minerals such as biotite,   Coarse sands in the narrow sea channels are particularly rich in K 2 O, CaO, Rb, Sr, Ba, and Tl, and have a distinctive positive Eu anomaly, but are extremely depleted in many elements, such as Sc and TiO 2 (Figure 8a).These features are indicative of the selective accumulation of biogenic calcareous materials, quartz, and alkali-feldspars.Quartz and alkali-feldspar are abundant in granitic and felsic volcanic rocks, and are highly resistant to weathering processes, so are preserved as coarse particles e.g., [41].In contrast, the coarse sediments in Iyo-Nada and the Bungo Channel are relatively abundant in MgO, P 2 O 5 , CaO, Sc, TiO 2 , V, Fe 2 O 3 , Co, and Zn in addition to the elements mentioned above (Figure 8a).These additional elements are found abundantly in mafic minerals such as biotite, hornblende, and pyroxene, as well as accessary minerals such as magnetite and ilmenite.We therefore assume that the fine particles of mafic and accessary minerals are swept away from the narrow sea channels by a strong tidal current, but this does not occur in Iyo-Nada and the Bungo Channel.In contrast, the coarse grains in all areas are rarely moved over any long distances by water power.

Denudation of the Basement Rocks of the Topographic Highs
Cluster analysis provides different results for coarse sediments from the Pacific Ocean (compare Figure 7a,b).This is because only three elements (Sc, TiO 2 , and Cr) abundant in sediments derived from mafic volcanic rocks and metamorphic rocks (dominantly greenschist) are used in the cluster analysis of the immobile element dataset (Figure 7b).Thus, cluster analysis using immobile elements may be poorly sensitive to such rock types.Consequently, we conclude that Figure 7a would provide a more plausible result: mafic volcanic and metamorphic rocks are the probable sources of coarse sediments in the Pacific Ocean.However, these rock types are located relatively distantly from the coarse sediment deposition sites (Figure 2).The high EFs of MnO, Fe 2 O 3 , Co, and As in coarse sediments from the Pacific Ocean (Figure 6b) can be explained by Fe hydroxide and Mn dioxide coatings on coarse particles and gravel (for example, T-Fe 2 O 3 and As; Figure 4).In this regard, relict quartz particles coated with Fe hydroxides are found in the Bungo Channel [26].However, this explanation does not fully account for the enrichment of MgO, P 2 O 5 , Sc, TiO 2 , and V.According to Okamura [19], the topographic highs in this region are formed from basement rocks of Miocene-Pliocene siltstone.These highs typically have a low sedimentation rate and are located in a region of erosion, due to the strong bottom flow associated with the Kuroshio Current [24].The spar and knoll sediments are therefore assumed to have formed by denudation of the basement rocks (siltstone), which have a different lithology to the southern part of Shikoku Island.Therefore, we assume that basement rocks of siltstone might be influenced by mafic volcanic activity at Miocene-Pliocene age. Figure 8b shows the abundance of elements in coarse or gravelly spar sediments and knoll sediments (mainly fine sands), normalized to the stream sediments in the Pacific Ocean region.Despite the fact that Ashizuri Spar is more than 100 km away from Muroto Spar and Tosa-Bae, the abundance patterns of the coarse spar sediments are very similar.Furthermore, the fine sands collected around Ashizuri Knoll and Muroto Knoll have similar abundance patterns to those of the coarse spar sediments.As such, the unexpected cluster analysis result for coarse sediments in this region can be explained by mineral segregation under a strong tidal current unique to the Seto Inland Sea, and by a denudation of old basement rocks (Miocene-Pliocene siltstone) in the Pacific Ocean.Silty and clayey sediments supplied from rivers take a long time to reach the sea floor.Therefore, their spatial distribution is a result of conveyance and dispersion by coastal sea currents.Ohta and Imai [7] and Ohta et al. [8] simply focused on the spatial distribution patterns of high Cr and Ni concentrations in silty sediments, which are indicative of transport of fine particles from land to sea or marine environments.Their spatial distribution can be easily discriminated because they are extremely enriched in silty marine sediments derived from ultramafic rocks in the adjacent terrestrial area.Ohta et al. [42] tried to visualize the particle transfer process using the spatial distribution patterns of the other elements.However, many elements are enriched in finer sediments due to the particle size effect, and so in order to highlight the particle transfer process, it is useful to detect outliers using the relationships between elemental concentrations and median diameter or mud content [42].
Water 2017, 9, 37 18 of 26 Unfortunately, in this study, data of median diameter or mud content are available for only half of the samples.ANOVA analyses suggest that Cs, Hg and Th concentrations in marine sediments are determined only by particle size and are not sensitive at all to either regional differences or interaction effects (η 2 < 0.01) (Table 4).The particle size effect of Cs has the largest η 2 score of the three elements, suggesting that the particle size effect has the maximum effect on Cs concentrations.Although concentrations of 13 elements, including Cr, Cu, and La, in marine sediments are also controlled dominantly by particle size effects, they are further influenced weakly by regional difference effects and interaction effects (0.01 ≤ η 2 < 0.06) (Table 4).In addition, Hg, Cu, and Cd concentrations in marine sediments might be influenced by metalliferous deposits and anthropogenic activity in the adjacent terrestrial area.For these reasons, we conveniently use Cs concentration as a proxy for the median diameter or mud content.
Figure 9 shows the relationships between Cr, Cu, Zn, Cd, Sb, and La concentrations and that of Cs.The concentrations of all six elements correlate positively and linearly with Cs, and it can be assumed that those samples plotting outside the linear trends are subject to influencing factors in addition to the particle size effect.Ohta et al. [42] detected outliers from a scatter diagram according to their best judgment.In this study, the outliers appearing in the scatter diagram are objectively obtained using Mahalanobis' generalized distance (D) [43].This method presupposes that the data follow a normal distribution similar to ANOVA.Thus, the data transformation is given in Table C1.The obtained D values gradually increase from the inner part of the data cluster to the outer part, giving an onion-like structure (Figure 9).D values of D ≤ 1.117, D ≤ 1.665, and D ≤ 2.146 indicate that 50%, 75%, and 90% of the total data points are included, respectively, and data points are defined as an outlier herein if D > 2.146, which includes the upper and lower 5% of multivariate outliers.The outliers were determined from all datasets, so silt and clay classified as outliers are expressed as cross symbols; while coarse sediments and fine sands classified as outliers are expressed as plus symbols in Figure 4.
Water 2017, 9, 37 18 of 26 suggesting that the particle size effect has the maximum effect on Cs concentrations.Although concentrations of 13 elements, including Cr, Cu, and La, in marine sediments are also controlled dominantly by particle size effects, they are further influenced weakly by regional difference effects and interaction effects (0.01 ≤ η 2 < 0.06) (Table 4).In addition, Hg, Cu, and Cd concentrations in marine sediments might be influenced by metalliferous deposits and anthropogenic activity in the adjacent terrestrial area.For these reasons, we conveniently use Cs concentration as a proxy for the median diameter or mud content.
Figure 9 shows the relationships between Cr, Cu, Zn, Cd, Sb, and La concentrations and that of Cs.The concentrations of all six elements correlate positively and linearly with Cs, and it can be assumed that those samples plotting outside the linear trends are subject to influencing factors in addition to the particle size effect.Ohta et al. [42] detected outliers from a scatter diagram according to their best judgment.In this study, the outliers appearing in the scatter diagram are objectively obtained using Mahalanobis' generalized distance (D) [43].This method presupposes that the data follow a normal distribution similar to ANOVA.Thus, the data transformation is given in Table C1.The obtained D values gradually increase from the inner part of the data cluster to the outer part, giving an onion-like structure (Figure 9).D values of D ≤ 1.117, D ≤ 1.665, and D ≤ 2.146 indicate that 50%, 75%, and 90% of the total data points are included, respectively, and data points are defined as an outlier herein if D > 2.146, which includes the upper and lower 5% of multivariate outliers.The outliers were determined from all datasets, so silt and clay classified as outliers are expressed as cross symbols; while coarse sediments and fine sands classified as outliers are expressed as plus symbols in Figure 4.

Transport of Elements in Silt Derived from the Parent Lithology
Figure 4 shows that the locations of the outliers calculated for Cr (cross symbols) correspond with a NNW-SSE distribution of silty-clayey sediments with high Cr concentrations in the Kii Channel.This distribution of Cr outliers indicates that silty particles that originated from ultramafic rocks are supplied through the Yoshino, Naka, and Kino Rivers to the Kii Channel, and are then widely dispersed by a periodic current, as Ohta and Imai [7] noted.Hiuchi-Nada also contains abundant silty sediments, however the distribution of outliers (corresponding to high concentration areas) for MgO, Cr, Co, and Ni is restricted to the vicinity of Shikoku Island (shown by Cr in Figure 4).This can be explained by the fact that the terrestrial sediment discharge yield to Hiuchi-Nada is minimal compared to the output of the Yoshino River, Naka River and Kino River systems (Table 1).
Y, Nb, Ln, Ta and Th are also enriched in the silty sediments of Hiuchi-Nada, especially in the eastern part (La in Figure 4).The differing distributions of La and Cr outliers suggest that Y, Nb, Ln, Ta and Th have a different origin from MgO, Cr, Co, and Ni.Minakawa et al. [44] reported that rare earth minerals, including allanite and monazite, occur within pegmatite veins of the Ryoke granitic rocks, which outcrop in northwest Shikoku Island, on the small islands in Aki-Nada and Bisan-Seto-Nada, and are partly exposed on the seafloor [12].Yanagi [10] suggested that fine particles transported by periodic currents accumulate in eastern Hiuchi-Nada where the amplitude of the tidal current is small and a counterclockwise circulation occurs.Consequently, we assume that silt enriched in rare earth minerals has been swept toward the eastern side of Hiuchi-Nada by periodic currents and accumulated there over a long period of time.
In Tosa Bay, tongue-shaped distribution patterns of high MgO, Cr, Co, Ni, Nb, Y, Ln, Ta, and Th concentrations and their outliers are apparent (shown by Cr and La in Figure 4).MgO, Cr, Co, and Ni originated from ultramafic rocks, while the other elements are derived from greenstones wedged between the Sambagawa metamorphic rocks and accretionary complexes e.g., [45].Ikehara [25] and Hoshino [46] reported that silty and clayey sediments are selectively deposited at the boundary between coastal waters with low salinity (due to River input) and the outer sea with higher salinity.This is supported by the fact that the water discharges from the Niyodo and Shimanto Rivers into Tosa Bay are among the largest in the region (Table A1).Thus, on the basis of this evidence, sedimentation at the boundary between water masses can explain this distinctive distribution [7].In addition, the countercurrent of the Kuroshio Current, which flows in a counterclockwise direction on the continental shelf (Figure 1b), may also contribute to the transport of fine particles enriched in the elements mentioned above [21].

Transfer of Materials Related to Mineral Deposits and Anthropogenic Activity
A number of silty sediment samples in Hiuchi-Nada are abundant in Cu (40.6-86.6 mg/kg) and are classified as outliers (Figure 4); these reflect the mining activity at the nearby Besshi Cu mines over the past 240 years [47].Similarly, the outlier observed within the Kii Channel sediments is related to activity at the Cu mine in the watershed of the Yoshino River.In addition, one outlier for Sb (1.89 mg/kg) is found near to the shore in southern Hiuchi-Nada.This enrichment of Sb suggests contamination from the Ichinose Sb mine, which is located close to the Besshi mine.
Osaka Bay is adjacent to Osaka City, which is the second largest city in Japan and has a long history of industrial development.Ohta et al. [32] found that Cu, Zn, Cd, Sn, Sb, Hg, and Pb are significantly enriched in stream sediments collected from the urban areas of Osaka Plain.Similarly, the silty sediments in Osaka Bay are highly enriched in Cr, Cu, Zn, As, Mo, Cd, Sn, Hg, Pb, and Bi (as shown by Cr, Cu, Zn, and As in Figure 4).Samples with the top 1% of concentrations for these elements, which are mostly classified as outliers, are found in Osaka Bay, and the area of high concentrations near the mouth of the Yodo and Yamato Rivers is clearly due to the discharge of contaminated sediments.Hoshika and Shiozawa [27], Hoshika et al. [48], and Manabe [49] reported the same results, showing large accumulations of Cr, Cu, Zn, Cd, Pb, and Hg in Osaka Bay.
A number of outliers of Cu, Zn, Cd, Sn, Pb, and Bi are also found in Harima-Nada, Hiroshima Bay and Suo-Nada (shown by Cu, Zn, and Cd in Figure 4).Their enrichment may be a result of the nearby metalliferous deposits and contamination by the coastal industrial zones [49,50].In particular, Ikuno mine is considered a possible source for the elements abundant in Harima-Nada because it is one of the largest hydrothermal-type polymetallic mines, bearing Cu, Zn, As, Cd, Sb, Sn, Pb, and Bi [31] (Figure 2).However, the spatial distributions of the outliers are different for each enriched element.A likely explanation for this is that the rivers flowing through the mineralized zones are small and have much smaller sediment discharges.Alternatively, the elements may have been released in the water during the oxidation of sulfide ores, and have subsequently absorbed into the surface sediment of the coastal seas [51].and Bi, as a result of early diagenetic processes in an oxic environment [7,42].In contrast, Cu, Sb, Cd, Hg, and U are enriched in deep-sea sediments during early diagenesis under anoxic conditions, as they are immobile in reducing waters [52].In addition, Cu and Hg are taken up by living organisms, transported to sediments, and bound within residual organic matter in bottom sediments [53][54][55].Results of this study reveal that Cu, Sb, Cd, Hg, and U are abundant in the silt and clay of the Tosa Basin, Hyuga Basin, and the Muroto Trough (shown by Cu, Sb and Cd in Figure 4).This suggests that these deep-sea basins are under anoxic conditions.However, outliers for these elements are scarce in the Pacific Ocean, in contrast to the Seto Inland Sea (Figure 4).Ikehara [25] reported that these deep-sea basins are hemi-pelagic and are partly associated with turbidites.Therefore, we assume that the enrichment of heavy metals is lost when gradual sedimentation is interrupted by turbidity flows, and that early diagenetic processes and the deposition of organic matter begin again after the turbidity flow event [56].As a result, no significant enrichment of heavy metals is visible in the Pacific Ocean.

Conclusions
The composition and transfer processes of coastal sea sediments are analyzed using a comprehensive geochemical database of many elements from the Seto Inland Sea and the Pacific Ocean in western Japan.The geochemical features of fine sands and silty-clayey sediments reflect those of stream sediments in the adjacent terrestrial areas.This pattern is comparable to the potential sediment yield from terrestrial areas and its supply to the adjacent marine environment.The spatial distribution of anomalous elemental concentrations in silty and clayey sediments suggests horizontal dispersion by a periodic current, as well as precipitation processes occurring at the boundary between water masses.However, the geochemical features of coarse sands and gravels differ significantly from those of stream sediments in the adjacent terrestrial areas, and fractionation of the mineralogical composition due to strong tidal currents, in addition to the denudation of Miocene-Pliocene basement rocks of the underwater topographic highs are possible reason for this inconsistency in geochemistry.As such, comprehensive geochemical maps of both the land and the sea are demonstrated as being useful for determining the source of fine marine sediments, as well as for tracking the transfer and deposition of silty-clayey sediments.

Figure 1 .
Figure 1.(a) Generalized location map of the study area; (b) Geographic map of the study area.Solid blue lines in terrestrial areas show major rivers.The abbreviations indicate the river names, as listed in TableA1.The blue colored area is the Seto Inland Sea.1: Kuroshio Current; 1': Kuroshio Counter-Current; 2: Tsushima-Current; 3: Oyashio-Current; 4: Liman-Current.Bathymetric depth contours are delineated using a dataset provided by the Japan Oceanographic Data Center.

Figure 1 .
Figure 1.(a) Generalized location map of the study area; (b) Geographic map of the study area.Solid blue lines in terrestrial areas show major rivers.The abbreviations indicate the river names, as listed in TableA1.The blue colored area is the Seto Inland Sea.1: Kuroshio Current; 1': Kuroshio Counter-Current; 2: Tsushima-Current; 3: Oyashio-Current; 4: Liman-Current.Bathymetric depth contours are delineated using a dataset provided by the Japan Oceanographic Data Center.

Figure 2 .
Figure 2. Geological map of the study area simplified from the Geological Map of Japan 1:1,000,000[17].

Figure 3 .
Figure 3. Sample locations for stream and coastal sea sediments.

Figure 3 .
Figure 3. Sample locations for stream and coastal sea sediments.

Figure 4 .
Figure 4. Spatial distributions of elemental concentrations in terrestrial and marine areas.The data are presented in wt % oxide (Al 2 O 3 , CaO, K 2 O, and T-Fe 2 O 3 ) and mg/kg element (Cr, Cu, Zn, As, Cd, Sb, Cs and La).Stars, circles, diamonds, squares and triangular symbols indicate major metalliferous deposits (as in Figure 1b).Cross and plus symbols indicate outliers obtained using Mahalanobis' generalized distances (details in text in Section 5.5.1).

Figure 5 .
Figure 5. Median concentrations of elements and oxides in sediments from the Seto Inland Sea and the Pacific Ocean.A, B, and C represent coarse sediment, fine sand, and silty-clayey sediments, respectively.η 2 (A), η 2 (B), and η 2 (A × B) are the size effects of regional differences, particle sizes, and the interaction effect, respectively, obtained from the two-way ANOVA (Table4).

Figure 5 .
Figure 5. Median concentrations of elements and oxides in sediments from the Seto Inland Sea and the Pacific Ocean.A, B, and C represent coarse sediment, fine sand, and silty-clayey sediments, respectively.

Figure 6 .
Figure 6.Enrichment factors of 49 elements and oxides in stream and marine sediments of (a) Seto Inland Sea and (b) Pacific Ocean, compared to the upper continental crust [36] (details in the text).
,b), which implies that they have different origins to the modern terrestrial

Figure 6 .
Figure 6.Enrichment factors of 49 elements and oxides in stream and marine sediments of (a) Seto Inland Sea and (b) Pacific Ocean, compared to the upper continental crust [36] (details in the text).

Figure 8 .
Figure 8.(a) Profiles of chemical compositions of coarse sediments in the narrow sea channels (around Aki-Nada, Bisan-Seto-Nada, and Awaji Island), Iyo-Nada, and the Bungo Channel, normalized to stream sediments in the Seto Inland Sea region; (b) Profiles of chemical compositions of coarse sediments in Ashizuri Spar, Muroto Spar and Tosa Bae, and fine sands in sea knolls normalized to stream sediments in the Pacific Ocean region.
sediments and fine sands in the Pacific Ocean Ashizuri Spar Muroto Spar and Tosa-Bae Fine sands in sea knolls a) Coarse sediments in the Seto Inland Sea Narrow sea channels Iyo-Nada Bungo-Channel

Figure 8 .
Figure 8.(a) Profiles of chemical compositions of coarse sediments in the narrow sea channels (around Aki-Nada, Bisan-Seto-Nada, and Awaji Island), Iyo-Nada, and the Bungo Channel, normalized to stream sediments in the Seto Inland Sea region; (b) Profiles of chemical compositions of coarse sediments in Ashizuri Spar, Muroto Spar and Tosa Bae, and fine sands in sea knolls normalized to stream sediments in the Pacific Ocean region.

5. 5 .
The Transfer of Silty and Clayey Sediments Influenced by Periodic Currents, Tidal Currents, and a Water Mass Boundary5.5.1.Mahalanobis' Generalized Distances for Finding Outliers

Figure 9 .Figure 9 .
Figure 9. Relationship between elemental concentrations and Cs concentration.Mahalanobis' generalized distances (D) increase from the central part to the outer part of the plot in an onion-like structure.Values of D ≤ 1.117, D ≤ 1.665, and D ≤ 2.146 indicate that 50%, 75%, and 90% of the total data points are included, respectively.5.5.2.Transport of Elements in Silt Derived from the Parent LithologyFigure4shows that the locations of the outliers calculated for Cr (cross symbols) correspond with a NNW-SSE distribution of silty-clayey sediments with high Cr concentrations in the Kii Channel.This distribution of Cr outliers indicates that silty particles that originated from ultramafic

5. 5 . 4 .
Geochemical Features of Silt and Clay in the Deep Sea Basins of the Pacific Ocean Silt and clay deposited in deep Pacific basins are typically enriched in MnO, V, Ni, Co, Mo, Sb, Pb,

Table 1 .
Sediment yield of each region and discharge to adjacent marine environment.

Table 3 .
Median elemental concentrations of marine sediments and stream sediments.

Table B1 .
Median elemental concentrations in stream sediments classified according to the dominant rock type in the drainage basin.