Next Article in Journal
Improved Planning of Energy Recovery in Water Systems Using a New Analytic Approach to PAT Performance Curves
Previous Article in Journal
Microbial Assisted Hexavalent Chromium Removal in Bioelectrochemical Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Possible Hydrochemical Processes Influencing Dissolved Solids in Surface Water and Groundwater of the Kaidu River Basin, Northwest China

1
Key Laboratory of Water Cycle and Related Land Surface Processes, Institute of Geographic Science and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
College of Geography and Environmental Science, Hainan Normal University, Haikou 571158, China
*
Author to whom correspondence should be addressed.
Water 2020, 12(2), 467; https://doi.org/10.3390/w12020467
Submission received: 2 December 2019 / Revised: 4 February 2020 / Accepted: 8 February 2020 / Published: 10 February 2020
(This article belongs to the Section Hydrology)

Abstract

:
Hydrochemical processes under intense human activities were explored on the basis of the hydrochemical characteristics of 109 surface water samples and 129 groundwater samples collected during August 2015 to September 2016, in the Kaidu River Basin. Results obtained in this study indicated that the water in the basin was neutral to slightly alkaline with low total dissolved solids. Rock weathering and evaporation controlled the natural hydrochemical mechanisms. Mountain groundwater and stream water were dominated by Ca2+-HCO3 type water, whereas the plains groundwater was dominated by mixed type water. The results of principal component analysis demonstrated that water-rock interaction and human activity explained 71.6% and 12.9% of surface water hydrochemical variations, respectively, and 75.1% and 14.2% of groundwater hydrochemical variations, respectively. Sulfate, chloride, and carbonate weathering were the major water-rock interaction processes. Livestock farming and agricultural activities were the primary human activities influencing the water hydrochemistry. In addition, cation exchange is another important process influencing the hydrochemical characteristics in the study area. This study would be helpful in forecasting of water quality in arid areas.

1. Introduction

Arid and semiarid areas have always suffered from water resource shortage. With increasing population and developing economy, overexploitation of water resources and water pollution have been increasing in severity, and therefore water resource shortage has become more serious [1]. Water quality in arid basins is influenced by a various of hydrochemical processes, such as water-rock interaction, land reclamation, wastewater infiltration, sewage exfiltration, and agricultural irrigation [2,3,4]. Due to the complex influences of multiple natural and anthropogenic processes in arid areas, identifying the major hydrochemical processes is difficult but highly useful [5,6].
Various studies have been conducted to evaluate the effect of natural and anthropogenic processes on hydrochemistry in different environments. Wang et al. [7] studied the hydrochemical characteristics of ground ice in permafrost regions on the Qinghai-Tibet plateau in China; the results indicated that soil moisture, air temperature, and the thickness of the active layer are the major factors controlling variations of soil water chemistry. In the Wiesent River Basin, southern Germany, bedrock geology remains the dominant controlling factor of the major ionic chemistry and agricultural influences are the strongest near the headwaters [8]. Human activities including groundwater exploitation, agricultural activities, industrial emission and even the construction of Three Gorges Dam were factors that led to the increase in pH, NO3-N, and Cl in groundwater from 1992 to 2010, in the Jianghan Plain in China [9]. Shi et al. [3] assessed the effect of seawater intrusion, precipitation infiltration, granite weathering, anthropogenic pollution, and concrete material dissolution to shallow groundwater in a highly urbanized coastal city, Shenzhen, China. The results of principal component analysis (PCA) demonstrated that water-rock interactions and concrete material dissolution are the most important factors affecting shallow groundwater hydrochemistry. Li et al. [10] claimed that halite and carbonate dissolution with strong cation exchange are the major sources of the high salinity in groundwater in the Northwest Namibia. These studies discussed hydrochemical processes under different environments, as well as the effect of these processes on water quality. However, in arid areas, scholars have focused more on water quantity rather than quality [11,12,13,14]. In fact, hydrochemical processes, particularly water quality, should receive the same level of attention in arid areas, as they are related to water resource usability.
The Kaidu River is the major recharge source of the Bosten Lake, which is the largest inland freshwater lake in China (with an area of 1210.5 km2 and storage of 73.03 × 108 km3 when the lake elevation is 1047 m above sea level (a.s.l.)). The river is the major water resource for life and economic development in the Bayingolin Mongol Autonomous Prefecture (Bazhou), as well as for the ecological construction and environmental protection in the Yanqi Basin, the Konqi River Basin, and the lower reaches of the Tarim River Basin. Security of water resources, in the Kaidu River Basin, is of great importance for regional sustainable development. Hydrochemical processes are key factors that influence the security of the water quality. Moreover, in recent years, with the rapid development of irrigation agriculture and animal husbandry, a large amount of residue of pesticide and chemical fertilizer is discharged into soil and water, threatening the security of water resources.
The hydrochemical data of surface water and groundwater, in the Kaidu River Basin, are presented in this study, with the objectives of exploring the predominant mechanisms and hydrochemical processes that control the hydrochemistry of both surface water and groundwater under the influence of intense human activities. This work should be useful for the understanding of water quality changes in arid areas, which is necessary for regional water resource management.

2. Study Site

