Isotope-Based Study on Nitrate Sources in a Karst Wetland Water, Southwest China

: Huixian karst wetland is the largest karst wetland in the subtropical low-altitude areas in China, which is known as the kidney of Lijiang River. This study was focused on researching the NO 3 − source and N cycling in Huixian karst wetland using multi-isotope, hydro-chemical data, and stable isotope analysis in R model, including 7 groundwater samples, 6 ﬁsh pond water samples, and 17 surface river water samples, respectively. Hydro-chemical data showed that the maximum concentrations of NO 3 − (49.3 mg · L − 1 ) and NH 4+ (0.6 mg · L − 1 ) in the groundwater exceeded the limit of Groundwater Quality Standard Class III Standard by 2.5 and 1.2 times, respectively; the maximum concentration of NH 4+ (8.2 mg · L − 1 ) in ﬁsh pond water exceeded the limit of the Groundwater Quality Standard Class III Standard by 16.3 times; the maximum concentrations of NO 3 − (24.5 mg · L − 1 ), NO 2 − (7.3 mg · L − 1 ), and NH 4+ (2.4 mg · L − 1 ) in surface river water were 1.2, 7.3, and 17.3 times higher than the limit of the limit of Groundwater Quality Standard Class III Standard, respectively. All water samples were in an aerobic environment with the range of dissolved oxygen (2.6–14.8 mg · L − 1 ). Therefore, nitriﬁcation mainly occurred in the process of nitrate transformation, signiﬁcant denitriﬁcation did not. Multi-isotope ( δ 15 N-NO 3 , δ 18 O-NO 3 , and 13 C DIC ) and stable isotope analysis in R model revealed that manure sewage, ammonium nitrogen fertilizer, and soil nitrogen were the three dominant NO 3 − sources, contributing to (39.8 ± 6.4)%, (33.4 ± 5.1)%, and (26.8 ± 3.9)%, respectively.


