Temporal Variations of Spring Water in Karst Areas: A Case Study of Jinan Spring Area, Northern China

: Jinan is known as “Spring City,” because of its famous 72 artesian springs. Spring water plays an important role in the social and economic development of Jinan. However, the accelerating process of urbanization and more intensive human activities have significantly affected the Jinan springs. Based on the data from four spring groups (2015–2018), the hydrochemical characteristics of spring water were analyzed and 14 parameters were selected to evaluate the quality of spring water. In addition, the main ions variation characteristics of spring water in Jinan from 1958 to 2018 were analyzed, based on the previous studies. Subsequently, the spring flow dynamics of the Jinan spring area from 1958 to 2012 were also discussed. Overall, the water quality of the four spring groups is good. The chemical composition of the spring water is mainly controlled by the dissolution of calcite, followed by dolomite, and de ‐ dolomitization. However, spring water is affected by human activities, and the hydrochemical type tends to evolve from HCO 3 ‐ Ca to HCO 3 ∙ SO 4 ‐ Ca or HCO 3 ‐ Ca ∙ Mg. From 1958 to 2013, the main ions in spring water increased, while a relatively steady and certain fluctuation trend was observed from 2015 to 2018. Since 1958, the spring water flow in Jinan has experienced four stages, from decline to recovery. The exploitation of karst groundwater, change of land use type, and decrease of atmospheric precipitation are the three most important factors, affecting spring discharge dynamics. The slowdown of the increase of the main ion concentration in Jinan spring water and the resumption of spring water flowing indicate that remarkable achievements have been made to protect spring water, but there is still a long way to go to fully protect Jinan spring water.


