Application of a Self-Organizing Map of Isotopic and Chemical Data for the Identiﬁcation of Groundwater Recharge Sources in Nasunogahara Alluvial Fan, Japan

: Paddy rice fields on an alluvial fan not only use groundwater for irrigation but also play an important role as groundwater recharge sources. In this study, we investigated the spatial distribution of isotopic and hydrochemical compositions of groundwater in the Nasunogahara alluvial fan in Japan and applied a self-organizing map (SOM) to characterize the groundwater. The SOM assisted with the hydrochemical and isotopic interpretation of the groundwater in the fan, and clearly classified the groundwater into four groups reflecting the different origins. Two groundwater groups with lower isotopic ratios of water than the mean precipitation values in the fan were influenced by the infiltration of river water flowing from higher areas in the catchments and were differentiated from each other by their Na + and Cl − concentrations. A groundwater group with higher isotopic ratios was influenced by the infiltration of paddy irrigation water that had experienced evaporative isotopic enrichment. Groundwater in the fourth group, which was distributed in the upstream area of the where dairy farms showed little influence of recharge waters from paddy rice fields. The findings of this study will contribute to proper management of the groundwater resources in the fan. These results show that multiple tracers could be e ﬀ ectively used to evaluate the inﬂuence of di ﬀ erent groundwater recharge sources, including paddy water in rice cultivation areas, and the application of an SOM can provide understandable and readily visualized results to assist in the interpretation of the characteristics of groundwater in the fan. The ﬁndings of this study can therefore contribute to proper planning for regional groundwater use and the preservation of spring waters in the fan.


Introduction
Alluvial fans are important lowland areas in Japan because approximately 70% of the total land area of Japan is mountainous or hilly. Vast paddy rice fields are situated on alluvial fans, which occupy more than half (54%) of Japan's lowland area [1]. River water is usually diverted at the apexes of such alluvial fans and then delivered through a network of irrigation canals to each paddy rice field by gravity [2]. Alluvial fans are also vital aquifer systems that support local agricultural, industrial, and socio-economic development worldwide, wherever groundwater is utilized as an important water resource [3]. Aquifers of permeable sand and gravel formed in alluvial fans are vulnerable to contamination from the surface, and therefore interactions between surface water and groundwater can even affect hydrochemical processes (e.g., [4][5][6][7]). In addition, the infiltration of irrigation water contributes considerably to groundwater recharge [8] and also produces complex groundwater flows in paddy-dominant alluvial fans. For the sustainable use of groundwater resources in fans, it is important to understand these groundwater flows and recharge sources.
Environmental isotopes and hydrochemical data are powerful tools for investigating groundwater flow conditions (e.g., [9,10]), including hydrogeological investigations of alluvial fans (e.g., [8,11,12]). A multifaceted analysis using a suite of tracers (multi-tracer approach) can be effectively used to Groundwater was not sampled at the sites denoted by white circles in March 2019 because the wells were dry owing to a drop in the groundwater level. The land use map was created by using National Land Numerical Information [29]. The gray topographic base map in (c) was from the Geospatial Information Authority of Japan [30].

Hydrogeological Settings
The Nasunogahara alluvial fan consists mostly of Quaternary deposits overlying a basement composed of pre-Quaternary green tuff and mudstone (Shioya and Arakawa groups). In order of decreasing depth, the Sakaibayashi gravel formation, Tatenokawa tuff formation, Kuroiso volcanic breccia, Nabekake conglomerate, Kanemarubara gravel formation, Nasuno gravel formation, the Samui and Okuzawa gravel formations, and a flood plain gravel bed overlie the basement rocks ( Figure 2).
An unconfined aquifer is in the flood plain gravel bed and the Samui and Okuzawa formations; in this aquifer, the groundwater table is higher during summer and disappears during winter. A second unconfined aquifer in the Nasuno gravel formation supplies most of the groundwater used for irrigation; it typically has a high hydraulic conductivity of 10 −4 to 10 −3 m/s [31,32]. The Tatenokawa tuff formation is impermeable and the gravels of overlying formations cover a paleotopography of hills and valleys on the upper surface of the Tatenokawa tuff. The Sakaibayashi gravel formation Groundwater was not sampled at the sites denoted by white circles in March 2019 because the wells were dry owing to a drop in the groundwater level. The land use map was created by using National Land Numerical Information [29]. The gray topographic base map in (c) was from the Geospatial Information Authority of Japan [30].

Hydrogeological Settings
The Nasunogahara alluvial fan consists mostly of Quaternary deposits overlying a basement composed of pre-Quaternary green tuff and mudstone (Shioya and Arakawa groups). In order of decreasing depth, the Sakaibayashi gravel formation, Tatenokawa tuff formation, Kuroiso volcanic breccia, Nabekake conglomerate, Kanemarubara gravel formation, Nasuno gravel formation, the Samui and Okuzawa gravel formations, and a flood plain gravel bed overlie the basement rocks ( Figure 2).
An unconfined aquifer is in the flood plain gravel bed and the Samui and Okuzawa formations; in this aquifer, the groundwater table is higher during summer and disappears during winter. A second unconfined aquifer in the Nasuno gravel formation supplies most of the groundwater used for irrigation; it typically has a high hydraulic conductivity of 10 −4 to 10 −3 m/s [31,32]. The Tatenokawa tuff formation is impermeable and the gravels of overlying formations cover a paleotopography of hills and valleys on the upper surface of the Tatenokawa tuff. The Sakaibayashi gravel formation underlies the Tatenokawa tuff formation and is a confined aquifer. In this study, groundwater samples were collected from depths down to 30.3 m in wells installed in the first and second unconfined aquifers, as Water 2020, 12, 278 4 of 21 described in Section 3.1. The groundwater level of the fan begins to rise in April, reaches a maximum in October, and then declines through March, echoing changes in precipitation [24]. Groundwater levels are thus high from April to September, the paddy rice irrigation period, suggesting that the groundwater may be recharged by paddy irrigation water.
Water 2020, 12, x FOR PEER REVIEW 4 of 20 underlies the Tatenokawa tuff formation and is a confined aquifer. In this study, groundwater samples were collected from depths down to 30.3 m in wells installed in the first and second unconfined aquifers, as described in Section 3.1. The groundwater level of the fan begins to rise in April, reaches a maximum in October, and then declines through March, echoing changes in precipitation [24]. Groundwater levels are thus high from April to September, the paddy rice irrigation period, suggesting that the groundwater may be recharged by paddy irrigation water.

Figure 2.
Geological cross section of the Nasunogahara alluvial fan (modified from the Kanto Regional Agricultural Administration Office [31]).