Introduction
Wetlands are important water resources and have the laudatory name of "the kidney of the earth". Similarly, karst wetlands are called "the kidneys of karst areas". In recent years, in order to pursue economic benefits, man-made unreasonable use of wetlands caused a sharp decline in wetland water environment quality: wetland water quality deteriorated and disappeared [1]. The eutrophication of wetland water has become an ecological problem of great concern in our country and even the whole world [2][3][4]. Nitrogen pollution in water is one the best representations of eutrophication [5]. Nitrogen is soluble and highly mobile in a water environment. Once exceeding the standard, it is very difficult to control nitrogen pollution [6]. There were many reports that karst water being polluted by NO 3 − due to man-made factors [7,8]. Drinking groundwater with a high concentration of nitrogen for a long time can cause cancer, teratogenicity, and mutagenesis, which is seriously harmful to health. Because of the sensitivity and fragility of karst water, when a karst aquifer is polluted by NO 3 − it not only seriously threatens the drinking water source in karst areas but also harms the surrounding ecological environment [9]. Therefore, it is necessary to effectively prevent and control NO 3 − pollution in karst areas, and the core of solving the problem is to study the transformation process and source of NO 3 − . In complex systems, however, where overlapping isotopic ranges of NO 3 − source end-numbers occur together with interference caused by nitrification, denitrification, ammonia volatilization, anammox, and/or other unknown processes, the apportionment of sources is complicated [10,11]. The identification of NO 3 − sources has been identified as the first step to mitigate NO 3 − pollution in aquatic environments. Given its inherent advantages, nitrate isotopes ( 15 N-NO 3 and 18 O-NO 3 ) have been widely used to trace NO 3 − sources in the past two decades [12][13][14][15]. In addition, previous studies have successfully used 15 N-NO 3 and 18 O-NO 3 combining with NO 3 − concentrations to quantify NO 3 − sources within several karst watersheds [16,17]. Therefore, a traceability study using nitrate nitrogen and oxygen isotopes is the most extensive research method to identify nitrate pollution sources in water. It has been confirmed that there was a positive correlation between man-made input of nitric acid and 13 C DIC in karst areas [18][19][20].
In recent years, the enrichment ratio of δ 15 N-NO 3 to δ 18 O-NO 3 and the relationship between δ 13 C DIC and NO 3 − have been successfully applied to identify the transformation process of NO 3 − in water [21,22]. Predecessors have successfully used the method of combining hydrochemistry with nitrate nitrogen and oxygen isotopes (δ 15 N-NO 3 and δ 18 O-NO 3 ) to trace the sources of NO 3 − in various types of water bodies, and many results have been achieved [23,24]. At present, the most effective method for tracing the source of nitrate is the stable isotope analysis in R model [25], which makes it possible for quantitative research on the NO 3 − source [26]. Therefore, multi-isotope qualitative analysis and stable isotope analysis in R model are the most effective means of analysis of NO 3 − source distribution characteristics, pollution status, and pollution causes in water [11,27].
Huixian karst wetland, known as "the kidney of the Lijiang River", is the largest subtropical low-altitude karst ecological wetland in China, which is located in Guilin, the most beautiful landscape in the world. In recent decades, in order to pursue economic benefits, human beings have increased the demand for Huixian wetland. Excessive reclamation of wetlands, excessive use of pesticides, fertilizers, and livestock breeding resulted in serious eutrophication of wetlands, which has attracted great attention of government departments and scholars. Peng et al. [28] determined the NO 3 − source and the linear relationship between groundwater runoff conditions and NO 3 − concentration in the core area of Huixian wetland by multi-isotopes means. Li et al. [29] analyzed the present situation of surface lake water quality in Huixian wetland and found that the concentration of TN and NH 3 -N decreased from upstream to downstream of the wetland. Wan et al. [30] analyzed the temporal and spatial differences in nitrogen and phosphorus emissions from the main surface water systems of wetlands by using ArcGIS software and a SWAT model (Environmental Systems Research Institute, Redlands, CA, USA). The results showed that total nitrogen (TN) was less affected by spatio-temporal factors, but greatly affected by land use types, industrial wastewater, aquaculture sewage, and domestic sewage. Xie [31] analyzed nitrate nitrogen and total nitrogen in surface water, soil, and groundwater of Huixian wetland, and revealed the temporal and spatial difference characteristics and influencing factors of nitrogen pollution. By using grey correlation analysis and correlation analysis, it was concluded that there was a close relationship between nitrogen concentration and pH, and there was a certain correlation between underlying surface properties and nitrogen and phosphorus emission concentration in the test area.
The previous studies have mainly focused on the water quality characteristics and influencing factors in a small area (about 20 km 2 ) in the core area of the wetland. However, the Huixian karst wetland is a wetland with karst groundwater recharge. It was not systematic enough to study only within the wetland core area. On the other hand, the area of fish ponds in the wetland is about 2 km 2 , accounting for about 1/10 of the wetland area. Changing part of the wetlands into fish ponds, in order to ensure that there is enough water in the fish ponds, fishermen cut off the original surface water, resulting in flooding in the wetlands near the fish ponds during the rainy season and serious water shortage in the dry season. As a result, the ecological environment of the wetland is seriously affected, and the discharge of water from fish ponds also has an important impact on the water quality of the wetland. However, most scholars have conducted sampling research on groundwater and surface river water, while ignoring fish pond water research in wetlands.
In addition, the dense population, cultivated land, and agricultural activities in wetlands have a huge influence on groundwater. Therefore, only by carrying out research on the hydrochemical characteristics of various types of water in the whole basin can we better and comprehensively grasp the causes of wetland hydrochemical formation. This study systematically sampled groundwater, surface rivers, and fish ponds in the entire catchment area of the wetland (about 90 km 2 ), and using hydrochemistry, multi-isotope (δ 15 N-NO 3 , δ 18 O-NO 3 ) methods, and δ 13 C DIC to analyze the NO 3 − transformation process in Huixian wetland water, and using the stable isotope analysis in R model to quantitatively calculate the contribution rates of different NO 3 − inputs. This study took the whole Huixian wetland basin as the research object, which has not been done before. Secondly, this study not only discussed the surface river water and groundwater in the study area, but also the artificial fish ponds which have a great impact on the wetlands. Thus, reasonable suggestions are put forward for water pollution control and water resources protection of Huixian wetland and similar wetlands. This study focused on: (1) revealing the characteristics of nitrogen pollution in Huixian wetland; (2) exploring the NO 3 − transformation of nitrogen in the study area; and (3) calculating the contribution ratio of different inputs of NO 3 − in Huixian wetland.