The Kaidu River Basin is located at the southern slope of the Tianshan Mountains and the northern rim of the Tarim River Basin (Figure 1). It covers an area of 44,147 km2, with an altitude ranging from 928 m a.s.l. to 4796 m a.s.l. The Kaidu River originates from glaciers in the Tianshan Mountains, runs through the small Yourdusi Basin, the large Yourdusi Basin and, then, along the mountain valley to the Dashankou hydrometric station, through the Yanqi Basin, and finally flows into the Bosten Lake. It is one of the largest rivers on the southern slope of the Tianshans (with a mean annual discharge of 108.42 m3/s) [15].
The study area is situated in central Eurasia, controlled by a typical temperate continental climate. Climate characteristics differ significantly between the mountains and piedmont plains. The average temperatures are −4.2 °C and 8.9 °C with average precipitation of 280.5 mm and 84.3 mm, in the mountains and the plains, respectively. Approximately 86% and 80% of annual precipitation in the mountains and the plains, respectively, is concentrated in the May to September period.
The Kaidu River Basin is located at the southern Tianshan geosynclinal folded zone; folds and faults have developed under this influence. The North-West-West (NWW) and Nort-East-East (NEE) trending faults are the major two groups of faults; the NWW trending fault is in accordance with the main tectonic line. The exposure strata are mainly Cenozoic-Quaternary sediments in the small Yourdusi Basin, large Yourdusi Basin, and Yanqi Basin; Paleozoic sediments in the alpine valley region; and Proterozoic era sediments in the valley area at the east of the small Yourdusi Basin. The outcropped lithologies in the Kaidu River Basin are mainly glutenite, argillaceous sandstone, mudstone, clay, and silt, which are rich in carbonate, sulfate, and halite [16,17]. The groundwater resource is richer in the mountainous areas than in the plains.
The annual river runoff at the mountain outlet station ranges from 24.61 × 10 8   m 3 to 57.26 × 10 8   m 3 , with a mean value of 35.32 × 10 8   m 3 . The maximum runoff occurred in July and August. The spatial and temporal variations of groundwater depth in the downstream of the Kaidu River Basin are very significant. The groundwater depth increases with the distance from the Kaidu River. The annual mean groundwater level ranges from 3 m to more than 28 m below the ground. In recent years, the groundwater level has decreased significantly, averaging a 5 m to 14 m decrease. Because of over pumping for irrigation, the groundwater level is the shallowest in winter and deepest during the growing season (June to August) and during the autumn irrigation season (October). In regions with large farmland areas, the groundwater level is deeper than in other areas.
The major land uses in the basins are grasslands (accounting for 61.60% of the total study area), deserts (27.72%), and farmland (6.07%) (Figure 1). In the mountains, the major soil types are tundra soil, meadow soil, steppe soil, chestnut soil, and brown desert soil. In the plains, the major soil types are moisture soil, brown desert soil, saline soil, and irrigation desert soil [16,17].

3. Data and Methods

3.1. Sampling

In order to avoid the effects of weather, clear weather (zero precipitation) for at least three consecutive days was required before sampling. From August 2015 to September 2016, five systematic sampling campaigns in four seasons (August 2015, January 2016, April 2016, August 2016, and September 2016) were conducted (Figure 1). The meteorological seasons are referred to as spring (March, April, and May), summer (June, July, and August), autumn (September, October, and November), and winter (December, January, and February). A total of 106 stream water samples, 121 groundwater samples, and 8 mountain groundwater samples were collected.
The stream water samples were collected from the middle of the river’s cross-section. Groundwater samples were collected from farm wells or domestic water wells along the river. Human activities have the greatest influence on shallow groundwater. Therefore, we only collected shallow groundwater from well depths less than 50 m. The well depth was recorded; however, if the wells were closed, we did not measure the groundwater table. The wells were purged at least half an hour before sampling. All water samples were collected in three 500 mL clean and dry high-density polyethylene plastic bottles. Prior to sampling, the bottles were soaked for 24 h in deionized water in the laboratory. After soaking, the bottles were rinsed three times with deionized water and, then, dried at room temperature. All bottles were washed three times again using the sample water in situ before sampling and, then, filled with sample water and sealed immediately using parafilm to prevent contamination. After collection, samples were transported immediately via a cooling box.
Analyses were conducted in the Central Laboratory of Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences immediately after returning from the field. The pH was measured using a pH meter (PHJS-4A, Inesa Inc., Shanghai, China). Major cationic and anionic analyses were conducted by ion chromatography systems (Dionex ICS-5000, Sunnyvale, CA, USA). The total dissolved solids (TDS) was measured using a weighing method, with a precision within 1 μg/L.

3.2. Saturation Indices (SIs)

As an important parameter in geochemistry and hydrogeology, the saturation index is used to identify whether a water sample tends to dissolve or precipitate a certain kind of mineral [18]. In this study, we calculated the saturation indices (SIs) of anhydrite, aragonite, calcite, dolomite, gypsum, and halite using the following equation [19]:
S I = log ( I A P k s ( T ) )
where IAP is the relevant activity product of the solution, k s ( T ) is the equilibrium constant of the reaction considered at the sample temperature (°C). The SIs were calculated using the geochemical computer model PHREEQC [20]. When the water is saturated with the dissolved mineral, the SI equals zero; positive values indicate oversaturation, while negative values indicate undersaturation [21,22].

3.3. Cation Exchange

The chloro-alkaline indices (CAI-I and CAI-II), proposed by Schoeller [23], are frequently used to explore the cation exchange between water and its host environment during water migration. The CAI-I and CAI-II are calculated by the following equations:
CAI I = Cl ( Na + + K + ) Cl
CAI II = Cl ( Na + + K + ) HCO 3 + SO 4 2 + CO 3 2
Negative Schoeller indices indicate Ca2+ or Mg2+ exchange in solution with Na+ or K+, i.e., Ca2+ or Mg2+ are removed from solution and Na+ or K+ are released into it. The unit of the ion content is denoted in meq/L. Negative values indicate chloro-alkaline disequilibrium. This reaction is known as the cation exchange reaction; the host rocks are the primary sources of dissolved solids in the water during this process, while the positive Scholler indices indicate the reverse ion exchange.

4. Results

4.1. Water Hydrochemistry

