Next Article in Journal
Maximization of Total Profit for Hybrid Hydro-Thermal-Wind-Solar Power Systems Considering Pumped Storage, Cascaded Systems, and Renewable Energy Uncertainty in a Real Zone, Vietnam
Previous Article in Journal
LCA Operational Carbon Reduction Based on Energy Strategies Analysis in a Mass Timber Building
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploring the Hydrogeochemical Formation and Evolution of the Karst Aquifer System in the Yufu River Based on Hydrochemistry and Isotopes

1
Water Resources Research Institute of Shandong Province, Jinan 250013, China
2
Shandong Provincial Key Laboratory of Water Resources and Environment, Jinan 250013, China
3
Key Laboratory of Groundwater Resources and Environment, Ministry of Education, Jilin University, Changchun 130021, China
4
College of New Energy and Environment, Jilin University, Changchun 130021, China
*
Authors to whom correspondence should be addressed.
Sustainability 2024, 16(15), 6580; https://doi.org/10.3390/su16156580
Submission received: 17 June 2024 / Revised: 12 July 2024 / Accepted: 29 July 2024 / Published: 1 August 2024

Abstract

:
Jinan, renowned as the “Spring City” in China, relies significantly on karst groundwater as an indispensable resource for socio-economic development, playing a crucial role in ecological regulation, tourism, and historical and cultural aspects. The Yufu River basin, situated within Jinan’s karst region, represents a vital riverine leakage zone. Therefore, investigating the evolutionary characteristics and causative mechanisms of surface water and groundwater at different aquifer levels in the Yufu River basin can provide a scientific foundation for the protection of Jinan’s springs. This study, based on hydrogeochemical and isotopic data from the river water, shallow groundwater, deep groundwater, and springs in the Yufu River basin, explored the hydrogeochemical evolution in this region. The findings revealed significant spatial variations in the hydrochemical parameters of the Yufu River basin. Groundwater received contributions from surface water, while springs represented a mixture from both surface water and various recharge aquifers. Dominant ions include Ca2+ and HCO3, with prevailing hydrochemical types being HCO3·SO4-Ca and HCO3-Ca. Atmospheric precipitation served as the primary source of recharge for surface water and groundwater in the Yufu River basin, albeit influenced by pronounced evaporation processes. The hydrochemical composition in the Yufu River basin was primarily attributed to water–rock interactions, mainly driven by the combined effects of carbonate rock, silicate rock, and gypsum weathering and dissolution. Among these, the weathering and dissolution of carbonate rocks played a dominant role, with human activities exerting a relatively minor influence on the hydrochemistry of the Yufu River basin.

1. Introduction

Karst water constitutes a crucial water source for urban supply, industrial, and agricultural activities, playing a pivotal role in safeguarding residents’ livelihoods and supporting sustainable economic and social development [1,2,3]. Jinan, renowned as China’s “City of Springs”, boasts abundant karst underground water resources, integral to Jinan’s socio-economic development [2,3,4]. These resources hold indispensable value in regional ecological regulation, tourism, and historical and cultural domains [5]. Since the 1970s, in tandem with population growth and rapid socio-economic development, the extraction of karst underground water in Jinan has steadily increased. This has resulted in a continuous decline in the groundwater table, water quality degradation, cessation of spring flow, and a substantial impact on the distinctive landscape of the “Spring City” [5,6,7]. It was not until the early 21st century, following the widespread implementation of ecological regulatory measures, that the four major spring groups experienced a resurgence [8]. The Yufu River basin, a crucial river-type leakage zone within the Jinan Spring area, has been targeted for a comprehensive multi-source recharge initiative [9]. Leveraging the natural advantages of the Yufu River’s channel, water conditions, and groundwater infiltration, the government has undertaken extensive recharge projects in the Yufu River basin [10]. Spring water protection has remained a focal and challenging aspect of hydrogeological efforts in Jinan’s spring area. Current research predominantly focuses on hydraulic connectivity [11,12,13]. Given the significant variations in water quality from diverse sources in the Yufu River, understanding the hydrochemical characteristics and formation mechanisms of river water and groundwater post-artificial recharge holds paramount significance for the sustainable utilization and conservation of karst water resources.
Groundwater chemical composition is a crucial component of the groundwater environment, influenced not only by precipitation composition and quantity but also by aquifer minerals, geological structures, groundwater dynamics, and human activities [14]. Specifically, the primary chemical characteristics of groundwater are closely associated with the processes of dissolution and precipitation of minerals in the aquifer, and the type of water–rock interaction predominantly governs the evolution of groundwater hydrochemistry [15,16]. These interactions serve as “fingerprints” for recording the geochemical history and current status of different aquifers, crucially providing key evidence for hydrogeochemical interactions occurring during groundwater mixing processes. Additionally, variations in groundwater flow conditions and hydrochemical sources and sinks induced by human activities significantly influence groundwater chemistry [15]. Therefore, investigating the chemical characteristics and formation mechanisms of groundwater is essential for analyzing the evolutionary processes of the groundwater environment, understanding the sources of ions, groundwater circulation, water–rock interactions, and identifying the main controlling factors of groundwater chemical components [17,18,19]. Additionally, hydrogen and oxygen stable isotopes are considered conservative tracers to describe the conceptual hydrogeological model features of karst aquifers, such as determining recharge areas, estimating recharge amounts, revealing responses to rainfall events, separating groundwater flow components, and estimating aquifer storage [20,21,22,23]. Hydrogen and oxygen isotopes have been successfully applied to reveal interactions between surface water and groundwater, identify sources of springs, and elucidate relationships among rainfall, surface water, and groundwater [16]. The integrated application of stable isotopes and hydrochemistry not only helps clarify the formation mechanisms of groundwater chemistry but also provides crucial information for assessing interactions between surface water and groundwater, particularly in understanding the evolution of karst hydrological systems with significant spatial and temporal variability.
Therefore, this study collected hydrochemistry and isotope samples of surface water, groundwater of different aquifers and spring water in Yufu River basin, and analyzed the characteristics and genetic mechanism of hydrochemistry evolution, in order to solve the following problems: (i) characteristics of hydrochemistry changes in surface water, shallow groundwater, deep groundwater and spring water in Yufu River basin; (ii) isotopic variation characteristics of surface water, shallow groundwater, deep groundwater and spring water in Yufu River basin; (iii) characteristics and genetic mechanism of hydrochemical evolution in Yufu River basin. This study can provide a scientific basis for better understanding the complexity and connectivity of karst groundwater in different layers in Jinan, and provide guidance for the protection and spring preservation of karst water resources in Jinan.

