Suitability Evaluation of River Bank Filtration along the Second Songhua River, China

The Second Songhua River is the biggest river system in Jilin Province, China. In recent years, the rapid economic development in this area has increased the prominence of water resources and water-related environmental problems; these include surface water pollution and the overexploitation of groundwater resources. Bank infiltration on the floodplains of the Second Songhua River is an important process of groundwater-surface water exchange under exploitation conditions. Understanding this process can help in the development of water resource management plans and strategies for the region. In this research, a multi-criteria evaluation index system was developed with which to evaluate the suitability of bank filtration along the Second Songhua River. The system was comprised of main suitability indexes for water quantity, water quality, the interaction intensity between surface water and groundwater, and the exploitation condition of groundwater resources. The index system was integrated into GIS (Geographic Information System) to complete the evaluation of the various indicators. According to the weighted sum of each index, the suitability of river bank filtration (RBF) in the study area was divided into five grades. Although the evaluation index system and evaluation method are applicable only to the Second Songhua River basin, the underlying principle and techniques it embodies can be applied elsewhere. For future generalization of the evaluation index system, the specific evaluation index and its scoring criteria should be modified appropriately based on local conditions.


Introduction
With its rich and stable water quantity, good water quality, and easy access and management, river bank filtration (RBF) is widely adopted in many countries and has become an important means of groundwater exploitation [1,2].In Holland, Germany, Slovakia and Hungary, the proportion of total drinking water supplied by RBF water supply has reached 5%, 16%, 50% and 45%, respectively [3].The northern regions of China have established a number of RBF waterworks since the early 1950s; these include the Songhua River waterworks in Heilongjiang Province, the Hunhe River waterworks in Liaoning Province, the Jiuwutan and Beijiao waterworks in Henan Province, and the Hanzhong waterworks in Shanxi Province [4].To date, more than 300 large riverside RBF waterworks [5][6][7][8] have been established in China.
The combined use of surface water and groundwater resources can be realized through RBF, which can stimulate and increase the recharge of river water to the aquifer [9].However, the unreasonable exploitation of groundwater resources will lead to a continually descending groundwater level, the disruption of the hydraulic connection between groundwater and surface water, saltwater intrusion and land subsidence [10].From the perspective of the surface water infiltration process and runoff infiltration process and runoff variation along the river, Qin (1995) proposed a vertical seepage method by combining Manning's Formula and Darcy's law to calculate the water recharge capacity of river water [11].The riverbed clogging issue was thought to be an important factor that affects the performance of the riverbank filtration system [12].Han (1996), according to the landform unit of different RBF waterworks, divided water sources into four types, including intermountain valley, intermountain basin, alluvial-proluvial fan and coastal plain [4]; the researchers concluded that the infiltration capacity calculation of surface water is the key link during the water resources evaluation of RBF waterworks.Grischek and Ray (2009) presented the surface-groundwater interaction issues at various geomorphologic settings, from the headwaters of a river to its confluence with oceans/lakes [13].
The quality of surface water can be improved by the physical, chemical and biological processes that take place in RBF process.As surface water infiltrates groundwater, the total organic carbon, turbidity, heavy metals, inorganic substances, viruses, parasites and some microorganisms in surface water can be effectively reduced.In RBF systems, 60% of the total organic carbon and 90% of the turbidity can be removed [3,14].Likewise, 90% of chromium and arsenic, and 50% of cadmium, zinc, lead, copper and nickel can be removed [15].The removal of ammonia nitrogen can reach 95% [16], and the removal of bacteria, parasites, and viruses can reach 100% [14].
In the past, a large number of studies have been conducted mainly for specific RBF waterworks.Different methods were developed to analyze the feasibility of RBF systems [17][18][19].For the feasibility research of RBF sites in a regional setting, Sandhu et al. (2011) provided the analysis results in India [14].So far, there are only few publications on the suitability evaluation of RBF on the scale of a river basin, because a generalization of site selection criteria is very difficult.
In this research, the downstream river basin of the Second Songhua River was set as the study area and the potential suitable area for RBF was evaluated.Based on the analysis of regional natural geography and hydrogeological conditions, an index system was established to evaluate the suitability of RBF along the main stem of the river in the study area.