Sample Acquisition
Many hand-dug wells for paddy irrigation (as well as for upland crop irrigation and domestic use) are distributed across the Nasunogahara alluvial fan, although there are a few wells in the upper part of the fan (Figure 1). We collected samples for the analyses of stable isotopic composition, radon radioisotope ( 222 Rn), and major ion concentrations of water from wells, river waters, and paddy waters, both during the irrigation period (IP, 15-18 May 2018) and during the non-irrigation period (NP, 4-7 March 2019). In the IP, groundwater was collected from 124 wells, but in the NP, some of the wells had dried up because of a drop in the groundwater level (white circles in Figure 1), so groundwater was collected from only 98 wells. The ground elevation of the wells ranged from 134 to 417 m a.s.l., and the well depth ranged from 3.3 to 30.3 m, with an average depth of 9.6 m below ground level. Toward the apex of the fan, the wells were deeper, and toward the toe of the fan, they were shallower. The mean depths of groundwater in the wells in the IP and NP were 3.8 and 2.5 m, respectively. Therefore, isotopic and chemical compositions of groundwaters were considered to be the same regardless of water depth. An acrylic well bailer was used to collect groundwater samples from the wells. Paddy waters were collected from the irrigation water inlet and outlet of seven paddy rice fields only in the IP. Spring waters were collected at six sites (SP1-6) in the IP and three sites (SP4-6) in the NP. Some springs dried up in the NP because of the drop in the groundwater level. A polyethylene cup or EVA (ethylene-vinyl acetate copolymer) bucket were used to sample the paddy and spring waters. River water samples were collected from each of the four rivers (Naka, Kuma, Sabi, and Houki rivers), either from the bank or from a bridge, using an EVA bucket. A relatively low rainfall was observed and the maximum daily rainfalls in and before the IP and NP were 8 and 14 mm, respectively. Therefore, the collected river waters were close to baseflow conditions. Rainwater was collected at site P1 (321 m a.s.l.; Figure 1) from 14 March 2018 to 21 August 2019 (527 days) and site P2 (179 m a.s.l.) from 9 August 2018 to 21 August 2019 (377 days) at a mean interval of 33 days in rain collectors designed, following recommendations described in International Atomic Energy Agency (IAEA)/Global Network of Isotopes in Precipitation (GNIP) [33], to prevent water reevaporation. The electrical conductivity (EC) of the sampled surface water and groundwater was measured by using a portable EC sensor (TOA-DKK, WM-32EP, Tokyo, Japan) immediately after the sampling. All water samples for chemical and isotopic analyses were transferred to clean  Figure 2. Geological cross section of the Nasunogahara alluvial fan (modified from the Kanto Regional Agricultural Administration Office [31]).

Sample Acquisition
Many hand-dug wells for paddy irrigation (as well as for upland crop irrigation and domestic use) are distributed across the Nasunogahara alluvial fan, although there are a few wells in the upper part of the fan (Figure 1). We collected samples for the analyses of stable isotopic composition, radon radioisotope ( 222 Rn), and major ion concentrations of water from wells, river waters, and paddy waters, both during the irrigation period (IP, 15-18 May 2018) and during the non-irrigation period (NP, 4-7 March 2019). In the IP, groundwater was collected from 124 wells, but in the NP, some of the wells had dried up because of a drop in the groundwater level (white circles in Figure 1), so groundwater was collected from only 98 wells. The ground elevation of the wells ranged from 134 to 417 m a.s.l., and the well depth ranged from 3.3 to 30.3 m, with an average depth of 9.6 m below ground level. Toward the apex of the fan, the wells were deeper, and toward the toe of the fan, they were shallower. The mean depths of groundwater in the wells in the IP and NP were 3.8 and 2.5 m, respectively. Therefore, isotopic and chemical compositions of groundwaters were considered to be the same regardless of water depth. An acrylic well bailer was used to collect groundwater samples from the wells. Paddy waters were collected from the irrigation water inlet and outlet of seven paddy rice fields only in the IP. Spring waters were collected at six sites (SP1-6) in the IP and three sites (SP4-6) in the NP. Some springs dried up in the NP because of the drop in the groundwater level. A polyethylene cup or EVA (ethylene-vinyl acetate copolymer) bucket were used to sample the paddy and spring waters. River water samples were collected from each of the four rivers (Naka, Kuma, Sabi, and Houki rivers), either from the bank or from a bridge, using an EVA bucket. A relatively low rainfall was observed and the maximum daily rainfalls in and before the IP and NP were 8 and 14 mm, respectively. Therefore, the collected river waters were close to baseflow conditions. Rainwater was collected at site P1 (321 m a.s.l.; Figure 1) from 14 March 2018 to 21 August 2019 (527 days) and site P2 (179 m a.s.l.) from 9 August 2018 to 21 August 2019 (377 days) at a mean interval of 33 days in rain collectors designed, following recommendations described in International Atomic Energy Agency (IAEA)/Global Network of Isotopes in Precipitation (GNIP) [33], to prevent water re-evaporation. The electrical conductivity (EC) of the sampled surface water and groundwater was measured by using a portable EC sensor (TOA-DKK, WM-32EP, Tokyo, Japan) immediately after the sampling. All water samples for chemical and isotopic analyses were transferred to clean polyethylene vials and brought to the laboratory of the Institute for Rural Engineering, National Agriculture, and Food Research Organization (Ibaraki, Japan), where they were filtered through syringe filters with a polytetrafluoroethylenemembrane (0.45 µm) and stored in a refrigerator until analysis.

Analysis of Environmental Isotopic and Hydrochemical Compositions
Major ions (Na + , K + , Mg 2+ , Ca 2+ , Cl − , NO 3 − , SO 4 2− ) were quantified using ion chromatography (TOA-DKK, ICA-2000, Tokyo, Japan) using multicomponent standard solutions. Nitrate (NO 3 − ) is reported as the nitrate-nitrogen (NO 3 -N) concentration. Alkalinity was determined using titration to pH 4.8 and is reported as HCO 3 − . For all water samples, errors in the cation-anion charge balance were less than 5%. 222 Rn, which is a radioactive gas generated by the decay of 226 Ra in geological strata, is soluble in water and has a half-life of 3.8 days. The 222 Rn concentration in groundwater finally reaches radioactive equilibrium. It reaches 92% of the value at secular equilibrium after two weeks underground. Using these characteristics, 222 Rn is often used as a target tracer to investigate the hydrological aspects of groundwater (e.g., estimation of infiltration from surface water to aquifers [34] and groundwater flow rate [35]), because far more 222 Rn is contained in groundwater than in surface water. The radioactivity of 222 Rn in samples was analyzed using a liquid scintillation counter (Packard, 2250 CA, CT, USA) after extraction with toluene within one week after collecting the water, and the concentrations of 222 Rn were calculated based on the count rate [36]. The stable isotopic compositions of the water were measured using a cavity ringdown spectrometer-based water isotope analyzer (Picarro, L2140-i, Santa Clara, CA, USA). They were reported with respect to Vienna Standard Mean Ocean Water (VSMOW) and are described here in "delta" notation as δ 18 O and δ 2 H (in parts per thousand, % ). Three working standards calibrated by international measurement standards were used for the analysis of δ 18 O and δ 2 H. The analytical precision in this study was ±0.05% for δ 18 O and ±0.50% for δ 2 H. The deuterium excess (d-excess, used to indicate the kinetic fractionation effect of evaporation [37]) of each water sample was calculated as d-excess = δ 2 H − 8δ 18 O [38]. The Mann-Whitney U test, which is applicable to non-normally distributed variables, was used to assess the independence of the measured isotopic and hydrochemical compositions between the IP and NP.