2. Materials and Methods

2.1. Study Area

Jinan Spring area is located in the transition zone between the middle Shandong mountains and the sloping plain in front of the mountains. The bottom layer is a slightly northward sloping monoclinal structure, which is a typical karst development area in northern China [10,11,24]. The strata, from top to bottom, consist of the Quaternary, Ordovician Beianzhuang Formation, Donghuangshan Formation, Sanshanzi Formation, Cambrian Chaomidian Formation, Zhangxia Formation, and the Yanshanian igneous rock body [8,11,25]. The aquifer system is composed of the Cambrian Zhangxia Formation limestone, Upper Fengshan Formation, Cambrian–Ordovician Sanshanzi Formation, Chaomidian Formation, and the Ordovician aquifer rock layers. The lithology mainly includes limestone, dolomite, dolomitic limestone, argillaceous limestone, and limy dolomite [8,24,25]. Karst fractures are well developed within the spring area, facilitating good connectivity, which is advantageous for the infiltration of atmospheric precipitation and the runoff, convergence, and discharge of groundwater [8,26]. The spring city is divided from south to north into an indirect recharge zone, a direct recharge zone, and a discharge zone [9].
The Yufu River is located in the western part of the Jinan Spring area, with the Ma Shan Fault on its west side. The riverbed is covered with cobblestones and gravel, serving as the primary conduit for surface leakage replenishing groundwater [27]. The exposed strata in the Yufu River basin mainly consist of the Lower, Middle, and Upper Cambrian series of the Paleozoic era, the Lower Ordovician series, and a 10-meter-thick Quaternary sediment layer composed of sand, gravel, and cobblestones [7,10,27]. A strong leakage zone has formed in the Xikema area, where the riverbank and valley are covered with Quaternary strata, and karst fissures are well developed [7,27].

2.2. Sampling and Analysis Methods

In the period from December 2021 to January 2022, a total of 19 samples were collected in the Yufu River basin. This comprised 5 samples of surface water, 9 samples of shallow groundwater, 4 samples of deep groundwater, and 1 sample of water from the Heihu Spring. The distribution of sampling points is illustrated in Figure 1. On-site procedures involved the use of a GPS device to record the latitude and longitude coordinates of sampling points. Water levels and well diameters were measured using a water level gauge and meter, respectively.
The DDS-307A conductivity meter was employed in the field to assess temperature (T), electrical conductivity (EC), and total dissolved solids (TDS) with resolutions of 0.1 °C, 1 μS/cm, and 1 ppm, respectively. pH and oxidation reduction potential (ORP) were gauged using a pHS-3G pH meter with resolutions of 0.01 pH and 1 mV, respectively. Alkalinity was determined via titration using the Aquamerck® universal titration kit (Newbury Township, OH, USA), with a resolution of 0.05 mmol/L, and each sample underwent titration 2–3 times for enhanced accuracy.
The analysis of cation and anion samples was conducted at the Shandong Provincial Institute of Water Science. The EDTA titration method was employed for testing Ca2+ and Mg2+; Na+ and K+ were measured using atomic absorption spectrometry; NO3, Cl, and SO42− were determined using ion chromatography, and HCO3 and CO32− were quantified through titration. Samples for dissolved organic carbon (DOC), particulate organic carbon (POC), δ18O, and δ2H were dispatched to the Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences. DOC analysis was performed using a total organic carbon analyzer (TOC-LCPH, Shimadzu, Japan), POC measurement utilized an elemental analyzer (Vario MACRO cube, Elementar, Germany), and δ18O and δ2H assessments were conducted using an LGR liquid water isotope analyzer (Los Gatos Research, San Jose, CA, USA).

3. Results and Discussion

3.1. Hydrochemical Characteristics of Yufu River Basin

The results revealed significant variations in hydrochemical parameters in the Yufu River basin (Figure 2). The pH values ranged from 7.02 to 8.19, with an average of 7.565, indicating a moderately alkaline to weakly alkaline water. Specifically, the river water’s pH ranged from 7.91 to 8.19, with an average of 8.096, the shallow groundwater pH varied from 7.06 to 7.63, with an average of 7.334, and the deep groundwater pH fluctuated from 7.02 to 7.60, with an average of 7.365. The pH of river water was significantly higher than that of groundwater, possibly due to the utilization of dissolved inorganic carbon (DIC) by aquatic photosynthetic organisms in river water, as well as the precipitation of calcite and the degassing of CO2, which reduced [HCO3] and led to higher dissolved carbon dioxide concentrations [28,29,30]. Similarly, variations in DO, DIC concentration, and ρCO2 further support these two reasons. The mean DO concentrations in surface water, shallow groundwater, and deep groundwater were 5.542 mg/L, 3.679 mg/L, and 3.27 mg/L, respectively. The mean DIC concentrations were 3.98 mmol/L, 5.4 mmol/L, and 4.475 mmol/L, and the mean ρCO2 concentrations were 1507.93 ppmv, 14,447.87 ppmv, and 10,366.35 ppmv, respectively. Surface water exhibited higher DO concentrations than groundwater, while ρCO2 and DIC concentrations were significantly lower in surface water than in groundwater, indicating the influence of both biological processes and water–rock interactions on the hydrochemistry of the study area. Additionally, the clear interaction between surface water and groundwater in the strong leakage zone results in elevated pH and DO levels in the groundwater compared to other regions, while ρCO2 and DIC concentrations were notably lower.
Additionally, the concentrations of DOC and POC in surface water were significantly higher than those in groundwater (Figure 2). Generally, the sources of TOC and POC can be classified into autochthonous and allochthonous components. (1) Allochthonous sources primarily originate from the surrounding land (such as plants or soil) or the upstream watershed ecosystem; (2) autochthonous sources stem from the degradation of organic organisms in the water column (including photodegradation and bacterial degradation) and the leakage of planktonic bacterial cells [31,32,33]. Typically, the origin of organic carbon can be determined by the C/N ratio of POC, where a C/N > 12 indicates allochthonous sources, and a C/N < 8 indicates autochthonous sources [34,35,36]. It can be observed that the C/N ratios of river water in the Yufu River basin range from 2.19 to 5.13, with an average of 3.49; shallow groundwater C/N ranges from 0.94 to 5.77, with an average of 3.18; deep groundwater varies from 1.03 to 2.22, with an average of 1.46; and spring water has a C/N of 5.09, all of which were less than 8 (Figure 2). Consequently, the organic carbon at the sampling points was predominantly of autochthonous sources, aligning with Zhang et al. [36]. This may be attributed to low rainfall and weak hydrodynamics during the winter, resulting in minimal allochthonous input and emphasizing autochthonous organic carbon sources in the study area [36]. Simultaneously, bacterial C/N ratios were typically around 2 to 3, while algal ratios were around 6 [35]. It is evident that only the C/N ratio of deep groundwater fell within the range of 2 to 3, indicating minimal influence from river water, while shallow groundwater and spring water were influenced by river water.
The variations in the ionic composition in the Yufu River basin are depicted in Figure 3. The average concentrations of cations and anions in the Yufu River water, groundwater, and spring water were consistent. The mean concentrations of major cations were Ca2+ > Na+ > Mg2+ > K+, while the mean concentrations of major anions were HCO3 > SO42− > Cl > NO3. The coefficient of variation was employed to indicate the complexity of factors influencing the formation and evolution of groundwater chemical components. The variation coefficients for surface water were relatively small, with all values except for Cl (0.66) being less than 0.5. For groundwater, the coefficients were larger, particularly for SO42− in deep groundwater, reaching a high value of 1.7. This indicated significant spatial variability and complexity in the sources of parameters in groundwater. In contrast, the spatial variability of parameters in surface water was lower, suggesting a more stable ion source and relatively minor overall spatial concentration changes.