Study Area
The Second Songhua River is the largest river in Jilin Province, and originates from Tianchi Lake in the Changbai Mountains.The catchment area is 7.3 × 10 4 km 2 , including the tributaries of the Yitong River, Yinma River, Mangniu River and others, and accounts for 38.2% of the total province area.The study area focuses on the plain area of the Second Songhua River catchment with a total area of 2.46 × 10 4 km 2 (Figure 1).The major cities in the study area include Changchun City, the provincial capital of Jilin Province, and Jilin City, Dehui City, Jiutai City, Nong'an County, and Fuyu County.The total length of the Second Songhua River is approximately 330 km in the study area.The Yinma River is the largest tributary of the Second Songhua River and stems from Hulanling Mountain in Panshi County, Jilin Province, finally entering the Second Songhua River in Nong'an County.The Yitong River is the tributary of the Yinma River and the second tributary of the Second Songhua River.The Yitong River stems from the north of Qingding Mountain in Yitong County, Jilin Province, and finally flows into the Second Songhua River.The average annual precipitation in the study area is 400-700 mm, decreasing from southeast to northwest (Figure 2).Both sides of the Second Songhua River are dominated by the valley plain geomorphic unit.The sedimentary plain is divided into multiple blocks between rivers with floodplain and terrace, and the surface elevation decreases gradually from the southeast to the northwest (Figure 3).Quaternary sediment with good permeability is widely distributed in the Second Songhua River Basin with a thickness of 10-30 m.The type of groundwater is mainly pore groundwater in the sandy gravel layer and the groundwater table is shallow, generally 1-3 m below the ground surface.The single well groundwater production rate is generally 1000-3000 m 3 /d.Surface water and groundwater are well connected, and groundwater recharges surface water in most areas along the river (Figure 4).The major cities in the study area include Changchun City, the provincial capital of Jilin Province, and Jilin City, Dehui City, Jiutai City, Nong'an County, and Fuyu County.The total length of the Second Songhua River is approximately 330 km in the study area.The Yinma River is the largest tributary of the Second Songhua River and stems from Hulanling Mountain in Panshi County, Jilin Province, finally entering the Second Songhua River in Nong'an County.The Yitong River is the tributary of the Yinma River and the second tributary of the Second Songhua River.The Yitong River stems from the north of Qingding Mountain in Yitong County, Jilin Province, and finally flows into the Second Songhua River.The average annual precipitation in the study area is 400-700 mm, decreasing from southeast to northwest (Figure 2).Both sides of the Second Songhua River are dominated by the valley plain geomorphic unit.The sedimentary plain is divided into multiple blocks between rivers with floodplain and terrace, and the surface elevation decreases gradually from the southeast to the northwest (Figure 3).Quaternary sediment with good permeability is widely distributed in the Second Songhua River Basin with a thickness of 10-30 m.The type of groundwater is mainly pore groundwater in the sandy gravel layer and the groundwater table is shallow, generally 1-3 m below the ground surface.The single well groundwater production rate is generally 1000-3000 m 3 /d.Surface water and groundwater are well connected, and groundwater recharges surface water in most areas along the river (Figure 4).The major cities in the study area include Changchun City, the provincial capital of Jilin Province, and Jilin City, Dehui City, Jiutai City, Nong'an County, and Fuyu County.The total length of the Second Songhua River is approximately 330 km in the study area.The Yinma River is the largest tributary of the Second Songhua River and stems from Hulanling Mountain in Panshi County, Jilin Province, finally entering the Second Songhua River in Nong'an County.The Yitong River is the tributary of the Yinma River and the second tributary of the Second Songhua River.The Yitong River stems from the north of Qingding Mountain in Yitong County, Jilin Province, and finally flows into the Second Songhua River.The average annual precipitation in the study area is 400-700 mm, decreasing from southeast to northwest (Figure 2).Both sides of the Second Songhua River are dominated by the valley plain geomorphic unit.The sedimentary plain is divided into multiple blocks between rivers with floodplain and terrace, and the surface elevation decreases gradually from the southeast to the northwest (Figure 3).Quaternary sediment with good permeability is widely distributed in the Second Songhua River Basin with a thickness of 10-30 m.The type of groundwater is mainly pore groundwater in the sandy gravel layer and the groundwater table is shallow, generally 1-3 m below the ground surface.The single well groundwater production rate is generally 1000-3000 m 3 /d.Surface water and groundwater are well connected, and groundwater recharges surface water in most areas along the river (Figure 4).Municipal and industrial wastewater seriously affects the water quality of Second Songhua River.Annual emissions of wastewater in this basin are about 645 million tons, of which 44% arises in Jilin Province.As total nitrogen, total phosphorus, ammonia nitrogen exceed the standard, water quality deterioration in the Yinma River is very serious; the river length in which water quality is inferior to Grade V water quality accounts for 42.9% of the total length [20].Changchun City is located in the middle reach of the Yitong River, and municipal and industrial sewage has had serious impacts on water quality of the river [21].The Second Songhua River in Jilin Province has abundant water resources, but now cannot meet the needs of economic development, especially the water requirement of large cities [22].