Study Area
Huixian karst wetland is located in the area of Sitang Town and Huixian Town, Lingui District, Guilin City, about 21 km away from the southwest of Guilin City. Geographical location: east longitude 110 • 09 23 -110 • 14 36 , north latitude 25 • 04 58 -25 • 08 25 , the altitude is between 150 and 170 m. It is located in the subtropical monsoon climate, with the average annual precipitation of 1886 mm. The rainy season is from April to August, and the precipitation accounts for 70% of the whole year. The dry season is from September to March of the next year, and the precipitation accounts for 30% of the whole year [32]. According to the hydrological observation data for many years, June to August is the wet season, and September to November is the monthly average water period, and the dry season is March to May and December to February of the following year.
The stratigraphic lithology in the study area is Donggangling (D 2 d) Formation limestone and dolomite, Dongcun (D 3 d) Formation limestone, Guilin (D 3 g) Formation limestone, Yanguan (C 1 y) Formation limestone and Datang (C 1 d) Formation limestone [32]. In addition, a small amount of Xindu (D 2 x) Formation siltstone is developed ( Figure 1). The water-bearing rock groups are mainly Quaternary loose rock pore aquifers, medium to thick pure carbonate rock aquifers and impure carbonate rock aquifers. The groundwater hydrological system of Huixian karst wetland is regulated by meteorology, hydrology, and underlying surface, and the recharge sources include atmospheric precipitation, swamp water infiltration recharge, and external water lateral recharge. Groundwater flow is comprehensively controlled by topography and structure, with fanshaped runoff from the north and south to the west, middle, and east. There are many forms of groundwater discharge, mainly undercurrent, springs, and molten pools. Surface water is the main source of recharge to groundwater in wet seasons, and surface water is mainly supplied by groundwater in dry seasons [24]. The wetland water system is developed under the effect of rivers including the Taiping River and Qingshui River in the west, which eventually merge into the Acacia River, and the tributaries in the east merge into the Liangfeng River. According to the survey, most of the original natural lakes and wasteland have been converted into fish farms, etc., which seriously damaged the original hydrological cycle conditions of the wetland. There are a large number of village residential areas in the study area, and the main sources of income for the villagers are crop cultivation on arable land, fish farming, and aquaculture in fish ponds. Therefore, arable land and fish ponds are the most widely distributed in the area, which are almost all over the study area.