SOM
An SOM is an unsupervised training algorithm of an artificial neural network model developed by Kohonen [39]. An SOM can project high-dimensional target data onto a low-dimensional format (typically a two-dimensional grid map called a "feature map") [40]. Therefore, it is an effective and powerful analysis tool for the detection of data characteristics using pattern classification and visualization [16]. An SOM consists of an input layer and an output layer. The output layer often consists of a rectangular grid map with a hexagonal lattice ("nodes"). Each node in the output layer has a vector of coefficients associated with the input data, the so-called reference vector (also referred to as the weight vector and the codebook vector). The vectors can be obtained using iterative updates during a training phase consisting of three main procedures: competition between nodes, selection of a winning node, and updating of the reference vectors [40,41].
The determination of the SOM structure (calculation of the total number of nodes and the ratio of the number of map rows to the number of map columns) is an important step in the application of an SOM. If the number of nodes (i.e., map size) is too small, it might not explain some important differences that should be detected. Conversely, if the map size is too big, the differences between adjoining nodes are too small [42]. The optimal number of nodes in an SOM can be selected by using the heuristic rule m = 5 √ n, where m denotes the total number of nodes and n is the number of samples in the data set [40,43]. The ratio of the number of rows to the number of columns in the feature map is determined by calculating the square root of the ratio of the largest eigenvalue of the correlation matrix of the input data to the second-largest eigenvalue [19]. In this study, the eigenvalues were obtained using principal component analysis of the input dataset. Normalization of the data was necessary prior to the application of the SOM to ensure that all values of the chemical parameters were given similar importance. The results of the SOM application are sensitive to the method used for data pre-processing because the SOM is trained to organize itself according to the Euclidean distances between input data [40]. In this study, all variables of the input data were transformed to have a range of 0 to 1 such that they would have equal importance in the formation of the SOM.
For clustering the trained SOM, several methods, each with advantages and disadvantages, are available. Agglomerative clustering and partition clustering result in clear groupings on the trained map, whereas application of a U-matrix (unified distance matrix), which is commonly used to create a rough visualization of the classification, may not define clear boundaries (e.g., [44]). Ward's agglomerative hierarchical clustering (Ward's method) is an easy way to cluster nodes on the trained map because it is not necessary in this method to determine the optimum number of groups before performing the clustering calculations [42]. Ward's method has been commonly used to subdivide the trained map into several groups according to the similarity of the reference vectors of the nodes (e.g., [43,45]). Thus, Ward's method was applied in this study to identify groups on the trained map.
In this study, the input layer of the SOM consisted of 11 chemical and isotopic parameters (δ 2 H, δ 18 O, 222 Rn, and the major ion concentrations) of groundwater samples collected at 124 and 98 wells in the IP and NP, respectively. In accordance with the heuristic rules described above, a feature map with 72 (9 × 8) nodes was selected. The "kohonen" package in the R statistical software environment (version 3.5.1, Vienna, Austria) [46] was used to implement the SOM. Default values of the package parameters (e.g., learning rate and size of the neighborhood) were used. The observed hydrochemical and isotopic data of the groundwater in the fan were preprocessed by the proposed SOM, and then the groundwater data were spatially clustered into different homogeneous groups on the trained map. To identify the groundwater recharge sources, the hydrochemical and isotopic compositions of the groups were plotted on various types of diagrams (i.e., hexa-diagrams, a trilinear diagram, a δ 18 O-δ 2 H scatter diagram, and spatial distribution maps) and compared.

Groundwater Level
We mapped groundwater levels, depth to the water table, and cross-sections of the water table in Figure 3. The groundwater table was almost parallel to the ground surface of the fan, which sloped toward the southeast, and the groundwater flowed in the same direction. were given similar importance. The results of the SOM application are sensitive to the method used for data pre-processing because the SOM is trained to organize itself according to the Euclidean distances between input data [40]. In this study, all variables of the input data were transformed to have a range of 0 to 1 such that they would have equal importance in the formation of the SOM. For clustering the trained SOM, several methods, each with advantages and disadvantages, are available. Agglomerative clustering and partition clustering result in clear groupings on the trained map, whereas application of a U-matrix (unified distance matrix), which is commonly used to create a rough visualization of the classification, may not define clear boundaries (e.g., [44]). Ward's agglomerative hierarchical clustering (Ward's method) is an easy way to cluster nodes on the trained map because it is not necessary in this method to determine the optimum number of groups before performing the clustering calculations [42]. Ward's method has been commonly used to subdivide the trained map into several groups according to the similarity of the reference vectors of the nodes (e.g., [43,45]). Thus, Ward's method was applied in this study to identify groups on the trained map.
In this study, the input layer of the SOM consisted of 11 chemical and isotopic parameters (δ 2 H, δ 18 O, 222 Rn, and the major ion concentrations) of groundwater samples collected at 124 and 98 wells in the IP and NP, respectively. In accordance with the heuristic rules described above, a feature map with 72 (9 × 8) nodes was selected. The "kohonen" package in the R statistical software environment (version 3.5.1, Vienna, Austria) [46] was used to implement the SOM. Default values of the package parameters (e.g., learning rate and size of the neighborhood) were used. The observed hydrochemical and isotopic data of the groundwater in the fan were preprocessed by the proposed SOM, and then the groundwater data were spatially clustered into different homogeneous groups on the trained map. To identify the groundwater recharge sources, the hydrochemical and isotopic compositions of the groups were plotted on various types of diagrams (i.e., hexa-diagrams, a trilinear diagram, a δ 18 Oδ 2 H scatter diagram, and spatial distribution maps) and compared.

Groundwater Level
We mapped groundwater levels, depth to the water table, and cross-sections of the water table in Figure 3. The groundwater table was almost parallel to the ground surface of the fan, which sloped toward the southeast, and the groundwater flowed in the same direction.  Houki R. Naka R.
Houki R. Naka R.  The contour line of the ground surface was modified from the Geospatial Information Authority of Japan [30]. The depth to the groundwater table in the upper to middle part of the fan (from the apex to about 260 m above sea level) was relatively large. In this area, the Sabi and Kuma rivers were losing streams and hence recharged the aquifer. The terrace cliffs on the right bank of the Naka River along the eastern edge of the fan varied from several meters to several tens of meters in height (Figure 3c), and therefore, the river bed of the Naka River was lower than the groundwater table. This means that there was no recharge from the Naka river to the aquifer. The riverbed of the Houki River was about 10 m higher than the groundwater surface near the middle part of the fan. The Houki River recharged the aquifer in this part and the infiltrated river water flowed in the downstream direction in accordance with hydraulic gradient, whereas the groundwater table became higher toward the toe of the fan. The lower part of the fan, where the water level became closer to the ground surface, corresponded to the groundwater discharge area.

