Variations of Groundwater Quality in the Multi-Layered Aquifer System near the Luanhe River, China

Nitrate pollution is an environmental problem in the North China Plain. This paper investigates the variation of groundwater levels and nitrate concentrations in an alluvial fan of the Luanhe river, northeast of the North China Plain. Three transects perpendicular to the riverbank were selected to investigate the exchange between river water and groundwater, and nitrate concentration with its isotopic composition (δN-NO3 and δO-NO3). The results showed that the groundwater level decreased slightly during the dry season, and increased regularly during the period of river stage rise. The groundwater is recharged by the river over 10 months each year. The nitrate concentration in the groundwater and river water varied with seasons. The nitrate concentration of groundwater in wells near the river is affected by the river water, which varied in basically the same way as the river. The nitrate concentrations in the zone of groundwater level depression cone were lower than those in the wells near the river, due to the long-term pumping of groundwater. However, the nitrate concentrations of river water have little influence on those of groundwater in wells far from the river. The values of δN-NO3 and the relationship between the two isotopes (δN-NO3 and δO-NO3) suggested that NO3-N was mainly attributable to sewage, livestock manure and natural soil organic matter. Due to the existence of a groundwater depression cone near the river, nitrate contamination can be transported into the aquifer with the flow. The average time lag of nitrate migration from the river to the zone of groundwater level depression cone is different in different sections, which shows an increasing trend from the upstream to downstream along the river, with an average of two to six months. It is mainly related to the stratigraphic structure, the migration distance, the hydraulic conductivities of the aquifer and the riverbed sediment. Compared with the case of considering the silt layer, the time lag of nitrate migration is greater than that of the case of ignoring the silt layer. The results will provide useful information for detecting nitrate concentrations in the alluvial fan area of the Luanhe river, northeast of the NCP (North China Plain).