Division of the Core Study Area of RBF
In order to focus on the optimum area, the core study area of RBF was firstly separated considering the following two criteria: 1. Natural geographical conditions.The plain area along the river was mainly considered as it is the main place of production and human activity, which has a strong demand for water resources.2. Hydraulic connection extent between river water and groundwater.The feasibility of RBF is getting worse and may even be impossible with increasing distance from river.In order to get more complete background information on potential RBF sites, a large study area was recommended.According to the runoff quantity, the study range of RBF was divided along the river (Table 1).Where there is a surface water divide in the study range, this natural feature is used as the boundary of the study area.
Table1.Core study range (distance from river) of river bank filtration (RBF) along river banks.The scope of the study area of RBF and the feasible region of RBF was determined by overlay analysis of areas that satisfied the above two criteria, such that any areas that did not satisfy both criteria were excluded.Municipal and industrial wastewater seriously affects the water quality of Second Songhua River.Annual emissions of wastewater in this basin are about 645 million tons, of which 44% arises in Jilin Province.As total nitrogen, total phosphorus, ammonia nitrogen exceed the standard, water quality deterioration in the Yinma River is very serious; the river length in which water quality is inferior to Grade V water quality accounts for 42.9% of the total length [20].Changchun City is located in the middle reach of the Yitong River, and municipal and industrial sewage has had serious impacts on water quality of the river [21].The Second Songhua River in Jilin Province has abundant water resources, but now cannot meet the needs of economic development, especially the water requirement of large cities [22].

Division of the Core Study Area of RBF
In order to focus on the optimum area, the core study area of RBF was firstly separated considering the following two criteria: Natural geographical conditions.The plain area along the river was mainly considered as it is the main place of production and human activity, which has a strong demand for water resources.

2.
Hydraulic connection extent between river water and groundwater.The feasibility of RBF is getting worse and may even be impossible with increasing distance from river.In order to get more complete background information on potential RBF sites, a large study area was recommended.According to the runoff quantity, the study range of RBF was divided along the river (Table 1).Where there is a surface water divide in the study range, this natural feature is used as the boundary of the study area.The scope of the study area of RBF and the feasible region of RBF was determined by overlay analysis of areas that satisfied the above two criteria, such that any areas that did not satisfy both criteria were excluded.

Evaluation Index System
The aim of this research is to find potential suitable areas for future water development plans and further detailed investigation of RBF systems.The specific goal of RBF is mainly to provide sufficient water resources with good water quality.A multi-criteria RBF evaluation index system was established based on water quantity, water quality, the development and utilization conditions of groundwater resources and the interaction intensity between surface water and groundwater.The index weight and detailed scoring criterion were all based on specialist marking methods, which is more effective for a specific site but may not be applicable to other areas (Table 2).The main reasons for index weight are explained as follows.
(1) The sum of evaluation index weight is equal to 1 by value assignment.The index weight of water quantity, water quality and interaction intensity between surface water and groundwater each accounts for 0.3 of total weight individually, and this allotment reflected the aim and decisive factors of RBF suitability evaluation along the Second Songhua River.The index weight of development and utilization conditions of groundwater resources only account for 0.1 of total weight, because the factor will influence groundwater development cost and help to confirm the priority, but it is not a decisive factor for RBF suitability evaluation.(2) For the index of water quantity, the groundwater index accounts for 0.20 and surface water index accounts for 0.1, because the water source of RBF waterworks comes from riverside groundwater and surface water via bank filtration.Hydraulic conductivity and aquifer thickness influence the infiltration rate and capacity of surface water to aquifer.(3) For water quality, the index weight of groundwater and surface water are equal to 0.15.The groundwater quality only indicates the current status and treatment effect of RBF.In order to avoid any possible groundwater contamination from surface water, the treatment effect of RBF should not be addressed too much and the surface water quality should not be neglected.(4) The interaction intensity between surface water and groundwater controls the water exchange efficiency.Only those places within the influence zone of surface water, the aquifer prone to receiving sufficient water quantity through RBF, the permeability of riverbed layer controls the real water exchange rate and the groundwater hydraulic gradient partly reflects the real condition of aquifer property.(4) The index weight of groundwater depth accounts for 0.1, because it will only influence the priority order of potential RBF areas.

Water Quantity Hydraulic Conductivity (K)
Aquifer hydraulic conductivity (K) reflects lithology and permeability of an aquifer.The scoring criterion of hydraulic conductivity is shown in Table 3. Aquifer Thickness (M) A suitable RBF site requires a certain scale of aquifer for adequate water production.The aquifer thickness scoring criterion is shown in Table 4. Runoff in Cross-Section (Q) The magnitude of runoff in a given cross-section of the landscape directly reflects the richness of the surface water quantity.The larger the annual average runoff, the greater the recharge potential of surface water to groundwater, and the greater the potential suitability of the area for RBF.The runoff in the cross-section scoring criterion is shown in Table 5.

Water Quality
Water quality status directly reflects whether surface water or groundwater is suitable for drinking, and indicates the degree of difficulty and cost of water treatment that will produce drinking water.In order to guarantee the water quality of potential waterworks, the water quality improvement effect of RBF was thought of as an additional assurance.According to China's Environmental quality standards for surface water (GB3838-2002) [23] and Quality standard for groundwater (GB/T 14848-93) [24] (see Appendix A), water of quality grades I, II and III can be used as drinking water directly.Water of quality grade IV must be treated before it is supplied for drinking.Because higher index values of the quality grade indicate lower quality water, the assigned index values decrease as water quality grades change from I to IV. Water of quality grade V is a negative factor for the suitability of RBF, and its index value is assigned a negative number.This can make the negative water quality factor play a decisive role in the process of suitability evaluation of RBF.The groundwater quality (G) and surface water quality (S) suitability scoring criterion is shown in Table 6.