Chemical Compositions
Analytical results of the chemical and isotopic parameters of the groundwater, along with those of spring water, paddy water, river water, and rainwater, are shown in Table 1. Four groups of the groundwater were classified using the SOM, as described later in Section 4.4. The chemical compositions of the groundwater, spring water, and river water are shown as a Piper (trilinear) diagram in Figure 4 and as hexa-diagrams in Figure 5. The Ca-HCO 3 type typically indicates shallow fresh groundwater, while Na-HCO 3 type typically indicates weathering and ion exchange processes. As shown by the key diagram and hexa-diagrams, some groundwater samples had a Ca-HCO 3 -to-Ca-SO 4 -Cl-type composition, whereas the Ca-HCO 3 -type composition was dominant. Groundwater samples were plotted near the center of the cation ternary diagram (Figure 4), suggesting little variability in the relative proportions of these cations. The anion ternary diagram indicated larger variability among samples: Groundwater samples were mostly HCO 3 − -type water, whereas some groundwater samples were scattered between HCO 3 − and SO 4 2− types.
Water 2020, 12, x FOR PEER REVIEW 7 of 20 The depth to the groundwater table in the upper to middle part of the fan (from the apex to about 260 m above sea level) was relatively large. In this area, the Sabi and Kuma rivers were losing streams and hence recharged the aquifer. The terrace cliffs on the right bank of the Naka River along the eastern edge of the fan varied from several meters to several tens of meters in height (Figure 3c), and therefore, the river bed of the Naka River was lower than the groundwater table. This means that there was no recharge from the Naka river to the aquifer. The riverbed of the Houki River was about 10 m higher than the groundwater surface near the middle part of the fan. The Houki River recharged the aquifer in this part and the infiltrated river water flowed in the downstream direction in accordance with hydraulic gradient, whereas the groundwater table became higher toward the toe of the fan. The lower part of the fan, where the water level became closer to the ground surface, corresponded to the groundwater discharge area.

Chemical Compositions
Analytical results of the chemical and isotopic parameters of the groundwater, along with those of spring water, paddy water, river water, and rainwater, are shown in Table 1. Four groups of the groundwater were classified using the SOM, as described later in Section 4.4. The chemical compositions of the groundwater, spring water, and river water are shown as a Piper (trilinear) diagram in Figure 4 and as hexa-diagrams in Figure 5. The Ca-HCO3 type typically indicates shallow fresh groundwater, while Na-HCO3 type typically indicates weathering and ion exchange processes. As shown by the key diagram and hexa-diagrams, some groundwater samples had a Ca-HCO3-to-Ca-SO4-Cl-type composition, whereas the Ca-HCO3-type composition was dominant. Groundwater samples were plotted near the center of the cation ternary diagram (Figure 4), suggesting little variability in the relative proportions of these cations. The anion ternary diagram indicated larger variability among samples: Groundwater samples were mostly HCO3 − -type water, whereas some groundwater samples were scattered between HCO3 − and SO4 2− types.
EC: electrical conductivity, IP: irrigation period, NP: non-irrigation period, -: not analyzed, 25th and 75th: 25th and 75th percentiles, mean values of rainwater were weighted using the precipitation amount.
Water 2020, 12, 278 9 of 21 Water 2020, 12, x FOR PEER REVIEW 10 of 20 Figure 5. Hexa-diagrams of groundwater, spring water, river water, and paddy water. Table 2 shows the matrix of the Pearson correlation coefficients for EC and major ions of all groundwater samples. Ca 2+ , which is highly correlated with EC, is a major controlling factor of EC. A strong positive correlation between Ca 2+ and HCO3 − indicates that a dominant water type of the groundwaters in this aquifer was Ca-HCO3 type. Concentrations of Ca 2+ and HCO3 − in aquifers of alluvial fans increased due to mineralization induced by a water-rock interaction during the groundwater flow process (e.g., [47]). Figure 6 shows the relationships between EC, Ca 2+ , and HCO3 − of the groundwater with groundwater level. These values did not increase with the downstream direction in the fan. The spatial distribution of the hexa-diagrams in Figure 5 indicates that the chemical compositions varied greatly depending on the location, especially the distances from the rivers.  Figure 5. Hexa-diagrams of groundwater, spring water, river water, and paddy water. Table 2 shows the matrix of the Pearson correlation coefficients for EC and major ions of all groundwater samples. Ca 2+ , which is highly correlated with EC, is a major controlling factor of EC. A strong positive correlation between Ca 2+ and HCO 3 − indicates that a dominant water type of the groundwaters in this aquifer was Ca-HCO 3 type. Concentrations of Ca 2+ and HCO 3 − in aquifers of alluvial fans increased due to mineralization induced by a water-rock interaction during the groundwater flow process (e.g., [47]). Figure 6 shows the relationships between EC, Ca 2+ , and HCO 3 − of the groundwater with groundwater level. These values did not increase with the downstream direction in the fan. The spatial distribution of the hexa-diagrams in Figure 5 indicates that the chemical compositions varied greatly depending on the location, especially the distances from the rivers.   Figure 7 illustrates the spatial distributions of EC and the concentrations of the three characteristic ions Ca 2+ , SO4 2− , and Cl − of groundwater in the IP and NP. The spatial distribution of EC, which was low along the Sabi and Kuma rivers in the central part of the fan, was very consistent with previous survey results [25,48]. Along these two rivers, Ca 2+ , like EC, was relatively low, whereas SO4 2− was relatively high. In contrast, Cl − concentrations were higher along the western edge of the fan near the Houki River than elsewhere in the fan. As is shown by the hexa-diagrams and ternary diagrams, the Naka, Sabi, and Kuma River waters were characterized by large SO4 2− contents, whereas the Houki River water samples had higher Na + and Cl − concentrations than the waters from the other rivers. Spring waters were plotted between the HCO3 − and SO4 2− types on the anion ternary diagram; SP1, SP2, and SP3 were characterized by larger SO4 2− contents, as indicated by their hexadiagram shapes. No significant seasonal difference between the IP and NP was observed in the chemical compositions of the river and spring waters.

Isotopic Compositions
The relationship between δ 18 O and δ 2 H (δ 18 O-δ 2 H scatter diagram) in all water samples is shown in Figure 8. In rainwater collected at sites P1 and P2, δ 18 O and δ 2 H varied over relatively wide ranges, from −12.9‰ to −3.1‰ and from −85‰ to −5‰, respectively. The local meteoric water line (LMWL) in this study area was determined to be δ 2 H = 8.48δ 18 O + 20.09. The slope of the LMWL was slightly   Figure 7 illustrates the spatial distributions of EC and the concentrations of the three characteristic ions Ca 2+ , SO 4 2− , and Cl − of groundwater in the IP and NP. The spatial distribution of EC, which was low along the Sabi and Kuma rivers in the central part of the fan, was very consistent with previous survey results [25,48]. Along these two rivers, Ca 2+ , like EC, was relatively low, whereas SO 4 2− was relatively high. In contrast, Cl − concentrations were higher along the western edge of the fan near the Houki River than elsewhere in the fan. As is shown by the hexa-diagrams and ternary diagrams, the Naka, Sabi, and Kuma River waters were characterized by large SO 4 2− contents, whereas the Houki River water samples had higher Na + and Cl − concentrations than the waters from the other rivers. Spring waters were plotted between the HCO 3 − and SO 4 2− types on the anion ternary diagram; SP1, SP2, and SP3 were characterized by larger SO 4 2− contents, as indicated by their hexa-diagram shapes.
No significant seasonal difference between the IP and NP was observed in the chemical compositions of the river and spring waters.
Water 2020, 12, x FOR PEER REVIEW 11 of 20 Figure 6. Relationships between EC, Ca 2+ , and HCO3 − with groundwater level. Figure 7 illustrates the spatial distributions of EC and the concentrations of the three characteristic ions Ca 2+ , SO4 2− , and Cl − of groundwater in the IP and NP. The spatial distribution of EC, which was low along the Sabi and Kuma rivers in the central part of the fan, was very consistent with previous survey results [25,48]. Along these two rivers, Ca 2+ , like EC, was relatively low, whereas SO4 2− was relatively high. In contrast, Cl − concentrations were higher along the western edge of the fan near the Houki River than elsewhere in the fan. As is shown by the hexa-diagrams and ternary diagrams, the Naka, Sabi, and Kuma River waters were characterized by large SO4 2− contents, whereas the Houki River water samples had higher Na + and Cl − concentrations than the waters from the other rivers. Spring waters were plotted between the HCO3 − and SO4 2− types on the anion ternary diagram; SP1, SP2, and SP3 were characterized by larger SO4 2− contents, as indicated by their hexadiagram shapes. No significant seasonal difference between the IP and NP was observed in the chemical compositions of the river and spring waters.