Introduction
Riparian zones are defined as a transition area between the terrestrial and aquatic ecosystems [1,2].They play a crucial role in ecological systems and biodiversity because of the interaction of groundwater from inland and the river water from stream.Therefore, it is essential to understand the dynamics of exchange between surface water and groundwater.Aquifer-stream interaction occurs continuously and can transit between gaining and losing depending on many natural factors, including stream water level, stream bank slope, hydraulic properties of stream bed, groundwater table depth, and aquifer's hydraulic properties [3][4][5][6].In addition, human activities also impact stream-aquifer interaction, such as groundwater pumping near a stream, which has resulted in a decrease of groundwater table, causing aquifer depletion and water quality degradation.Therefore, understanding the hydraulic connectivity between streams and alluvial aquifers is a requirement for sustainable utilization of water resources, both for groundwater and surface water management.
The nitrate contamination in groundwater has become an environmental problem worldwide, especially in countries where groundwater is the main source of potable water.The research shows that nitrate is the most frequent pollutant in groundwaters [7][8][9].Since the end of the twentieth century, many countries and regions in Europe and America have carried out research on the influencing factors and control counter measures of nitrate pollution in groundwater [10][11][12].In many parts of China, groundwater is also polluted by nitrate to varying degrees [13].The northeast of the North China Plain (NCP) is known as a hotspot of severe groundwater depletion caused by irrigation, and the intensive use of nitrogen fertilizers has led to increasing groundwater nitrate pollution which, in places, exceeds the national maximum concentration for drinking water of 45 mg/L [14].Nitrate pollution is a serious problem affecting drinking water for a vast population.Nitrate leaching to groundwater has been found to be related to Nitrogen surpluses, hydrological conditions, land use and soil type [15].In response to this concern, there is a need for field observations and examination of the dispersal of nitrate pollutants in groundwater systems.Recently, an increasing number of studies on this specific topic and in this area have been carried out, including the field observations, data analysis and numerical modelling [13,14,[16][17][18].However, the data are still scarce with respect to the nitrate concentration affected by the overexploitation and fertilizers.
Understanding the nitrate sources and transformation processes is critical to detecting nitrate pollution [19,20].In contaminated zones, nitrate isotopes (δ 15 N-NO 3 , δ 18 O-NO 3 ) can be used to identify the source of nitrogen as mineral or organic fertilizer, sewage, or atmospheric deposition [21][22][23].The isotope mass-balance mixing models based on the δ 15 N-NO 3 and δ 18 O-NO 3 approach have been used to quantify nitrate sources [21,24,25].The method of combining isotopes and multivariate statistical analysis was also used to analyze the nitrate source in groundwater in multiple land-use areas [26].The δ 15 N-NO 3 composition of most nitrate sources falls between −10‰ and +25‰ [27].The δ 18 O-NO 3 of atmospheric nitrate shows a wide range between+25‰ and +70‰, and that of organic fertilizer has a value ranging from +17 to +25‰ [20,28,29].
In the vicinity of rivers, nitrate may be retained, transformed, or degraded by microbial communities [30].Here, the interaction between surface water and groundwater has an obvious influence on the migration process of nitrate of groundwater.Time series of contamination concentrations have been used to determine the travel time of river water into aquifers induced by pumping [31][32][33].However, there is a time lag for the nitrate movement from the surface water into the groundwater which depends on the aquifer properties [34].Therefore, determining the duration of time lag is of critical importance for mastering the mechanism of nitrate transport.
In the alluvial fan area of the Luanhe river, northeast of the NCP, groundwater is an important resource to be used for domestic, agricultural and industrial purposes, due to the lack of quantity of surface water resources.The previous research focused on the environmental geologic problem in the alluvial fan of the piedmont plain zone, such as the decrease of groundwater table and land subsidence [35].Huang et al. [36] analyzed the characteristics of the shallow groundwater quality in the alluvial fan area of the Luanhe river, North China.However, there is no report on nitrate pollution in this area, including nitrate distribution in groundwater, identification of pollution sources, and influencing factors.In this paper, we present field measurements of water table and nitrate concentration in the wells of the alluvial fan area of Luanhe river, North China.The environmental isotopic tracers δ 15 N-NO 3 and δ 18 O-NO 3 were used to determine the nitrate sources.The propagation time of nitrate into the aquifers following the river rise was estimated.The primary objectives of the study are to (1) quantify the nitrate concentration of groundwater and its isotopic composition, (2) identify the sources of nitrate pollution, and (3) estimate the time lag of nitrate in the Luanhe River transporting from the river to adjacent alluvial aquifer.The results will provide useful information for detecting nitrate concentrations in the alluvial fan area of the Luanhe river, northeast of the NCP.

Study Area
Our study area is located in the northeast of Luan county, Hebei Province, China.It is situated in the west of Luanhe River, the south of Yanshan Mountains, with a latitude 39 • 42 51" to 39 • 45 03" and a longitude 118 • 44 56" to 118 • 45 54" (Figure 1).Luanhe alluvial plain is a compound alluvial fan located at the front of the Yanshan Mountains, covering the northeastern part of the North China Plain.The study region is located in the front edge of the alluvial plain formed by Luanhe River.The elevation of the relief topography is between 20 m and 30 m.About 20% of the study area has been covered by forest, and 80% of it is covered by crops (wheat and corn).The area has a continental, semi-humid climate with an average annual temperature of 10.5 • C. Annual rainfall in this area is around 636.7 mm, and 86% of precipitation occurs from July to September.The mean annual pan evaporation is 1030.9mm.The major surface water resource in the area is the Luanhe River.In Luan county, the river has a length of 41 km from North to South, a drainage area of 91 km 2 , and a streambed width of 700 m to 1100 m. transporting from the river to adjacent alluvial aquifer.The results will provide useful information for detecting nitrate concentrations in the alluvial fan area of the Luanhe river, northeast of the NCP.

Study Area
Our study area is located in the northeast of Luan county, Hebei Province, China.It is situated in the west of Luanhe River, the south of Yanshan Mountains, with a latitude 39°42'51"to 39°45'03"and a longitude 118°44'56" to 118°45'54" (Figure 1).Luanhe alluvial plain is a compound alluvial fan located at the front of the Yanshan Mountains, covering the northeastern part of the North China Plain.The study region is located in the front edge of the alluvial plain formed by Luanhe River.The elevation of the relief topography is between 20 m and 30 m.About 20% of the study area has been covered by forest, and 80% of it is covered by crops (wheat and corn).The area has a continental, semi-humid climate with an average annual temperature of 10.5 °C.Annual rainfall in this area is around 636.7 mm, and 86% of precipitation occurs from July to September.The mean annual pan evaporation is 1030.9mm.The major surface water resource in the area is the Luanhe River.In Luan county, the river has a length of 41 km from North to South, a drainage area of 91 km 2 , and a streambed width of 700 m to 1100 m.

Field Site
The field survey showed that groundwater in the alluvial aquifer was contaminated on both sides of the river.Therefore, three hydrogeological cross sections were selected from north to south of the study area along Luan River, as shown in Figure 1.The monitored transects were arranged in the right bank side of the river, taking into account the symmetry.The distance between the upstream cross section A-A' and downstream cross section B-B' is about 5.5 km.There are five monitoring wells ('W') placed on the cross sections A-A' and B-B', respectively.W1 from the A-A′ and W6 from the B-B′ sections are situated in the river flood plain, and the distance of both W1 and W6 from the Luanhe river is 10 m.The middle cross section C-C' consisted of three wells W11-W13, which are 20 m, 50 m and 800 m away from the river bank, respectively.Therefore, there are a total of thirteen monitoring wells in the study area, as shown in Figure 1, which are suitable for local study.More details about the field arrangement can be found in Huang et al. [36].Some parameters related to the wells are listed in Table 1.

Field Site
The field survey showed that groundwater in the alluvial aquifer was contaminated on both sides of the river.Therefore, three hydrogeological cross sections were selected from north to south of the study area along Luan River, as shown in Figure 1.The monitored transects were arranged in the right bank side of the river, taking into account the symmetry.The distance between the upstream cross section A-A' and downstream cross section B-B' is about 5.5 km.There are five monitoring wells ('W') placed on the cross sections A-A' and B-B', respectively.W1 from the A-A and W6 from the B-B sections are situated in the river flood plain, and the distance of both W1 and W6 from the Luanhe river is 10 m.The middle cross section C-C' consisted of three wells W11-W13, which are 20 m, 50 m and 800 m away from the river bank, respectively.Therefore, there are a total of thirteen monitoring wells in the study area, as shown in Figure 1, which are suitable for local study.More details about the field arrangement can be found in Huang et al. [36].Some parameters related to the wells are listed in Table 1.

Geology and Hydrogeology
Based on the boring logs and high-density resistivity method, the stratigraphic structure at different cross section was identified.Sediments types include silty sand, clay, and gravel from the top to bottom (Figure 2).The sand layer and gravel layer are distributed continuously; however, the intermediate clay layer is discontinuous.The thicknesses of the silty sand, clay and gravel are about 1-7 m, 0-15 m and 8-23 m, respectively.In the vertical profile from the top to bottom, the particle size of various components is from fine to coarse, with a typical dual structure.The bottom of the sediments is bedrock, which is about 30 m below the surface based on the high-density resistivity method.The hydraulic conductivities of the layers are estimated based on the pumping test and slug test data.The permeability of the upper sand layer and lower gravel layer is high, and the hydraulic conductivities of these two layers are 1-6 and 100-350 m/day, respectively.The permeability of middle clay layer is relatively low, and the hydraulic conductivity of which is 0.01-0.1 m/day.The well W5 is located in a second-grade terrace, while the other wells are distributed in a first-grade terrace.The permeability of the layers increases with depth gradually.In the direction perpendicular to the river, the thickness of the silty sand and gravel become thin gradually, and that of clay increases.Therefore, the permeability of the aquifer near the river bank is greater, whereas it decreases far away from the river bank.Along the river longitudinally, from the upstream section A-A' to downstream section B-B', the thickness of sand and gravel becomes thin gradually, and the thickness of clay increases.Thus, the permeability of the layers has a decreasing trend from the upstream to downstream of the river.

Groundwater Level and Quality
Measurements of groundwater level and the river stage were conducted manually at time interval of five days from 11 November 2010 to 11 September 2011.Groundwater and river water samples were collected at time interval of two months from November 2010 to September 2011, and then they were transported to the laboratory for chemical analysis of major ions and nitrate concentration.In order to avoid the air entering groundwater during the sample, the observed wells were sampled for water quality at a depth 1 m below the groundwater level.Samples were collected using 1 L bottles and labeled.In brief, all wells were pumped for about 10 min before sampling using a 3-inch water pump.There were only twelve samples from wells at each time, because well W4 was destroyed.
In the current study, the focus was on nitrate in the wells sampled once every two months from November 2010 to September 2011.Huang et al. [36] provided the details of the sampling and analysis procedures for the major ion concentrations.The dissolved nitrate was analyzed spectrophotometrically on an automated ion flow injection analyzer (QuickChem Methods 10-107-04-1-C).The other physical parameters such as electrical conductivity (EC), total dissolved

Groundwater Level and Quality
Measurements of groundwater level and the river stage were conducted manually at time interval of five days from 11 November 2010 to 11 September 2011.Groundwater and river water samples were collected at time interval of two months from November 2010 to September 2011, and then they were transported to the laboratory for chemical analysis of major ions and nitrate concentration.In order to avoid the air entering groundwater during the sample, the observed wells were sampled for water quality at a depth 1 m below the groundwater level.Samples were collected using 1 L bottles and labeled.In brief, all wells were pumped for about 10 min before sampling using a 3-inch water pump.There were only twelve samples from wells at each time, because well W4 was destroyed.
In the current study, the focus was on nitrate in the wells sampled once every two months from November 2010 to September 2011.Huang et al. [36] provided the details of the sampling and analysis procedures for the major ion concentrations.The dissolved nitrate was analyzed spectrophotometrically on an automated ion flow injection analyzer (QuickChem Methods 10-107-04-1-C).The other physical parameters such as electrical conductivity (EC), total dissolved solids (TDS), and temperature and pH of the samples were measured in the field using portable EC and pH meters.
Isotopic compositions of 18 O and 15 N are reported using the standard per mil (δ) notation, which are calculated in terms of the ‰ deviation of the isotope ratio of the sample relative to the ratio in the V-SMOW(Vienna Standard Mean Ocean Water) standard.δ 18 O was prepared using the CO 2 equilibration method [37].Reproducibility of the samples calculated from standards systematically interspersed in the analytical batches was ±0.20‰ for δ 18 O.The dissolved NO 3 − is first transformed to AgNO 3 and then combusted to N 2 in sealed quartz tubes for δ 15 N analysis through a gas chromatograph into a Thermo-Finnigan Delta Plus isotope-ratio mass spectrometer.

Analysis of Surface Water-Groundwater Exchange
The water exchange between the river and aquifer can be quantitatively expressed by Darcy's law, namely where L is the stream length [L], W is the stream width [L], M is the thickness of riverbed sediment [L], K is the hydraulic conductivity [L/T] of the stream sediments, H R is the river level [L], H is the groundwater level near river [L], and Q is the quantity of recharge and discharge of groundwater [L 3 /T].When the river level is higher than the groundwater level of the aquifer, the value of Q is larger than 0, indicating that the ground water is recharged by the surface water.In contrast, the groundwater discharges from the aquifer to the river.

Analysis of the Migration Time of Pollution
When the groundwater velocity is calculated, the migration time of pollution in the groundwater can be estimated, according to the distance between the monitoring wells and the river.The equation for calculating the groundwater velocity V [L/T] is expressed as where L 1 is the distance between the river and well [L].
Then, the actual flow velocity u [L/T] can be calculated as where the effective porosity n e [-] can be expressed as where n is the porosity [-].
Then, the migration time of pollution from the river water to groundwater t [T] can be estimated

Groundwater and Surface Water Interactions
The drilling data of the thirteen monitoring wells show that there is a low permeable layer at the bottom of Luanhe river.The stream and multi-layered aquifer system can exchange water through the thin layer with low permeable material.Figure 3 shows the variations of the groundwater level observed in the wells and the stream level on the cross sections A-A', B-B' and C-C', respectively.One can see that the observed groundwater level in the wells and stream level varied with seasons.During the dry season from November 2010 to April 2011, the observed groundwater level decreased slightly and the groundwater discharged to the river.Then, the groundwater level exhibited a gradually increasing trend from May to September 2011 during the period of river stage rise.In addition, from Figure 2, one can see that there are cones of groundwater depression in the vicinity of wells W2, W7 and W12, because there are a large number of wells for civil use with an amount of pumping more than 10 thousand m 3 per day.Figure 4 shows the contour map of the average groundwater level in the study area.The groundwater level decreases from the section A-A' to section B-B'.This indicates that the groundwater flow direction is from north to south, which is basically consistent with that of the river water.
section A-A' to section B-B'.This indicates that the groundwater flow direction is from north to south, which is basically consistent with that of the river water.In order to investigate the interaction between river water and groundwater, the quantity of recharge from the river water to groundwater and that of discharge from groundwater to river water were calculated, based on the equation (1).Although groundwater flow paths in the riparian areas can be significantly more complex than simple vertical fluxes [38,39], it is still useful to compare river level with nearby groundwater level to determine if the reach is likely to be gaining or losing water [40,41].The hydraulic conductivity of riverbed sediment was estimated to be 0.032 m/d based on the field test.The total length of the calculated river reach is 5.5 km.The length of each section was determined according to the distribution of seepage field.The calculated width is equal to half of the width of river, considering the symmetry.Based on the field investigation, the thickness of the bottom sediment is 0.6 m, 0.7 m and 0.8 m for the sections A-A', B-B' and C-C',respectively.The river level is the level corresponding to each section.The groundwater level is equal to nearby groundwater level in each section (well 1 for section A-A', well 11 for section C-C' and well 6 for section B-B').Therefore, the amount of recharge and discharge is estimated based on the river level and groundwater level according to the Darcy's law.The amount of recharge and discharge in different sections are shown in Figure 5. From the figure, one can see that the river water recharges groundwater during most of time from November 2010 to September 2011.The amount of groundwater supplied by the river is 346.3 × 10 4 m 3 /a, and the discharge of groundwater to the river is about 14.08 × 10 4 m 3 /a during this period (Table 2).