Water Sampling and Testing
The sampling campaign of the entire watershed in Huixian wetland was performed from 13 to 15 January 2021, including 7 groundwater samples (GW01-GW07), 6 fish pond water samples (PW01-PW06), and 17 surface river water samples (SW01-SW17). Moreover, the groundwater is collected from the outlet of the underground river and the karst pool water. The distribution of sampling points is shown in Figure 1.
The pH, dissolved oxygen (DO), temperature (T), oxidation reduction potential (Eh), and electrical conductivity of the water body were tested with a hand-held water quality meter (WTW3430, GER) on site, and the resolutions was 0.01, 0.01 mg·L −1 , and 0.01 • C, respectively. Ca 2+ and HCO 3 − were tested with a hardness tester and alkalinity meter on site, respectively. The samples for testing conventional ions of water chemistry were collected with a 500 mL polyethylene plastic bottle, and the samples for testing cations, heavy metals, and trace elements were acidified with high-grade pure nitric acid (HNO 3 ) to pH < 2 on site. The δ 15 N-NO 3 , δ 18 O-NO 3 , and δ 13 C DIC samples were collected in 50 mL brown polyethylene plastic bottles, and 3-5 drops of saturated HgCl 2 were added to the δ 13 C DIC samples. The sampling bottle was rinsed 3-4 times with filtered (0.45 µm) water before loading the sample. All water samples were refrigerated at 4 • C for testing. The NO 3 − and Cl − were measured by ion chromatograph (DionexICS-1100) with the resolutions of 0.05 and 0.1 mg·L −1 , respectively. The NO 2 − and NH 4 + were measured by T6 Xinyue visible spectrophotometer with the accuracy of 0.002 mg·L −1 and 0.02 mg·L −1 , respectively. The carbon isotope of DIC (δ 13 C DIC ) was measured by GasBench II (Thermo Fisher Scientific Corp., Waltham, MA, USA) after being acidified by phosphoric acid. The value of δ 13 C DIC was calibrated using three international standards (NBS-18, NBS-19, and NBS-30) and four internal standards (GBW04405, GBW04406, GBW04416, and GBW04417), and the analytical precision was generally better than 0.2‰. The above sample determination work was completed in the Karst Geological Resources and Environmental Supervision and Testing Center of the Ministry of Natural Resources.
Nitrate isotopes ( 15 N-NO 3 and 18 O-NO 3 ) were determined using the bacterial denitrifier method after NO 2 removal with sulfamic acid [33,34]. The isotopic compositions of produced N 2 O (NO 3 − → N 2 O) were determined on Delta-plus (Thermo-Finnigan Corp, Sili-con Valley, CA, USA) fitted to auto-sampler equipment (CTC Analytics AG, Zwingen, Basel, Switzerland) and a Trace Gas Pre-concentrator unit. The analytical quality was controlled by blanks and four international standards (USGS-32, USGS-34, USGS-35, and IAEA-N3) and the analytical precision was generally better than 0.1‰ for δ 15 N-NO 3 and 0.5‰ for 18 O-NO 3 , respectively. The relative ratio of N and O isotopes in nitrate is expressed by δ: In the Equation (1)

Stable Isotope Analysis in R Mixing Model
The Stable Isotope Analysis in R (SIAR) model is a stable isotope mixture model based on R statistical software developed by Parnell et al. Its application makes the source analysis of pollutants change from qualitative to quantitative [11]. The premise of the stable isotope analysis in R model is that there is no obvious denitrification in the water body. By defining j isotopes of n mixtures from k sources, and considering the influence of isotopic fractionation, it is expressed as Equation (2): In the Equation, X ij is the δ value of the mixture isotope, where: I = 1, 2, 3. . . n, j = 1, 2, 3. . . J; S jk is the j isotope value of the kth endmember (k = 1, 2, 3. . . ), µ jk is the mean, ω jk is the standard deviation; P k is the contribution ratio of endmember k, which needs to be predicted according to the stable isotope analysis in R model; C jk is the fractionation factor of endmember k on the j isotope, λ jk is the mean of the fractionation factor, τ jk is the standard deviation; ε jk is the residual error, representing the undetermined variable between different individual mixtures, the mean is 0, and the standard deviation is σ i . The R Programming Language was used to calculate the contribution ratio of each pollution source by running the package of Stable Isotope Analysis in R.
Therefore, the steps of the methodology in study are shown in Figure 2.