3.2. Stable Isotope Characterization

The δ18O values of Yufu River water ranged from −8.12‰ to −7.38‰, with an average of −7.66‰ and a coefficient of variation of 0.04. The δ18O values of shallow groundwater vary from −8.43‰ to −6.49‰, with an average of −7.75‰ and a coefficient of variation of 0.10. In deep groundwater, the δ18O values ranged from −9.55‰ to −7.61‰, with an average of −8.71‰ and a coefficient of variation of 0.10. The δ18O value for spring water is −7.55‰. The δ2H values of Yufu River water ranged from −63.37‰ to −60.23‰, with an average of −62.92‰ and a coefficient of variation of 0.02. Shallow groundwater δ2H values varied from −67.71‰ to −57.96‰, with an average of −61.51‰ and a coefficient of variation of 0.06. In deep groundwater, δ2H values ranged from −74.53‰ to −63.70‰, with an average of −69.39‰ and a coefficient of variation of 0.08. The δ2H value for spring water was −62.24‰. Surface water showed more positive δ18O and δ2H values, and spatially, there was no significant difference in the δ18O and δ2H values of surface water. This was attributed to surface water being an open system, and most surface water has undergone intense evaporation [37]. Additionally, the δ18O and δ2H values of groundwater near surface water sampling points were more positive than in other areas, while deep groundwater had more negative values. This suggested that in the strong leakage zone, groundwater may undergo evaporation along the river channel or undergo a transformation between surface water and groundwater (Figure 2).
In general, δ18O and δ2H were effective tracers for identifying groundwater recharge sources [22,23,25]. As depicted in Figure 4, most samples fall below the atmospheric precipitation line (GMWL: δ2H = 8δ18O + 10) [38], indicating that atmospheric precipitation is the source of recharge for both groundwater and surface water. According to the local atmospheric precipitation line obtained by Guo et al. [39], which is δ2H = 8.26δ18O + 13.38, the stable isotope regression lines for surface water, shallow groundwater, and deep groundwater, based on the least squares method, were δ2H = 3.53δ18O − 34.44, δ2H = 4.19δ18O − 30.59, and δ2H = 5.81δ18O − 18.76, respectively (Figure 4). Surface water exhibited a smaller slope, indicating exposure to a highly evaporative environment, while the slopes of shallow and deep groundwater successively increase, suggesting a mixture with surface water or an evaporation process [37,40] (Figure 4). The intersection of the evaporation line and the atmospheric precipitation line represents the initial stable isotope composition of precipitation. Most samples have isotopes more positive than the precipitation’s isotopic composition, further confirming the occurrence of evaporation in precipitation. Springwater samples were located near surface water samples but fell on the stable isotope regression lines of shallow and deep groundwater, indicating that spring water was a mixture of surface water and different recharge aquifers (Figure 4).

3.3. Hydrochemical Origin and Evolution Process

3.3.1. Hydrogeochemical Type

The Piper trilinear diagram is widely employed in the study of groundwater hydrochemical characteristics, commonly used for the classification of groundwater hydrochemical components and the analysis of water quality evolution characteristics [41,42]. As depicted in Figure 5, sampling points for surface water, groundwater, and springs in the Yufu River basin were all situated in the upper-left region of the diamond, indicating that the concentration of alkaline earth metal ions (Ca2+ and Mg2+) was greater than that of alkali metal ions (Na+ and K+). This suggested that the chemical nature of groundwater was dominated by alkaline earth metals and weak acids. Additionally, the distribution of spring water was in the middle of the sampling points, indicating that the spring water receives contributions from multiple sources. From Figure 3 and Figure 5, the main cations and anions in the water of the Yufu River basin were Ca2+ and HCO3, respectively. Using the Schöeller classification method to categorize the hydrochemical types, three types were identified in the study area: HCO3·SO4-Ca type, HCO3-Ca type, and SO4-Ca type. Surface water, G4, G6, G5, G13, and the spring were all of the HCO3·SO4-Ca types, while the rest of the sampling points exhibited the HCO3-Ca type. Additionally, from the distribution map of the sampling points, it is evident that G4, G6, and G5 are located near surface water, indicating that the hydrochemical type of groundwater was influenced by surface water (Figure 5).

3.3.2. Water–Rock Interaction

