Next Article in Journal
A Review on the Prevalence of Arcobacter in Aquatic Environments
Previous Article in Journal
SpillwayPro: Integrated Water Surface Profile, Cavitation, and Aerated Flow Analysis for Smooth and Stepped Chutes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrochemical Evolution and Quality Assessment of Groundwater in the Sanjiang Plain, China

1
Key Laboratory of Groundwater Resources and Environment, Jilin University, Ministry of Education, Changchun 130021, China
2
Jilin Provincial Key Laboratory of Water Resources and Environment, Jilin University, Changchun 130021, China
3
College of New Energy and Environment, Jilin University, Changchun 130021, China
*
Author to whom correspondence should be addressed.
Water 2022, 14(8), 1265; https://doi.org/10.3390/w14081265
Submission received: 4 March 2022 / Revised: 19 March 2022 / Accepted: 8 April 2022 / Published: 13 April 2022
(This article belongs to the Section Hydrogeology)

Abstract

:
Groundwater is subjected to contamination threats from human activities, such as agriculture, especially long-term farming in the Sanjiang Plain, China. Identifying the sources and distribution of pollution is essential for its reasonable prevention and control. In this study, we analysed the chemical characteristics of 389 samples at 60 shallow groundwater monitoring points from 2011 to 2015 in the Sanjiang Plain using traditional hydrochemical methods, water quality assessment, Pearson’s correlation, and principal component analysis (PCA). Although groundwater type in this area was predominantly HCO3-Ca·Mg, three forms of nitrogen pollution (ammonia, nitrate, and nitrite) were all detected in this area. The interaction of natural geochemical and anthropogenic factors during hydrochemical formation is confirmed by the high coefficients of variation and Gibbs plots of the main ions in the water. The overall shallow groundwater situation was described as good, with more than 40% and 90% of groundwater samples suitable for drinking and irrigation according to the quality assessment, respectively. The proportion of poor water quality in the wet season was higher than that in the dry season. NO3-N and NH3-N were identified as the major anthropogenic pollutants in the study area. Results from Pearson’s correlation and principal component analysis shows two main pollutants fall into two chemical controlling factors together with natural chemical parameters, which implies that the migration and transformation of pollutants may have affected the overall hydrochemical characteristics of the regional groundwater. Therefore, findings from this paper can provide insight into the chemical evolution of groundwater in response to long-term agricultural activities and can help contribute to better management of groundwater resources and agricultural sustainable development.

1. Introduction

The Sanjiang Plain has the largest area of natural marsh wetland distribution in China. Due to the special climatic conditions and geographical location, its soil has a remarkably high humus content, making it one of the four most precious native black soil belts globally [1]. Fertile land has also greatly promoted agricultural development in this area. Large wetlands areas have developed into farmlands since large-scale reclamation in the 1950s. In 2005, the arable land area was approximately 55,688 km2 more than tripled from 17,134 km2 in 1954 [2,3]. From 2010 to 2018, the arable land area increased by 3828.17 km2. This area has become an important commodity grain base in China, with a per capita grain output three times higher than the national average [4]. However, long-term agricultural development has had a significant influence on the ecological environment. Specifically, the use of chemical fertilizers in agriculture has been increasing yearly, from less than 0.75 g/m2 in the 1950s to more than 15 g/m2 after 2000s [5]. The excessive accumulation of chemical fertilisers changed the soil’s properties and seriously damages the health of the underground ecosystem [6,7].
Groundwater is an important carrier of underground environment materials and energy, and its chemical composition reflects the underground environment characteristics. Soil pollution caused by long-term agricultural activities indirectly generates signals on the groundwater chemical composition [8]. Sun et al. declare that the spatial variability of nitrate concentration in groundwater, total soil nitrogen, and shallow soil soluble nitrate-nitrogen have a moderately strong spatial autocorrelation, indicating that the vertical migration of pollutants is the primary cause of groundwater pollution [9]. Moreover, the relationship between nitrate concentration in groundwater and land use dominated by agriculture has been widely discussed. By comparing nitrate concentrations in shallow groundwater at different locations, Nolan et al. found that the most polluted wells were under farmland in the United States [10]. Yu et al. calculated that the average nitrogen concentration of rice fields and other agricultural land was about 5.1 mg/L and found that intensive agricultural land is the main cause of groundwater nitrogen pollution [11]. Therefore, high nitrogen concentrations in groundwater have been regarded as typical pollutants induced by agricultural activities. Based on the traceability analysis of stable isotopes of NO3-δN, Cui et al. and Xu Chunying et al. verified that the main source of nitrate pollution in groundwater is chemical fertilisers [12,13]. In recent years, the rapid development of agriculture in the Sanjiang Plain has already caused groundwater pollution. The continuous decline of groundwater quality has led to water shortages in some farms and residential areas in the region [14]. Large-scale changes in land use are believed to cause nitrogen pollution in groundwater. Many studies have shown that high nitrate concentrations were found in shallow groundwater in this area [6,15,16]. However, groundwater is a dynamic mobile system, and nitrate is a highly migrated component in water. To identify and evaluate the impact of human activities such as agriculture on an underground environment, an in-depth multi-components analysis and discussion of groundwater chemical evolution over several years is necessary [17,18,19].
Hydrogeochemical methods are used to explore the evolution of groundwater quality and to analyse the origin of groundwater [20,21,22]. Generally, traditional hydrogeochemical methods can display an obvious overview of groundwater chemical characteristics but often fail to reflect the subtle and critical environmental and geological phenomena in local areas [23,24]. Besides, multivariate statistical methods are nowadays widely applied into groundwater chemistry research in Bangladesh and the Middle East, such as factor analysis, cluster analysis, and principal component analysis (PCA) [25,26,27,28]. Based on multivariate statistical methods, suitability of groundwater quality were assessed for irrigation and drinking, and the chemical origins of groundwater can be identified as well [29,30].
Therefore, based on multi-year monitoring data of shallow groundwater chemistry in the Sanjiang Plain, this study: (1) presents the chemical characteristics and evolution of shallow groundwater (within 60 m below the surface); (2) identifies pollution factors in the chemical composition of groundwater meant for drinking and irrigation; and (3) applies multivariate statistical methods to identify the agricultural impacts on groundwater chemistry throughout the study area. This research will provide more insight into the influence of agriculture on groundwater chemical evolution and provide a theoretical basis for the rational development, utilisation, and protection of groundwater resources in regional agricultural water resources.