Physicochemical Characteristics
All sample water chemistry and isotopic test results are shown in Table 1. The DO concentration of groundwater in the study area was in the range of 2.6-14.8 mg·L −1 with an average value of 6.8 mg·L −1 , the DO concentration of fish pond water was in the range of 2.2-10.1 mg·L −1 with an average value of 6.5 mg·L −1 , the DO concentration in surface river water ranged from 6.3 to 14.2 mg·L −1 with an average value of 10.0 mg·L −1 . It showed that the three water bodies in the study area were all in an aerobic environment. Compared with groundwater and fish pond water, surface water has a higher DO concentration, which was mainly because surface river water has better fluidity and openness, and was closely related to atmospheric exchange.
The Cl − concentration in groundwater was at a low concentration level (1.1-7.4 mg·L −1 ), with an average value of 5.4 mg·L −1 ; the concentration range in fish ponds was 3.4-18.6 mg·L −1 , with an average value of 7.6 mg·L −1 ; except for the higher concentrations of SW02, SW05, and SW06 (147.2 1 , 198.8, and 104.4 mg·L −1 , respectively), the surface river water was generally at a low level. The concentration level was 2.0-24.4 mg·L −1 , with an average value of 9.3 mg·L −1 . There were two main reasons for the increase in Cl − concentration at the SW02, SW05, and SW06 sampling points: firstly, the three sampling points are located in residential cultivated areas and are subject to the use of artificial chlorine-containing fertilizers. It has been confirmed that the concentration is in direct proportion; the second is the dissolution of salt rock. It has been confirmed that the dissolution of rock salt is one of the main sources of Cl − in karst water [35]. Through statistical analysis of data, it was concluded that the order of Cl − concentration in three water bodies was surface river water > fish pond water > groundwater. This was mainly because the use of man-made chlorine-containing fertilizers increased the concentration of Cl − into surface river water and fish ponds, while the concentration of Cl − into groundwater is relatively small. The DIC concentration in groundwater was in the range of 18.8-35.2 mg·L −1 , with an average value of 27.3 mg·L −1 ; the concentration range in fish ponds was 31.6-55.0 mg·L −1 , with an average value of 47.9 mg·L −1 ; the concentration in river water was generally was in the range of 17.8-51.8 mg·L −1 , with an average value of 35.7 mg·L −1 . It was concluded that the order of DIC concentration in three water bodies was fish pond water > surface river water > groundwater.

Nitrogen Characteristics
According to the concentration data of NO 3 − , NO 2 − , and NH 4 + , the concentration box of three nitrogen elements in water is shown in Figure 3. The maximum concentrations of NO 3 − and NH 4 + in underground water were 49.3 and 0.6 mg·L −1 , which were 2.5 and 1.2 times the limit.
Among the six fish pond samples, a higher concentration of NH 4 + was detected. The maximum mass concentration of NH 4 + in fish pond water was 8.2 mg·L −1 , which was 16.3 times the limit. The maximum concentrations of NO 3 − , NO 2 − , and NH 4 + were 24.5, 7.3, and 8.6 mg·L −1 in surface river water, respectively, which were 1.2, 7.3, and 17.3 times the limit.
To sum up, the order of the average concentration of nitrate was groundwater > surface river water > fish pond water, and the order of the average concentration of nitrite and ammonium ion was fish pond > surface river water > groundwater. In the three types of water bodies, the nitrogen concentrations in some sampling sites exceeded the limits of Class III water in the groundwater quality Standard (GB/T14848-2017) and the Sanitary Standard for drinking Water (GB5749-2006). Among them, the maximum concentration of NO 3 − and NH 4 + in groundwater exceeded the limit standard, the concentration of NH 4 + in fish pond water exceeded the limit by nearly one order of magnitude, and the concentration of NO 3 − , NO 2 − , and NH 4 + in surface river water exceeded the limit standard, among which the concentration of NH 4 + has exceeded the limit standard by more than 10 times. In the study area, the long-term use of these three types of water by humans and organisms will pose potential health hazards. The over-standard points of water samples were mainly located in residential, cultivated land, and orchards. On the one hand, the topography was flat, the hydrodynamic condition was poor, and the concentration of pollutants was easy to gather, which led to the increase in concentration. On the other hand, the orchard and cultivated land used a large number of chemical fertilizers, resulting in the increase in NO 3 − through infiltration and surface runoff. The high concentrations of NO 2 − and NH 4 + were mainly in surface river water and fish pond water, while the concentration of NO 3 − was relatively low. According to the two stages of nitrification, it can be concluded that the first stage of nitrification mainly occurs in fish pond water and surface river water (the initial stage of nitrification), while the second stage (nitrification) has not occurred obviously, indicating that fish pond water and surface river water were polluted by nitrogen. Fishermen may have fed too much nitrogenous feed to fish and nearby villagers used nitrogen fertilizer on fruit trees [36]. Surface river water and fish pond water play an important role in the transformation of the three waters. As important water resources in the study area, relevant departments should strengthen prevention and control. Among them, the water area of fish ponds accounts for about 10% of the water area of Huixian wetlands. During the high and low water periods, the polluted fish pond water seeps into the groundwater, which seriously affects the quality of groundwater, harms the health, and hinders the sustainable development of the ecosystem. Fishermen should reduce the use of nitrogenous feed, feed it reasonably, and apply fertilizer scientifically in agricultural activities to ensure that the water in fish ponds is not polluted.
It can be concluded from Figure 4 that, in three different types of water bodies, the order of δ 15 N-NO 3 was groundwater > surface river water > fish pond water; the order of δ 18 O-NO 3 was groundwater > fish pond water > surface river water; and the order of δ 13 C DIC was fish pond water > surface river water > groundwater.