(1)
The Gibbs model
The Gibbs model provides a clearer understanding of the origin and evolution patterns of groundwater chemical components, characterizing whether various ions in groundwater have undergone processes such as precipitation, water–rock interaction, and evaporation concentration [42,43,44]. It analyzes the main controlling factors in the hydrochemical evolution process. The Gibbs model classifies the controlling factors of hydrochemical components into three types: precipitation control, water–rock interaction, and evaporation concentration [42,45]. According to the Gibbs model diagram (Figure 6), the water bodies in the Yufu River basin were concentrated in the water–rock control zone, with TDS concentrations ranging from 100 to 1000 mg/L. The Na+/Na++Ca2+ and Cl+/Cl++HCO3 ratios were all less than 0.5, indicating that the hydrochemical composition in the Yufu River basin was primarily influenced by water–rock interactions. This suggested that while atmospheric precipitation was the main source of surface water and groundwater in the Yufu River basin (Figure 4), water–rock interactions, ion exchange, adsorption, and other processes during runoff and infiltration reduced the impact of atmospheric precipitation on hydrochemical components.
(2)
The mineral saturation index and ion ratio
The characteristics of groundwater chemical components result from the combined influence of various factors. For instance, water–rock reactions are a crucial factor affecting the majority of groundwater characteristics and are controlled by the hydrodynamics of groundwater flow, the physicochemical environment, and the mineral types within the aquifer [14,25,46]. Therefore, the compounds and their relationships in groundwater can explain the origin of chemical components and the hydrogeochemical evolution process. The mineral saturation index describes the dissolution/precipitation tendency of minerals in groundwater [14,39]. The saturation indices (SI) for calcite, dolomite, and gypsum are shown in Figure 7. The calcite saturation index (SIc) ranged from 0.47 to 0.68 for surface water, 0.16 to 0.75 for shallow groundwater, 0.19 to 0.31 for deep groundwater, and 0.47 for spring water. Overall, calcite minerals in the Yufu River basin were generally oversaturated, indicating that groundwater has undergone sufficient reactions with carbonate rocks during the infiltration of atmospheric precipitation and runoff processes, leading to oversaturation of SIc (Figure 7a). The dolomite saturation index (SId) ranged from 0.36 to 0.78 for surface water, −1.02 to 0.94 for shallow groundwater, −0.32 to 0.20 for deep groundwater, and was 0.37 for spring water. Some groundwater samples did not reach saturation for dolomite minerals (Figure 7b). The gypsum saturation index (SIg) ranged from −1.65 to −1.28 for surface water, −1.69 to −1.38 for shallow groundwater, −1.92 to −0.02 for deep groundwater, and was −1.32 for spring water. All gypsum minerals were undersaturated (Figure 7c), which indicated that gypsum dissolution was not the main source of ions in the surface and groundwater of the Yufu River basin.
The impact of weathering from three different rock types, carbonate rocks, evaporites, and silicate rocks, on the chemical composition of groundwater can be discerned through the analysis of relationships such as HCO3/Na+ vs. Ca2+/Na+ and Mg2+/Na+ vs. Ca2+/Na+. As depicted in Figure 7d,e, surface water and groundwater were distributed between silicate and carbonate rocks, indicating that the ion composition in the Yufu River basin is primarily influenced by the joint weathering of carbonate and silicate rock minerals. The predominant dissolved components in rock dissolution in groundwater are Ca2+, Mg2+, HCO3, and SO42−. The weathering of calcite and dolomite is depicted in Figure 7f as 1:1 and 1:2, respectively, from which it can be seen that most of the water samples were located between 1:1 and 1:2, indicating that the water of the Yufu River was affected by a combination of both calcite and dolomite weathering, and that the water samples are closer to 1:1, and therefore, the groundwater ions are more affected by calcite. As the dissolution of gypsum and salts (e.g., NaCl, CaCl2) does not affect alkalinity, an increase in alkalinity and the total concentration of major cations in water results mainly from the dissolution of carbonate and silicate [37,46]. However, Figure 7g shows that the ratio of total cations to bicarbonate exceeds 1:1, indicating the dissolution of gypsum in the study area. Moreover, if all Ca2+, Mg2+, and HCO3 are derived solely from the dissolution of carbonate, the relationship between Ca2++Mg2+ and HCO3 should be close to 1:1. However, samples in the study area lie above the 1:1 line, indicating the presence of gypsum dissolution. Considering the thin layer of gypsum in the geological formation, the dissolution of gypsum needs to be considered. To further assess the influence of calcite, dolomite, and gypsum dissolution on groundwater chemistry, Ca2+-SO42− was used to represent Ca2+ unaffected by gypsum, and Ca2+-0.33HCO3 represented Ca2+ unaffected by carbonate rocks. The effects of calcite and dolomite are indicated by 1:2 and 1:4, respectively (Figure 7j), with the majority of the samples distributed in the middle of the 1:2 line and the 1:4 line and closer to 1:2, suggesting that calcite dissolution has a greater effect on the chemical composition of the groundwater. In Figure 7k, the 1:1 line represents gypsum dissolution, and most samples deviated from this line, suggesting a minimal impact of gypsum on the hydrochemistry of the study area. Additionally, the Ca2+/Mg2+ ratio can be used to analyze the proportion of calcite and dolomite dissolution. Lines denoted as 10:1, 3:1, and 1:1 represent calcite dissolution, mixed dissolution of calcite and dolomite, and dolomite dissolution, respectively. Most water samples fell within the 10:1 and 1:1 range, closer to 3:1, indicating that the dissolution of calcite and dolomite is the primary source of Ca2+ in the study area.

3.4. Impact of Human Activities