Water in the basin was neutral to slightly alkaline with the pH varying from 6.82 to 8.30 (mean 7.95). The TDS values (Figure 2a) increased along the water flow direction, indicating more soluble matter downstream of the basin [24,25]. The mean TDS ranged from 130.92 mg/L (mountain groundwater) to 422.58 mg/L (plain groundwater) (Table 1). It is lower than the drinking water standard of WHO [26,27].
The mean ion concentrations were the lowest in mountain groundwater and the highest in plain groundwater. HCO3 and Ca2+ were the dominant ions; the mean concentration of HCO3 ranges from 116.33 to 203.85 mg/L. The mean concentration of Ca2+ ranges from 36.09 to 51.56 mg/L, while the mean concentration of SO42− ranges from 11.92 to 118.93 mg/L. The mean concentration of Mg2+ ranges from 7.14 to 27.46 mg/L, whereas those of NO3 and K+ were the lowest. The mean concentration of NO3 ranges from 0.43 to 2.25 mg/L, while the mean concentration of K+ ranges from 1.37 to 3.81 mg/L.
The contributions of different ions to total ions varied with water types. In stream water, Ca2+ accounted for 65% of total cations (based on meq L−1 values), followed by Mg2+ (22%), Na+ (12%), and K+ (1%), and anions by HCO3(71%) > SO42− (21%) > Cl (8%). They were very similar to the major ions of the mountain groundwater, with cations as Ca2+ (63%) > Mg2+ (19%) > Na+ (17%) > K+ (1%) and anions as HCO3 (76%) > SO42− (13%) > Cl (11%). The order of cation concentration is in accordance with that of standard crustal materials (Ca2+ > Mg2+ > Na+ > K+) [28]. Cations followed the order of Ca2+ (35%) = Na+ (35%) > Mg 2+ (29%) > K+ (1%) and anions of HCO3 (46%) > SO42− (33%) > Cl (21%) for the plain groundwater (Table 1). The contribution of NO3 to total anions was less than 1% in the Kaidu River Basin (Table 1).
The predominant hydrochemical types of stream water are Ca2+-HCO3 type water followed by mixed type water. The hydrochemical types of mountain groundwater are very similar to stream water. The predominant hydrochemical types of the plain groundwater are also Ca2+-HCO3 and mixed type water. The Na+, Cl, and SO42− counted as a much larger proportion for groundwater hydrochemistry than for surface water. The water types of both stream water and groundwater in the plain did not change significantly with seasons (Figure 3).
Though there are some differences, the seasonal variations of the major chemical ions in the stream water and the plain groundwater are very similar, except for Na+ and Cl (Figure 4). Owing to the wide distribution of halite, the concentration of Na+ and Cl in groundwater is high all year round [29].
Figure 2 displayed the spatial distributions of the major ions of groundwater in the plains. The SO42−, Na+, and Cl concentrations were much higher in the southern and the eastern part of the plain. The concentrations of HCO3 and Ca2+ decreased along the water flow path. This indicates that the water in the mountain has a stronger dissolving ability for carbonate as compared with water in the plain in the downstream of the Kaidu River. The influence of evaporation is stronger on water in plains as compared with water in the mountain areas.

4.2. Water-Rock Interaction

The dominant hydrochemical formation mechanism (including precipitation dominance, rock weathering dominance, or evaporation dominance) is determined according to the ratios of Na+/(Na+ + Ca2+) against TDS [30]. In stream water and mountain groundwater, the average values of Na+/(Na++Ca2+) is lower than 0.5, and the TDS values are 188.85 and 130.91 mg/L, respectively (Figure 5), showing that rock weathering is the dominant hydrochemical mechanism. In plain groundwater, the average value of Na+/(Na+ + Ca2+) is 0.64, which is much higher than 0.5, and the TDS is 461.30 mg/L, both values are much higher than that of stream water and mountain groundwater (Figure 5). This indicates that the hydrochemistry of groundwater is not only controlled by rock weathering but also by evaporation processes.
Bivariate plots are effective tools to reflect the major water–rock interaction processes [3,25]. In this study, all water samples were plotted around the 1:1 line in the (Ca2+ + Mg2+) vs. (HCO3 + SO42−) plot (Figure 6a). Moreover, the carbonate minerals and sulfate minerals are widely distributed in the area. This indicates that the carbonate and sulfate weathering, similar to dolomite and gypsum weathering, are important geochemical processes [31].
Correlations among the major ions in stream water and groundwater are shown in Table 2. Most of the major cations and anions in stream water and plain groundwater were significantly correlated with TDS, indicating that these ions are dissolved continuously, resulting in the rise of TDS [24]. However, for mountain groundwater, only Ca2+ and HCO3 were significantly correlated with TDS. This is due to the high proportion of Ca2+ + HCO3 (71%) to the total ions, indicating that carbonate weathering is the major hydrochemical process in mountain groundwater.
The high positive correlation between Ca2+ and HCO3 in stream water and mountain groundwater indicates calcite weathering, as expressed in Reaction (1) [24]. The high correlations between Mg2+ and HCO3, Ca2+ and HCO3, and Ca2+ and Mg2+ in stream water indicate the weathering of dolomite, as expressed in Reaction (2) [32]. The high correlations between Mg2+ and SO42−, and Na+ and SO42− in stream water and mountain groundwater indicate the dissolution of sulfate [33]. The high correlations between Na+ and Cl in stream water and plain groundwater, K+ and Cl, and K+ and Na+ in stream water indicate the dissolution of evaporites [34]. The high correlation between Mg2+ and Cl in stream water and plain groundwater indicates dissolution of magnesium chloride. The wide distribution of saline soil in the plains offers the source of magnesium chloride [17]. The high correlation between Mg2+ and Na+ in stream water indicates the weathering of certain kinds of minerals that contain these two ions, such as anorthose, as shown in Reaction (3) [35]. The high correlation between K+ and Mg2+ in stream water indicates biotite weathering, as expressed in Reaction (4) [31].
CaCO 3 + CO 2 + H 2 O Ca 2 + + 2 HCO 3   Reaction   1
CaMg ( CO 3 ) 2 + 2 CO 2 + 2 H 2 O Ca 2 + + Mg 2 + + 4 HCO 3
Dolomite Reaction 2
Na 6 Ca 4 Al 14 Si 26 O 80 + 31 H 2 O + 14 H 2 CO 3 7 AlSi 2 O 5 ( OH ) 4 + 6 Na + + 4 Ca 2 + + 14 HCO 3 + 12 H 4 SiO 4
Anorthose kaolinite orthosilicic acid Reaction 3
14 KAlMg 3 Si 3 O 10 ( OH ) 2 + 98 H 2 CO 3 + 7 H 2 O 7 Al 2 Si 2 O 5 ( OH ) 4 + 42 Mg 2 + + 14 K + + 98 HCO 3 + 28 H 4 SiO 4
Biotite kaolinite orthosilicic acid Reaction 4
The calculated SI values of aragonite, calcite, and dolomite ranged from −0.97 to 0.72, −0.81 to 0.87, and −2.53 to 1.64, with an average of −0.04, 0.12, and −0.18, respectively (Figure 7), indicating that the aragonite, calcite, and dolomite were saturated. They would precipitate under conditions of carbon dioxide content decrease or temperature decrease. However, the SI values of anhydrite, gypsum, and halite ranged from −3.80 to −1.20, −3.55 to −0.95, and −5.94 to −10.67, with an average of −2.23, −1.98, and −7.91 (Figure 7), respectively, indicating the undersaturation of these minerals. They dissolve continually under suitable conditions. In basins in the watershed, including the small Yourdusi Basin, the large Yourdusi Basin, and the Yanqi Basin, lithology is mainly controlled by characteristics of the quaternary sediment. There is a wide distribution of carbonate minerals and sulphate minerals, such as aragonite, calcite, dolomite, anhydrite, and gypsum [36]. Moreover, in the Yanqi Basin, controlled by dry climate and extensive evaporation, there is widespread halite in the topsoil [34].
According to the Na+ vs. Cl plot and (Na+ + K+) vs. Cl plot, though stream water and mountain groundwater samples were plotted very close to the 1:1 line, most samples were scattered above the 1:1 line (Figure 6), indicating other Na+ sources, except for halite, dissolution [34]. The cation exchange between Ca2+ or Mg2+ and Na+ or K+ helps explain the relatively high Na+ concentration in water [33]. In the Kaidu River Basin, the average values of Scholler indices of stream water and groundwater were negative (Table 1), indicating that ion exchange is an important water–rock interaction process [24]. In this process, Ca2+ or Mg2+ was removed from solution and Na+ was released into it (Reaction (5)).
Ca 2 + / Mg 2 + + Na X Na + Ca / Mg + X   Reaction   5