Qualitative Analysis of the Source of Nitrate
The characteristics of the concentration of Cl − and TN can identify the source of groundwater pollution [37]. According to the concentration of Cl − and TN of three kinds of water samples in the study area, the relation diagram was drawn ( Figure 5). The maximum concentration of TN among samples in the study area did not exceed 60 mg·L −1 and did not exceed the concentration limit of groundwater seriously polluted by human/animal excreta [38]. However, the concentrations of Cl − and TN in some samples were higher than 7 and 10 mg·L −1 , respectively, and the concentrations exceed the limits of Cl − and TN concentrations of man-made pollution sources of groundwater in karst areas [39,40], indicating that the three types of water bodies in the study area were affected by human activities, and the sources of nitrogen and Cl − were diversified, including nitrogen fertilizer, mixed source, and manure sewage. Individual samples had the characteristics of high Cl − and low TN, and Cl − may come from chlorine-containing fertilizers, such as KCl used by crops in the watershed.
In order to further identify the source of NO 3 − , the relationship between δ 15 N-NO 3~N O 3 − and δ 15 N-NO 3~δ 18 O-NO 3 was interpreted [23,41]. It is pointed out that when the water body has the characteristics of high δ 15 N-NO 3 value and high NO 3 − concentration, NO 3 − comes from manure, when there is low δ 15 N-NO 3 value and high NO 3− concentration, NO 3 − comes from chemical fertilizer, and when there is medium δ 15 N-NO 3 value and medium NO 3 − concentration, NO 3 − comes from domestic sewage [42]. The variation ranges of δ 15 N-NO 3 and NO 3 − in different water bodies are large ( Figure 6), indicating that they did not have the characteristics of unique sources, and the sources of NO 3 − in water bodies were complex, which was consistent with the sources of nitrogen interpreted by the characteristics of the concentration of Cl − and TN.
The physical, chemical, and biological processes of nitrogen in nature will lead to the fractionation of nitrogen isotopes. The main processes of nitrogen isotope fractionation are ammonia fixation, assimilation, mineralization, nitrification, denitrification, and ammonia volatilization. Among them, the fractionation effects of nitrification, denitrification, and assimilation are obvious, and the fractionation coefficients are between −12 and 29‰, −40-5‰, and −27-0‰, respectively [43,44]. However, the fractionation effects caused by biogeochemical processes of nitrogen and oxygen are also different under different conditions, and the fractionation effects of different forms of nitrogen in the process of assimilation are also different. Conditions such as light and dissolved oxygen can also affect the isotope fractionation effect by affecting biogeochemical reactions [45]. Because of the overlap of δ 15 N-NO 3 value distribution in atmospheric nitrogen deposition and soil, chemical fertilizer, and organic fertilizer, the source of nitrogen cannot be accurately determined by nitrogen isotopes alone. Therefore, the source of nitrogen can be determined more accurately by using δ 15 N and δ 18 O bistable isotope tracer technique [46]. Previous studies have shown that the δ 15 N and δ 18 O values of nitrate produced by inorganic chemical fertilizers are −6-+6‰ and +17-+25‰, respectively [47]; the δ 15 N and δ 18 O values of nitrate produced by animal manure are +5-+25‰ and +5-+7‰, respectively [48]; the δ 15 N and δ 18 O values of nitrate produced by domestic sewage are +4-+19‰ and −5-+10‰, respectively [49]; the δ 15 N and δ 18 O values of nitrate produced by organic nitrogen mineralization in soil are 0-+8‰ and −10-+10‰, respectively [50]. According to the relationship between δ 15 N-NO 3 and δ 18 O-NO 3 , it can be seen from Figure 7 that almost all the groundwater, fish pond water, and surface river water samples in the study area fall within the range of ammonium nitrogen fertilizer, soil nitrogen, and manure sewage, indicating that ammonium nitrogen fertilizer, soil nitrogen, and manure sewage were the main input sources of NO 3 − in the study area, and the effects of rainfall and nitrate fertilizer were not obvious. In addition, it can also be seen from the distribution areas of different types of water points that groundwater was more distributed in animal manure and domestic sewage areas, fish pond water was more in chemical fertilizer areas, and surface river water was more in soil nitrogen, animal manure, and domestic sewage areas. Generally speaking, the groundwater, fish pond water, and surface river water in the study area were affected by chemical fertilizer, animal manure, and domestic sewage.