The dissolution of calcite, dolomite, and gypsum, along with human activities, constitutes the primary sources of Ca2+ and Mg2+ in groundwater. SO42− is derived from the dissolution of gypsum and human activities. Typically, the correlation between Ca2+ + Mg2+ and HCO3 + SO42− is used to assess the impact of human activities on ions in water. If all Ca2+ and Mg2+ originate from the dissolution of calcite, dolomite, and gypsum, all samples should fall on the 1:1 line. As depicted in Figure 8a, the content of Ca2+ + Mg2+ in the study area exhibited the trend deep groundwater < surface water < shallow groundwater < spring. Most Ca2+ + Mg2+ and HCO3 + SO42− ion relationships were distributed along the 1:1 line, indicating that Ca2+, Mg2+, HCO3, and SO42− in the Yufu River basin primarily originated from the dissolution of carbonate minerals such as calcite and dolomite and evaporite minerals such as gypsum. Additionally, some samples deviated from the 1:1 line, suggesting human activity influence in the study area. The Na+/Cl ratio characterizes the enrichment of Na+ in groundwater. Under natural conditions, Na+ in groundwater primarily originates from atmospheric precipitation, the dissolution of evaporite minerals, and silicate mineral dissolution. Human activities such as industrial production and domestic sewage infiltration can also influence its content. The (K+ + Na+)/Cl ratio in the Yufu River basin was distributed both above and below the 1:1 line, indicating the influence of multiple factors such as mineral dissolution, atmospheric precipitation, and human activities. Human activities, such as those from daily life, industrial and agricultural activities, introduced SO42− and Cl into the groundwater system. Under the strong influence of human activities, the total cation value with Cl + SO42− approaches 1:1. However, in this study, all sampling points fell on the 1:1 line, and the correlation between NO3/Cl is not significant (r2 = 0.077), indicating a minimal impact of human activities on SO42− and Cl. Studies suggested that Cl was a conservative ion in groundwater circulation, unaffected by physical, chemical, and biological processes, making it a reliable indicator of human activities [37]. The assessment of human activity impact on groundwater can be further determined by the relationships NO3/Cl and NO3/Ca2+ vs. SO42−/Ca2+. From Figure 8e,f, it can be observed that despite the relatively minor overall impact of human activities on the hydrochemistry of the study area, the infiltration of domestic sewage, among these activities, had a more significant and pronounced effect on the chemistry of the groundwater.

4. Conclusions

This study investigated the hydrochemical and isotopic characteristics of water in the Yufu River basin, Jinan, focusing on understanding the dynamics of karst groundwater across varying depths. The analysis encompassed water samples from river water, shallow groundwater, deep groundwater, and springs, revealing substantial variability in hydrochemical parameters and isotopic compositions. The study aimed to infer factors influencing hydrochemistry, providing crucial insights for a better understanding of the complexity and connectivity of karst groundwater at different depths in Jinan. The results revealed significant variations in hydrochemical parameters among river water, shallow groundwater, deep groundwater, and spring in the Yufu River basin. There was evident interaction between shallow groundwater and surface water, with the spring representing a mixture from both surface water and various recharge aquifers. The predominant ions in the hydrochemical composition of the study area were Ca2+ and HCO3, with the main hydrochemical types being HCO3·SO4-Ca and HCO3-Ca. Stable isotopes (18O and 2H) served as effective tracers for identifying groundwater recharge sources. Atmospheric precipitation served as the primary source of recharge for surface water and groundwater in the Yufu River basin, undergoing substantial evaporation. Springwater samples were located near surface water samples but fell on the stable isotope regression lines of shallow and deep groundwater, indicating that spring water was a mixture of surface water and different recharge aquifers. Water–rock interactions, primarily dominated by carbonate weathering and dissolution, with secondary effects from silicate and gypsum weathering, controlled the hydrochemical composition. Although human activities had a relatively minor overall impact on groundwater quality, domestic sewage posed a significant potential threat.
Based on the results of this study, we recommend that future research focus on several key aspects. Further isotope studies are necessary to precisely delineate the recharge sources and mixing processes of different aquifers. The seasonal variability of hydrochemical and isotopic parameters should be investigated to understand temporal dynamics and develop effective water management strategies. Additionally, advanced isotopic analysis techniques and numerical models should be utilized to explore the interaction mechanisms between groundwater and surface water and their responses to environmental changes. These insights provide valuable guidance for the sustainable management and protection of karst water resources in the region.

Author Contributions