4.3. Effect of Human Activity

NO3 concentration is an indicator of intensity of human activity [37,38]. The NO3 concentration in plain groundwater ranged from 0 to 5.4 mg/L and decreased from west to east, in the spring of 2016 (Figure 2b). The NO3 concentration was lower than 2 mg/L in surface water, while it was higher in pastoral areas with altitudes ranging from 1400 to 3100 m a.s.l. as compared with the high mountains where the altitude is higher than 3100 m or the arable land where the altitude is lower than 1400 m (Figure 8). This indicates that NO3 is mainly obtained from livestock farming in the mountainous area and irrigation in the plains. A large amount of nitric oxides discharge into the environment due to animal excrement and fertilization. The mountain area is the major water gathering area, as well as the major pastoral region in the basin. The plain in the downstream of the river is dominated by irrigation agriculture, which is heavily dependent on drip irrigation. The soil water content is very low in the study area, therefore agricultural waste is not be immediately discharged into groundwater or river. Thus, livestock farming has a larger impact on NO3 concentration than irrigation agriculture.
The population of Yanqi oasis had increased before 2010 and then decreased. The population was 432,369, in 2005, increased to 475,178, in 2010, then decreased to 443,998, in 2015. The farmland area, orchard area, number of livestock, and the use of pesticides and chemical fertilizer have increased significantly. The sown area increased from 93,520 ha, in 2005, to 121,246.7 ha, in 2015. The use of chemical fertilizer increased from 24,072 ton, in 2005, to 52,051 ton, in 2015. The pesticide use increased from 39,088 kg, in 2005, to 405,604 kg, in 2015. The population of livestock increased from 325,640 sheep units, in 2005, to 404,843,0 sheep units, in 2015 (Table 3). Pesticides and fertilizer residuals in soil did not discharge into groundwater or rivers in arid areas without extreme events such as heavy rainstorms; however, they were a potential menace for regional water resource security [4,39]. In addition, animal husbandry which was mainly distributed in water source areas directly contaminated the water.

4.4. Estimation of the Contribution Sources

PCA was conducted to evaluate the contributions of the above discussed factors to water hydrochemistry, in the Kaidu River Basin. Table 4 shows the PCA results of stream water and groundwater.
For surface water, the first three major factors accounted for 84.5% of the total hydrochemical variance. The first factor explained 46.7% of the total variance with high loadings of SO42−, Cl, Na+, K+, and Mg2+, attributed to sulfate and chloride weathering, such as anhydrite, gypsum, and halite. The second factor accounted for 25.0% of the total variance with high loadings of HCO3, Ca2+, and Mg2+, owing to carbonate weathering, such as aragonite, calcite, and dolomite. The third factor accounted for 12.9% of the total variance with high loadings of NO3. It represented an anthropogenic pollution source such as the widely distributed livestock farming in the mountain.
For the plain groundwater, the first three major factors accounted for 89.3% of the total variance. The first factor accounted for 55.3% of the total variance with high loadings of SO42−, Cl, Na+, and Mg2+, indicating the weathering of sulfate and chloride, such as anhydrite, gypsum and halite. The second factor explained 19.9% of the total variance with high loadings of HCO3, Ca2+, and Mg2+, indicating carbonate weathering, such as aragonite, calcite, and dolomite. The third factor accounted for 14.2% of the total variance with high loadings of NO3 and K+, indicating anthropogenic pollution sources including agriculture activities and livestock farming in the plains.

5. Discussion

Hydrochemical processes have a direct effect on regional water quality. Factors that influence hydrochemical processes include both natural and human factors. Natural factors include bed rock properties, geomorphological features, and climatic features. Studies at the Hongfeng Lake in Guizhou province, China found that non-point source pollution caused by heavy rainfall results in increasing the total nitrogen content of the lake [40]. Studies in mountain catchments in northern Taiwan found that the hydrochemical processes during typhoon period are significantly different from other times, the high predictability of ion input and output budgets using stream discharge during the non-typhoon period largely disappeared during the typhoon periods [41]. Chiogna et al. [42] found that hydrochemical processes are significantly affected by drought. For groundwater in the Northwest Namibia, Li et al. [10] found that the TDS increase in groundwater is primarily due to mineral dissolution, rather than evaporation.
The Kaidu River Basin, which is located in an arid area, is characterized by poor precipitation, frequent extreme climatic events, intense evaporation, and a wide distribution of different types of dissoluble minerals. The hydrochemical processes in the basin are affected by bedrock interaction, evaporation, ion exchange, and extreme climatic events. For example, the TDS in surface water and shallow groundwater increased significantly after heavy rainfall.
Human factors that influence hydrochemical processes are mainly related to the types and intensities of human activities. Approximately 27% of the TDS in the Taizi River (Northeast China) resulted from anthropogenic inputs [43]. Base flow solute concentrations (particularly sulfate, chloride, bicarbonate alkalinity, and sodium) in the Chattahoochee River Basin (Georgia, USA) increased with the degree of urbanization [44]. Overuse of chemical fertilizers, industrial sewage discharging, and the vulnerability of groundwater systems to contamination resulted in the NO3 contamination of stream water and groundwater in the Hutuo River Basin, northern China [45]. With widely distributed irrigation agriculture, on the one hand, a large amount of pesticide and chemical fertilizer discharges into soil; on the other hand, irrigation, particularly drip irrigation, results in salts concentration in the topsoil [39,46]. Livestock farming, concentrated farming, or nomadism discharge large amounts of nitrogen, phosphorus, potassium, bacteria, methane, ammonia, and other materials into environment [47,48].
Nomadic herding is highly developed in the mountainous areas of the study site, whereas irrigation agriculture is widely distributed downstream of the catchment. In recent years, the population of livestock, irrigation area, and amount of pesticide and fertilizer used have increased significantly (Table 3). A significant number of chemical elements such as nitrogen, phosphorous, and potassium are discharged into soil and water. These change the hydrochemical processes, threaten regional water source quality, and enhance the spatial and temporal variance of hydrochemical characters. For example, the TDS of surface water and shallow groundwater increased significantly after autumn irrigation.