Introduction
Hosting large reservoirs of groundwater, karst aquifers are of great significance in coordinating the balance of water environment and human needs [1,2]. Jinan is famous for numerous karst springs. Spring water is a good "city card" for the city and served as a history and culture carrier [3,4]. The four most famous spring groups are including Five Dragon Pool (FDP), Pearl Spring (PS), Spouting Spring (SS), and Black Tiger Spring (BTS). The development of karst in the Jinan spring area is a typical representative of the karst spring system in northern China. Karst groundwater resources are abundant, due to their unique geologicalhydrogeological conditions and monoclinic geological structure characteristics [5] and have always been the main source of water supply for industry, agriculture, and domestic use in Jinan area. Therefore, these groundwater resources are considered as an important resource for social and economic development. However, with the construction of urban modernization, spring water resources have been seriously damaged and polluted [6][7][8][9][10][11]. With the acceleration of industrialization and urbanization, the exploitation of karst water has been increasing, since the 1980s. Furthermore, human activities have polluted groundwater to a certain extent, resulting in the detection of industrial and living pollutants in groundwater, with a gradual increase in the degree of pollution. Current literature research suggests the chemical composition of spring water has a continuous upward trend and the quality of spring water becomes worse year by year [6,[10][11][12]. During recent years, extensive work has been done to protect Jinan spring water, including studies on hydrochemistry, water quality, and dynamics of spring water, which shows the importance of Jinan spring water as well as the attention required for its protection [13]. However, relatively few comprehensive studies exist on the hydrochemistry and flow dynamics of Jinan spring water. Lu et al. [10] determined the evolution of karst water quality in the Jinan spring area and its response to human activities. Yang et al. [8] discussed the impact of major engineering constructions on Jinan spring water and surrounded areas. Based on the groundwater table dynamics and precipitation data of the Jinan spring area, Qi et al. [14] analyzed the multi-time scale characteristics, correlation, and response of groundwater table dynamics of Jinan spring area to precipitation, using wavelet transform method. Zhou et al. [15] analyzed the chemical characteristics of karst water in the Jinan area and evaluated the water quality. Yin et al. [5] investigated the hydrochemistry and environmental isotopes of spring water in the Jinan spring area, through combining the actual hydrogeological conditions of the karst spring area. Xing et al. [3] studied the mixing ratios of spring water recharge sources in the Jinan, through the analysis of spring water level, tracer test, hydrochemical data, and karst development degree and suggested the presence of seasonal differences in spring water recharge sources.
Spring water is the soul of Jinan and its protection has attracted great attention from all walks of life [3]. To protect springs, a lot of work and efforts have been done in the city, which included putting forward the spring protection plans and promulgating the spring protection regulations [16,17]. Despite achieving remarkable results in the regard, due to the expansion of the city and the rapid increase of the urban population, more concerns have been raised, hence the protection of spring water in Jinan has a long way to go. Compared to the 1950s, the urban area has increased from dozens of square kilometers to more than 500 square kilometers. At present, reasonable exploitation of groundwater is required to reduce the impact of spring water on spring flow and control the pollution sources, to protect Jinan spring water. Therefore, based on the previous studies, this study comprehensively analyzed the temporal variation characteristics of spring water, using the long-term data of water quality (four spring groups), spring discharge, precipitation, groundwater exploitation, and urban land area in Jinan spring area. This study aims to (1) determine the temporal variations of spring water quality in the Jinan spring area; (2) characterize the temporal variations of spring water flow; (3) identify the major influence factors, causing spring water changes. This study contributes to further understanding of the temporal variation characteristics of Jinan spring water and can provide useful information for policymakers to protect the Jinan spring water.

Study Area
The Jinan spring area is located at the junction of the northern margin of the mountain area of Shandong province and the piedmont inclined plain, with the south high and the low north and the east high and the low west ( Figure 1). The area is surrounded by the Taishan Mountains from the south and the Piedmont inclined plain and the Yellow River alluvial plain from the north. The central part of the area is dominated by low hills. Jinan is located in the mid-latitude inland zone, with a warm temperate continental climate. The annual average temperature is 14.30 °C, precipitation is 641.68 mm, evaporation (measured on the free water surface) is 1500-1900 mm , with uneven distribution of precipitation during the year. June-September is the flood season, accounting for about 73% of the annual precipitation. The great disparity in annual precipitation and periodic changes are the main characteristics of the weather in the area, which keeps the dry and wet change every 10 years. Figure 1. Location of four spring groups and regional hydrogeological map [16,17].
From the oldest to the newest in Jinan, the exposed sediments in the area are including the Archean Taishan rock group, the Paleozoic Cambrian, Ordovician, Carboniferous, and Permian, and Cenozoic Tertiary and Quaternary. Jinan spring water mainly occurs in the limestone aquifer of Cambrian Zhangxia Formation and Jiulong Group Jiaomidian Formation, Jiulong Group Sanshanzi Formation, and Ordovician Majiagou Group. The former is mainly distributed in the southern mountainous areas, mostly depression springs, while the latter is mostly distributed in the urban area of Jinan City, which is an ascending spring. Among them, the most famous are the four spring groups, whose spring waters flow into the Daming Lake along the moat and then into the Xiaoqing River.
Jinan spring area is located in the Jinan monoclinic tectonic area, on the north wing of Taishan anticline and the strata incline northward as a whole. The north-inclined monoclinic structure is consistent with the south-high-north-low terrain, providing a geological and geomorphological foundation for the formation of spring water. The huge limestone layer with karst fissures provides a good place and channel for the recharge, storage, and migration of karst groundwater, as the source of spring water in Jinan. Jinan magmatic rock mass extends in the East-West direction, constituting a natural barrier for karst groundwater of monoclinic structure in Jinan. Indeed, the barrier of the magmatic rock mass is the key to the formation of spring water in the area. The aquifers of the four spring groups in Jinan mainly consist of dolomite and dolomitic limestone of Ordovician Cambrian Sanshanzi formation, which are cut by Qianfoshan fault and Wenhuaqiao fault, and the karst fissures are relatively developed. The shallow karst development form is mainly honeycomb-shaped, and a few large karst caves are developed, with a diameter of less than 2.0 m for a single karst cave. Karsts are very well developed in the area of BTS-SS-FDP. From 2010 to 2016, the SS group (37.9%) and BTS group (32.5%) accounted for 70.4% of the four spring groups, accounting for a large proportion, while, FDP group and PS group accounted for 20.4% and 9.2%, respectively.
According to the characteristics of aquifer medium and the movement and storage of groundwater in the aquifer, Jinan area can be divided into four major water-bearing formations, including loose rock pore water-bearing formation, carbonate rock fissure karst water-bearing formation, clastic rock with carbonate rock-fissure karst water-bearing formation, and massive rock fracture water-bearing formation. Meanwhile, the infiltration of atmospheric precipitation surface water and leakage of irrigation water are identified as sources of karst groundwater recharge in the study area. Among them, atmospheric precipitation infiltration recharge is the main source of karst groundwater in the Jinan spring area, with a recharge area in the southern mountainous area ( Figure 1). The runoff direction and runoff path of karst groundwater are controlled by topography, landform, lithology, and structure [18]. The movement direction of karst groundwater is consistent with the inclination direction of topography and rock strata. The runoff is generally from the south to the north, with a larger hydraulic gradient in the southern mountainous area. Groundwater drainage in the spring area mainly includes spring water drainage, artificial mining drainage, subsurface drainage, surface drainage, and evaporation drainage. Groundwater exploitation and spring drainage are the main ways of groundwater drainage in the Jinan spring area [3].

Sampling and Analyses
A total of 80 groundwater samples were collected from the four spring groups from March 2015 to July 2018. The sampling was performed every two months. A dry and clean polyethylene plastic bottle was used to hold spring water. Before sampling, the sampling bottle was cleaned and rinsed 2-3 times with the spring water to be taken. All spring water samples were refrigerated and sent to the laboratory of Shandong Geological Environmental Monitoring Station, as soon as possible for water quality analysis. The pH value of spring water was determined in the laboratory, using pH acidity meter (PHS-3C, Lei magnetic, Shanghai, China). Total dissolved solids (TDS), K + , Na + , Ca 2+ , Mg 2+ , H2SiO3, and total hardness (TH) were analyzed, using inductively coupled plasma-optical emission spectrometer (optima7000DV, Perkinelmer, Waltham, MA, U.S.). SO4 2− , Cl − , NO3 − , and F − were measured using ion chromatography (ICS-600 model, Thermo, Waltham, MA, U.S.). HCO3 − and chemical oxygen demand (COD) were determined by titration, while, NO2 − was measured by colorimetry. In addition, the normalized inorganic charge balance (NICB) was calculated, as follows:

Methods
In this study, hydrochemical types were presented by the Piper diagram and Durov diagram, as effective tools for determining water types. Descriptive statistics, correlation analysis, and principal component analysis (PCA) were performed using SPSS 22.0 software. The main controlling factors of the spring water chemistry of four spring groups were analyzed using the ratio of major ions. The equilibrium state of the minerals in spring water was described by the saturation index (SI), which was calculated using Phreeqc interactive software (version: 3.4.0-12927, USGS). In addition, the spring water quality of four spring groups for drinking was characterized by a water quality index (WQI), as a widely-used tool to determine water quality [19][20][21]. The procedures of calculating WQI are illustrated in Figure 2. In this study, a concept of discharge precipitation ratio (DPR) is introduced to determine the recharge conditions of precipitation to groundwater in the spring area. Consequently, the annual groundwater withdrawal (Qk) in the spring area is added to the total annual spring water flow (Qq) to form the annual total discharge (Qz) of the groundwater in the spring area. We refer to the ratio of the total discharge of groundwater to the annual precipitation (P) in the spring area as DPR, as follows:

Spring Water Chemistry
Statistical results of the main chemical components (March 2015 to July 2018) in the spring water of the four groups in Jinan are presented in Table 1. Overall, the spring water of the four springs in Jinan consists of weak alkalinity water, with the TDS lower than 1000 mg/L and the TH lower than 500 mg/L. The dominant cations are including Ca 2+ and Mg 2+ , while, HCO3 − and SO4 2− are the main anions in the spring water of the four groups. Except for NO2 − , Na + , and COD, the coefficient of variation (Cv) of other chemical components are small, indicating their relatively stable conditions in the spring water. As NO2 − and COD are one of the signs of pollution, the large coefficient of variation for these two components indicates the spring water is affected by human activities, although their content in spring water is very low.
TH and TDS are considered as significant parameters for evaluating water quality suitability [22]. The four spring groups of FDP, PS, SS, and BTS has an average TDS concentration of 566.60 mg/L, 561.14 mg/L, 622.11 mg/L, and 701.01 mg/L, respectively, and the average TH concentrations of 331.47 mg/L, 320.47 mg/L, 361.51 mg/L, and 404.55 mg/L, respectively (Table 1). Based on the TH and TDS, the water quality of the four spring groups for the duration of March 2015 to July 2018 is excellent. NO3 -can usually reflect the impact of human activities, such as industrial and agricultural production, on groundwater [23]. The mean NO3 − concentrations of FDP, PS, SS, and BTS are 37.47 mg/L, 33.80 mg/L, 43.44 mg/L, and 52.60 mg/L, respectively (Table 1). Among them, only the mean NO3 − concentration of the BTS exceeds the WHO drinking water standards (i.e., 50 mg/L), however, the value does not exceed China's national standards for drinking water (i.e., 88.57 mg/L) [24]. Other major ions, such as Cl − , SO4 2− , Na + , and F − are lower than the WHO standards [25].
Based on the comparison between the four spring groups, except for the similar value of Fand pH, the chemical components of the four spring groups are different. Among them, the average concentrations of TDS, TH, Ca 2+ , HCO3 − , and NO3 − present an order of BTS > SS > FDP > PS, while Mg 2+ , Cl − , SO4 2− , and COD demonstrate an order of BTS > SS > PS > FDP, while K + follows the order of PS > SS > BTS > FDP (Figures 3 and 4).  Table 1. Statistical results of hydrochemical components and SI of four spring groups. Schoeller diagram is a commonly used method for hydrochemical characterization, which can clearly show the changes in the concentration of major ions, in water samples [26]. According to the Schoeller diagram, obvious differences are observed in the main ion concentration of each spring group. Meanwhile, the same trend is observed for the water sample curves of the four groups, with a slight difference at individual water sample points, indicating the same source of recharge for these four groups.

Water Types
Piper diagram [27] and Durov diagram [28] are identified as effective tools for understanding the hydrochemical facies and hydrogeochemical regimes [29,30]. It can be seen from both diagrams in Figure 5 that the cations are mainly concentrated at the Ca 2+ end (i.e., region A in the figure), while, the anions are mainly distributed at the HCO3 − end (i.e., region E in Figure 5). In the diamond part of the Piper diagram, the spring water sample points are mainly located in the region 5, while in the Durov diagram, the water sample points are mainly distributed in the small square on the upper right-hand side, which indicates the water chemistry of the spring water is mainly HCO3-Ca type. This can be explained by a large number of carbonate rocks, distributed in the area, therefore, the weathering of carbonates is the main source of Ca 2+ and HCO3 − in the spring water. In addition, the evolution trend of hydrochemical types of spring water can be observed in the Piper diagram and Durov diagram ( Figure 5). In the Piper diagram, the water samples gradually change from zone 5 to zone 9 over time, that is, over time, the spring hydrochemical type gradually evolves from the HCO3-Ca type to the mixed type, due to the increase of SO4 2− . In the Durov diagram, in addition to the change in the water type, the TDS values gradually increase over time, and the pH values of many water samples in 2015, 2017, and 2018 are higher than those in other years.

Spring Water Chemistry Formation Mechanisms
Generally speaking, rock weathering, atmospheric precipitation, and evaporation are three important factors, affecting the groundwater hydrochemistry in non-coastal areas. The Gibbs's model diagram [31] can identify the dominant controlling factors of hydrochemistry. In this study, all the spring water samples, collected from the four spring groups, fall into the rock weathering dominance zone (Figure 6), indicating that water-rock interaction plays a significant role in the spring water chemistry.  The ratios of major ion can help to further understand the types of rock weathering sources related to the hydrochemical characteristics of groundwater [32][33][34]. As presented in Figure 7a,b, most water samples are concentrated in the zones below the y = x line, which indicates that the water chemistry of the samples is affected by the dissolution of carbonate minerals, which provides the main source forCa 2+ , Mg 2+ , and HCO3 − in the spring water. In addition, most of spring water samples are plotted in the zones between the y = x line and y = 2x line, but in Figure 7c the water samples are close to y = 2x line, while in Figure 7d, samples are inclined to y = x line, indicating that the spring water chemistry is mainly affected by the dissolution of calcite, followed by the dissolution of dolomite. Most spring water samples fall above the dissolution line (y = x line) of the halite (Figure  7e), indicating the dissolution of halite is not the source of Na + and Cl − . Except for a few water samples, Cl − contents in most water samples of the four spring groups increase with the increase of NO3 − (Figure 7f), which suggests a certain impact of human activities on spring water, in the study area. The saturation index (SI) is an indicator that is often used in hydrogeochemistry to determine the saturation state of mineral phases in groundwater [35][36][37]. In this regard, the states of the mineral, in the supersaturated state, the equilibrium state, and the unsaturated state are indicated by SI > 0, SI = 0, and SI < 0, respectively. Table 1  In addition, the amount of Ca 2+ and SO4 2− increases with the increase of the SI of gypsum, which indicates that the dissolution of gypsum has a certain influence on the hydrochemistry of spring water. Simultaneously, the content of Mg 2+ in most spring water samples increases as the SI of gypsum increases, which may indicate the presence of dedolomitization. As gypsum is not saturated and continues to dissolve, the concentrations of Ca 2+ and SO4 2− in water increase, and Mg 2+ in dolomite is replaced by Ca 2+ to produce calcite precipitation, which leads to dedolomitization [38,39].

Principal Component Analysis
In this study, the PCA method is used to extract the principal components, affecting the hydrochemical components through dimensionality reduction of hydrochemical data (from March 2015 to July 2018) of the four spring groups to determine the impact of different factors on the spring water. The results show that KMO (Kaiser-Meyer-Olkin value) = 0.647 > 0.6 and p = 0.000 < 0.05, suggesting the factor model has good adaptability and can be used for PCA. As presented in Table 2, four principal components with an eigenvalue >1 are calculated, accounting for 79.7% of the total variance. PC1 accounted for 33.8% of the total variance of spring water data, which is significantly correlated with TH, TDS, Ca 2+ , HCO3 − , and NO3 − . The Ca 2+ and HCO3 − of spring water are mainly governed by weathering of carbonate rocks. NO3 − represents the impact of anthropogenic inputs (such as fertilization) on spring water. Therefore, PC1 mainly reflects the water-rock interaction of carbonate and the influence of human activities on the chemical composition of spring water. In addition, NO3 − has a good correlation with Ca 2+ and HCO3 − , which indicates that the increase of NO3 − concentration, caused by human activities, affects the dissolution of carbonate minerals [40].
PC2, accounted for 26% of the total variance, shows positive correlations with Na + , K + , Cl − , and SO4 2− , Na + and Cl − may originate from the dissolution of halite, however, previous analysis on the ratios of major ion (see the previous section) shows that the dissolution of halite is very weak, therefore, PC2 does not represent the dissolution of halite. Weathering of silicate rocks can provide Na + and K + , however, PC2 is negatively correlated with H2SiO3, hence, PC2 is not the product of silicates weathering. Yin et al. [5] suggested that sulfur in the spring water of the four spring groups mainly derives from fire coal, showing a growing trend. Therefore, PC2 may indicate that spring water is affected by precipitation inputs, influenced by fire coal. PC3 accounted for 11.7% of the total variance and has significant correlations with Mg 2+ and COD. Mg 2+ may be represented as a product of the dedolomitization. COD is a pollution index, indicating the content of organic matter in water. The main sources of organic pollution in groundwater are including industrial, agricultural, and domestic pollutions. Thus, PC3 may reflect the effects of dedolomitization and human activities.
PC4 accounted for 8.2% of the total variance and is positively correlated with pH and NO2 − . PC4 likely represents the anthropogenic inputs. However, the pH values and concentrations of NO2 − in the spring water remains below the WHO standards [25], with relatively small variance, which indicated that the human activities, represented by PC4 have less influence on the quality of spring water, in the study area.

Temporal Variations of Spring Water Quality
Based on the hydrochemical data of previous studies on the four spring groups and the data from this study, the temporal variations of spring water quality are discussed. Overall, significant changes are observed in the main ions in the spring water of the four spring groups from 1984 to 2013, with a continuously increasing trend. However, from 2015 to 2018, a relatively stable trend, with obvious fluctuations is observed in the main ions ( Figure 9).
Regarding Mg 2+ , the concentration of Mg 2+ in spring water of the four spring groups gradually increases, before 2006, while the growth rate slows down significantly, between 2006 and 2013. As for Cl − , overall, the Cl − content shows a gradual increase before 2010 and increase slowly after 2010 (Figure 9). The concentration of SO4 2− in spring water shows the characteristics of continuous increase. In addition, the concentration of SO4 2− in PS spring water has obvious fluctuations. Ca 2+ content presents a gradual increase before 2005, followed by a slow increase. The contents of Na + and Cl − in the four spring groups also show a slow upward trend, between 2015 and 2018. In addition, we gathered data on TDS, TH, SO4 2− , and Cl − in SS from 1958 to 2012, through previous studies [4]. Combined with the sampling data, collected in this study, we suggest the changes in the content of these chemical components in spring water from 1958 to 2012 can be divided into two stages ( Figure 10). Before 1990, the content of each component is low, with a slow upward trend. From 1990 to 2012, the contents of each component in spring water increase, rapidly. However, the chemical composition of SS spring water is stable, from 2015 to 2018. The four spring groups are located in the urban area of Jinan, with a relatively close distance between the spring groups, suggesting the same source of recharge. Therefore, considering the chemical composition analysis of the SS, the chemical composition of the spring water of four spring groups from the 1960s to the present has experienced three phases of slow increase, increase, and stability. To better understand the water quality of the four spring groups in Jinan, considering the hydrochemical data from March 2015 to July 2018, 14 parameters (including K + , Na + , Ca 2+ , HCO3 − , Mg 2+ , Cl − , SO4 2− , F − , NO2 − , NO3 − , TH, COD, TDS, and pH) are selected to evaluate the water quality of spring water, using WQI method [19][20][21]. The WQI values of spring water are calculated according to the flow chart, shown in Figure 2. The weight of each water quality parameter is determined, according to its relative importance (Table 3), before the evaluation.
Only one sample from the BTS site, collected in July 2017, shows a WQI value of more than 50 ( Figure 11)    The ratio of karst water exploitation to spring water self-discharge is less than 1:2. In this period, the average discharge of spring is 3.67 × 10 5 m 3 /d, and the karst water level is between 28.75 m and 31.86 m. From 1962 to 1964, with high precipitation, the exploitation of karst water is less than 1.65 × 10 5 m 3 /d, the water level is stable at 31.54-31.86 m, and the total discharge of spring group is maintained at 4.66-5.02 × 10 5 m 3 /d. (1968)(1969)(1970)(1971)(1972)(1973)(1974)(1975)(1976)(1977)(1978)(1979)(1980) The exploitation of karst water in spring area has increased from 2.99 × 10 5 to 6.49 × 10 5 m 3 /d, with a decrease in the flow of spring water from 1.43 × 10 5 to 9.39 × 10 4 m 3 /d. In this period, some springs are temporarily cut off.

Flowing decline stage (1981-2002)
After 1980, the exploitation of karst groundwater in Jinan continues to increase, with a gradual decrease in the spring water flow. The annual average spring water flow is less than 1.

Influencing Factors of Spring Water Quality
As the Jinan spring groups are located in the urban area of Jinan, their runoff area and direct supply area are mostly urban built-up areas [8]. The rapid development of the economy has resulted in stronger and more intensive human activities, which has put tremendous pressure on the groundwater environment. Therefore, the variation of chemical characteristics of spring water in Jinan is closely related to the intensity of human activities [6]. As shown in Figure 13, since the 1980s, the gross domestic product (GDP) and the prime operating revenue of industrial enterprises, above designated size in Jinan, have been divided into three stages including slow growth, rapid growth, and growth. This corresponds to the temporal variation of chemical components in spring water, indicating that human activities are the main factor, affecting the chemical variation characteristics of spring water in Jinan. In addition, excessive exploitation of groundwater and changes in land-use patterns have changed the dynamic conditions of groundwater, accelerating the transportation of pollutions from contaminated rivers and sewage to groundwater, hence affecting the quality of spring water.

Influencing Factors of Spring Water Flow
The influencing factors of spring water flow in Jinan are including the recharge conditions, precipitation, and groundwater exploitation.
Restricting by the Yellow River, from the north, hinders the urban development in Jinan City in the northern area. Since the founding of New China, urban construction has rapidly developed and the urban area has grown from 23 km 2 in the 1950s to more than 500 km 2 in 2015 ( Figure 14). The development direction of Jinan is to the east and south. Considering that the eastern and southern mountainous areas are the recharge areas of Jinan karst water and Jinan spring groups [8], urban development has resulted in a serious reduction of the effective recharge area of the spring water, which on the one hand reduces the amount of precipitation infiltration recharge spring water, and on the other hand increases the peak flow in the flood season (take the Xiaoqing River Basin as an example, Figure 15). With the continuous development of the social economy and the increase of population, the karst water production increases year by year and the groundwater level decreases year by year, followed by a decrease in the spring water flow. Since 2003, the groundwater exploitation in the urban area has been close to zero ( Figure 16). Since then, the springs have been continuously flowing up to now, which suggests karst water exploitation is the main cause of spring water cut-off in Jinan. The decrease of precipitation in the Jinan area is also considered an important factor, affecting spring discharge. The low average precipitation in the 1980s (i.e., less than 600 mm) and droughts, occurred in 1988-1989 and 1999-2002 [8], result in the reduction of the spring water recharge ( Figure  16). In addition, the decrease of river seepage, caused by reservoir construction and soil erosion, caused by overexploitation in southern mountainous areas, are also identified as deriving factors for the decrease of spring water flow, in Jinan.
The conditions of spring area are complex, as the total discharge of karst water of spring area can only be calculated by considering the karst water exploitation and spring water flow, without considering evaporation. However, the recharge of the whole spring area cannot be accurately calculated for various reasons. Therefore, the concept of DPR is introduced in this study.  1955 1959 1963 1967 1971 1975 1979 1983 1987 1991 1995 1999 2003 2007 Urban groundwater discharge includes the groundwater withdrawal and spring water flow. If the recharge conditions remain the same, the DPR value should also be stable, however, it is not the case ( Figure  17). As shown in Figure 17, the DPR value fluctuates between 0.05 and 0.10 before 2000, with a maximum of 0.14, and maintains a relatively stable trend. When the annual precipitation is high, the DPR value is small, and vice versa. This coincides with the fact that the surface runoff of the spring area is high when the precipitation is high, with a non-linear relationship, which is due to the limited storage space of the groundwater. However, after 2000, smaller values of DPR observed, i.e., less than 0.05, with a gradually decreasing trend. It indicates that the recharge conditions of precipitation to groundwater have deteriorated, indeed, the ground hardening in the recharge area has reduced the recharge of groundwater in the spring area, under the same precipitation conditions. In addition, in this study, we also add the karst water exploitation in the urban and spring areas to the spring water flow, as the total discharge of groundwater in urban and spring areas, respectively, and compare them with the precipitation of the current year (P1), the last two years (P2), the last three years (P3), the last four years (P4), and the past five years (P5). As shown in Figure 18, a trend graph is made and its polynomial trend equation and corresponding correlation coefficients are calculated. Based on the results: (1) The correlation between discharge and precipitation becomes better with the increase of precipitation years, which indicates that the precipitation in the last year, the last two years, and even the last four years affect the total discharge of karst water in that year. Thus, it can be concluded that annual precipitation can affect spring water flow in the next five years at least.
(2) According to the dynamic changes of DPR in the spring area, before 2000, the DPR value shows an initial increase, followed by a decrease, which indicates that the exploitation of groundwater can increase the recharge of precipitation to groundwater, i.e., "recharge increment" [41]. However, after 2000, a sudden decrease trend is observed, causing the ground hardening in the recharge area, which results in the decrease of groundwater recharge by precipitation.  16 1958 1961 1964 1967 1970 1973 1976 1979 1982 1985 1988 1991 1994 1997

Conclusions
Spring water has an irreplaceable effect on the city of Jinan. To better protect the spring water in the city, the data on the water quality and spring discharge of the four major spring groups are used to discuss the dynamic change of spring water time and its influencing factors. The main findings of our study are as follows: (1) From 2015 to 2018, the TDS and TH of the four spring groups in Jinan are less than 1000 mg/L and 500 mg/L, respectively, and WQI indicates good water quality in the area. Compared to the other three spring groups, the contents of chemical components in the water samples collected from the BTS site, are higher. The hydrochemistry of spring water is mainly controlled by the dissolution of calcite, dolomite, and gypsum, as well as the occurrence of dedolomitization in the study area. However, the quality of spring water is affected by anthropogenic inputs. Since 1958, the main ions in spring water have experienced three stages of slow increase, increase, and slow increase, in response to intensive urban development in the area. (2) From 1958 to 2012, Jinan spring flow has experienced four stages, including the natural outflow stage, flowing transitional stage, flowing decline stage, flowing recovery stage. Groundwater exploitation, urban construction land area, and precipitation are the three most important factors, affecting the spring flow dynamics in Jinan. The changes inland use not only affects the spring water flowing in Jinan, but also increases the risk of flooding of the urban area. Based on the results from DPR, annual precipitation can affect spring flow in the next five years, at least.
Finally, Jinan has achieved remarkable results in spring water protection, however, there is still a long way to go. In the future, the protection of spring water should be carried out simultaneously from two aspects: spring water quality and spring water flow. To restore and protect spring water in Jinan, the dynamic monitoring of spring water, the control of pollution sources, and the prohibition of groundwater overexploitation should be strengthened. In addition, the ecological environment of spring water recharge areas should be strongly protected.