Table 2.
Total recharge and discharge of groundwater in the study area(10 4 m 3 /a).

Recharge and Discharge
The Section of A-A'

The Section of B-B'
The Section of C-C'  In order to investigate the interaction between river water and groundwater, the quantity of recharge from the river water to groundwater and that of discharge from groundwater to river water were calculated, based on the equation (1).Although groundwater flow paths in the riparian areas can be significantly more complex than simple vertical fluxes [38,39], it is still useful to compare river level with nearby groundwater level to determine if the reach is likely to be gaining or losing water [40,41].The hydraulic conductivity of riverbed sediment was estimated to be 0.032 m/d based on the field test.The total length of the calculated river reach is 5.5 km.The length of each section was determined according to the distribution of seepage field.The calculated width is equal to half of the width of river, considering the symmetry.Based on the field investigation, the thickness of the bottom sediment is 0.6 m, 0.7 m and 0.8 m for the sections A-A', B-B' and C-C', respectively.The river level is the level corresponding to each section.The groundwater level is equal to nearby groundwater level in each section (well 1 for section A-A', well 11 for section C-C' and well 6 for section B-B').Therefore, the amount of recharge and discharge is estimated based on the river level and groundwater level according to the Darcy's law.The amount of recharge and discharge in different sections are shown in Figure 5. From the figure, one can see that the river water recharges groundwater during most of time from November 2010 to September 2011.The amount of groundwater supplied by the river is 346.3 × 10 4 m 3 /a, and the discharge of groundwater to the river is about 14.08 × 10 4 m 3 /a during this period (Table 2).

Chemistry and Nitrate of Groundwater
The total dissolved solids (TDS), pH, electrical conductivity, temperature and nitrate concentration measured in the groundwater and surface water are shown in Table 3.The electrical conductivity (EC), which is an indirect measure of TDS, varies between 477 and 1844 µS/cm with an average of 972 µS/cm (Table 3).In general, the samples are slightly alkaline with pH values ranging between 7.0 and 8.5.The temperature of surface water samples ranges from 2.9 • C to 26.7 • C, and that of groundwater samples is between 9.9 • C and 18.6 • C. TDS values range from as low as 294 mg/L to a maximum of 1500 mg/L with an average of 629 mg/L.Table 3 shows that nitrate concentration is spatially and temporally variable.Three typical two-dimensional profiles of A-A', B-B' and C-C' were selected to study the nitrate variation in groundwater and surface water in different months.Figure 6a shows the variations in nitrate concentrations of groundwater and surface water in the transect A-A'.The nitrate concentration of surface water varied through time from 3.34 mg/L to 6.88 mg/L.The nitrate concentrations of surface water were higher from November 2010 to March 2011 than those sampled during the months from May to September 2011.The river level rises in May 2011, because the upstream sluice was open and the river water was released.This suggests that the nitrate concentration was diluted by the river water flow from May 2011.For groundwater samples at wells W1-3, nitrate concentration changes very large with time ranging from 4.18 mg/L to 41.6 mg/L.Compared with the nitrate concentration of groundwater in wells far from the river, the nitrate concentration of river water was lower, indicating that the pollution of river water has little influence on groundwater far from the river.The nitrate concentration in well W1 was affected by the river level, because it is near the river.The observed nitrate concentration increased from 17.3 mg/L to 29.3 mg/L when the groundwater level was lower than the river level, and then decreaseed to 4.18 mg/L as the groundwater level increased from November 2010 to July 2011.Then, it tended to increase slowly from July to September 2011.The nitrate concentration in well W2 as lower than that of well W1.This may be due to the long-term pumping of groundwater in well W2.For W3, the nitrate concentration decreased from 47 mg/L to 17 mg/L over the period from November 2010-March 2011, and then increased to 28.2 mg/L in September 2011.Generally, the nitrate concentrations of groundwater in the transect A-A' were low from May to July of 2011, because the groundwater level rises as the rainfall increases during this period, which make the nitrate in the aquifer diluted.
Compared with Figure 6a, in Figure 6b, one can see that the variation in nitrate concentrations of surface water in the transect A-A' was similar to those of surface water in the transect B-B'.Nitrate concentrations of surface water varied little over time, ranging from 0.07 mg/L to 6.48 mg/L.However, the nitrate concentrations of surface water in the transect B-B' were lower than those of the surface water in transect A-A'.This can be attributed to the different structure of aquifer of these two transects.This was mainly due to the large particle size of soil in the transect A-A' compared with that of transect B-B' (Figure 2).Therefore, the velocity of groundwater flow in the transect A-A' is faster than that in the transect B-B'.As a result, the nitrate in groundwater is easily discharged to the river water in the transect A-A'. Figure 6b shows that the nitrate concentration of groundwater undergoes large changes through time ranging from 0.07 mg/L to 60.3 mg/L in the transect B-B'.Similar to the well W1, the well W6 near the river is also affected by the river water, and the nitrate concentration variation is basically the same as that of the river.The nitrate concentration increased from 0.97 mg/L to 7.36 mg/L, and then decreased to 0.09 mg/L as the groundwater level increased from November 2010 to July 2011.It increased to 1 mg/L slowly from July to September 2011.The variation in nitrate concentration of well W7 is basically consistent with that of the river water in the transect B-B'; however, it was lower than that of the river water during the period of November 2010-March 2011.This was mainly due to the groundwater pumping for drinking near well W7.In particular, the nitrate concentration in well W7 was lower than the detection limit in July 2011, which was assumed to be 0 mg/L.For the sampling well W8, the nitrate concentration increased from 23.7 mg/L to 24.6 mg/L, and decreased to be 15.0 mg/L from November 2010 to March2011.Then, it continued to increase, reaching a maximum value of 30.5 mg/L in September 2011.The variation trend of nitrate concentration in well W9 was similar to that in well W8.It decreased from 42.6 mg/L to 6.83 mg/L, and then increased to 30.4 mg/L from November 2010 to September 2011.The nitrate concentration in well W10 showed considerable variation, ranging from 20.2 mg/L to 60.3 mg/L.During the period of January-July 2011, the nitrate concentration as above the national drinking water standards (45 mg/L) in China.The nitrate concentrations in well W8-10 were greater than those of well W6-7.This was mainly due to the fact that the observation well W8-9 are located in farmland, and the surface fertilization results in higher nitrate concentrations in the wells.For well W10, it might be due to the discharge of domestic sewage, because it is located near a village.The nitrate concentration in the river water of the C-C' section is the average value of nitrate concentration in the A-A' section and that in the B-B' section, because the C-C' section is located between the A-A' section and B-B' section.Figure 6c shows that the variation in nitrate concentrations of surface water in the transect C-C' changes from 1.67 mg/L to 6.43 mg/L.The nitrate concentration is higher when the river level is low; however, it is lower when the river level is high.This indicates that the nitrate concentration is diluted as the increasing of river water, and the nitrate source of river water is related to the pollution of surface water.From Figure 6c, one can see that the nitrate concentration of groundwater changes considerably over time, with values ranging between 1.83 mg/L and 34.9 mg/L in the transect C-C'.The nitrate concentration in well W11 increased from 28.7 mg/L to 34.9 mg/L, and then decreased to be 9.6 mg/L from November 2010 to September 2011.However, there is no observed data in well W11 after May 2011, because the well was scrapped.For well W12, the nitrate concentration firstly increased from 17.7 mg/L to 24.9 mg/L, and then decreased to 20.4 mg/L from November 2010 to September 2011.The nitrate concentration in well W13 decreased from 17.2 mg/L to 13 mg/L, and then increased to a maximum value of 24.6 mg/L in March 2011.Finally, it decreased to 1.83 mg/L from March to September 2011.

Nitrate Sources
The nitrate in river water and groundwater mainly comes from atmospheric deposition, inorganic fertilizers, soil, livestock manure and sewage.The nitrate from different sources has different isotopic ranges, and the nitrate source can be judged according to the value of δ 15 N-NO 3 in groundwater.The δ 15 N-NO 3 values for atmospheric nitrogen deposition ranges from −13‰ to 13‰ [20,22].Nitrate originating from synthetic fertilizers and natural soil organic matter is characterized by δ 15 N-NO 3 values ranging from −4‰ to 4‰ and from 4‰ to 8‰, respectively, which is distinct from nitrogen in sewage (8‰~15‰) and livestock manure(10‰~22‰) [42][43][44].Combined with the hydrogeological condition of the study area, the nitrate source of groundwater and river water can be deduced based on the values of δ 15 N-NO 3 .However, despite its effectiveness, nitrate sources cannot be identified by the δ 15 N signature alone due to the δ 15 N overlapping ranges and the occurrence of nitrate isotopic fractionation [19,45].The relationship between the two isotopes (δ 15 N-NO 3 and δ 18 O-NO 3 ) can be further used to judge the nitrate source of groundwater and river water [46,47].
The values of 18 O and 15 N are listed in Table 3.The average standard deviation of normalized data is 0.3‰ for δ 15 N. Isotopic composition of nitrate in groundwater and river water showed seasonal variation.The range of δ 15 N-NO 3 values was from 3.8 to 139‰ for groundwater and from 7.2 to 120.1‰ for surface water, showing a large variance.The maximum δ 15 N-NO 3 value of 120.1‰ occurred at surface water in the transect A-A' in July 2011, when the river NO 3 -N concentration was at its minimum.The δ 18 O-NO 3 values ranged from −7.66‰ to −6.34‰ for groundwater and from −7.4‰ to −6.5‰ for surface water, with the mean values of −7.15 ‰ for groundwater and −6.95 ‰ for surface water.The highest δ 18 O-NO 3 value occurred at well W1 in September 2011 (−6.34‰), while the lowest value (−7.66‰) of δ 18 O-NO 3 was obtained at well W7.As can be seen from Figure 7, there was a tendency for both δ 15 N-NO 3 and δ 18 O-NO 3 to increase with N-NO 3 concentration.According to the bi-plot of δ 15 N-NO 3 versus δ 18 O-NO 3 for potential nitrate sources, the isotopic signatures of nitrate in groundwater and river water fall within the window for sewage, livestock manure and natural soil organic matter (Figure 8).This indicates that a significant positive relationship was discernable for the whole dataset.from nitrogen in sewage (8‰~15‰) and livestock manure(10‰~22‰) [42][43][44].Combined with the hydrogeological condition of the study area, the nitrate source of groundwater and river water can be deduced based on the values of δ 15 N-NO3.However, despite its effectiveness, nitrate sources cannot be identified by the δ 15 N signature alone due to the δ 15 N overlapping ranges and the occurrence of nitrate isotopic fractionation [19,45].The relationship between the two isotopes (δ 15 N-NO3 and δ 18 O-NO3) can be further used to judge the nitrate source of groundwater and river water [46,47].The values of 18 O and 15 N are listed in Table 3.The average standard deviation of normalized data is 0.3‰ for δ 15 N. Isotopic composition of nitrate in groundwater and river water showed seasonal variation.The range of δ 15 N-NO3 values was from 3.8 to 139‰ for groundwater and from 7.2 to 120.1‰ for surface water, showing a large variance.The maximum δ 15 N-NO3 value of 120.1‰ occurred at surface water in the transect A-A' in July 2011, when the river NO3-N concentration was at its minimum.The δ 18 O-NO3 values ranged from −7.66‰ to −6.34‰ for groundwater and from −7.4‰ to −6.5‰ for surface water, with the mean values of −7.15 ‰ for groundwater and −6.95 ‰ for surface water.The highest δ 18 O-NO3 value occurred at well W1 in September 2011 (−6.34‰), while the lowest value (−7.66‰) of δ 18 O-NO3 was obtained at well W7.As can be seen from Figure 7, there was a tendency for both δ 15 N-NO3 and δ 18 O-NO3 to increase with N-NO3 concentration.According to the bi-plot of δ 15 N-NO3 versus δ 18 O-NO3 for potential nitrate sources, the isotopic signatures of nitrate in groundwater and river water fall within the window for sewage, livestock manure and natural soil organic matter (Figure 8).This indicates that a significant positive relationship was discernable for the whole dataset.

The Time of Nitrate Transport from River Water to Groundwater
Recently, an increasing number of studies on solute transport have been published, including the field observations, data analysis and numerical modelling [16][17][18]48,49].However, data are still scarce with respect to the nitrate concentration affected by the over exploitation.Generally, some time is needed for nitrate transport from the river to groundwater.The migration time of groundwater flow near the bank is affected by the hydraulic conductivities of the riverbed sediments and adjacent aquifers, the distance between river and inland wells, and the thickness of riverbed sediments.The nitrate in river water movement into the adjacent aquifers should be considered in zone with groundwater level depression cones.The pollution of river water has little