Interaction Intensity between Surface Water and Groundwater
A large proportion of RBF waterworks primarily capture surface water; the interaction intensity between surface water and groundwater plays an essential role.

Hydraulic Gradient (I) under the Current Condition
The hydraulic gradient directly reflects recharge-discharge relationship between surface water and groundwater, and it also indicates recharge-discharge conditions of the interaction zone.A positive hydraulic gradient was defined to represent scenarios in which river water recharges groundwater; similarly, a negative hydraulic gradient indicates that river water is recharged by groundwater.A positive value of hydraulic gradient is beneficial to RBF in most situations; however, a very large value for the hydraulic gradient might mean that riverside groundwater is being exploited intensively or that the interaction between surface water and groundwater is poor, because the very low permeability of aquifer or riverbed will obviously enlarge hydraulic gradient and decrease water exchange intensity.A negative hydraulic gradient indicates that the surface water cannot be induced into the ambient aquifer at present.Nevertheless, a negative gradient might be reversed and to be beneficial for RBF under groundwater exploitation condition.If the present negative hydraulic gradient is excessively small, it may indicate that the permeability of aquifer or riverbed is very poor, and the large groundwater recharge potential from river water should not be expected under any exploitation conditions.The hydraulic gradient suitability scoring criterion is shown in Table 7.The possible influence zone width by surface water indicates the range of the hydraulic connection between surface water and groundwater.The greater the hydraulic connection range between surface water and groundwater, the superior the RBF site.The possible influence zone width (L) was assumed to be equal to the ratio of the aquifer thickness (M) to the hydraulic gradient (I) in a specified cross section.The suitability scoring criterion for the possible influence zone width by surface water is shown in Table 8.Permeability of Riverbed Layer (R) Riverbed permeability indicates the exchange capacity between surface water and groundwater.Riverbed permeability is calculated using the hydrogeological cross section of the riverbed obtained from a field investigation as standard.Then, the lithology under the riverbed is analyzed to obtain its permeability, either by measurement or practical experience.The riverbed permeability suitability scoring criterion is shown in Table 9.The Exploitation Condition of the Groundwater Resource When considering the development and utilization of groundwater resources, groundwater depth is a primary concern because it indicates the cost of building RBF wells.If the groundwater depth is too far below the ground surface, it may be in a highlands area or in a cone of a groundwater depression.In those places, the extraction of groundwater requires more energy.The groundwater depth scoring criterion is shown in Table 10.Because groundwater depth is strongly influenced by regional characteristics, this standard should be adapted to different hydrogeology conditions.

Suitability Index
According to the evaluation index system created above, the complex suitability index for a potential RBF site can be calculated by weighted summation using Equation (1).
A " X K ˆWK `XM ˆWM `XQ ˆWQ `XG ˆWG `XS ˆWS `XI ˆWI `XL ˆWL `XR ˆWR `XD ˆWD (1) In Equation (1), A is the suitability index of a potential RBF site, X is the score of individual indices as defined in Tables 3-10 and W is the weight of each corresponding index as defined in Table 2.The RBF suitability was then classified into five grades according to the suitability index value (Table 11).In areas for which RBF suitability grades are between I and III, water quality satisfies human drinking water requirements and does not need treatment, and water quantity is adequate.Furthermore, the exchange between river water and groundwater is intensive, the water quantity and water quality of RBF is guaranteed, and there are good development and utilization conditions of groundwater resources.In short, areas in grades I-III are suitable for RBF.In areas for which the RBF suitability grade is IV, either water quality does not satisfy drinking water requirements directly or the water recharge is insufficient.After water treatment and limited exploration, such areas may be suitable for limited-scale RBF waterworks.Areas in grade V are unsuitable for RBF works because of unacceptable water quality, high water treatment costs or insufficient water recharge to sustain extraction at the required rate and volume.

Spatial Analysis Based on GIS
GIS (Geographic Information System) has many powerful functions such as data editing, data management, data conversion, projection transformation, geographic analysis, metadata management, spatial analysis, and overlay analysis [25].Applying GIS to a regional suitability evaluation of potential RBF sites can improve the accuracy and efficiency of evaluation.Spatial analysis based on GIS was adopted to facilitate the evaluation of the various indicators.