6. Conclusions

The major hydrochemical processes, in the Kaidu River Basin, were assessed based on hydrochemical data of stream water and groundwater samples. The predominant hydrochemical water type in stream water and mountain groundwater was Ca2+-HCO3, while mixed type water predominated the basin groundwater. Groundwater showed the highest TDS values. Owing to the frequent exchange between surface water and groundwater, there was no significant difference between the solute concentrations in stream water and mountain groundwater. Rock weathering and evaporation were the dominant natural hydrochemical mechanisms. On the one hand, water–rock interaction accounted for 71.6% and 75.1% of the total hydrochemical variations of surface water and groundwater, respectively. Of these, sulfate and chloride weathering explained 46.7% and 55.3%, respectively, while carbonate weathering explained 25.0% and 19.9%, respectively, of the total hydrochemical variations in surface water and groundwater. On the other hand, human activities explained 12.9% and 14.2% of these hydrochemical variations, respectively. Ion exchange was also determined as an important hydrochemical process. For stream water, livestock farming was found to be the most important artificial pollution source, whereas for groundwater in the plains, irrigation was the major contributor.

Author Contributions

Conceptualization, H.C.; data curation, D.L.; formal analysis, D.L.; funding acquisition, S.J.; investigation, H.C.; methodology, D.L.; project administration, A.L.; resources, A.L.; software, H.C.; supervision, S.J.; validation, D.L.; visualization, H.C.; writing–original draft, H.C.; writing–review and editing, D.L. All authors have read and agreed to the published version of the manuscript.

Funding