Isotopic Compositions
The relationship between δ 18 O and δ 2 H (δ 18 O-δ 2 H scatter diagram) in all water samples is shown in Figure 8. In rainwater collected at sites P1 and P2, δ 18 O and δ 2 H varied over relatively wide ranges, from −12.9‰ to −3.1‰ and from −85‰ to −5‰, respectively. The local meteoric water line (LMWL) in this study area was determined to be δ 2 H = 8.48δ 18 O + 20.09. The slope of the LMWL was slightly   The upper and lower maps show the data in the IP and NP, respectively.

Isotopic Compositions
The relationship between δ 18 O and δ 2 H (δ 18 O-δ 2 H scatter diagram) in all water samples is shown in Figure 8. In rainwater collected at sites P1 and P2, δ 18 O and δ 2 H varied over relatively wide ranges, from −12.9% to −3.1% and from −85% to −5% , respectively. The local meteoric water line (LMWL) in this study area was determined to be δ 2 H = 8.48δ 18 O + 20.09. The slope of the LMWL was slightly larger than the global meteoric water line (GMWL) slope of 8 [49] and the slope of the meteoric water line on the Pacific Ocean side in Japan (7.8) [50]. The weighted means of δ 18 O and δ 2 H (weighted using the precipitation amount during the sampling time interval) were −8.3% and −54% at P1 and P2, respectively; these values, which were plotted below the LMWL, reflect the low d-excess value of rainwater in summer in Japan [51] because the region of the Nasunogahara fan receives more precipitation in summer than in winter. compositions in groundwater [54]. In this study, paddy waters collected at paddy field outlets had higher isotopic ratios than those collected at the inlets because evaporative isotopic enrichment occurred during the residence time of the water in the paddy rice fields. The slope of the regression line fitted to all groundwater data was 5.04, between those of the LMWL and the PWRL (Figure 8b). The isotopic ratios of the river water samples were smaller than those of groundwater samples and plotted near the LMWL.
The spatial distributions of δ 18 O, d-excess, and the 222 Rn concentration in groundwater in the IP and NP are shown in Figure 9. The groundwater along the Sabi and Kuma rivers had comparatively low δ 18 O (less than −8.5‰) and high d-excess values (more than 12.0‰), whereas groundwater in the southeastern part of the fan had high δ 18 O (more than −8.0‰) and low d-excess values (less than 11.0‰). Compared with the IP, the domain with a high δ 18 O and low d-excess values in the southeastern fan shrank in the NP. Spatial differences in the 222 Rn concentration in groundwater were indistinct, although concentrations in the NP were higher than those in the IP. The 222 Rn concentration in waters of the four rivers was close to zero (Table 1).    The slope of the regression line between δ 18 O and δ 2 H of paddy waters, denoted as PWRL in Figure 8, was 3.6. The plots of the paddy water data deviated from the LMWL and the d-excess value of the paddy water (−2.2% ) was extremely low compared to those of the groundwater, rainwater, and river water (Table 1). Generally, kinetic fractionation of oxygen and hydrogen isotopes during evaporation causes the remaining water to become enriched in the heavier isotopes ( 18 O and 2 H); as a result, trends with slopes ranging from 4 to 7 are seen on the δ 18 O-δ 2 H scatter diagram, whereas the equilibrium isotopic fractionation line lies near the LMWL and has a slope close to 8 [52]. If there is no rain, irrigation water is supplied to the paddy fields almost every day and is affected by evaporative enrichment, even in just one day. Paddy waters affected by evaporative isotopic enrichment would have higher isotopic ratios and lower d-excess values [53] and effects the isotopic compositions in groundwater [54]. In this study, paddy waters collected at paddy field outlets had higher isotopic ratios than those collected at the inlets because evaporative isotopic enrichment occurred during the residence time of the water in the paddy rice fields. The slope of the regression line fitted to all groundwater data was 5.04, between those of the LMWL and the PWRL (Figure 8b). The isotopic ratios of the river water samples were smaller than those of groundwater samples and plotted near the LMWL.
The spatial distributions of δ 18 O, d-excess, and the 222 Rn concentration in groundwater in the IP and NP are shown in Figure 9. The groundwater along the Sabi and Kuma rivers had comparatively low δ 18 O (less than −8.5% ) and high d-excess values (more than 12.0% ), whereas groundwater in the southeastern part of the fan had high δ 18 O (more than −8.0% ) and low d-excess values (less than 11.0% ). Compared with the IP, the domain with a high δ 18 O and low d-excess values in the southeastern fan shrank in the NP. Spatial differences in the 222 Rn concentration in groundwater were indistinct, although concentrations in the NP were higher than those in the IP. The 222 Rn concentration in waters of the four rivers was close to zero (Table 1). Water 2020, 12, x FOR PEER REVIEW 13 of 20        The SOM nodes were subdivided into four groups (groups 1 to 4, comprising 34, 73, 72, and 43 groundwater samples, respectively) according to the dendrogram produced by the hierarchical cluster analysis (Figure 11). Figure 12 shows the spatial distributions of the color-coded hexa-diagrams of the four groups in the IP and NP. Group 1 was mostly distributed near the Houki River, which flows along the western edge of the fan. Group 2 was distributed in the central and lower parts of the fan. Group 3 was distributed along the Sabi and Kuma rivers. Group 4 was distributed across the middle part of the fan. The SOM nodes were subdivided into four groups (groups 1 to 4, comprising 34, 73, 72, and 43 groundwater samples, respectively) according to the dendrogram produced by the hierarchical cluster analysis (Figure 11). Figure 12 shows the spatial distributions of the color-coded hexadiagrams of the four groups in the IP and NP. Group 1 was mostly distributed near the Houki River, which flows along the western edge of the fan. Group 2 was distributed in the central and lower parts of the fan. Group 3 was distributed along the Sabi and Kuma rivers. Group 4 was distributed across the middle part of the fan.   The SOM nodes were subdivided into four groups (groups 1 to 4, comprising 34, 73, 72, and 43 groundwater samples, respectively) according to the dendrogram produced by the hierarchical cluster analysis (Figure 11). Figure 12 shows the spatial distributions of the color-coded hexadiagrams of the four groups in the IP and NP. Group 1 was mostly distributed near the Houki River, which flows along the western edge of the fan. Group 2 was distributed in the central and lower parts of the fan. Group 3 was distributed along the Sabi and Kuma rivers. Group 4 was distributed across the middle part of the fan.   The EC of group 3 groundwater was lower than that of the other groups (Table 1), indicating that it contained lower amounts of dissolved components. Groundwater in groups 1 and 4 had higher EC than groundwater in groups 2 and 3. In particular, concentrations of Na + and Cl − were notably higher in group 1, whereas those of Ca 2+ and HCO 3 − were higher in group 4. The mean NO 3 -N concentration was the highest in group 4, but all groundwater samples in this study had nitrate concentrations below the standard level of 10 mg/L NO 3 -N in Japan. The width of each hexa-diagram is an approximate indication of the total ionic content. Thus, the hexa-diagram width of group 3, which had a low EC, was small, whereas the hexa-diagram widths of groups 1 and 4, which had higher ECs, were large ( Figure 12). Chemical compositions of groundwaters in the four groups, spring water, and river water are shown as a trilinear diagram in Figure 13. Groundwater samples of group 4 were mostly HCO 3 − -type water, whereas those of groups 1 to 3 tended to be scattered between HCO 3 − and SO 4 2−