The Time of Nitrate Transport from River Water to Groundwater
Recently, an increasing number of studies on solute transport have been published, including the field observations, data analysis and numerical modelling [16][17][18]48,49].However, data are still scarce with respect to the nitrate concentration affected by the over exploitation.Generally, some time is needed for nitrate transport from the river to groundwater.The migration time of groundwater flow near the bank is affected by the hydraulic conductivities of the riverbed sediments and adjacent aquifers, the distance between river and inland wells, and the thickness of riverbed sediments.The nitrate in river water movement into the adjacent aquifers should be considered in zone with groundwater level depression cones.The pollution of river water has little

The Time of Nitrate Transport from River Water to Groundwater
Recently, an increasing number of studies on solute transport have been published, including the field observations, data analysis and numerical modelling [16][17][18]48,49].However, data are still scarce with respect to the nitrate concentration affected by the over exploitation.Generally, some time is needed for nitrate transport from the river to groundwater.The migration time of groundwater flow near the bank is affected by the hydraulic conductivities of the riverbed sediments and adjacent aquifers, the distance between river and inland wells, and the thickness of riverbed sediments.The nitrate in river water movement into the adjacent aquifers should be considered in zone with groundwater level depression cones.The pollution of river water has little influence on the zone in the west of the groundwater level depression cone.According to the hydraulic conductivity of the river bed and the aquifer, the distance between the wells and river and the thickness of riverbed, the migration time of groundwater flow in the zone of groundwater level depression cone may be calculated.In order to estimate the time of nitrate transport from the river to groundwater, two cases have been considered to compute the time lag of nitrate, i.e., the cases of considering the silt of the riverbed and ignoring it.
The average time of river water passing through the riverbed is 45 days, considering the influence of bottom sediments of the A-A' section.However, the average migration time of nitrate from the river to the zone of groundwater level depression cone ranges from 56 to 77 days.This indicates that there is obvious hysteresis in the nitrate transporting from the river water to groundwater.As can be seen from Figure 6a, the peak of nitrate concentration in the river water occurred in January 2011; however, that in the well W1 was in March 2011.The peak nitrate concentration in well W2 was not clear.The nitrate concentration in well W2 can be determined to have peaked in April 2011, according to the estimation of the migration time of nitrate contamination in the river.The peak value of nitrate concentration in the river water in January 2011 corresponds to that in well W1 in March 2011 and that in well W2 in April 2011; therefore, the time lag is between two and three months.The average migration time of nitrate from the river to the zone of groundwater level depression cone is 12 to 23 days when the effect of riverbed sediment is not considered.One can see that the time lag of nitrate migration in the nearest well is less than that in the well far away from the river for the A-A' section under the condition of two cases.Taking the case of ignoring the silt of the riverbed as an example, the time lag of nitrate contamination in well W1 is 12 days, while that in well W2 is 23 days.This indicates that the time of contamination transporting to groundwater increases with the distance.
In the transect B-B', the average time of river water passing through the riverbed is 59 days, which is larger than that in the transect A-A', because the hydraulic conductivity of riverbed sediment in the transect B-B' is lower than that in the transect A-A'.Similar to the transect A-A', there is a time lag in the nitrate transporting from the river water to groundwater in the the transect B-B'.The average migration time of nitrate from the river to the zone of groundwater level depression cone ranges from 72 to 85 days.In addition, Figure 6b shows that the peak value of nitrate concentration in the river appeared in November 2010.However, the peaks of nitrate concentration in the wells W6 and W7 were in March 2011, which indicates that the time lag is four months.The migration time of nitrate from the river to the zone of groundwater level depression cone is ranges from 14 to 27 days when the effect of riverbed sediment is ignored.
Similar to the transects A-A' and B-B', the average time of river water passing through the riverbed is 53 days in the transect C-C'.The average migration time of nitrate from the river to the zone of groundwater level depression cone ranges from 65 to 82 days when the effect of bottom sediments of the river is considered.Additionally, it can be seen from Figure 6c that there is obvious hysteresis in the nitrate transporting from the river to groundwater.The peak value of nitrate concentration in the river water in November 2010 corresponds to that of well W11 in January 2011 and that of well W12 in May 2011, respectively.Thus, the time lag for the migration of nitrate from the river to groundwater is two to six months, which is longer than that obtained by the analytical method.This is because the migration of nitrate is considered in the analytical method without considering the dispersion effect.However, without consideration of the effect of riverbed sediment, the migration time of nitrate from the river to the zone of groundwater level depression cone ranges from 13 to 20 days.