Nitrification
The nitrification process is mainly completed by two stages: the first stage is nitrosation that is NH 4 + oxidation to NO 2 − (Equation (3)), and the second stage is nitrification that is NO 2 − oxidation to NO 3 − (Equation (4)). Nitrification will lead to an increase in NO 3 − in groundwater (Equation (5)). In the process of nitrification, nitrifying bacteria absorb NH 4 + , and release protons (H + ) in order to maintain the charge balance, and the pH of groundwater decreases. The acid (H + ) produced during nitrification will accelerate the weathering of carbonate rocks (Equation (4)), making more old carbon from carbonate rocks enter into groundwater, resulting in a deviation of δ 13 C DIC value in groundwater [12]. If nitrification occurs in the karst aquifer, there should be a negative correlation between δ 15 N-NO 3 and δ 13 C DIC and a positive correlation between NO 3 − and δ 13 C DIC in groundwater.
There was a significant negative correlation between δ 15 N-NO 3 and δ 13 C DIC in the three types of water bodies in the study area (Figure 8a), which proved that nitrification was the main controlling factor of the nitrogen cycle in the whole basin. The positive correlation between NO 3 − and δ 13 C DIC in the three types of water bodies was not very obvious (Figure 8b), but the overall indication was that nitrification occurred obviously.   In the NO 3 − produced by nitrification, about 2/3 of the oxygen comes from the surrounding water, and the remaining 1/3 of the oxygen comes from the oxygen dissolved in the water (23.5‰) [12]. Therefore, the nitrification process in water can be indicated by the measured δ 18 O-NO 3 value [51]: The δ 18 O-NO 3 values of most sampling sites in the study area ( Figure 7) fall in the range of typical δ 18 O-NO 3 values produced by nitrification (−10-10%) [12], which shows that nitrification in study area played a leading role in the source of nitrate. Only one point in groundwater has δ 18 O-NO 3 values greater than 10‰, which may be caused by higher δ 18 O-H 2 O atmospheric precipitation input, isotope fractionation, or denitrification of δ 18 O-H 2 O in the surrounding water during nitrification.

Denitrification
The process of denitrification causes significant isotope fractionation, which leads to the enrichment of heavy δ 15 N-NO 3 and δ 18 O-NO 3 in groundwater residues. The ratio of δ 15 N-NO 3 to δ 18 O-NO 3 is usually between 1.3 and 2.1 [52][53][54]. As can be seen from Figure 9, only two spots in all samples were within the range of denitrification and the fitting formula of δ 15 N-NO 3 and δ 18 O-NO 3 was y = 0.857x + 3.26. In addition, it is generally believed that denitrification will occur when the concentration of DO is less than 2 mg·L −1 [44], while the range of DO of each sampling in the study area was 2.24~14.80 mg·L −1 , and the average value was greater than 6 mg·L −1 . It was an aerobic environment, which also showed that significant denitrification does not easily occur. It showed that nitrification is the dominant process of the whole study area.