SOM and Clustering Results
types. In addition, the relative abundance of Cl − in group 1 samples was somewhat large.
Water 2020, 12, x FOR PEER REVIEW 15 of 20 The EC of group 3 groundwater was lower than that of the other groups (Table 1), indicating that it contained lower amounts of dissolved components. Groundwater in groups 1 and 4 had higher EC than groundwater in groups 2 and 3. In particular, concentrations of Na + and Cl − were notably higher in group 1, whereas those of Ca 2+ and HCO3 − were higher in group 4. The mean NO3-N concentration was the highest in group 4, but all groundwater samples in this study had nitrate concentrations below the standard level of 10 mg/L NO3-N in Japan. The width of each hexa-diagram is an approximate indication of the total ionic content. Thus, the hexa-diagram width of group 3, which had a low EC, was small, whereas the hexa-diagram widths of groups 1 and 4, which had higher ECs, were large ( Figure 12). Chemical compositions of groundwaters in the four groups, spring water, and river water are shown as a trilinear diagram in Figure 13. Groundwater samples of group 4 were mostly HCO3 − -type water, whereas those of groups 1 to 3 tended to be scattered between HCO3 − and SO4 2− types. In addition, the relative abundance of Cl − in group 1 samples was somewhat large. The δ 18 O-δ 2 H scatter diagram of groundwater in the four groups is shown in Figure 14. The groundwater in each of the four groups had characteristic δ 18 O and δ 2 H values: group 2 waters had the highest values, followed by waters of groups 4 and 1, and group 3 waters had the lowest values ( Figure 14, Table 1). Group 2 waters had the lowest d-excess value and group 3 waters had the highest d-excess value.

Influence of Recharge Sources on Groundwater Hydrochemical and Isotopic Compositions
It has been pointed out that shallow groundwater in the Nasunogahara alluvial fan is recharged by multiple sources, including the infiltration of rainwater, river water, and paddy irrigation water (e.g., [24]). The isotopic ratios of groundwater in this study, which were plotted as an elongated cluster between the LMWL and the PWRL in the δ 18 O-δ 2 H diagram (Figure 8), also indicate a mixture of groundwaters recharged from three different sources. The fact that the isotopic ratios of the groundwater were plotted between the LMWL and PWRL is strong evidence that groundwater in the fan was influenced by the infiltration of not only rainwater but also paddy irrigation water ( Figure  8). In addition, the comparison between groundwater levels and ground surface elevations indicated that there was infiltration from the Houki, Sabi, and Kuma rivers to the aquifer (Figure 3). The isotopic ratios of river waters (mean values of four rivers of −9.9‰ for δ 18 O and −63‰ for δ 2 H) were much less than the weighted means of rainwater (see Section 4.3) (Figure 8). The reason for the lower ratios of these river waters was that the rivers gathered precipitation that falls at high altitude in the mountains, which had lower isotopic ratios because of the isotopic altitude effect. The isotopic ratios of some groundwater samples were lower than the weighted mean rainwater δ 18 O and δ 2 H values, suggesting that these samples represented a mixture of waters with low isotopic ratios. We therefore inferred that the groundwaters in the western part and upper part of the fan were affected by the infiltration of waters from the Houki River, and the Kuma and Sabi rivers, respectively. The river waters contained baseflows supplied from the mountainous watersheds, but 222 Rn concentrations of river waters were low. This indicates that dissolved 222 Rn in river waters was lost through dispersion to the atmosphere and radioactive decay before reaching the fan. It had been expected that the 222 Rn concentration of the groundwater near the rivers decreased due to river water infiltration, but the concentrations were not lower than those far away from the rivers (Figure 9). This means that the residence time of infiltrated waters to reaching sampling wells was more than two weeks, which were required to reach radioactive equilibrium.
The hydrochemical data also implied that infiltration of river water recharged some groundwaters. The groundwater around the Sabi and Kuma rivers, the EC, and some dissolved ion concentrations, such as that of Ca 2+ , were relatively low, but SO4 2-was relatively high; thus, we can

Influence of Recharge Sources on Groundwater Hydrochemical and Isotopic Compositions
It has been pointed out that shallow groundwater in the Nasunogahara alluvial fan is recharged by multiple sources, including the infiltration of rainwater, river water, and paddy irrigation water (e.g., [24]). The isotopic ratios of groundwater in this study, which were plotted as an elongated cluster between the LMWL and the PWRL in the δ 18 O-δ 2 H diagram (Figure 8), also indicate a mixture of groundwaters recharged from three different sources. The fact that the isotopic ratios of the groundwater were plotted between the LMWL and PWRL is strong evidence that groundwater in the fan was influenced by the infiltration of not only rainwater but also paddy irrigation water (Figure 8). In addition, the comparison between groundwater levels and ground surface elevations indicated that there was infiltration from the Houki, Sabi, and Kuma rivers to the aquifer ( Figure 3). The isotopic ratios of river waters (mean values of four rivers of −9.9% for δ 18 O and −63% for δ 2 H) were much less than the weighted means of rainwater (see Section 4.3) (Figure 8). The reason for the lower ratios of these river waters was that the rivers gathered precipitation that falls at high altitude in the mountains, which had lower isotopic ratios because of the isotopic altitude effect. The isotopic ratios of some groundwater samples were lower than the weighted mean rainwater δ 18 O and δ 2 H values, suggesting that these samples represented a mixture of waters with low isotopic ratios. We therefore inferred that the groundwaters in the western part and upper part of the fan were affected by the infiltration of waters from the Houki River, and the Kuma and Sabi rivers, respectively. The river waters contained baseflows supplied from the mountainous watersheds, but 222 Rn concentrations of river waters were low. This indicates that dissolved 222 Rn in river waters was lost through dispersion to the atmosphere and radioactive decay before reaching the fan. It had been expected that the 222 Rn concentration of the groundwater near the rivers decreased due to river water infiltration, but the concentrations were not lower than those far away from the rivers (Figure 9). This means that the residence time of infiltrated waters to reaching sampling wells was more than two weeks, which were required to reach radioactive equilibrium.
The hydrochemical data also implied that infiltration of river water recharged some groundwaters. The groundwater around the Sabi and Kuma rivers, the EC, and some dissolved ion concentrations, such as that of Ca 2+ , were relatively low, but SO 4 2− was relatively high; thus, we can infer that the groundwaters were influenced by the infiltration of waters from these rivers, which had similar compositions, and this inference was supported by the low isotopic ratios and high d-excess values of the groundwaters. The groundwater in the western part of the fan had relatively high Na + and Cl − concentrations, reflecting the infiltration of Houki River water, which was also characterized by relatively high Na + and Cl − concentrations ( Figures 5 and 7). Similarly, the relatively low isotopic ratios and high EC of the groundwaters suggested a large contribution of infiltration from the Houki River (Figures 7 and 9). Thus, groundwater in the fan reflected three different recharge sources and the fraction of the contribution from each source varied by location.