X.C.: Investigation, Visualization, Resources, Writing—original draft. C.H.: Investigation, Visualization, Writing—original draft, Writing—review and editing. S.L.: Investigation, Visualization. Z.W.: Investigation, Visualization. D.L.: Conceptualization, Visualization, Methodology. Q.G.: Investigation, Visualization. W.Z.: Conceptualization, Resources, Methodology, Investigation, Data curation, Supervision, Writing—original draft, Writing—review and editing, Funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by Shandong Key Laboratory of Water Resources and Environment (Project “Analysis of groundwater microecological characteristics and response elements in the seepage zone of the YuFu River”) (PAN: 2024-07) and Optional Subjects of the Water Resources Research Institute of Shandong Province (SDSKYZX202121-1).

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding authors.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Liang, Y.; Gao, X.; Zhao, C.; Tang, C.; Shen, H.; Wang, Z.; Wang, Y. Review: Characterization, evolution, and environmental issues of karst water systems in Northern China. Hydrogeol. J. 2018, 26, 1371–1385. [Google Scholar] [CrossRef]
  2. Shang, H.; Qi, X.; Zhang, M.; Li, H.; Li, G.; Yang, L. Characteristics, Distribution, and Source Analysis of the Main Persistent Toxic Substances in Karst Groundwater at Jinan in North China. J. Chem. 2020, 2020, 4217294. [Google Scholar] [CrossRef]
  3. Zhu, H.H.; Xing, L.T.; Meng, Q.H.; Xing, X.R.; Peng, Y.M.; Li, C.S.; Li, H.; Yang, L.Z. Water Recharge of Jinan Karst Springs, Shandong, China. Water 2020, 12, 694. [Google Scholar] [CrossRef]
  4. Cao, H.; Dong, W.H.; Chen, H.L.; Wang, R.Q. Groundwater vulnerability assessment of typical covered karst areas in northern China based on an improved COPK method. J. Hydrol. 2023, 624, 129904. [Google Scholar] [CrossRef]
  5. Luo, Q.; Yang, Y.; Qian, J.; Wang, X.; Chang, X.; Ma, L.; Li, F.; Wu, J. Spring protection and sustainable management of groundwater resources in a spring field. J. Hydrol. 2020, 582, 124498. [Google Scholar] [CrossRef]
  6. Liu, X.; Hu, L.; Sun, K. In Analysis of Spring Flow Change in the Jinan City under Influences of Recent Human Activities. In Proceedings of the 8th International Water Resources Management Conference of ICWRS, Beijing, China, 13–15 June 2018; pp. 263–268. [Google Scholar]
  7. Zhang, J.; Jin, M.; Cao, M.; Huang, X.; Zhang, Z.; Zhang, L. Sources and behaviors of dissolved sulfate in the Jinan karst spring catchment in northern China identified by using environmental stable isotopes and a Bayesian isotope-mixing model. Appl. Geochem. 2021, 134, 105109. [Google Scholar] [CrossRef]
  8. Meng, Q.H.; Xing, L.T.; Liu, L.; Xing, X.R.; Zhao, Z.H.; Zhang, F.J.; Li, C.L. Time-lag characteristics of the response of karst springs to precipitation in the northern China. Environ. Earth Sci. 2021, 80, 348. [Google Scholar] [CrossRef]
  9. Chen, X.; Guan, Q.; Li, F.; Liu, D.; Han, C.; Zhang, W. Study on the ecological control line in the major leakage area of Baotu spring in Shandong province, eastern China. Ecol. Indic. 2021, 133, 108467. [Google Scholar] [CrossRef]
  10. Zhang, Z.; Wang, W. Managing aquifer recharge with multi-source water to realize sustainable management of groundwater resources in Jinan, China. Environ. Sci. Pollut. Res. 2021, 28, 10872–10888. [Google Scholar] [CrossRef] [PubMed]
  11. Qian, J.; Zhan, H.; Wu, Y.; Li, F.; Wang, J. Fractured-karst spring-flow protections: A case study in Jinan, China. Hydrogeol. J. 2006, 14, 1192–1205. [Google Scholar] [CrossRef]
  12. Kang, F.; Jin, M.; Qin, P. Sustainable yield of a karst aquifer system: A case study of Jinan springs in northern China. Hydrogeol. J. 2011, 19, 851–863. [Google Scholar] [CrossRef]
  13. Zhang, Z.; Wang, W.; Qu, S.; Huang, Q.; Liu, S.; Xu, Q.; Ni, L. A New Perspective to Explore the Hydraulic Connectivity of Karst Aquifer System in Jinan Spring Catchment, China. Water 2018, 10, 1368. [Google Scholar] [CrossRef]
  14. Deng, X.; Xing, L.T.; Zhang, F.J.; Xing, X.R.; Zhang, Y.F.; Yu, M.; Liu, S.Z.; Pan, W.Y. Mapping regional and nested flow systems in the karst aquifers of Jinan spring using hydrochemical and isotope data. Water Supply 2023, 23, 3323–3344. [Google Scholar] [CrossRef]
  15. Chen, Z.; Zhou, P.; Wang, G.; Mao, H.; Feng, X.; Huang, H. Groundwater chemistry and isotope for interpreting the hydrogeological conditions and hydrochemical evolution of multilayer aquifer system of Donghai island, China. Appl. Geochem. 2023, 159, 105833. [Google Scholar] [CrossRef]
  16. Ke, Y.L.; Song, X.F.; Yang, L.H.; Yang, S.T. River-Spring Connectivity and Hydrogeochemical Processes in a Karst Water System of Northern China: A Case Study of Jinan Spring Catchment. Water 2024, 16, 829. [Google Scholar] [CrossRef]
  17. de Carvalho Filho, C.A.; Bomtempo, V.L.; Santos Cota, S.D.; Pelogia Minardi, P.S.; Passos, R.G. Use of major ions to evaluate groundwater chemistry and identify hydrochemical processes in a shallow coastal aquifer in southeast Brazil. Environ. Earth Sci. 2022, 81, 423. [Google Scholar] [CrossRef]
  18. Wang, X.; Chen, R.; Li, H.; Li, K.; Liu, J.; Liu, G. Detection and attribution of trends in flood frequency under climate change in the Qilian Mountains, Northwest China. J. Hydrol.-Reg. Stud. 2022, 42, 101153. [Google Scholar] [CrossRef]
  19. Xie, C.; Liu, H.; Li, X.; Zhao, H.; Dong, X.; Ma, K.; Wang, N.; Zhao, L. Spatial characteristics of hydrochemistry and stable isotopes in river and groundwater, and runoff components in the Shule River Basin, Northeastern of Tibet Plateau. J. Environ. Manag. 2024, 349, 119512. [Google Scholar] [CrossRef]
  20. Mao, H.; Wang, G.; Shi, Z.; Liao, F.; Xue, Y. Spatiotemporal Variation of Groundwater Recharge in the Lower Reaches of the Poyang Lake Basin, China: Insights From Stable Hydrogen and Oxygen Isotopes. J. Geophys. Res.-Atmos. 2021, 126, e2020JD033760. [Google Scholar] [CrossRef]
  21. Wu, H.; Wu, J.; Li, J.; Fu, C. Spatial variations of hydrochemistry and stable isotopes in mountainous river water from the Central Asian headwaters of the Tajikistan Pamirs. Catena 2020, 193, 104639. [Google Scholar] [CrossRef]
  22. Fang, L.; Gao, R.; Wang, X.; Liu, T. Isotopes-based characterization of precipitation compositions and atmospheric water vapor sources over typical Eurasian steppes in south mongolian Plateau. J. Hydrol. 2022, 615, 128724. [Google Scholar] [CrossRef]
  23. Li, Z.; Gui, J.; Feng, Q.; Zhang, B.; Cui, Q. Contribution of the cryosphere to runoff in “Chinese water tower” based on environmental isotopes. Geosci. Front. 2023, 14, 101613. [Google Scholar] [CrossRef]
  24. Wang, J.L.; Jin, M.G.; Lu, G.P.; Zhang, D.; Kang, F.X.; Jia, B.J. Investigation of discharge-area groundwaters for recharge source characterization on different scales: The case of Jinan in northern China. Hydrogeol. J. 2016, 24, 1723–1737. [Google Scholar] [CrossRef]
  25. Wu, X.C.; Li, C.S.; Sun, B.; Geng, F.Q.; Gao, S.; Lv, M.H.; Ma, X.Y.; Li, H.; Xing, L.T. Groundwater hydrogeochemical formation and evolution in a karst aquifer system affected by anthropogenic impacts. Environ. Geochem. Health 2020, 42, 2609–2626. [Google Scholar] [CrossRef]
  26. Sivelle, V.; Labat, D. Short-term variations in tracer-test responses in a highly karstified watershed. Hydrogeol. J. 2019, 27, 2061–2075. [Google Scholar] [CrossRef]
  27. Zhang, J.; Cao, M.; Jin, M.; Huang, X.; Zhang, Z.; Kang, F. Identifying the source and transformation of riverine nitrates in a karst watershed, North China: Comprehensive use of major ions, multiple isotopes and a Bayesian model. J. Contam. Hydrol. 2022, 246, 103957. [Google Scholar] [CrossRef]
  28. Sun, H.L.; Han, C.H.; Liu, Z.H.; Wei, Y.; Ma, S.; Bao, Q.; Zhang, Y.; Yan, H. Nutrient limitations on primary productivity and phosphorus removal by biological carbon pumps in dammed karst rivers: Implications for eutrophication control. J. Hydrol. 2022, 607, 127480. [Google Scholar] [CrossRef]
  29. Macpherson, G.L. CO2 distribution in groundwater and the impact of groundwater extraction on the global C cycle. Chem. Geol. 2009, 264, 328–336. [Google Scholar] [CrossRef]
  30. Yan, H.; Liu, Z.; Sun, H. Large degrees of carbon isotope disequilibrium during precipitation-associated degassing of CO2 in a mountain stream. Geochim. Cosmochim. Acta 2020, 273, 244–256. [Google Scholar] [CrossRef]
  31. Jin, B.C.; Lin, Z.H.; Liu, W.Y.; Xiao, Y.; Meng, Y.; Yao, X.L.; Zhang, T.T. Spatiotemporal variations of dissolved organic matter in a typical multi-source watershed in northern China: A fluorescent evidence. Environ. Sci. Pollut. Res. 2022, 29, 20517–20529. [Google Scholar] [CrossRef]
  32. Liu, S.; He, Z.; Tang, Z.; Liu, L.; Hou, J.; Li, T.; Zhang, Y.; Shi, Q.; Giesy, J.P.; Wu, F. Linking the molecular composition of autochthonous dissolved organic matter to source identification for freshwater lake ecosystems by combination of optical spectroscopy and FT-ICR-MS analysis. Sci. Total Environ. 2020, 703, 134764. [Google Scholar] [CrossRef]
  33. Zhou, Y.; Davidson, T.A.; Yao, X.; Zhang, Y.; Jeppesen, E.; de Souza, J.G.; Wu, H.; Shi, K.; Qin, B. How autochthonous dissolved organic matter responds to eutrophication and climate warming: Evidence from a cross-continental data analysis and experiments. Earth-Sci. Rev. 2018, 185, 928–937. [Google Scholar] [CrossRef]
  34. Jiang, Y.B.; Ji, H.B. Isotopic indicators of source and fate of particulate organic carbon in a karstic watershed on the Yunnan-Guizhou Plateau. Appl. Geochem. 2013, 36, 153–167. [Google Scholar] [CrossRef]
  35. Lv, Y.C.; Ren, Y.X.; Li, X.Y. Influence of reservoir management on the source and transport of particulate organic carbon in surface waters of the lower Yellow River. J. Soils Sediments 2022, 22, 2548–2556. [Google Scholar] [CrossRef]
  36. Zhang, S.; Lu, X.X.; Sun, H.Q.; Han, J.T.; Higgitt, D.L. Geochemical characteristics and fluxes of organic carbon in a human-disturbed mountainous river (the Luodingjiang River) of the Zhujiang (Pearl River), China. Sci. Total Environ. 2009, 407, 815–825. [Google Scholar] [CrossRef]
  37. Gao, S.; Li, C.; Jia, C.; Zhang, H.; Lv, M.; Sun, B.; Chen, H.; Gang, S.; Meng, F. Hydrochemical and stable isotope (δ2H and δ18O) characteristics and hydrogeochemical processes in the Baotu Spring Basin, Eastern China. Arab. J. Geosci. 2021, 14, 2084. [Google Scholar] [CrossRef]
  38. Craig, H. Standard for Reporting Concentrations of Deuterium and Oxygen-18 in Natural Waters. Science 1961, 133, 1833–1834. [Google Scholar] [CrossRef]
  39. Guo, Y.; Qin, D.J.; Li, L.; Sun, J.; Li, F.L.; Huang, J. A Complicated Karst Spring System: Identified by Karst Springs Using Water Level, Hydrogeochemical, and Isotopic Data in Jinan, China. Water 2019, 11, 947. [Google Scholar] [CrossRef]
  40. Brooks, J.R.; Gibson, J.J.; Birks, S.J.; Weber, M.H.; Rodecap, K.D.; Stoddard, J.L. Stable isotope estimates of evaporation: Inflow and water residence time for lakes across the United States as a tool for national lake water quality assessments. Limnol. Oceanogr. 2014, 59, 2150–2165. [Google Scholar] [CrossRef]
  41. Qu, S.; Liao, F.; Wang, G.C.; Wang, X.B.; Shi, Z.M.; Liang, X.Y.; Duan, L.M.; Liu, T.X. Hydrochemical evolution of groundwater in overburden aquifers under the influence of mining activity: Combining hydrochemistry and groundwater dynamics analysis. Environ. Earth Sci. 2023, 82, 135. [Google Scholar] [CrossRef]
  42. Xiao, Y.; Hao, Q.; Zhang, Y.; Zhu, Y.; Yin, S.; Qin, L.; Li, X. Investigating sources, driving forces and potential health risks of nitrate and fluoride in groundwater of a typical alluvial fan plain. Sci. Total Environ. 2022, 802, 149909. [Google Scholar] [CrossRef]
  43. Singh, G.; Rishi, M.S.; Herojeet, R.; Kaur, L.; Sharma, K. Evaluation of groundwater quality and human health risks from fluoride and nitrate in semi-arid region of northern India. Environ. Geochem. Health 2020, 42, 1833–1862. [Google Scholar] [CrossRef]
  44. Marghade, D.; Malpe, D.B.; Rao, N.S.; Sunitha, B. Geochemical assessment of fluoride enriched groundwater and health implications from a part of Yavtmal District, India. Hum. Ecol. Risk Assess. 2020, 26, 673–694. [Google Scholar] [CrossRef]
  45. Sun, J.; Yan, B.Z.; Li, Y.; Sun, H.X.; Wang, Y.H.; Chen, J.Q. Characterization and Cause Analysis of Shallow Groundwater Hydrochemistry in the Plains of Henan Province, China. Sustainability 2021, 13, 12586. [Google Scholar] [CrossRef]
  46. Gao, S.; Li, C.S.; Liu, Y.Y.; Sun, B.; Zhao, Z.Q.; Lv, M.H.; Gang, S.T. Hydrogeochemical Characteristics and Evolution Processes of Karst Groundwater Affected by Multiple Influencing Factors in a Karst Spring Basin, Eastern China. Water 2023, 15, 3899. [Google Scholar] [CrossRef]