Conclusions
In the present study, the variation of groundwater level and nitrate distribution in the alluvial fan area of Luanhe river, North China Plain, was investigated.Based on the field observations, a three-layered structure was identified which is characterized by a high permeable sand layer, a sand gravel layer and a clay aquitard between them.The measured results showed that the groundwater level varied with seasons.During the dry season from November 2010 to April 2011, the observed groundwater level decreased slightly and the groundwater discharged to the river.Then, the groundwater level increased gradually from May to September 2011 during the period of river stage rise.Typically, the river water recharged groundwater from November 2010 to September 2011.The amount of groundwater supplied by the river is 346.3 × 10 4 m 3 /a, and the discharge of groundwater to the river is about 14.08 × 10 4 m 3 /a during this period.
The nitrate concentration of groundwater and river water varied with seasons.The nitrate of groundwater in wells near the river is affected by the river water, and the nitrate concentration variation is basically the same as that of the river.The nitrate concentration decreases as the groundwater level increased from November 2010 to July 2011.The nitrate concentrations in the zone of groundwater level depression cone were lower than those in the wells near the river, due to the long-term pumping of groundwater.However, the nitrate concentration of river water has little influence on that of groundwater in wells far from the river.
The values of δ 15 N-NO 3 and the relationship between the two isotopes (δ 15 N-NO 3 and δ 18 O-NO 3 ) were used to judge the nitrate source of groundwater and river water.The results showed that the main sources of nitrate in groundwater are sewage, livestock manure and natural soil organic matter.Due to the existence of a groundwater depression cone near the wells W2, W7 and W12, the nitrate contamination can be transport into the aquifer near the river with the river flow.There is obvious hysteresis in the nitrate transporting from the river water to groundwater.The average time lag of nitrate migration from the river to the zone of groundwater level depression cone is different in different sections, which shows an increasing trend from the upstream to downstream along the river, with an average of two to six months.This is mainly related to the stratigraphic structure, the migration distance, the hydraulic conductivities of the aquifer and the riverbed sediment.Compared with the case of considering the silt layer, the time lag of nitrate migration is greater than that of the case in which the silt layer was ignored.The results will provide useful information for detecting nitrate concentrations in the alluvial fan area of the Luanhe river, northeast of the NCP.