The Core Study Area of RBF
According to the terrain classification standard, the altitude range of 0 m-200 m belongs to the plain area.The boundary of the right bank of the Second Songhua River is mainly determined by the surface water divide.The annual average runoff of the Yinma River is 9.44 ˆ10 8 m 3 (Dehui Hydrological Station) and it belongs to river grade V.The annual average runoff of the Yitong River is 3.7 ˆ10 8 m 3 (Nong'an Hydrological Station) and this river also belongs to river grade V. So, the part of Dehui City that is on the left bank of the Yinma River and Nong'an County on the left bank of the Yitong River should take a hydraulic connection between river water and groundwater as standard.Thus, using the Yitong River and Yinma River as the center, an area extending 10 km from both sides of the rivers was defined as the core study area near the rivers (Figure 5).

Water Quantity Suitability
The aquifer lithology of the Second Songhua River plain mainly is gravel and coarse sand, and the aquifer hydraulic conductivity is approximately 60 m/day.The aquifer lithology of the Yinma River and the Yitong River plain is mainly medium sand and sandy clay, and the aquifer hydraulic conductivity is only 12.5 m/day.The distribution of aquifer hydraulic conductivity in the core study area is shown in Figure 6.

Water Quantity Suitability
The aquifer lithology of the Second Songhua River plain mainly is gravel and coarse sand, and the aquifer hydraulic conductivity is approximately 60 m/day.The aquifer lithology of the Yinma River and the Yitong River plain is mainly medium sand and sandy clay, and the aquifer hydraulic conductivity is only 12.5 m/day.The distribution of aquifer hydraulic conductivity in the core study area is shown in Figure 6.
The aquifer thickness distribution in the core study area is shown in Figure 7.In the area traversed by the Second Songhua River downstream of Jilin city and in the reach from the Yinma River estuary to Fuyu City upstream of the Second Songhua River, the aquifer thickness is usually 10 m-30 m, and as much as 50 m in some sections.In the area traversed by the Yinma River, the aquifer thickness of the middle reach from Jiutai to Dehui is 10 m-30 m, whereas the aquifer thickness of other regions is only 5 m-10 m.In the area traversed by the Yitong River, the aquifer thickness of the Changchun City reach and of the Nong'an County reach is relatively thin, but in other locations the thickness is 10 m-30 m.
The aquifer lithology of the Second Songhua River plain mainly is gravel and coarse sand, and the aquifer hydraulic conductivity is approximately 60 m/day.The aquifer lithology of the Yinma River and the Yitong River plain is mainly medium sand and sandy clay, and the aquifer hydraulic conductivity is only 12.5 m/day.The distribution of aquifer hydraulic conductivity in the core study area is shown in Figure 6.The aquifer thickness distribution in the core study area is shown in Figure 7.In the area traversed by the Second Songhua River downstream of Jilin city and in the reach from the Yinma River estuary to Fuyu City upstream of the Second Songhua River, the aquifer thickness is usually 10 m-30 m, and as much as 50 m in some sections.In the area traversed by the Yinma River, the aquifer thickness of the middle reach from Jiutai to Dehui is 10 m-30 m, whereas the aquifer thickness of other regions is only 5 m-10 m.In the area traversed by the Yitong River, the aquifer thickness of the Changchun City reach and of the Nong'an County reach is relatively thin, but in other locations the thickness is 10 m-30 m.The runoff of the Second Songhua River is very large (approximately 150 × 10 8 m 3 /a).The runoff of the Yinma River and the Yitong River is 9 × 10 8 m 3 /a and 3 × 10 8 m 3 /a, respectively.The annual average runoff of rivers in the study area is shown in Figure 8.
According to the distribution of aquifer hydraulic conductivity, aquifer thickness and runoff quantity in the region, the water quantity suitability scores assessed by using Equation ( 1) is shown in Figure 9.The score distribution of all indexes is shown in Appendix B. The runoff of the Second Songhua River is very large (approximately 150 ˆ10 8 m 3 /a).The runoff of the Yinma River and the Yitong River is 9 ˆ10 8 m 3 /a and 3 ˆ10 8 m 3 /a, respectively.The annual average runoff of rivers in the study area is shown in Figure 8.
According to the distribution of aquifer hydraulic conductivity, aquifer thickness and runoff quantity in the region, the water quantity suitability scores assessed by using Equation ( 1) is shown in Figure 9.The score distribution of all indexes is shown in Appendix B.
As seen in Figure 9, the suitability of water quantity in most regions of the core study area is better than grade III.Only in the middle reach of the Second Songhua River and in the Songhua Lake upstream reach in Jilin City is the water quantity suitability relatively poor (grade IV).Overall, the suitability of water quantity in the study area is good and can basically ensure an adequate water supply for RBF wells.
of the Yinma River and the Yitong River is 9 × 10 m /a and 3 × 10 m /a, respectively.The annual average runoff of rivers in the study area is shown in Figure 8.
According to the distribution of aquifer hydraulic conductivity, aquifer thickness and runoff quantity in the region, the water quantity suitability scores assessed by using Equation ( 1) is shown in Figure 9.The score distribution of all indexes is shown in Appendix B.  As seen in Figure 9, the suitability of water quantity in most regions of the core study area is better than grade III.Only in the middle reach of the Second Songhua River and in the Songhua Lake upstream reach in Jilin City is the water quantity suitability relatively poor (grade IV).Overall, the suitability of water quantity in the study area is good and can basically ensure an adequate water supply for RBF wells.

Water Quality Suitability
The groundwater quality assessment was based on the Quality standard for groundwater (GB/T 14848-93) [24] using the monitoring results for chloride, fluoride, sulfate, ammonia, total hardness, total iron, nitrite and nitrate as the evaluation index.The single-factor evaluation, which contrasts each water quality index to the national standard, was adopted to assess the water quality grade.The status of groundwater quality in the study area is shown in Figure 10.

Water Quality Suitability
The groundwater quality assessment was based on the Quality standard for groundwater (GB/T 14848-93) [24] using the monitoring results for chloride, fluoride, sulfate, ammonia, total hardness, total iron, nitrite and nitrate as the evaluation index.The single-factor evaluation, which contrasts each water quality index to the national standard, was adopted to assess the water quality grade.The status of groundwater quality in the study area is shown in Figure 10.
The surface water quality assessment was based on the Environmental quality standards for surface water (GB 3838-2002) [23] as a standard; the Water Quality Index was adopted for the assessment.The status of surface water quality in the study area is shown in Figure 11.Using Equation ( 1), the scoring distributions of groundwater and surface water quality in the study area were integrated to produce the water quality suitability distribution shown in Figure 12.
The groundwater quality assessment was based on the Quality standard for groundwater (GB/T 14848-93) [24] using the monitoring results for chloride, fluoride, sulfate, ammonia, total hardness, total iron, nitrite and nitrate as the evaluation index.The single-factor which contrasts each water quality index to the national standard, was adopted to assess the water quality grade.The status of groundwater quality in the study area is shown in Figure 10.The surface water quality assessment was based on the Environmental quality standards for surface water (GB 3838-2002) [23] as a standard; the Water Quality Index was adopted for the assessment.The status of surface water quality in the study area is shown in Figure 11.Using Equation ( 1), the scoring distributions of groundwater and surface water quality in the study area were integrated to produce the water quality suitability distribution shown in Figure 12.As seen in Figure 12, water quality is distributed across the core study in discrete areas.In the Fuyu City reach and in the downstream reach of the Yinma River estuary of the Second Songhua River, the water quality suitability is grade III.Other parts of the Second Songhua River have grade I water quality suitability.In the Jiutai City reach of the Yinma River, the water quality suitability is grade I, but other parts of the Yinma River are grade III.In the Yitong River basin, the water quality suitability is mainly grade III and grade V, distributed sporadically, except for the Changchun City reach where the water quality suitability is grade I.In general, the water quality suitability is of an acceptable standard in most of the study area.As seen in Figure 12, water quality is distributed across the core study in discrete areas.In the Fuyu City reach and in the downstream reach of the Yinma River estuary of the Second Songhua River, the water quality suitability is grade III.Other parts of the Second Songhua River have grade I water quality suitability.In the Jiutai City reach of the Yinma River, the water quality suitability is grade I, but other parts of the Yinma River are grade III.In the Yitong River basin, the water quality suitability is mainly grade III and grade V, distributed sporadically, except for the Changchun City reach where the water quality suitability is grade I.In general, the water quality suitability is of an acceptable standard in most of the study area.As seen in Figure 12, water quality is distributed across the core study in discrete areas.In the Fuyu City reach and in the downstream reach of the Yinma River estuary of the Second Songhua River, the water quality suitability is grade III.Other parts of the Second Songhua River have grade I water quality suitability.In the Jiutai City reach of the Yinma River, the water quality suitability is grade I, Water 2016, 8, 176 13 of 24 but other parts of the Yinma River are grade III.In the Yitong River basin, the water quality suitability is mainly grade III and grade V, distributed sporadically, except for the Changchun City reach where water quality suitability is grade I.In general, the water quality suitability is of an acceptable standard in most of the study area.

Interaction Intensity between Surface Water and Groundwater
The interaction intensity between surface water and groundwater was assessed using the river level (H) and the groundwater level (h) at 1 km (L) distant from the river bank to calculate the hydraulic gradient (I = [(H -h)/L]) at intervals along the river of not more than 10 km.The cross sections at which hydraulic gradients were calculated are shown in Figure 13.According to these calculations, the hydraulic gradient scoring index was assigned and distributed as shown in Figure 14.Then, using the calculated cross sections of the hydraulic gradient, the possible influence zone widths of surface water were determined and assigned a score (from Table 8), as shown in Figure 15.Then, using the calculated cross sections of the hydraulic gradient, the possible influence zone widths of surface water were determined and assigned a score (from Table 8), as shown in Figure 15.Then, using the calculated cross sections of the hydraulic gradient, the possible influence zone widths of surface water were determined and assigned a score (from Table 8), as shown in Figure 15.
The locations of hydrogeology profiles based on drilling results are shown in Figure 16.For each cross section, the lithology under the riverbed was analyzed and, using the Handbook of Hydrogeology [26] to obtain a value of the lithology permeability, the lithology permeability was multiplied by 0.1 to determine riverbed permeability.For example, Figure 17 shows the No. 3 hydrogeology profile from Xinlitun to Shaokou; the lithology layer incised by the river is comprised of sand gravel and its hydraulic conductivity (based on experience) is 50 m/day-100 m/day.In this assessment, the hydraulic conductivity of the sand and gravel (75 m/day) was adopted; thus, the riverbed permeability at this cross section was 7.5 m/day.The distribution of riverbed permeability is shown in Figure 18.According to the riverbed permeability suitability evaluation criterion (Table 9), the riverbed permeability was assigned a score.Then, using the calculated cross sections of the hydraulic gradient, the possible influence zone widths of surface water were determined and assigned a score (from Table 8), as shown in Figure 15.The locations of hydrogeology profiles based on drilling results are shown in Figure 16.For each cross section, the lithology under the riverbed was analyzed and, using the Handbook of Hydrogeology [26] to obtain a value of the lithology permeability, the lithology permeability was multiplied by 0.1 to determine riverbed permeability.For example, Figure 17 shows the No. 3 hydrogeology profile from Xinlitun to Shaokou; the lithology layer incised by the river is comprised of sand gravel and its hydraulic conductivity (based on experience) is 50 m/day-100 m/day.In this assessment, the hydraulic conductivity of the sand and gravel (75 m/day) was adopted; thus, the riverbed permeability at this cross section was 7.5 m/day.The distribution of riverbed permeability is shown in Figure 18.According to the riverbed permeability suitability evaluation criterion (Table 9), the riverbed permeability was assigned a score.Using Equation (1), the distribution of scores for hydraulic gradient, possible interaction zone widths of surface water, and riverbed permeability in the core study area were combined to determine the suitability of interaction intensity between surface water and groundwater, as shown in Figure 19.As shown in Figure 19, the suitability of the interaction intensity between surface water and groundwater is grade IV and grade V at the Fuyu part and the Jilin part, respectively, of the Second Songhua River; in other parts of the river the suitability of interaction is superior to grade III.The suitability of the interaction intensity between surface water and groundwater at most parts of the Yitong River and Yinma River is inferior to grade III, except at the Changchun part of the Yitong River and the middle part between Jiutai and Dehui of the Yinma River, where the suitability of the interaction intensity between surface water and groundwater is grade III.Using Equation (1), the distribution of scores for hydraulic gradient, possible interaction zone widths of surface water, and riverbed permeability in the core study area were combined to determine the suitability of interaction intensity between surface water and groundwater, as shown in Figure 19.Using Equation (1), the distribution of scores for hydraulic gradient, possible interaction zone widths of surface water, and riverbed permeability in the core study area were combined to determine the suitability of interaction intensity between surface water and groundwater, as shown in Figure 19.As shown in Figure 19, the suitability of the interaction intensity between surface water and groundwater is grade IV and grade V at the Fuyu part and the Jilin part, respectively, of the Second Songhua River; in other parts of the river the suitability of interaction is superior to grade III.The suitability of the interaction intensity between surface water and groundwater at most parts of the Yitong River and Yinma River is inferior to grade III, except at the Changchun part of the Yitong River and the middle part between Jiutai and Dehui of the Yinma River, where the suitability of the interaction intensity between surface water and groundwater is grade III.As shown in Figure 19, the suitability of the interaction intensity between surface water and groundwater is grade IV and grade V at the Fuyu part and the Jilin part, respectively, of the Second Songhua River; in other parts of the river the suitability of interaction is superior to grade III.The suitability of the interaction intensity between surface water and groundwater at most parts of the Yitong River and Yinma River is inferior to grade III, except at the Changchun part of the Yitong River and the middle part between Jiutai and Dehui of the Yinma River, where the suitability of the interaction intensity between surface water and groundwater is grade III.

The Exploitation Condition of Groundwater Resources
The variation of groundwater depth is shown as Figure 20, respectively.The groundwater depth in the core study area is generally less than 15 m, which meets the criteria for groundwater exploitation and is beneficial to the construction of RBF works.

The Comprehensive Suitability Evaluation of RBF Sites
The various index scores were integrated using Equation (1) to determine the overall suitability of locations in the core study area for RBF works (Figure 21).Five typical points were selected to show the characteristics of each grade of suitability (Table 12).

The Comprehensive Suitability Evaluation of RBF Sites
The various index scores were integrated using Equation (1) to determine the overall suitability of locations in the core study area for RBF works (Figure 21).Five typical points were selected to show the characteristics of each grade of suitability (Table 12).

The Comprehensive Suitability Evaluation of RBF Sites
The various index scores were integrated using Equation (1) to determine the overall suitability of locations in the core study area for RBF works (Figure 21).Five typical points were selected to show the characteristics of each grade of suitability (Table 12).* The meaning of indexes is explained in Table 2.
Based on the multi-criteria analysis, the suitability of many locations in the study area was classified as grade V.These locations were the Jilin City reach, the Yinma River estuary and the Fuyu City reach of the Second Songhua River, as well as the whole Yitong River and the Dehui City reach of the Yinma River.These locations are unsuitable as potential RBF sites.The remainder of the study area was classified with a better than grade III suitability, and should be suitable for locating RBF works.

Conclusions
Many important factors affect the site selection for RBF works such as groundwater and surface water quantity, current water quality situation, hydraulic interaction degree and exchange relationship between groundwater and surface water.In this research, a multi-criteria index system was developed for the regional assessment of RBF site suitability in the Second Songhua River and Yinma River.The system was based on a detailed analysis of physical geography and geological and hydrogeological conditions (which were considered to be the main influential factors of RBF), as well as on the development and utilization of water resources and water demand.Scoring criteria based on specialist marking methods were used to determine weighting coefficients and weighted scores.The evaluation method was integrated into the spatial analysis features of GIS to determine the distribution of suitable areas for RBF works.By identifying the suitability grade of areas along the Second Songhua River, the regional evaluation system highlights locations that should be targeted for RBF works.The suitability evaluation system developed in this research thus provides a scientific basis for making decisions for regional industry distribution and relevant groundwater development plans.
However, this research still has some limitations that must be mentioned here and will be studied further.(1) The evaluation index system to assess RBF site suitability that was established in this research is suitable for use only in the Second Songhua River catchment.Unfortunately, this catchment is not universally representative of other catchments because of its very specific hydrogeological conditions.(2) The index system is not very comprehensive; some factors such as riverbed thickness and river width were missing because of the limited field data in the study area.(3)Some index value such as riverbed permeability was inferred by limited information, and this would bring some uncertainty to the final results.(4) The specialist marking method was adopted to determine the preference weights; the applicability was limited within specific areas.The AHP (Analytic Hierarchy Process) method and Fuzzy theory would be better to adopt in a general index system [27,28].
(5) For the specialist marking method, it is important to perform a sensitivity analysis on the preference weights and provide measures for assessing the sensitivity due to changes in these weights [28,29].Nevertheless, the primary evaluation principle and the index system are reasonable, and can be used as a scientific reference for any other regional site selection or feasibility research of RBF.

Figure 2 .
Figure 2. Rainfall isoline in the study area.

Figure 3 .
Figure 3. Surface topography of the study area.

Figure 2 .
Figure 2. Rainfall isoline in the study area.

Figure 2 .
Figure 2. Rainfall isoline in the study area.

Figure 3 .
Figure 3. Surface topography of the study area.

Figure 3 .
Figure 3. Surface topography of the study area.

Figure 4 .
Figure 4. Contour of the water table.

Figure 4 .
Figure 4. Contour of the water table.

Table 5 .
Runoff in Cross Section (Q) Scoring Criterion.

Figure 5 .
Figure 5.The core study area for evaluating RBF site suitability.

Figure 5 .
Figure 5.The core study area for evaluating RBF site suitability.

Figure 9 .
Figure 9. Distribution of water quantity suitability.

Figure 12 .
Figure 12.Distribution of water quality suitability in the core study area.

Figure 12 .
Figure 12.Distribution of water quality suitability in the core study area.

Figure 12 .
Figure 12.Distribution of water quality suitability in the core study area.

Figure 15 .
Figure 15.Distribution of the possible influence zone width of surface water.

Figure 15 .
Figure 15.Distribution of the possible influence zone width of surface water.

Figure 19 .
Figure 19.Suitability of interaction intensity between surface water and groundwater.

Figure 19 .
Figure 19.Suitability of interaction intensity between surface water and groundwater.

Figure 19 .
Figure 19.Suitability of interaction intensity between surface water and groundwater.

Figure 21 .
Figure 21.Overall suitability of locations within the core study area for RBF works.

Figure 21 .
Figure 21.Overall suitability of locations within the core study area for RBF works.

Figure 21 .
Figure 21.Overall suitability of locations within the core study area for RBF works.

Figure A4 .
Figure A4.Distribution of groundwater quality scores.

Figure A5 .
Figure A5.Distribution of surface water quality scores.

Figure A5 .
Figure A5.Distribution of surface water quality scores.

Figure A5 .
Figure A5.Distribution of surface water quality scores.

Figure A9 .
Figure A9.Distribution of the groundwater depth scores influence zone width of surface water.

Figure A7 . 24 Figure A7 .
Figure A7.Distribution of scores for the possible.

Figure A9 .
Figure A9.Distribution of the groundwater depth scores influence zone width of surface water.

Figure A7 .
Figure A7.Distribution of scores for the possible.

Figure A9 .
Figure A9.Distribution of the groundwater depth scores influence zone width of surface water.

Figure A9 .
Figure A9.Distribution of the groundwater depth scores influence zone width of surface water.

Table 1 .
Core study range (distance from river) of river bank filtration (RBF) along river banks.

Table 2 .
RBF suitability evaluation index system.

Table 6 .
Groundwater quality (G) and surface water quality (S) scoring criterion.
Possible Influence Zone Width of Surface Water under the Condition of Groundwater Exploitation (L)

Table 8 .
Scoring criterion for the possible influence zone width of surface water under the condition of groundwater exploitation (L).

Table 12 .
All single index values and scores for five typical points.

Table 12 .
All single index values and scores for five typical points.

Table 12 .
All single index values and scores for five typical points.