Characterization of Groundwater Using SOM
The SOM classified the groundwater in the fan into four groups with different hydrochemical and isotopic compositions. The groundwater of each group was inferred to be affected by different recharge sources. Group 1 groundwaters, distributed in the western part of the fan, had low isotopic ratios and high Na + and Cl − , suggesting the infiltration of the Houki River. The relatively high isotopic ratios and low d-excess values of group 2 groundwaters indicated a large recharge contribution due to infiltration from paddy rice fields. Group 3 groundwaters around the Sabi and Kuma rivers were inferred to be influenced by infiltration of these rivers, which had similar chemical and isotopic compositions: low EC, low dissolved ion concentrations excluding SO 4 2− , and low stable isotopic ratios. Group 4 groundwaters, distributed at higher elevations than group 2 groundwaters, had lower isotopic ratios than the weighted means of rainwater, suggesting that these samples were affected by a mixture of Sabi and Kuma river waters. However, the relatively large EC and ionic contents indicated a smaller contribution of infiltration from these river waters compared to group 3 groundwaters. In the upper part of the fan, there are few paddy rice fields but many livestock farms (Figure 1). Reflecting this difference in land use, NO 3 -N concentrations in group 4 groundwaters were higher than those in the other groups (Table 1). We inferred that these higher concentrations reflected the high nitrogen load associated with livestock farming (e.g., [27]). It can be readily seen that the contribution of infiltration from paddy rice field to groundwater recharge was smaller for group 4 than for group 2 groundwaters. In contrast, smaller ionic contents, including NO 3 -N, in group 2 groundwaters indicate a dilution due to infiltration of paddy waters. The application and results of the SOM assisted the interpretation of the difference in hydrochemical and isotopic compositions depending on the location and the influence of multiple recharge sources on the groundwater. Data that produce a scattered distribution in the trilinear diagram, hexa-diagrams, and δ 18 O-δ 2 H diagram make it difficult to resolve ambiguous boundaries and to produce reasonable classification of the groundwaters influenced by multiple recharge sources, especially for the shallow groundwater with similar chemical and isotopic compositions. The SOM application could automatically establish the four different groups, independent of human-subjective criteria, and summarized the spatial distribution of each group with the readily understandable visualization (Figure 15).
Water 2020, 12, x FOR PEER REVIEW 17 of 20 infer that the groundwaters were influenced by the infiltration of waters from these rivers, which had similar compositions, and this inference was supported by the low isotopic ratios and high d-excess values of the groundwaters. The groundwater in the western part of the fan had relatively high Na + and Cl − concentrations, reflecting the infiltration of Houki River water, which was also characterized by relatively high Na + and Cl − concentrations ( Figures 5 and 7). Similarly, the relatively low isotopic ratios and high EC of the groundwaters suggested a large contribution of infiltration from the Houki River (Figures 7 and 9). Thus, groundwater in the fan reflected three different recharge sources and the fraction of the contribution from each source varied by location.