2. Materials and Methods

2.1. Study Area

The study area falls in the flooded area of the Songhua River and Naoli River basins in the Sanjiang Plain, northeast China. It covers an area of 26,164 km2, located between the latitudes 129°5731″ N and 133°5644″ N, and longitudes 46°637″ E to 47°4543″ E (Figure 1). With a temperate continental climate, the average annual precipitation and potential evaporation are 440–650 mm and 550–840 mm, respectively. The terrain gradually decreased from southwest to northeast. Black soil covers an area of 11,101 km2, accounting for 72.8% of the total farmland area. The aquifer media distributed in the floodplain of the Songhua River, Naoli River, and its main tributaries comprise Holocene fine silt sand, fine sand, and gravel with a general thickness of 10–25 m. The types of groundwater mainly include floodplain pore water in valleys, terrace pore water-weakly confined water, and micro-pore fissure water in mountain terraces. The main recharge source of groundwater is the infiltration of atmospheric precipitation, and discharge from mainly lateral runoff, artificial exploitation, and evaporation. The groundwater depth along the Songhua River and the Naoli River basins is shallower than 2.5 m and gradually increases with the distance from the rivers, whereas the depth could increase to more than 10 m downstream of the Naoli River because of artificial exploitation. The agriculture irrigation area covers an area of 16,173 km2, accounting for about 62% of the total area of the Sanjiang Plain.

2.2. Data Collection