The research is supported by the National Natural Science Foundation of China (41630859, 41671026, 41701024) and Major Project of the Science and Technology of Qinghai, China (2019-SF-A4-1).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Munia, H.; Guillaume, J.H.A.; Mirumachi, N.; Porkka, M.; Wada, Y.; Kummu, M. Water stress in global transboundary river basins: Significance of upstream water use on downstream stress. Environ. Res. Lett. 2016, 11, 014002. [Google Scholar] [CrossRef]
  2. Gil-Marquez, J.M.; Barbera, J.A.; Andreo, B.; Mudarra, M. Hydrological and geochemical processes constraining groundwater salinity in wetland areas related to evaporitic (karst) systems. A case study from Southern Spain. J. Hydrol. 2017, 544, 538–554. [Google Scholar] [CrossRef]
  3. Shi, X.; Wang, Y.; Jiao, J.J.; Zhong, J.; Wen, H.; Dong, R. Assessing major factors affecting shallow groundwater geochemical evolution in a highly urbanized coastal area of Shenzhen City, China. J. Geochem. Explor. 2018, 184, 17–27. [Google Scholar] [CrossRef]
  4. Solgi, E.; Sheikhzadeh, H.; Solgi, M. Role of irrigation water, inorganic and organic fertilizers in soil and crop contamination by potentially hazardous elements in intensive farming systems: Case study from Moghan agro-industry, Iran. J. Geochem. Explor. 2018, 185, 74–80. [Google Scholar] [CrossRef]
  5. Abou Zakhem, B.; Al-Charideh, A.; Kattaa, B. Using principal component analysis in the investigation of groundwater hydrochemistry of Upper Jezireh Basin, Syria. Hydrol. Sci. J.-J. Des Sci. Hydrol. 2017, 62, 2266–2279. [Google Scholar] [CrossRef]
  6. Kim, D.-M.; Yun, S.-T.; Cho, Y.; Hong, J.-H.; Batsaikhan, B.; Oh, J. Hydrochemical assessment of environmental status of surface and ground water in mine areas in South Korea: Emphasis on geochemical behaviors of metals and sulfate in ground water. J. Geochem. Explor. 2017, 183, 33–45. [Google Scholar] [CrossRef]
  7. Wang, W.H.; Wu, T.H.; Zhao, L.; Li, R.; Xie, C.W.; Qiao, Y.P.; Zhang, H.W.; Zhu, X.F.; Yang, S.H.; Qin, Y.H.; et al. Hydrochemical characteristics of ground ice in permafrost regions of the Qinghai-Tibet Plateau. Sci. Total Environ. 2018, 626, 366–376. [Google Scholar] [CrossRef]
  8. van Geldern, R.; Schulte, P.; Mader, M.; Baier, A.; Barth, J.A.C.; Juhlke, T.R.; Lee, K. Insights into agricultural influences and weathering processes from major ion patterns. Hydrol. Process. 2018, 32, 891–903. [Google Scholar] [CrossRef] [Green Version]
  9. Niu, B.B.; Wang, H.H.; Loaiciga, H.A.; Hong, S.; Shao, W. Temporal variations of groundwater quality in the Western Jianghan Plain, China. Sci. Total Environ. 2017, 578, 542–550. [Google Scholar] [CrossRef]
  10. Li, Z.; Wang, G.; Wang, X.; Li, W.; Shi, Z.; Heike, W.; Shoopala, U.; Collen-Issia, U. Groundwater quality and associated hydrogeochemical processes in Northwest Namibia. J. Geochem. Explor. 2018, 186, 202–214. [Google Scholar] [CrossRef]
  11. Chen, Y.N.; Li, W.H.; Deng, H.J.; Fang, G.H.; Li, Z. Changes in Central Asia’s Water Tower: Past, Present and Future. Sci. Rep. 2016, 6, 35458. [Google Scholar] [CrossRef] [PubMed]
  12. Bolch, T. HYDROLOGY Asian glaciers are a reliable water source. Nature 2017, 545, 161–162. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Pritchard, H.D. Asia’s glaciers are a regionally important buffer against drought. Nature 2017, 545, 169. [Google Scholar] [CrossRef] [PubMed]
  14. Huss, M.; Hock, R. Global-scale hydrological response to future glacier mass loss. Nat. Clim. Chang. 2018, 8, 135. [Google Scholar] [CrossRef] [Green Version]
  15. Deng, H.J.; Chen, Y.N.; Wang, H.J.; Zhang, S.H. Climate change with elevation and its potential impact on water resources in the Tianshan Mountains, Central Asia. Glob. Planet. Chang. 2015, 135, 28–37. [Google Scholar] [CrossRef]
  16. Wang, S.H. Hydrochemical and isotopic characteristics of groundwater in the Yanqi Basin of Xinjiang province, northwest China. Environ. Earth Sci. 2013, 3, 1–14. [Google Scholar] [CrossRef]
  17. Zhang, Z.Y.; Abuduwaili, J.; Yimit, H. The Occurrence, Sources and Spatial Characteristics of Soil Salt and Assessment of Soil Salinization Risk in Yanqi Basin, Northwest China. PLoS ONE 2014, 9, e106079. [Google Scholar] [CrossRef]
  18. Deutsch William, J.; Randy, S. Groundwater Geochemistry: Fundamentals and Applications to Contamination; CRC Press: Boca Raton, FL, USA, 1997; pp. 392–393. [Google Scholar]
  19. Lloyd, J.W.; Heathcote, J.A. Natural Inorganic Hydrochemistry in Relation to Groundwater: An Introduction; Clarendon Press: Oxford, UK, 1985. [Google Scholar]
  20. Parkhurst, D.L.; Appelo, C.A.J. User’s Guide to PHREEQC—A Computer Program for Speciation, Batch-Reaction, One-Dimensional Transport, and Inverse Geochemical Calculations; U U.S. Geological Survey: Denver, CO, USA, 1999. [Google Scholar]
  21. Marion, G.M. The geochemistry of natural waters: Surface and groundwater environments. J. Environ. Qual. 1998, 27, 245. [Google Scholar] [CrossRef]
  22. Appelo, C.A.J.; Postma, D. Geochemistry, Groundwater and Pollution; CRC Press: Boca Raton, FL, USA, 2004. [Google Scholar]
  23. Schoeller, H. Qualitative evaluation of groundwater resources, In Methods and Techniques of Groundwater Investigations and Devement; UNESCO: Paris, France, 1965; pp. 54–83. [Google Scholar]
  24. Liu, F.; Song, X.; Yang, L.; Zhang, Y.; Han, D.; Ma, Y.; Bu, H. Identifying the origin and geochemical evolution of groundwater using hydrochemistry and stable isotopes in the Subei Lake basin, Ordos energy base, Northwestern China. Hydrol. Earth Syst. Sci. 2015, 19, 551–565. [Google Scholar] [CrossRef] [Green Version]
  25. Weynell, M.; Wiechert, U.; Zhang, C.J. Chemical and isotopic (O, H, C) composition of surface waters in the catchment of Lake Donggi Cona (NW China) and implications for paleoenvironmental reconstructions. Chem. Geol. 2016, 435, 92–107. [Google Scholar] [CrossRef]
  26. World Health Organization. Guidelines for Drinking-Water Quality, 4th ed.; World Health Organization: Geneva, Switzerland, 2011; pp. 104–108. [Google Scholar]
  27. Ma, J.; He, J.; Qi, S.; Zhu, G.; Zhao, W.; Edmunds, W.M.; Zhao, Y. Groundwater recharge and evolution in the Dunhuang Basin, northwestern China. Appl. Geochem. 2013, 28, 19–31. [Google Scholar] [CrossRef]
  28. Taylor S, R. Abundance of chemical elements in the continental crust: A new table. Geochim. Cosmochim. Acta 1964, 28, 1273–1285. [Google Scholar] [CrossRef]
  29. Reimann, C.; de Caritat, P. Establishing geochemical background variation and threshold values for 59 elements in Australian surface soil. Sci. Total Environ. 2017, 578, 633–648. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Gibbs, R.J. Mechanisms controlling world water chemistry. Science 1970, 170, 1088–1090. [Google Scholar] [CrossRef]
  31. Wang, C.; Dong, Z.; Qin, X.; Zhang, J.; Du, W.; Wu, J. Glacier meltwater runoff process analysis using delta D and delta O-18 isotope and chemistry at the remote Laohugou glacier basin in western Qilian Mountains, China. J. Geogr. Sci. 2016, 26, 722–734. [Google Scholar] [CrossRef] [Green Version]
  32. Rao, N.S. Geochemistry of groundwater in parts of Guntur District, Andhra Pradesh, India. Environ. Geol. 2002, 41, 552–562. [Google Scholar]
  33. Yin, L.; Hou, G.; Dou, Y.; Tao, Z.; Li, Y. Hydrogeochemical and isotopic study of groundwater in the Habor Lake Basin of the Ordos Plateau, NW China. Environ. Earth Sci. 2011, 64, 1575–1584. [Google Scholar] [CrossRef]
  34. Xu, H.; Hou, Z.H.; An, Z.S.; Liu, X.Y.; Dong, J.B. Major ion chemistry of waters in Lake Qinghai catchments, NE Qinghai-Tibet plateau, China. Quat. Int. 2010, 212, 35–43. [Google Scholar] [CrossRef]
  35. Kumar, K.; Miral, M.S.; Joshi, S.; Pant, N.; Joshi, V.; Joshi, L.M. Solute dynamics of meltwater of Gangotri glacier, Garhwal Himalaya, India. Environ. Geol. 2009, 58, 1151–1159. [Google Scholar] [CrossRef]
  36. Chen, J. Original Appearances and its Evolution and Reformation in Mesozoie Yanqi Basin in Xinjiang; Northwest University: Xian, China, 2007. [Google Scholar]
  37. Jakobczyk-Karpierz, S.; Sitek, S.; Jakobsen, R.; Kowalczyk, A. Geochemical and isotopic study to determine sources and processes affecting nitrate and sulphate in groundwater influenced by intensive human activity—Carbonate aquifer Gliwice (southern Poland). Appl. Geochem. 2017, 76, 168–181. [Google Scholar] [CrossRef]
  38. Soana, E.; Balestrini, R.; Vincenzi, F.; Bartoli, M.; Castaldelli, G. Mitigation of nitrogen pollution in vegetated ditches fed by nitrate-rich spring waters. Agric. Ecosyst. Environ. 2017, 243, 74–82. [Google Scholar] [CrossRef]
  39. Kandler, M.; Blechinger, K.; Seidler, C.; Pavlu, V.; Sanda, M.; Dostal, T.; Krasa, J.; Vitvar, T.; Stich, M. Impact of land use on water quality in the upper Nisa catchment in the Czech Republic and in Germany. Sci. Total Environ. 2017, 586, 1316–1325. [Google Scholar] [CrossRef]
  40. Longyang, Q. Assessing the effects of climate change on water quality of plateau deep-water lake-A study case of Hongfeng Lake. Sci. Total Environ. 2019, 647, 1518–1530. [Google Scholar] [CrossRef] [PubMed]
  41. Chang, C.T.; Huang, J.C.; Wang, L.X.; Shih, Y.T.; Lin, T.C. Shifts in stream hydrochemistry in responses to typhoon and non-typhoon precipitation. Biogeosciences 2018, 15, 2379–2391. [Google Scholar] [CrossRef] [Green Version]
  42. Chiogna, G.; Skrobanek, P.; Narany, T.S.; Ludwig, R.; Stumpp, C. Effects of the 2017 drought on isotopic and geochemical gradients in the Adige catchment, Italy. Sci. Total Environ. 2018, 645, 924–936. [Google Scholar] [CrossRef] [PubMed]
  43. Cao, Y.; Tang, C.; Cao, G.; Wang, X. Hydrochemical zoning: Natural and anthropogenic origins of the major elements in the surface water of Taizi River Basin, Northeast China. Environ. Earth Sci. 2016, 75, 811. [Google Scholar] [CrossRef]
  44. Rose, S. The effects of urbanization on the hydrochemistry of base flow within the Chattahoochee River Basin (Georgia, USA). J. Hydrol. 2007, 341, 42–54. [Google Scholar] [CrossRef]
  45. Zhang, Q.; Wang, H.; Wang, Y.; Yang, M.; Zhu, L. Groundwater quality assessment and pollution source apportionment in an intensely exploited region of northern China. Environ. Sci. Pollut. Res. 2017, 24, 16639–16650. [Google Scholar] [CrossRef]
  46. Joimel, S.; Cortet, J.; Jolivet, C.C.; Saby, N.P.A.; Chenot, E.D.; Branchu, P.; Consales, J.N.; Lefort, C.; Morel, J.L.; Schwartz, C. Physico-chemical characteristics of topsoil for contrasted forest, agricultural, urban and industrial land uses in France. Sci. Total Environ. 2016, 545, 40–47. [Google Scholar] [CrossRef] [Green Version]
  47. Makris, K.C.; Andra, S.S.; Hardy, M.; Sarkar, D.; Datta, R.; Bach, S.B.H.; Mullens, C.P. Nitrous oxide supersaturation at the liquid/air interface of animal waste. Environ. Pollut. 2009, 157, 3508–3513. [Google Scholar] [CrossRef]
  48. Afazeli, H.; Jafari, A.; Rafiee, S.; Nosrati, M. An investigation of biogas production potential from livestock and slaughterhouse wastes. Renew. Sustain. Energy Rev. 2014, 34, 380–386. [Google Scholar] [CrossRef]