Quantitative Analysis of Nitrate Sources
The previous analysis showed that nitrification dominated the nitrogen cycle in the water body, and denitrification was not obvious. The analysis confirmed that the main pollution of wetland water nitrate came from ammonia fertilizer, soil nitrogen, animal manure, and domestic sewage. The stable isotope analysis in R model can be used to quantitatively analyze the contribution ratio of different NO 3 − inputs in the water body. Combined with the hydrogeological conditions and human activities in the study area, the typical isotopic values of ammonium nitrogen fertilizer, soil nitrogen, and manure sewage were taken from the literature [5], and the fractionation factor C jk was taken as 0.
According to the calculation results of stable isotope analysis in R model (Figure 10), the order of the contribution of different pollution sources to NO 3 − in the study area was manure sewage > ammonium nitrogen fertilizer > soil nitrogen, the average contribution rate of animal manure and domestic sewage was (39.8 ± 6.4)%, the average contribution rate of ammonia fertilizer was (33.4 ± 5.1)%, while soil nitrogen was relatively low, the average contribution rate was (26.8 ± 3.9)%. The main reason was that the Huixian wetland mainly consisted of cultivated land, mainly vegetable cultivation and aquaculture, producing a large amount of excess fertilizer and manure, which flowed into the surrounding surface water through rainfall and seeped into the ground through leaching; thus affecting the quality of surface water and groundwater. Figure 10. The contribution ratio of three types of pollution sources to water bodies in the study area.

1.
The maximum concentration of NO 3 − and NH 4 + in groundwater exceeded the limit standard by 2.5 and 1.2 times, the maximum concentration of NH 4 + in fish pond water exceeded the limit by 16.3 times, and the maximum concentration of NO 3 − , NO 2 − , and NH 4 + in surface river water exceeded the limit standard 1.2, 7.3, and 17.3 times, respectively. It showed that groundwater, fish pond, and surface river water were polluted by nitrogen to varying degrees. The δ 15 N-NO 3 values in the water bodies of the study area ranged from −15.7 to 20.8‰, with an average of 5‰; the δ 18 O-NO 3 values ranged from −6.9 to 14.3‰, with an average of 1.7‰; the δ 13 C DIC values ranged from −13.5 to 0.1‰, with an average of −8.3‰.The concentrations of Cl − and TN in some samples of three kinds of water bodies in the study area were higher than 7 and 10 mg·L −1 , respectively, indicating that the three types of bodies in the study area were affected by human activities. It suggested ammonium nitrogen fertilizer, soil nitrogen, and manure sewage were the main input sources of NO 3 − in the study area by the relationship of Cl −~T N, the relationship of δ 15 N-NO 3~N O 3 − , and the relationship of δ 15 N-NO 3~δ 18 O-NO 3 .

2.
The significant negative correlation between δ 15 N-NO 3 and δ 13 C DIC proved that nitrification was the main controlling factor of the nitrogen cycle in the study area. Meanwhile, the δ 18 O-NO 3 values of most sampling sites in the study area fall in the range of typical δ 18 O-NO 3 values produced by nitrification (−10-10%), which also showed that nitrification played a leading role in the source of nitrate. Only two points in groundwater have δ 18 O-NO 3 values greater than 10‰, which may be caused by higher δ 18 O-H 2 O atmospheric precipitation input, isotope fractionation, or denitrification of δ 18 O-H 2 O in the surrounding water during nitrification. In addition, the range of dissolved oxygen at each sampling point in the study area was 2.3~14.8 mg·L −1 , and the average value was greater than 6 mg·L −1 . It was an aerobic environment. Therefore, nitrification mainly occurred in the process of nitrate conversion in the study area, and significant denitrification did not occur. 3.
The order of the contribution of NO 3 − sources to the water body in the study area: the contribution rate of manure and sewage was (39.8 ± 6.4)%, the contribution rate of ammonium nitrogen fertilizer was (33.4 ± 5.1)%, and the contribution rate of soil nitrogen was (26.8 ± 3.9)%. In the prevention and control of nitrogen pollution in Huixian wetland, we should focus on three kinds of pollution sources: manure sewage, ammonium nitrogen fertilizer, and soil nitrogen, and reduce the discharge of fecal sewage and apply fertilizer scientifically and quantitatively, so as to reduce the supply of nitrogen sources in wetlands in the study area and cause good results for the control of nitrogen pollution in the study area.