Characterization of Groundwater Using SOM
The SOM classified the groundwater in the fan into four groups with different hydrochemical and isotopic compositions. The groundwater of each group was inferred to be affected by different recharge sources. Group 1 groundwaters, distributed in the western part of the fan, had low isotopic ratios and high Na + and Cl − , suggesting the infiltration of the Houki River. The relatively high isotopic ratios and low d-excess values of group 2 groundwaters indicated a large recharge contribution due to infiltration from paddy rice fields. Group 3 groundwaters around the Sabi and Kuma rivers were inferred to be influenced by infiltration of these rivers, which had similar chemical and isotopic compositions: low EC, low dissolved ion concentrations excluding SO4 2-, and low stable isotopic ratios. Group 4 groundwaters, distributed at higher elevations than group 2 groundwaters, had lower isotopic ratios than the weighted means of rainwater, suggesting that these samples were affected by a mixture of Sabi and Kuma river waters. However, the relatively large EC and ionic contents indicated a smaller contribution of infiltration from these river waters compared to group 3 groundwaters. In the upper part of the fan, there are few paddy rice fields but many livestock farms ( Figure 1). Reflecting this difference in land use, NO3-N concentrations in group 4 groundwaters were higher than those in the other groups (Table 1). We inferred that these higher concentrations reflected the high nitrogen load associated with livestock farming (e.g., [27]). It can be readily seen that the contribution of infiltration from paddy rice field to groundwater recharge was smaller for group 4 than for group 2 groundwaters. In contrast, smaller ionic contents, including NO3-N, in group 2 groundwaters indicate a dilution due to infiltration of paddy waters.
The application and results of the SOM assisted the interpretation of the difference in hydrochemical and isotopic compositions depending on the location and the influence of multiple recharge sources on the groundwater. Data that produce a scattered distribution in the trilinear diagram, hexa-diagrams, and δ 18 O-δ 2 H diagram make it difficult to resolve ambiguous boundaries and to produce reasonable classification of the groundwaters influenced by multiple recharge sources, especially for the shallow groundwater with similar chemical and isotopic compositions. The SOM application could automatically establish the four different groups, independent of humansubjective criteria, and summarized the spatial distribution of each group with the readily understandable visualization (Figure 15).   • EC and stable isotopic ratios were low.
• Groundwater was affected by infiltration from the Sabi and Kuma Rivers.
• The area expanded during non-irrigation period.
• Stable isotopic ratios were higher.
• Contribution of recharge from paddy rice field was large.
• The area shrank during non-irrigation period.
• Na + and Cl − concentrations were high.
• Groundwater was affected by infiltration from the Houki Rivers.
• Nitrate nitrogen concentration was high.
• Contribution of recharge from paddy rice field was relatively small. Figure 15. Characteristics of the Nasunogahara alluvial fan groundwaters in each of the four groups.
Medians of all variables plotted on the radar charts were transformed to range from 0 to 1. Figure 12 illustrates the differences in the spatial distribution of the four groundwater groups classified by the SOM between the IP and NP. We focus here on the central and eastern parts of the fan because we could not collect groundwater samples in the western part of the fan in the NP. One large seasonal change was that the distribution area of group 2 shrank in the NP, whereas that of group 3 expanded in the downstream direction. EC, δ 18 O, δ 2 H, and d-excess values and 222 Rn concentrations of groundwater at the observation sites classified into group 2 in the IP differed significantly between the IP and NP (p < 0.05) (Figure 16). Although the EC apparently decreased in the NP in the southeastern part of the fan (where only group 2 samples were distributed in the IP) (Figure 7), the EC in group 2 groundwaters did not differ significantly between the IP and NP (p = 0.12). The major ion concentrations in group 2 groundwaters also did not differ significantly between the two periods (p > 0.10). The median δ 18 O and δ 2 H values of group 2 groundwaters in the IP were −8.0% and −53% , respectively, and the median values at the same sites in the NP decreased to −8.4% and −55% , respectively, whereas the median d-excess value increased from 11.0% in the IP to 12.3% in the NP. These results indicated that, with respect to the relation between δ 18 O and δ 2 H, the group 2 groundwaters deviated more from the LMWL in the IP than in the NP. The high isotopic ratios of the group 2 groundwaters indicated that they were likely more greatly affected by recharge from paddy rice fields than groundwaters of the other groups. This inference is consistent with the higher isotopic ratios that were observed in the IP than in the NP. As a result of this seasonal difference, some sampling sites shifted from group 2 in the IP to group 3 in the NP, reflecting a reduced influence of recharge from paddy rice fields in the NP and a corresponding relative increase in the influence of infiltration of river waters at these sites. In addition, 222 Rn concentrations of groundwaters were lower in the IP than in the NP ( Figure 16). Generally, 222 Rn concentrations in groundwater decrease during the irrigation period because of downward flow of soil water pushed out by irrigation water [55]. Accordingly, the difference in 222 Rn concentrations between the IP and NP also indicated a change in the relative contributions of infiltration from paddy rice fields and rivers to groundwater recharge.
Water 2020, 12, x FOR PEER REVIEW 18 of 20 Figure 12 illustrates the differences in the spatial distribution of the four groundwater groups classified by the SOM between the IP and NP. We focus here on the central and eastern parts of the fan because we could not collect groundwater samples in the western part of the fan in the NP. One large seasonal change was that the distribution area of group 2 shrank in the NP, whereas that of group 3 expanded in the downstream direction. EC, δ 18 O, δ 2 H, and d-excess values and 222 Rn concentrations of groundwater at the observation sites classified into group 2 in the IP differed significantly between the IP and NP (p < 0.05) (Figure 16). Although the EC apparently decreased in the NP in the southeastern part of the fan (where only group 2 samples were distributed in the IP) (Figure 7), the EC in group 2 groundwaters did not differ significantly between the IP and NP (p = 0.12). The major ion concentrations in group 2 groundwaters also did not differ significantly between the two periods (p > 0.10). The median δ 18 O and δ 2 H values of group 2 groundwaters in the IP were −8.0‰ and −53‰, respectively, and the median values at the same sites in the NP decreased to −8.4‰ and −55‰, respectively, whereas the median d-excess value increased from 11.0‰ in the IP to 12.3‰ in the NP. These results indicated that, with respect to the relation between δ 18 O and δ 2 H, the group 2 groundwaters deviated more from the LMWL in the IP than in the NP. The high isotopic ratios of the group 2 groundwaters indicated that they were likely more greatly affected by recharge from paddy rice fields than groundwaters of the other groups. This inference is consistent with the higher isotopic ratios that were observed in the IP than in the NP. As a result of this seasonal difference, some sampling sites shifted from group 2 in the IP to group 3 in the NP, reflecting a reduced influence of recharge from paddy rice fields in the NP and a corresponding relative increase in the influence of infiltration of river waters at these sites. In addition, 222 Rn concentrations of groundwaters were lower in the IP than in the NP ( Figure 16). Generally, 222 Rn concentrations in groundwater decrease during the irrigation period because of downward flow of soil water pushed out by irrigation water [55]. Accordingly, the difference in 222 Rn concentrations between the IP and NP also indicated a change in the relative contributions of infiltration from paddy rice fields and rivers to groundwater recharge.

Conclusions
The groundwater in the Nasunogahara alluvial fan, which is mostly covered by paddy rice fields, was investigated by using environmental isotopes and hydrochemical investigations and the application of an SOM. The obtained isotopic and hydrochemical data indicated that the groundwater in the fan was affected by three different recharge sources: precipitation, river waters, and paddy rice field irrigation water. Furthermore, the SOM clearly classified the groundwater in the fan into four groups (groups 1 to 4), reflecting different recharge sources. The characteristics of the groundwater of each group can be summarized as follows: 1. Group 1 groundwater, distributed around the Houki River, which flows along the western edge of the fan, had relatively low isotopic ratios, but high EC values and high Na + and Cl − concentrations. Group 1 groundwater was thus inferred to have been greatly affected by infiltration from the Houki River, which had water of a different chemical composition compared to the other rivers.

Conclusions
The groundwater in the Nasunogahara alluvial fan, which is mostly covered by paddy rice fields, was investigated by using environmental isotopes and hydrochemical investigations and the application of an SOM. The obtained isotopic and hydrochemical data indicated that the groundwater in the fan was affected by three different recharge sources: precipitation, river waters, and paddy rice field irrigation water. Furthermore, the SOM clearly classified the groundwater in the fan into four groups (groups 1 to 4), reflecting different recharge sources. The characteristics of the groundwater of each group can be summarized as follows: 1.
Group 1 groundwater, distributed around the Houki River, which flows along the western edge of the fan, had relatively low isotopic ratios, but high EC values and high Na + and Cl − concentrations. Group 1 groundwater was thus inferred to have been greatly affected by infiltration from the Houki River, which had water of a different chemical composition compared to the other rivers.

2.
Group 2 groundwater, distributed in the central and lower part of the fan, had high isotopic ratios and was inferred to be recharged mainly by infiltration of paddy waters that had been affected by evaporative isotopic enrichment. However, the area of group 2 shrank during the non-irrigation period when the infiltration of paddy water did not occur.

3.
Group 3 groundwater, distributed around the Sabi and Kuma Rivers, which flows down the center of the fan, had a low EC and low isotopic ratios, indicating a greater influence of infiltration of water from the two rivers compared to other recharge sources. 4.
Group 4 groundwater, distributed on the upstream side of group 2 groundwater, had lower isotopic ratios than group 2 groundwater. In the upper part of the fan, few paddy rice fields but many livestock farms are distributed throughout this region; thus, recharge from paddy rice fields was relatively small. Furthermore, NO 3 -N concentrations in group 4 groundwaters were higher than those in the other groups.
These results show that multiple tracers could be effectively used to evaluate the influence of different groundwater recharge sources, including paddy water in rice cultivation areas, and the application of an SOM can provide understandable and readily visualized results to assist in the interpretation of the characteristics of groundwater in the fan. The findings of this study can therefore contribute to proper planning for regional groundwater use and the preservation of spring waters in the fan.
Author Contributions: Design of the study framework, T.T. and S.I.; director of the field survey and analysis T.T.; water sample collection, all authors; writing, T.T.; reviewing and editing, K.S., S.I., and S.Y. All authors have read and agreed to the published version of the manuscript.
Funding: This study was partly supported by JSPS KAKENHI (grant number JP19K06301).