Figure 1 .
Figure 1.Location of the observation wells and cross-sections in the study area.

Figure 1 .
Figure 1.Location of the observation wells and cross-sections in the study area.

Figure 2 .
Figure 2. Schematic hydrogeological cross section of the Luanhe alluvial fan and groundwater level change along (a) transect A-A', (b) transect B-B' and (c) transect C-C'.

Figure 2 .
Figure 2. Schematic hydrogeological cross section of the Luanhe alluvial fan and groundwater level change along (a) transect A-A', (b) transect B-B' and (c) transect C-C'.

Figure 3 .
Figure 3. Temporal variations of the groundwater level observed in the wells along (a) transect A-A', (b) transect B-B' and (c) transect C-C', respectively.

Figure 3 .
Figure 3. Temporal variations of the groundwater level observed in the wells along (a) transect A-A', (b) transect B-B' and (c) transect C-C', respectively.

Figure 4 .
Figure 4. Contour map of the average groundwater level in the study area, where the red lines represent the zone of the study area and the cross-sections, and the black lines represent the average groundwater level.

Figure 4 .
Figure 4. Contour map of the average groundwater level in the study area, where the red lines represent the zone of the study area and the cross-sections, and the black lines represent the average groundwater level.

Figure 5 .
Figure 5.The amount of recharge and discharge in different sections (a) transect A-A', (b) transect B-B', and (c) transect C-C'.

Figure 5 .
Figure 5.The amount of recharge and discharge in different sections (a) transect A-A', (b) transect B-B', and (c) transect C-C'.

Figure 6 .
Figure 6.Nitrate concentration of groundwater and surface water for different transects, (a) A-A', (b)B-B' and (c) C-C'.

Figure 6 .
Figure 6.Nitrate concentration of groundwater and surface water for different transects, (a) A-A', (b) B-B' and (c) C-C'.

Table 2 .
Total recharge and discharge of groundwater in the study area (10 4 m 3 /a).

Table 3 .
Summary of physical and chemical data obtained for surface and groundwater samples in the study area.