Figure 1. Geographic location and land cover of study area and sampling stations. (A) small Yourdusi Basin; (B) large Yourdusi Basin; (C) Yanqi Basin; (D) Bosten Lake; (a) Dashankou station.
Figure 1. Geographic location and land cover of study area and sampling stations. (A) small Yourdusi Basin; (B) large Yourdusi Basin; (C) Yanqi Basin; (D) Bosten Lake; (a) Dashankou station.
Water 12 00467 g001
Figure 2. Spacial distributions of TDS and the major ions in plain groundwater. The unit is mg/L.
Figure 2. Spacial distributions of TDS and the major ions in plain groundwater. The unit is mg/L.
Water 12 00467 g002
Figure 3. Piper diagrams of (a) stream water and (b) groundwater in different seasons.
Figure 3. Piper diagrams of (a) stream water and (b) groundwater in different seasons.
Water 12 00467 g003
Figure 4. Seasonal variation of (ae) major chemistry ions and (h) TDS of stream water and groundwater.
Figure 4. Seasonal variation of (ae) major chemistry ions and (h) TDS of stream water and groundwater.
Water 12 00467 g004
Figure 5. Gibbs diagram of (a) stream water and (b) groundwater.
Figure 5. Gibbs diagram of (a) stream water and (b) groundwater.
Water 12 00467 g005
Figure 6. Scatter diagrams showing relationship between (a) Ca2+ + Mg2+ vs. HCO3+SO42−; (b) HCO3 vs. Ca2+; (c) HCO3 vs. Mg2+; (d) SO42− vs. Ca2+; (e) SO42− vs. Mg2+; (f) SO42− vs. Na+; (g) Cl vs. Ca2+; (h) Cl vs. Mg2+; and (i) Cl vs. (Na+ + K+) of stream water and groundwater in different seasons.
Figure 6. Scatter diagrams showing relationship between (a) Ca2+ + Mg2+ vs. HCO3+SO42−; (b) HCO3 vs. Ca2+; (c) HCO3 vs. Mg2+; (d) SO42− vs. Ca2+; (e) SO42− vs. Mg2+; (f) SO42− vs. Na+; (g) Cl vs. Ca2+; (h) Cl vs. Mg2+; and (i) Cl vs. (Na+ + K+) of stream water and groundwater in different seasons.
Water 12 00467 g006
Figure 7. Saturation indices (SI) of some minerals versus TDS (mg/L) in stream water and groundwater.
Figure 7. Saturation indices (SI) of some minerals versus TDS (mg/L) in stream water and groundwater.
Water 12 00467 g007
Figure 8. Altitude versus NO3 concentrations of the surface water.
Figure 8. Altitude versus NO3 concentrations of the surface water.
Water 12 00467 g008
Table 1. Statistics of the stream water and groundwater hydrochemistry and Schoeller indices.
Table 1. Statistics of the stream water and groundwater hydrochemistry and Schoeller indices.
LocationWater Type HCO3 (mg/L)Cl
(mg/L)
SO42−
(mg/L)
NO3
(mg/L)
Ca2+
(mg/L)
Mg2+
(mg/L)
K+
(mg/L)
Na+
(mg/L)
TDS
(mg/L)
pHCAI-ICAI-II
mountainstream watermean136.29 7.33 28.68 0.70 44.72 8.44 1.37 8.01 166.69 8.04 −1.22 −0.06
min42.29 0.73 0.00 0.03 13.79 0.94 0.21 0.74 44.99 7.56 −4.96 −0.39
max315.78 36.27 174.69 1.60 92.17 41.67 4.59 36.10 431.53 8.29 0.53 0.11
groundwatermean116.33 7.96 11.92 0.43 36.09 7.14 1.37 8.27 130.91 8.05 −1.22 −0.10
min65.10 2.56 0.00 0.17 24.06 4.98 0.10 4.12 84.09 7.93 −2.29 −0.37
max160.17 16.63 23.09 0.64 48.11 8.80 4.07 21.29 158.39 8.25 0.40 0.08
plainstream watermean176.87 14.43 47.26 0.80 51.18 15.70 2.33 19.36 238.69 8.07 −1.07 −0.13
min117.19 6.26 31.44 0.16 43.30 10.79 0.15 6.90 193.37 7.77 −2.29 −0.37
max249.34 35.34 92.38 1.80 65.44 23.48 4.31 53.06 335.35 8.30 −0.59 −0.04
groundwatermean203.85 57.52 118.93 2.25 51.56 27.46 3.81 61.39 422.58 7.91 −1.15 −0.26
min67.71 7.65 19.74 0.00 12.52 4.56 0.59 7.27 194.89 7.27 −3.04 -0.75
max656.37 185.97 346.41 10.13 153.01 117.39 25.03 201.24 948.19 8.28 0.24 0.04
Table 2. Correlation matrix for the major parameters in stream water and groundwater.
Table 2. Correlation matrix for the major parameters in stream water and groundwater.
HCO3ClSO42−Ca2+Mg2+K+Na+TDS
Stream water
HCO310.403 **0.217 *0.735 **0.578 **0.1680.261 *0.720 **
Cl 10.523 **0.425 **0.699 **0.717 **0.871 **0.775 **
SO42− 10.469 **0.731 **0.410 **0.581 **0.775 **
Ca2+ 10.500 **0.359 **0.275 **0.770 **
Mg2+ 10.582 **0.716 **0.883 **
K+ 10.662 **0.558 **
Na+ 10.725 **
TDS 1
Plain groundwater
HCO31−0.051−0.0120.470 **0.447 **0.348 **0.0120.307 **
Cl 10.913 **0.310 **0.695 **0.327 **0.773 **0.903 **
SO42− 10.399 **0.681 **0.373 **0.724 **0.931 **
Ca2+ 10.382 **0.311 **−0.0270.516 **
Mg2+ 10.521 **0.459 **0.811 **
K+ 10.184 *0.461 **
Na+ 10.754 **
TDS 1
Mountain groundwater
HCO31−0.4420.1430.780 *−0.0290.197−0.7010.754 *
Cl 1−0.0640.1940.712 *0.6900.2390.097
SO42− 10.1920.447−0.0840.4560.544
Ca2+ 10.4770.713 *−0.5920.912 **
Mg2+ 10.5390.3410.589
K+ 1−0.4510.493
Na+ 1−0.291
TDS 1
* Indicates that the correlation is significant at 0.05 level; ** indicates that the correlation is significant at 0.01 level.
Table 3. Population and husbandry development during the recent 10 years, in the Yanqi Oasis.
Table 3. Population and husbandry development during the recent 10 years, in the Yanqi Oasis.
YearPopulationSown Area
(ha)
Irrigated Area (ha)Orchard
(ha)
Net Content(t)Pesticide(kg)Livestock
(Sheep Unit)
Nitrogenp FertilizerPotassic FertilizerCompound Fertilizer
2005432,36993,52069,26011,42015,5075929450218639,0883,256,460
2010475,178116,213.392,60018,161.220,77212,55626805194293,7292,282,150
2015443,998121,246.7104,76023,36426,10915,86039986084405,6044,048,430
The data were obtained from the Xinjiang statistical yearbooks of 2006, 2011, and 2015.
Table 4. Results of principal component analysis for major ions in stream water and groundwater.
Table 4. Results of principal component analysis for major ions in stream water and groundwater.
Factor 1Factor 2Factor 3
SWGWSWGWSWGW
Ca2+0.290.320.870.79−0.060.08
Mg2+0.820.850.450.460.020.01
Na+0.90.840.25−0.30.170.05
K+0.80.270.190.30.050.81
HCO30.23−0.030.930.950.020.22
SO42−0.80.960.12−0.03−0.080.12
Cl0.910.990.26−0.030.110.01
NO30.08−0.09−0.030.030.990.93
% Variance46.6655.2924.9819.8512.8714.16
% Cumulative46.6655.2971.6475.1484.5189.3
SourceSulfate and chlorideCarbonateHuman activity
SW, stream water and mountain groundwater; GW, plain groundwater.

Share and Cite

MDPI and ACS Style

Li, D.; Chen, H.; Jia, S.; Lv, A. Possible Hydrochemical Processes Influencing Dissolved Solids in Surface Water and Groundwater of the Kaidu River Basin, Northwest China. Water 2020, 12, 467. https://doi.org/10.3390/w12020467

AMA Style

Li D, Chen H, Jia S, Lv A. Possible Hydrochemical Processes Influencing Dissolved Solids in Surface Water and Groundwater of the Kaidu River Basin, Northwest China. Water. 2020; 12(2):467. https://doi.org/10.3390/w12020467

Chicago/Turabian Style

Li, Dalong, Haiyan Chen, Shaofeng Jia, and Aifeng Lv. 2020. "Possible Hydrochemical Processes Influencing Dissolved Solids in Surface Water and Groundwater of the Kaidu River Basin, Northwest China" Water 12, no. 2: 467. https://doi.org/10.3390/w12020467

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