Figure 1. Map showing the location of sampling sites along the Yufu River basin. Sampling points G1, G2, G4–G11 are shallow groundwater. G3, G11–G13 are deep groundwater.
Figure 1. Map showing the location of sampling sites along the Yufu River basin. Sampling points G1, G2, G4–G11 are shallow groundwater. G3, G11–G13 are deep groundwater.
Sustainability 16 06580 g001
Figure 2. Spatial variation in hydrochemical parameters in the Yufu River basin (the orange rectangles represent surface water sampling sites, and the gray rectangles represent deep groundwater sampling sites).
Figure 2. Spatial variation in hydrochemical parameters in the Yufu River basin (the orange rectangles represent surface water sampling sites, and the gray rectangles represent deep groundwater sampling sites).
Sustainability 16 06580 g002
Figure 3. Changes in anions and ions in Yufu River basin. (a) K+; (b) Na+; (c) Ca2+; (d) Mg2+; (e) HCO3; (f) SO42−; (g) Cl; (h) NO3.
Figure 3. Changes in anions and ions in Yufu River basin. (a) K+; (b) Na+; (c) Ca2+; (d) Mg2+; (e) HCO3; (f) SO42−; (g) Cl; (h) NO3.
Sustainability 16 06580 g003
Figure 4. Scatter diagram of stable isotopes of surface water, shallow groundwater, deep groundwater, and spring water samples.
Figure 4. Scatter diagram of stable isotopes of surface water, shallow groundwater, deep groundwater, and spring water samples.
Sustainability 16 06580 g004
Figure 5. Piper diagram of surface water and groundwater in Yufu River basin.
Figure 5. Piper diagram of surface water and groundwater in Yufu River basin.
Sustainability 16 06580 g005
Figure 6. The Gibbs graph of surface water and groundwater in Yufu River basin.
Figure 6. The Gibbs graph of surface water and groundwater in Yufu River basin.
Sustainability 16 06580 g006
Figure 7. Scatter diagrams of mineral saturation index and ionic relationships. (ac) The saturation indices for calcite, dolomite, and gypsum. (d) HCO3/Na+ (meq/L) vs. Ca2+/Na+ (meq/L). (e) Mg2+/Na+ (meq/L) vs. Ca2+/Na+ (meq/L). (f) Ca2++Mg2+ (mmol/L) vs. HCO3(mmol/L). (g) Total cations vs. HCO3. (h) Ca2++Mg2+ (meq/L) vs. HCO3(meq/L). (i) Ca2+—SO42−(mmol/L) vs. HCO3(mmol/L). (j) Ca2+—0.33HCO3(mmol/L)/ vs. SO42−(mmol/L). (k) Ca2+ (mmol/L) vs. Mg2+(mmol/L). (l) Ca2+ (meq/L) + vs. SO42−(meq/L).
Figure 7. Scatter diagrams of mineral saturation index and ionic relationships. (ac) The saturation indices for calcite, dolomite, and gypsum. (d) HCO3/Na+ (meq/L) vs. Ca2+/Na+ (meq/L). (e) Mg2+/Na+ (meq/L) vs. Ca2+/Na+ (meq/L). (f) Ca2++Mg2+ (mmol/L) vs. HCO3(mmol/L). (g) Total cations vs. HCO3. (h) Ca2++Mg2+ (meq/L) vs. HCO3(meq/L). (i) Ca2+—SO42−(mmol/L) vs. HCO3(mmol/L). (j) Ca2+—0.33HCO3(mmol/L)/ vs. SO42−(mmol/L). (k) Ca2+ (mmol/L) vs. Mg2+(mmol/L). (l) Ca2+ (meq/L) + vs. SO42−(meq/L).
Sustainability 16 06580 g007
Figure 8. Scatter diagrams of ionic relationships. (a) Ca2++Mg2+ (meq/L) vs. HCO3+SO42− (meq/L). (b) K++Na+ (meq/L) vs. Cl(meq/L). (c) Total cations (meq/L) vs. Cl+ SO42− (meq/L). (d) NO3 (meq/L) vs. Cl (meq/L). (e) NO3/Cl (meq/L) vs. Cl (meq/L). (f) NO3/Ca2+ (meq/L) vs. SO42−/Ca2+ (meq/L).
Figure 8. Scatter diagrams of ionic relationships. (a) Ca2++Mg2+ (meq/L) vs. HCO3+SO42− (meq/L). (b) K++Na+ (meq/L) vs. Cl(meq/L). (c) Total cations (meq/L) vs. Cl+ SO42− (meq/L). (d) NO3 (meq/L) vs. Cl (meq/L). (e) NO3/Cl (meq/L) vs. Cl (meq/L). (f) NO3/Ca2+ (meq/L) vs. SO42−/Ca2+ (meq/L).
Sustainability 16 06580 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, X.; Han, C.; Li, S.; Wang, Z.; Liu, D.; Guan, Q.; Zhang, W. Exploring the Hydrogeochemical Formation and Evolution of the Karst Aquifer System in the Yufu River Based on Hydrochemistry and Isotopes. Sustainability 2024, 16, 6580. https://doi.org/10.3390/su16156580

AMA Style

Chen X, Han C, Li S, Wang Z, Liu D, Guan Q, Zhang W. Exploring the Hydrogeochemical Formation and Evolution of the Karst Aquifer System in the Yufu River Based on Hydrochemistry and Isotopes. Sustainability. 2024; 16(15):6580. https://doi.org/10.3390/su16156580

Chicago/Turabian Style

Chen, Xuequn, Cuihong Han, Shuxin Li, Zezheng Wang, Dan Liu, Qinghua Guan, and Wenjing Zhang. 2024. "Exploring the Hydrogeochemical Formation and Evolution of the Karst Aquifer System in the Yufu River Based on Hydrochemistry and Isotopes" Sustainability 16, no. 15: 6580. https://doi.org/10.3390/su16156580

APA Style

Chen, X., Han, C., Li, S., Wang, Z., Liu, D., Guan, Q., & Zhang, W. (2024). Exploring the Hydrogeochemical Formation and Evolution of the Karst Aquifer System in the Yufu River Based on Hydrochemistry and Isotopes. Sustainability, 16(15), 6580. https://doi.org/10.3390/su16156580

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