Groundwater quality data were collected from the Heilongjiang Provincial Water Conservancy and Hydroelectric Power Investigation, Design, and Research Institute in China. However, considering the comparability of different source data, only the records of several regional shallow groundwater investigations were used in this study, which included: (1) 33 samples in May 2011 (wet season, 2011-W); (2) 33 samples in October 2011 (dry season, 2011-D); (3) 52 samples in May 2012 (wet season, 2012-W); (4) 52 samples in October 2012 (dry season, 2012-D); (5) 51 samples in May 2013 (wet season, 2013-W); (6) 51 samples in October 2013 (dry season, 2013-D); (7) 57 samples in May 2015 (wet season, 2015-W); and (8) 60 samples in October 2015 (dry season, 2015-D). In summary, 389 samples from shallow groundwater from 2011–2015 were included.
These samples were analysed with water chemical analysis. Among the water index, general indicators (pH, EC, total dissolved solids (TDS), and total hardness (TH); concentrations of major ions (Na+, K+, Ca2+, Mg2+, Cl, SO42−, HCO3, CO32−), and regional enriched metal elements (Fe, Mn) were used to analyse the evolution of regional groundwater characteristics. In addition, the nitrogen pollution index (including NH3-N, NO3-N, and NO2-N) was used to provide general information of groundwater quality contaminated by agricultural activities.

2.3. Traditional Hydrochemical Analysis

Traditional hydrogeochemical methods are widely applied to explore groundwater quality evolution and groundwater recharge sources [22]. First, Piper Trilinear diagram of major cations and anions were drawn to analyse the chemical components of groundwater samples. A Gibbs plot was then employed to identify the sources of water compositions. It can directly reflect the dominant factors, such as precipitation, rock weathering, or evaporation, and provide a qualitative reference for their impacts on water chemistry [31].

2.4. Water Quality Index

The water quality index (WQI) was computed according to drinking requirements to identify the pollution components in groundwater and its suitability to drinking purposes [32]. Five steps were involved in the WQI calculation: (i) according to the World Health Organization Water Quality Standard (WHO2011), 11 available parameters including K+, Na+, TH, Cl, SO42−, TDS, pH, NO3-N, NO2-N, Fe, and Mn were selected for calculation; (ii) each parameter was given a weight (wi) in terms of its relative importance in the overall water quality (Table 1); (iii) the relative weights (Wi) of each parameter involved were calculated based on (1); (iv) quality rating scale (qi) for each parameter was calculated following (2); (v) WQI is the weighted sum of all the parameters’ quality rating scale (3). The computed WQI values were categorised into five classes: excellent, good, poor, very poor, and unsuitable for drinking.
W i = w i / i = 1 n w i
q i = C i / S i * 100
WQI = W i · q i
where Wi is the relative weight, wi is the weight of each parameter, n is the number of parameters, qi is the quality rating scale, Ci is the concentration of each chemical parameter in each sample in mg/L, and Si is the World Health Organization standard for each chemical parameter in mg/L according to the guidelines of the WHO (2011).

2.5. Risk Assessment for Irrigation

With the expansion of farming areas, groundwater exploitation has become a primary way of agricultural irrigation due to its availability and low cost. However, poor quality irrigation water adversely affects farmland soil. Therefore, to evaluate the hazard risk of groundwater used to irrigate black soil, the following irrigational quality parameters were estimated: percent sodium (%Na) (4), sodium adsorption ratio (SAR) [33] (5), and the permeability index (PI) (6) were used, which were calculated as follows (all ionic concentrations are in meq/L):
% Na = Na + + K + Ca 2 + + Mg 2 + + Na + + K + × 100
SAR = Na + Ca 2 + + Mg 2 + 2
PI = Na + + HCO 3 Ca 2 + + Mg 2 + × 100
The salinity hazard was determined by the total concentration of soluble salts in water. It was classified into four levels: low (EC < 250 μs/cm), medium (250–750 μs/cm), high (750–2250 μs/cm), and very high (2250–5000 μs/cm) [34]. Sodium hazard is expressed in terms of the percent sodium (Na%) and sodium adsorption ratio (SAR). Water was classified as low (<10), medium (10–18), high (18–26), and very high (>26) sodium hazards based on the sodium adsorption ratio [34,35]. In addition, the classification of water-based on percent sodium was classified as excellent (<20%), good (20–40%), permissible (40–60%), doubtful (60–80%), and unsuitable (>80%) [36]. Lastly, water was classified based on PI into three classes: class I (>75%), class II (25–75%), and class III (<25%). (WHO1989)

2.6. Statistical Analysis

Descriptive statistics provide a general view of the hydrochemical characteristics [37]. In addition, interrelationships in a set of chemical variables were determined using factor analysis. The underlying but unobservable or latent processes which the raw dataset cannot directly reveal can be more clearly interpreted by factor analysis [38]. In this study, the R-mode was used for factor analysis to identify the hydrochemical controlling factors involved. Hypothesis testing showed that TDS was significantly associated with major ions (p < 0.01), and ammonia nitrogen was significantly associated with heavy metals. The dimensionality of the parameters was then reduced by PCA to reveal the association patterns between the hydrochemical variables. The measurements of pH and the concentrations of TDS, TH, and the major ions in the groundwater were analysed to identify the main factors controlling the chemical components of groundwater [39,40].

3. Results and Discussion

3.1. Chemical Characteristics of Groundwater

3.1.1. General Hydrochemistry

The statistical analysis of groundwater hydrochemical parameters helps understand water chemical composition’s enrichment and change law. Table 2 shows the maximum, minimum, average, variance, and coefficient of variation for each parameter. The results indicate that 98.47% of the samples were associated with neutral water (pH 6.2–8.3) except for individual extreme points of pH. As an important indicator of water palatability, TDS is between 109 and 2950 mg/L, and 51.57% of the samples exceeded the acceptable limit for drinking (500 mg/L), indicating about half of the water samples was distorted in taste due to high TDS. The TH of the groundwater reflects the lithological formation of the aquifer. Generally, water is classified based on TH value as soft water (<150 mg/L) and hard water (>150 mg/L). Among the samples, 76.68% were hard water, and 23.32% were soft water. The TH value was between 40–2430 mg/L, with an average of 284 mg/L. According to the statistical analysis, the cation abundance in the water samples was Ca2+ > Mg2+ > Na+ > K+. The calcium and magnesium ions concentrations ranged from 9.94–579 mg/L and 2.31–238 mg/L, with average values of 74.62 mg/L and 23.79 mg/L, respectively. A total of 2.16% of the sample’s calcium ion concentration was relatively high, exceeding the World Health Organization guideline value (200 mg/L). The sodium and potassium ions concentrations ranged from 4.88–179 mg/L and 0.29–218 mg/L, with the average values being 35.78 mg/L and 5.89 mg/L, respectively. The order of the anion concentration was HCO3 > Cl > SO42−, the observed bicarbonate and chlorine concentrations ranged from 22.8–1430 and 0.19–1570 mg/L, and the average values were 171.89 and 69.70 mg/L, respectively. Nitrogen pollution is mainly caused by ammonia nitrogen, nitrate, and nitrite pollution. The NH3-N, NO3-N, and NO2-N concentrations in groundwater samples were 0–4.51, 0–148, and 0–4.2 mg/L, with average values of 0.36, 24.03, and 0.04 mg/L, respectively (Table 2).
The degree of dispersion of the samples was a necessary reference for classification. The coefficients of variation of K, NO2-N, Mn, and Fe were >200%. The wetland part of the Sanjiang Plain experiences slow sedimentation leading to the accumulation of iron in the underlying groundwater and sediments. In addition, the complexation of natural organic compounds can prevent or delay the precipitation of oxyhydroxides. Therefore, the swamp system in the Sanjiang Plain is considered to play a vital role in species control and iron output [41]. The coefficients of variation of Cl, SO4, NO3-N, and NH3-N were >100%, indicating that the groundwater may have been affected by human activities. The coefficients of variation of Na+, TH, TDS, Ca2+, Mg2+, HCO3, and EC are >60%, and the index values had certain discrete distribution characteristics, which may be indicators that exceed the standard.

3.1.2. Hydrochemical Facies

Groundwater hydrochemical facies shows the chemical composition categories of groundwater, and reflects the variability of hydrochemical parameters in response to chemical processes during transport in the lithologic framework [42,43].
Piper diagram can provide a comprehensive view of the hydrochemical phases and types of groundwater [44,45,46]. According to the Piper diagram (Figure 2), the cation plot shows that most samples are concentrated in Zone A, indicating that the groundwater is mainly of the “calcium” type. About a quarter of the samples fall in the “no dominant’’ type (zone B) and “sodium” type (zone D) due to the increasing proportion of sodium in water. In the anion plot, sample points are more dispersedly distributed in zones E and B, suggesting the dominance of bicarbonate or mixed anions. The central diamond field shows the overall features of groundwater chemistry [45]. As shown in Figure 2, almost all groundwater samples are plotted in zones I and IV, illustrating that HCO3-Ca∙Mg and Cl∙ SO4-Ca∙Mg were the dominant water types. Only a few samples occurred in zones II and III. Similar results were also reported by Lu et al. [47] for the Abujiao River Basin, which is part of the Sanjiang Plain. Additionally, samples in zone I were mainly characterised by freshwater (TDS < 1000 mg/L), and samples in zone 4 (Cl-Ca type) and near zone 5 were all high in TDS.
According to the distribution of varied hydrochemical facies over the region (Figure 1), we observed that the pore phreatic water and weak confined water along the river valley and terrace were mostly of HCO3-Ca·Mg type. On the contrary, the chemical type of micro-pore fissure phreatic water in the mountain front of the evaluation area was mostly of HCO3·SO4-Ca·Mg type, and gradually changed into HCO3·Cl-Ca·Mg type in the lower reaches. A high density of Cl·HCO3 or Cl·SO4 water mainly appeared in the vicinity of Jiamusi, and the proportion of Na among the cations was also significantly increased.

3.1.3. Mechanisms Controlling Groundwater Chemistry

The mechanisms governing the major dissolved chemical ions constituents were summarized into three dominant processes of precipitation, rock weathering, and evaporation (Gibbs, 1970). Two diagrams of the weight ratios of Na+/(Na+ + Ca2+) and Cl/(Cl + HCO3) (x-axis) against TDS (y-axis) were presented by Gibbs to identify the functional sources of dissolved chemical compositions. The Gibbs diagram (Gibbs, 1970) was primordially used to represent the mechanism controlling the evolution of the world’s surface water, and it has been widely used in groundwater research.
According to the Gibbs diagram (Figure 3), 99.13% of the sampling points in the relationship chart between Na+/(Na+ + Ca2+) and TDS fell within the rock weathering zone, while 96.07% was observed for Cl/(Cl + HCO3). This indicates that the cation content is mainly affected by rock weathering, and anion content is controlled by rock weathering and evaporation. Moreover, the “bicarbonate type” samples were mainly plotted in the rock dominance zone, suggesting that rock weathering and leaching are the major mechanisms controlling groundwater chemistry in these areas [48]. The other types of distributions all showed a slight increasing trend towards the evaporative dominance band. Some samples still fell outside the Gibbs field, which may have been affected by anthropogenic sources [49]. However, anthropogenic activities, especially irrigation, can also impact evaporation [50] since large-scale irrigation increases the groundwater level. Our results indicate that rock weathering is the most important factor influencing groundwater evolution; evaporation and anthropogenic activities may contribute to water chemistry.

3.2. Groundwater Quality Assessment

3.2.1. Groundwater Quality for Drinking

Water Quality Assessment Index (WQI) is an aggregated communication tool for monitoring water quality [51], developed to summarise a huge amount of water quality data in a format that is easy to present and understand, with less information than the original data. In this study, WQI was employed to evaluate groundwater quality status for drinking water purposes. According to the WHO (2011) standard, the distribution weight (Wi) and relative weight (wi) of each parameter are listed in Table 1, and the calculated WQI values of samples in this area ranges from 5.65–3268.19. The distinguished classifications based on WQI scores are: (I) <50, excellent water; (II) 50–100, good water; (III) 100–200, poor water; (IV) 200–300, very poor water; and (V) >300, unfit for drinking water. As shown in Figure 4a, the proportion of poor-quality water in the wet season was higher than that in the dry season, and the water quality deteriorated first and then improved. The over-standard ratio was 27.03% in 2011, increasing to 59.65% in 2013 and dropping to 23.73% in 2015. The proportion of available drinking water in each period was higher than 40%, and the proportion of poor-quality water was higher than 23.73%.
Furthermore, statistics of each chemical parameter compared to those in the WHO standards are presented in Figure 4b. Five chemical parameters, distinguished by more than 40% of the samples, were above the recommended WHO limit in the following order: TH, TDS, NO3-, Mn, and Fe. Hardness in water is mainly caused by magnesium and calcium [52], the main cations in the groundwater of this region. Thus, the TH and TDS of more than 50% of the water samples exceeded those of the WHO’s drinking water standards. This high TH and TDS can be mainly attributed to the dissolution of soluble salts and minerals, and evaporation because most of the study area is located downstream of the watershed, as well as wastewater infiltration. NO3- is one of the main pollutants in groundwater in the world’s cultural regions [53] In the study area, the nitrate concentration of 49.76% of water samples were beyond the permissible limits of WHO standards, and are associated increased with soil cover nitrate fertilisers from agricultural practices. In addition, NO3-N, NO2-N, and NH3-N in water can be interconverted under certain favourable conditions [6]. Although only 0.72% water samples had NO2 content above the permissible limit, the over-limit ratio of NH3-N is up to 25.48%, according to the Chinese drinking water standard. However, groundwater’s Fe concentrations in this area reached high levels because of their geochemical nature [54]. Therefore, NO3-N and NH3-N were identified as the major anthropogenic pollutants in the study area.

3.2.2. Groundwater Quality for Irrigation

In this basin, groundwater is mainly extracted for drinking and irrigation, while water quality requirements vary for different purposes. To assess the suitability of groundwater for irrigation purposes, three indices relevant to sodium were used in the study: %Na, SAR, and PI.
The %Na values of the study area varied from 9.78–67.72%, with a mean value of 29.89%. In terms of the classification suggested by Wilcox (1955), 94.69% of the samples indicated %Na values in good and permissible limits, whereas 3.86% and 1.45% of samples were doubtful and unsuitable for irrigation purposes respectively (Figure 5a). SAR is the main factors that determine the suitability of water for agricultural purposes. Naseh et al. evaluated the agricultural water suitability of the Ghaen Plain according to SAR. [55]. As shown in Figure 5b, approximately 97.10% of the samples with low SAR pose almost no risk of exchangeable sodium and are of excellent water quality for irrigation. Only 2.90% of samples of high and very high sodium water (in C4S1 and C5S1) were regarded as unfavourable, as they contained detrimental levels of exchangeable sodium in soils. Figure 5c presents the PI value distribution of the samples that ranges from 13.42–112.17, with an average of 52.87. According to the PI value, approximately 69.31% of the groundwater samples belong to Category I, indicating that the water was good or moderate for irrigation purposes. Of the groundwater samples, 26.73% belonged to Type II (barely suitable), and the remaining (3.96%) belonged to Type III (not suitable for irrigation). Therefore, assessment of all the indices shows that most groundwater in this area is safe for irrigation and pose no risk of soil damage. By contrast, the proportions in each quality category assessment varied assessment with different indices because these indices considered different reactions of the sodium in water with soil, according to the equations of each index.

3.3. Statistical Analysis of Hydrochemistry

3.3.1. Ion Correlation Matrix Analysis

Pearson’s correlation expressed the relationships among the different parameters [56]. According to the correlation coefficient (R), the relevance degree can be divided into three levels: strong (R > 0.7), moderate (0.7 > R > 0.3), and weak (R < 0.3) [57].
Figure 6 shows the correlation curves between all pairs of the 15 parameters. In the study area, TDS exhibited a strong positive relationship with TH (R = 0.93), Ca2+ (R = 0.92), Mg2+ (R = 0.84), Cl (R = 0.77), and Na+ (R = 0.71), which is a common understanding of natural geochemistry of groundwater in alluvial plains. Regarding nitrogen as a typical agricultural pollutant, there was no strong correlation among the three forms of nitrogen (NO3-N, NO2-N, and NH3-N), whereas the NO3-N concentration in groundwater was proportional to TDS. The moderate positive correlation between NH3-N and Fe and Mn suggests that reduction conditions are one of the reasons for the high NH3-N concentration. The weak correlation between NO2-N and other parameters was due to its low concentration.
The annual average histogram of each parameter is shown in Figure 5. During 2011–2015, TDS and the concentration of its closely related ions fluctuated in a certain value range, whereas NO3-N and NO2-N concentration initially increased and then decreased after implementing the “three reductions” in agriculture since 2013 [58]. However, NH3-N did not show an obvious decrease during this period.

3.3.2. Principal Component Analysis

PCA was used to track the contribution of natural and human factors to the groundwater hydrochemistry components. Twenty-five indicators were analyzed by PCA for dimensionality reduction in the southern Tehran aquifer [59]. For the Sanjiang Plain, four principal components (PCs) with eigenvalues above 1 were retained, which explained 71.75% of the total variance (Figure 7). PC1 with the greatest variance (40.66%) decides the hydrochemical characteristics and is affected by positive loadings of TDS (0.964), TH (0.961), Ca2+ (0.948), Mg2+ (0.902), Cl (0.845), Na+ (0.748), SO42− (0.668), and NO3-N (0.645) (Table 3). The loadings on PC1 indicate that multiple processes, including carbonate weathering, anthropogenic factors, and evaporite dissolution, may be the main sources of these ions [60]. PC2 accounted for 12.90% of the total variance, with strong loadings for NH3-N (0.811), Mn (0.673), and Fe (0.528). The positive loadings of these three parameters imply an anoxic and reducing environment and the influence of human activities because of the inclusion of NH3-N. PC3 accounted for 10.09% of the total variance and mainly comprised pH (0.644) and HCO3 (0.640). The loadings on PC3 further indicate that mineral weathering affects water chemistry. PC4 explained 8.10% of the variance and had positive NO2-N (0.529) and K+ (0.52) loadings. The high positive loading of NO2-N indicates the contribution of anthropogenic factors on water chemistry.

4. Conclusions

The Sanjiang Plain plays a pivotal role in food production in China. However, large-scale reclamation and agricultural development, human activities, and excessive fertiliser application have led to the deterioration of the groundwater environment in the Sanjiang Plain, and the nitrogen content of the groundwater also shows complex changes. Therefore, chemical investigation of the groundwater was carried out to rationalise the development of groundwater resource utilisation and protection. The results are as follows:
(1)
In the study area, the cation abundance in shallow groundwater was Ca2+ > Na+ > Mg2+ > K+, and the order of the anion concentration was HCO3 > Cl > SO42−. Three forms of nitrogen pollution, ammonia, nitrate, and nitrite were detected in this area. In addition, the coefficients of variation of Cl, SO4, NO3-N, and NH3-N were higher than 1, indicating that the groundwater may have been affected by human activities.
(2)
There are four types of hydrogeochemical groundwater, and the main HCO3-type water accounted for 63.09% of the groundwater samples. The hydrochemical type, indicated by the Piper diagram, transits from HCO3–Ca·Mg or HCO3·SO4–Ca·Mg to HCO3·Cl–Ca·Mg type along the flow path in nature, but some abnormal points should be the result of human interference.
(3)
Water assessment by WQI and irrigation index showed that the proportion of available drinking water in each period was higher than 40%, and more than 90% of the groundwater in the study area was suitable for irrigation. In addition, nitrogen was identified as a typical anthropogenic pollutant in this area. Comparing the seasonal and annual variations in groundwater quality, the proportion of poor water quality in the wet season was higher than that in the dry season, and the water quality deteriorated first, but improved during the final year of the 2011–2015 period.
(4)
Although groundwater chemistry is still of rock dominance, it is increasingly influenced by agricultural activities. Three of the four principal components extracted by PCA are relevant to anthropogenic factors, explaining 71.75% of the variance in groundwater chemistry, and the first two PCs show the interaction of both natural and anthropogenic factors. The typical pollutants of the three nitrogen forms (NO3, NH3, and NO2) were determined by PC1, PC2, and PC4, respectively.

Author Contributions

Y.L.: Conceptualization, Methodology, Writing—review & editing. X.Y.: Resources, Data curation, Formal analysis, Supervision, Validation. Y.Z.: Writing—original draft, Visualization. X.D.: Funding acquisition, Project administration. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China, grant number 41972247.

Acknowledgments

The completion of this article is inseparable from the contributions of all the authors. Their support is gratefully acknowledged. I would also like to thank the anonymous reviewers for their constructive and valuable suggestions in bringing the paper to its present form.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Han, X.; Li, N. Research progress and prospect of black land in Northeast China. J. Geogr. Sci. 2018, 38, 10. [Google Scholar] [CrossRef]
  2. Liu, J.; Sheng, L.; Lu, X.; Liu, Y. A dynamic change map of marshes in the Small Sanjiang Plain, Heilongjiang, China, from 1955 to 2005. Wetl. Ecol. Manag. 2015, 23, 419–437. [Google Scholar] [CrossRef]
  3. Song, K.; Liu, D.; Wang, Z.; Zhang, B.; Jin, C.; Li, F.; Liu, H. Land use changes and driving forces in the Sanjiang Plain since 1954. Acta Geogr. Sinica 2008, 63, 93–104. [Google Scholar] [CrossRef]
  4. Shu, L.; Wang, Z.; Yuan, Y.; Zhang, F.; Liu, W.; Lu, C. Land use changes in typical areas of the Sanjiang Plain in the past 40 years and their effects on groundwater. J. Hydraul. Eng. 2021, 52, 11. [Google Scholar] [CrossRef]
  5. Zhang, F.; Qi, F.; Yang, X.; Wang, P.; Wang, C. Groundwater Environmental Quality Assessment of Sanjiang Plain in Heilongjiang Province; National Seminar on Groundwater Pollution; China Geological Survey: Beijing, China, 2008. [Google Scholar]
  6. Cao, Y.; Tang, C.; Song, X.; Liu, C.; Zhang, Y. Characteristics of nitrate in major rivers and aquifers of the Sanjiang Plain, China. J. Environ. Monit. 2012, 14, 2624–2633. [Google Scholar] [CrossRef]
  7. Wang, X.; Zhang, G.; Xu, Y.J.; Sun, G. Identifying the regional-scale groundwater-surface water interaction on the Sanjiang Plain, Northeast China. Environ. Sci. Pollut. Res. 2015, 22, 16951–16961. [Google Scholar] [CrossRef]
  8. Li, T.T.; Zhong, M.S.; Jiang, L.; Fan, Y.L.; Yao, J.J.; Xia, T.X.; Jia, X.Y. Comparison of two methods for assessing impact of contaminated soil on groundwater quality. Res. Environ. Sci. 2013, 26, 793–799. [Google Scholar]
  9. Sun, H.-Y.; Wei, X.-F.; Jia, F.-C.; Li, D.-J.; Li, J.; Li, X.; Yin, Z.-Q. Source of Groundwater Nitrate in Luanping Basin Based on Multi-environment Media Nitrogen Cycle and Isotopes. Huanjing Kexue Environ. Sci. 2020, 41, 4936–4947. [Google Scholar] [CrossRef]
  10. Nolan, B.T.; Stoner, J.D. Nutrients in Groundwaters of the Conterminous United States, 1992−1995. Environ. Sci. Technol. 2000, 34, 1156–1165. [Google Scholar] [CrossRef] [Green Version]
  11. Yu, Z.-Q.; Nakagawa, K.; Berndtsson, R.; Hiraoka, T.; Suzuki, Y. Groundwater nitrogen response to regional land-use management in South Japan. Environ. Earth Sci. 2021, 80, 634. [Google Scholar] [CrossRef]
  12. Choi, W.-J.; Han, G.-H.; Lee, S.-M.; Lee, G.-T.; Yoon, K.-S.; Choi, S.-M.; Ro, H.-M. Impact of land-use types on nitrate concentration and δ15N in unconfined groundwater in rural areas of Korea. Agric. Ecosyst. Environ. 2007, 120, 259–268. [Google Scholar] [CrossRef]
  13. Xu, C.; Li, Y.; Li, Q.; Wang, L.; Dong, Y.; Jia, X. Current status of groundwater nitrate pollution in Weifang, Shandong and its δ~(15)n traceability. Acta Ecol. Sin. 2011, 31, 78–83. [Google Scholar] [CrossRef]
  14. Zhao, W.; Shang, Z.; Cui, F. Problem existing in agricultural development and utilizaton of groundwater in the Sanjing Plain and countermeasures for water resources protection. J. Heilongjiang Water Conserv. Sci. Technol. 2010, 2, 145–146. [Google Scholar] [CrossRef]
  15. Cao, Y.; Tang, C.; Song, X.; Liu, C.; Zhang, Y. Residence time as a key for comprehensive assessment of the relationship between changing land use and nitrates in regional groundwater systems. Environ. Sci. Process. Impacts 2013, 15, 876–885. [Google Scholar] [CrossRef]
  16. Lu, L.; Cheng, H.; Pu, X.; Liu, X.; Cheng, Q. Nitrate behaviors and source apportionment in an aquatic system from a watershed with intensive agricultural activities. Environ. Sci. Process. Impacts 2015, 17, 131–144. [Google Scholar] [CrossRef]
  17. Fang, J.-J.; Zhou, A.-G.; Ma, C.-M.; Liu, C.-F.; Cai, H.-S.; Gan, Y.-Q.; Liu, Y.-D. Evaluation of nitrate source in groundwater of southern part of North China Plain based on multi-isotope. J. Central South Univ. 2015, 22, 610–618. [Google Scholar] [CrossRef]
  18. Guo, Q.; Zhou, Z.; Huang, G.; Dou, Z. Variations of Groundwater Quality in the Multi-Layered Aquifer System near the Luanhe River, China. Sustainability 2019, 11, 994. [Google Scholar] [CrossRef] [Green Version]
  19. Zhang, Y.; Lee, D.; Ding, J.; Lu, J. Environmental Impact of High Concentration Nitrate Migration in Soil System Using HYDRUS Simulation. Int. J. Environ. Res. Public Health 2020, 17, 3147. [Google Scholar] [CrossRef]
  20. Bu, J.; Sun, Z.; Ma, R.; Liu, Y.; Gong, X.; Pan, Z.; Wei, W. Shallow Groundwater Quality and Its Controlling Factors in the Su-Xi-Chang Region, Eastern China. Int. J. Environ. Res. Public Health 2020, 17, 1267. [Google Scholar] [CrossRef] [Green Version]
  21. Roy, P.D.; Selvam, S.; Gopinath, S.; Logesh, N.; Sánchez-Zavala, J.L.; Lakshumanan, C. Geochemical evolution and seasonality of groundwater recharge at water-scarce southeast margin of the Chihuahuan Desert in Mexico. Environ. Res. 2022, 203, 111847. [Google Scholar] [CrossRef] [PubMed]
  22. Wu, J.; Li, P.; Qian, H. Hydrochemical characterization of drinking groundwater with special reference to fluoride in an arid area of China and the control of aquifer leakage on its concentrations. Environ. Earth Sci. 2015, 73, 8575–8588. [Google Scholar] [CrossRef]
  23. Gao, Y.; Qian, H.; Ren, W.; Wang, H.; Liu, F.; Yang, F. Hydrogeochemical characterization and quality assessment of groundwater based on integrated-weight water quality index in a concentrated urban area. J. Clean. Prod. 2020, 260, 121006. [Google Scholar] [CrossRef]
  24. Redwan, M.; Moneim, A.A.A.; Amra, M.A. Effect of water–rock interaction processes on the hydrogeochemistry of groundwater west of Sohag area, Egypt. Arab. J. Geosci. 2016, 9, 111. [Google Scholar] [CrossRef]
  25. Howladar, M.F.; Rahman, M. Characterization of underground tunnel water hydrochemical system and uses through multivariate statistical methods: A case study from Maddhapara Granite Mine, Dinajpur, Bangladesh. Environ. Earth Sci. 2016, 75, 1501. [Google Scholar] [CrossRef]
  26. Chen, M.; Wu, Y.; Gao, D.; Chang, M. Identification of coal mine water-bursting source using multivariate statistical analysis and tracing test. Arab. J. Geosci. 2017, 10, 28. [Google Scholar] [CrossRef]
  27. Ferati, F.; Kerolli-Mustafa, M.; Kraja-Ylli, A. Assessment of heavy metal contamination in water and sediments of Trepça and Sitnica rivers, Kosovo, using pollution indicators and multivariate cluster analysis. Environ. Monit. Assess. 2015, 187, 338. [Google Scholar] [CrossRef] [PubMed]
  28. Sunkari, E.D.; Seidu, J.; Ewusi, A. Hydrogeochemical evolution and assessment of groundwater quality in the Togo and Dahomeyan aquifers, Greater Accra Region, Ghana. Environ. Res. 2022, 208, 112679. [Google Scholar] [CrossRef]
  29. Khelif, S.; Boudoukha, A. Multivariate statistical characterization of groundwater quality in Fesdis, East of Algeria. J. Water Land Dev. 2018, 37, 65–74. [Google Scholar] [CrossRef] [Green Version]
  30. Rahman, M.A.T.M.T.; Saadat, A.H.M.; Islam, S.; Al-Mansur, A.; Ahmed, S. Groundwater characterization and selection of suitable water type for irrigation in the western region of Bangladesh. Appl. Water Sci. 2017, 7, 233–243. [Google Scholar] [CrossRef] [Green Version]
  31. Gibbs, R.J. Mechanisms Controlling World Water Chemistry. Science 1970, 170, 1088–1090. [Google Scholar] [CrossRef] [PubMed]
  32. World Health Organization. Guidelines for Drinking-Water Quality, 4th ed.; World Health Organization: Geneva, Switzerland, 2011. [Google Scholar]
  33. Alrajhi, A.; Beecham, S.; Bolan, N.S.; Hassanli, A. Evaluation of soil chemical properties irrigated with recycled wastewater under partial root-zone drying irrigation for sustainable tomato production. Agric. Water Manag. 2015, 161, 127–135. [Google Scholar] [CrossRef]
  34. Richards, L.A. Determination of the properties of saline and alkali soils. In Diagnosis and Improvement of Saline and Alkali Soils, Agriculture Handbook No. 60; United States Department of Agriculture: Washington, DC, USA, 1954; pp. 7–53. [Google Scholar]
  35. Wilcox, L. Classification and Use of Irrigation Waters; United States Department of Agriculture: Washington, DC, USA, 1955; p. 969. [Google Scholar]
  36. Nishanthiny, S.C.; Thushyanthy, M.; Barathithasan, T.; Saravanan, S. Irrigation water quality based on hydro chemical analysis, Jaffna, Sri Lanka. Am. J. Agric. Environ. Sci. 2010, 7, 100–102. [Google Scholar]
  37. He, X.-L.; Wu, Y.-H.; Zhou, J.; Bing, H.-J. Hydro-chemical Characteristics and Quality Assessment of Surface Water in Gongga Mountain Region. Environ. Sci. 2016, 37, 3798–3805. [Google Scholar]
  38. Wang, F. Factor Analysis and Principal-Components Analysis. Int. Encycl. Hum. Geogr. 2009, 12, 1–7. [Google Scholar] [CrossRef]
  39. Rao, N.S.; Rao, J.P.; Subrahmanyam, A. Principal component analysis in groundwater quality in a developing urban area of Andhra Pradesh. J. Geol. Soc. India. 2007, 69, 959–969. [Google Scholar]
  40. Wu, T.-N.; Su, C.-S. Application of Principal Component Analysis and Clustering to Spatial Allocation of Groundwater Contamination. In Proceedings of the 2008 Fifth International Conference on Fuzzy Systems and Knowledge Discovery, Jinan, China, 18–20 October 2008; Volume 4, pp. 236–240. [Google Scholar]
  41. Zou, Y.; Jiang, M.; Yu, X.; Lu, X.; David, J.L.; Wu, H. Distribution and biological cycle of iron in freshwater peatlands of Sanjiang Plain, Northeast China. Geoderma 2011, 164, 238–248. [Google Scholar] [CrossRef]
  42. Kshetrimayum, K. Hydrochemical evaluation of shallow groundwater aquifers: A case study from a semiarid Himalayan foothill river basin, northwest India. Environ. Earth Sci. 2015, 74, 7187–7200. [Google Scholar] [CrossRef]
  43. Qian, J.; Peng, Y.; Zhao, W.; Ma, L.; He, X.; Lu, Y. Hydrochemical processes and evolution of karst groundwater in the northeastern Huaibei Plain, China. Appl. Hydrogeol. 2018, 26, 1721–1729. [Google Scholar] [CrossRef]
  44. Khalili, Z.; Asadi, N. Groundwater quality assessment using Aq.QA and GIS technique in Azna city, Lorestan province, Iran. Sustain. Water Resour. Manag. 2021, 7, 21. [Google Scholar] [CrossRef]
  45. Pant, R.R.; Zhang, F.; Rehman, F.U.; Wang, G.; Ye, M.; Zeng, C.; Tang, H. Spatiotemporal variations of hydrogeochemistry and its controlling factors in the Gandaki River Basin, Central Himalaya Nepal. Sci. Total Environ. 2018, 622, 770–782. [Google Scholar] [CrossRef] [PubMed]
  46. Piper, A.M. A graphic procedure in the geochemical interpretation of water-analyses. Eos Trans. Am. Geophys. Union 1944, 25, 914–928. [Google Scholar] [CrossRef]
  47. Lu, L.U.; Dai, E.; Cheng, Q. The sources and fate of nitrogen in groundwater under different land use types: Stable isotope combined with a hydrochemical approach. Acta Geogr. Sin. 2019, 74, 1878–1889. [Google Scholar]
  48. Wang, B.; Lee, X.-Q.; Yuan, H.-L.; Zhou, H.; Cheng, H.-G.; Cheng, J.-Z.; Zhou, Z.-H.; Xing, Y.; Fang, B.; Zhang, L.-K.; et al. Distinct patterns of chemical weathering in the drainage basins of the Huanghe and Xijiang River, China: Evidence from chemical and Sr-isotopic compositions. J. Southeast Asian Earth Sci. 2012, 59, 219–230. [Google Scholar] [CrossRef]
  49. Selvakumar, S.; Chandrasekar, N.; Kumar, G. Hydrogeochemical characteristics and groundwater contamination in the rapid urban development areas of Coimbatore, India. Water Resour. Ind. 2017, 17, 26–33. [Google Scholar] [CrossRef]
  50. Ravikumar, P.; Somashekar, R. Environmental Tritium (3H) and hydrochemical investigations to evaluate groundwater in Varahi and Markandeya River basins, Karnataka, India. J. Environ. Radioact. 2011, 102, 153–162. [Google Scholar] [CrossRef] [Green Version]
  51. Vasanthavigar, M.; Srinivasamoorthy, K.; Vijayaragavan, K.; Ganthi, R.R.; Chidambaram, S.; Anandhan, P.; Manivannan, R.; Vasudevan, S. Application of water quality index for groundwater quality assessment: Thirumanimuttar sub-basin, Tamilnadu, India. Environ. Monit. Assess. 2010, 171, 595–609. [Google Scholar] [CrossRef]
  52. Adimalla, N.; Li, P.; Venkatayogi, S. Hydrogeochemical Evaluation of Groundwater Quality for Drinking and Irrigation Purposes and Integrated Interpretation with Water Quality Index Studies. Environ. Process. 2018, 5, 363–383. [Google Scholar] [CrossRef]
  53. Adimalla, N.; Li, P. Occurrence, health risks, and geochemical mechanisms of fluoride and nitrate in groundwater of the rock-dominant semi-arid region, Telangana State, India. Hum. Ecol. Risk Assess. Int. J. 2019, 25, 81–103. [Google Scholar] [CrossRef]
  54. Chi, G.; Zhu, B.; Huang, B.; Chen, X.; Shi, Y. Spatiotemporal dynamics in soil iron affected by wetland conversion on Sanjiang Plain. Land Degrad. Dev. 2021, 32, 4669–4679. [Google Scholar] [CrossRef]
  55. Naseh, M.R.V.; Noori, R.; Berndtsson, R.; Adamowski, J.; Sadatipour, E. Groundwater Pollution Sources Apportionment in the Ghaen Plain, Iran. Int. J. Environ. Res. Public Health 2018, 15, 172. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Sreedevi, P.D.; Sreekanth, P.D.; Ahmed, S.; Reddy, D.V. Appraisal of groundwater quality in a crystalline aquifer: A chemometric approach. Arab. J. Geosci. 2018, 11, 211. [Google Scholar] [CrossRef]
  57. Emenike, C.P.; Tenebe, I.T.; Jarvis, P. Fluoride contamination in groundwater sources in Southwestern Nigeria: Assessment using multivariate statistical approach and human health risk. Ecotoxicol. Environ. Saf. 2018, 156, 391–402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Lu, H.; Jing, Y.; Zhang, H. Beidahuang’s “three cuts” led to the rapid development of green and organic agriculture. China State Farm 2015, 10, 39. [Google Scholar]
  59. Ghahremanzadeh, H.; Noori, R.; Baghvand, A.; Nasrabadi, T. Evaluating the main sources of groundwater pollution in the southern Tehran aquifer using principal component factor analysis. Environ. Geochem. Health 2017, 40, 1317–1328. [Google Scholar] [CrossRef]
  60. Elumalai, V.; Nwabisa, D.P.; Rajmohan, N. Evaluation of high fluoride contaminated fractured rock aquifer in South Africa–Geochemical and chemometric approaches. Chemosphere 2019, 235, 1–11. [Google Scholar] [CrossRef]
Figure 1. Study area and location of water sampling points.
Figure 1. Study area and location of water sampling points.
Water 14 01265 g001
Figure 2. Piper diagram for the groundwater.
Figure 2. Piper diagram for the groundwater.
Water 14 01265 g002
Figure 3. Gibbs diagrams showing the mechanism controlling groundwater chemistry of (a) TDS vs. cations and (b) TDS vs. anions.
Figure 3. Gibbs diagrams showing the mechanism controlling groundwater chemistry of (a) TDS vs. cations and (b) TDS vs. anions.
Water 14 01265 g003
Figure 4. Grades histogram of groundwater samples for drinking purposes in seasons (a) and parameters (b).
Figure 4. Grades histogram of groundwater samples for drinking purposes in seasons (a) and parameters (b).
Water 14 01265 g004
Figure 5. Irrigations hazards classifications by sodium percentage (a), sodium adsorption ratio, SAR (b), and permeability index, PI (c).
Figure 5. Irrigations hazards classifications by sodium percentage (a), sodium adsorption ratio, SAR (b), and permeability index, PI (c).
Water 14 01265 g005
Figure 6. Relationships between different chemical parameters and histogram of each parameter’s annual average change.
Figure 6. Relationships between different chemical parameters and histogram of each parameter’s annual average change.
Water 14 01265 g006
Figure 7. Principal components load diagrams.
Figure 7. Principal components load diagrams.
Water 14 01265 g007
Table 1. The relative chemical weight of the physicochemical parameters (all values are in mg/L, except for pH).
Table 1. The relative chemical weight of the physicochemical parameters (all values are in mg/L, except for pH).
Chemical ParametersWHO Standards (2011)Weight (wi)Relative Weight (Wi)
K20020.048
Na20020.048
TH (as CaCO3)100–20020.048
Cl25030.071
SO425040.095
TDS50050.119
pH (on scale)6.5–8.540.095
NO35050.119
NO2350.119
Mn0.150.119
Fe0.350.119
Total 421
Table 2. Maximum, minimum and average chemical composition and variance of groundwater samples (all values are in mg/L, except for pH and EC).
Table 2. Maximum, minimum and average chemical composition and variance of groundwater samples (all values are in mg/L, except for pH and EC).
maxminavesdcv
pH8.36.26.830.315%
EC (μs/cm)457017192064370%
TDS295010957637064%
Na1794.8835.7826.5174%
K2180.295.8914.76250%
Ca5799.9474.626790%
Mg2382.3123.7921.4190%
HCO3143022.8171.89113.4566%
SO44680.2165.9367.74103%
Cl15700.1969.7126.51181%
TH24304028424887%
NO3-N148024.0330.78128%
NH3-N4.5100.360.52143%
NO2-N4.200.040.26603%
Fe79.802.468.26336%
Mn10.300.621.28205%
Table 3. Element matrix of PCA.
Table 3. Element matrix of PCA.
PC1PC2PC3PC4
Cl0.845−0.031−0.08−0.385
SO40.668−0.027−0.240.362
HCO30.330.460.640.078
K0.2960.1430.4740.52
Na0.748−0.123−0.0570.365
Ca0.9480.008−0.005−0.211
Mg0.9020.010.152−0.223
TH0.9610.0050.051−0.222
TDS0.9640.0370.0780.022
pH−0.1810.2920.644−0.143
NH3-N0.0330.811−0.2420.011
NO3-N0.645−0.417−0.1690.092
NO2-N0.1410.189−0.2430.529
Fe−0.1540.528−0.324−0.261
Mn0.350.673−0.3390.068
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ye, X.; Zhou, Y.; Lu, Y.; Du, X. Hydrochemical Evolution and Quality Assessment of Groundwater in the Sanjiang Plain, China. Water 2022, 14, 1265. https://doi.org/10.3390/w14081265

AMA Style

Ye X, Zhou Y, Lu Y, Du X. Hydrochemical Evolution and Quality Assessment of Groundwater in the Sanjiang Plain, China. Water. 2022; 14(8):1265. https://doi.org/10.3390/w14081265

Chicago/Turabian Style

Ye, Xueyan, Yan Zhou, Ying Lu, and Xinqiang Du. 2022. "Hydrochemical Evolution and Quality Assessment of Groundwater in the Sanjiang Plain, China" Water 14, no. 8: 1265. https://doi.org/10.3390/w14081265

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop