Climate Change and Its Influence on the Karst Groundwater Recharge in the Jinci Spring Region, Northern China

Due to climate change and human activities over the last fifty years, the spring flow volume of karst groundwater has sharply diminished in China. Climate change is one of the critical factors that initiates a series of karst hydrogeologic and water ecological environmental problems, because the precipitation shows a decreasing trend while the temperature shows an increasing trend. The Jinci Spring is one of the largest, most famous springs in northern China. This study employed data from the Taiyuan Meteorological Station and ten precipitation stations in and around the Jinci Spring region as well as the runoff data gathered from two hydrological monitoring stations during 1960–2012. The sliding average method and the Mann-Kendall test were used to analyze the variation tendency of precipitation, temperature, and land evaporation in this area. Finally, the following were calculated: the varying pattern of the karst groundwater recharge amount and the response of the recharge amount to precipitation, land evaporation, and river runoff by quantitative analysis. The results indicated that the precipitation and land evaporation amount decreased at first and then subsequently increased. Likewise, the variation trend of the karst groundwater recharge amount in the spring region was roughly consistent with the precipitation variation pattern. In contrast, the temperature displayed an increasing trend. The climate change resulted in a reduction of the karst groundwater recharge amount, and it had the greatest influence in the 1990s, which caused the karst groundwater recharge amount to decrease 26.75 mm as compared to that of the 1960s (about 39.68% lower than that of the 1960s). The Jinci Spring had zero flow during this period. The reduction in precipitation was one of main factors that caused the cutoff of the Jinci Spring.


Introduction
The karst groundwater in the Jinci Spring region is endowed with rich groundwater resources, stabilized dynamic state, and good water quality.Additionally, it has become the water supply source for urban water supply and industrial water.Due to both climate change and human activity, the Jinci Spring has been cutoff, and the karst groundwater level in the spring source area has declined, which has greatly influenced the urban water supply security and water ecological environment.This was a notable event marking a water environment degradation in Taiyuan over the last 50 years.The recharge of karst groundwater in the Jinci Spring region is primarily affected by local climate change; thus, the influence of climate change on the recharge of groundwater in the spring region has particular significance to the optimization of utilization, the implementation of scientific scheduling, and the management of regional water resources.In general, water management is a strong key to use water resources in a correct way [1].
Since the 21st century, countries in Europe and America have begun to research the quantitative calculation methods of groundwater recharge on different temporal-spatial scales [2,3].The water balance model includes surface runoff as well as groundwater recharge, and it also takes into consideration the influence of climate change and human activities on groundwater recharge, while simultaneously researching the influence of climate change on water resources [4,5].Doveri, M. et al. utilized water isotopes as environmental tracers for the conceptual understanding of groundwater flow [6].Cervi, F. et al. studied 12 springs from the northern Apennines of Italy by means of a comprehensive hydrogeological investigation to unravel recharge processes taking place in a highly fractured slab of flysch rock hosting the corresponding aquifers [7].At present, the development of a groundwater recharge model and a sensitivity analysis of groundwater to climate change are still in the research and exploration stages in China [8].Based on precipitation data obtained from over 600 precipitation stations and the annual average temperature data obtained for over 100 meteorological stations in northern China, Weitai Wang et al. analyzed the change tendency and the change relationship among the sunspot quantity, the precipitation, and climate change as well as estimated the influence of climate change on regional karst groundwater recharge in 2012 [9].
There has been considerable research focused on the Jinci Spring; however, much of the study has been concentrated on analyzing the cause of the spring flow exhaustion, analyzing the influence of human activities on the spring flow, and analyzing the conditions of the Jinci Spring reflow [10,11].In contrast, there has been considerably less research conducted on the influence of climate change on karst groundwater recharge in the spring region [12].For example, only Qinghai Guo et al. analyzed the guiding significance of the spring flow process in the context of climate change [13].
In this paper, the sliding average method and the Mann-Kendall test method have been adopted to analyze the annual variations of climatic factors [14], such as precipitation, temperature, and land evaporation [15,16].The purpose of this study is to reveal the influence of atmospheric precipitation and Fenhe River leakage on karst groundwater recharge in the spring region and to provide a scientific basis for the reasonable exploitation and utilization of karst groundwater resources in the spring region and the establishment of karst groundwater management and protection policies [17,18].

Research Area Overview
The Jinci Spring is located 25 km southwest of Taiyuan, Shanxi Province, and it is one of the largest karst springs in northern China.The spring region is 2030 km 2 (Figure 1).The formation of the Jinci Spring was primarily a result of the karst groundwater of Taiyuan Xishan becoming blocked by the stratum of quaternary aquitard on the eastern side of the Taiyuan-Jiaocheng fault zone when there was basin movement.As a result, it overflowed through the ground and formed the spring.It is a classic piedmont fault overflow spring.
The Jinci Spring region has a typical temperate semi-arid continental monsoon climate, and it is typically characterized by droughts and windy conditions, a high precipitation concentration, intensive evaporation, four distinct seasons, a large temperature difference between day and night, and a short frost-free season.The multi-year (1960-2012) average precipitation of this region was 457.8 mm.The inter-annual, annual, and geographical distributions are extremely uneven.The precipitation of the highest flow year is over three times greater than that of the lowest flow year, and 60% of the annual precipitation is concentrated in the flood season, which extends from June to September.The multi-year average evaporation amount in this region is 1871.8mm (an observed value of 20 cm from an observation vessel), and the multi-year average temperature is 10 • C. The Jinci Spring region has a typical temperate semi-arid continental monsoon climate, and it is typically characterized by droughts and windy conditions, a high precipitation concentration, intensive evaporation, four distinct seasons, a large temperature difference between day and night, and a short frost-free season.The multi-year (1960-2012) average precipitation of this region was 457.8 mm.The inter-annual, annual, and geographical distributions are extremely uneven.The The Jinci Spring region is part of the Fenhe River system of the Yellow River basin.The river system within the area is relatively developed, and the Fenhe River is the largest river in the area; it runs through this area from west to east and then turns south after it reaches the mountainous outlet in Lancun.The Tianchi River, Shizi River, Tunlanchuan, Yuanpinchuan, Dachuanhe, Liulin River and others are the first-grade tributaries of the Fenhe River, but they are seasonal rivers.The spring region is composed of two parts, mountain area and plain area.It is bound by the Taiyuan-Jiaocheng fault.The northwestern part of the fault is a middle-low mountain area, and the terrain is higher in the northwest than in the southeast.The area is about 1771 km 2 , and the elevation is 900-1500 m.The eastern fault zone is a basin plain area, with an area of 259 km 2 and an elevation of 780-850 m.The mountainous area and the basin plain area connect directly with terrain change.
The main aquifer of the Jinci Spring karst groundwater is the limestone of the upper and lower Majiagou Formation and the Fengfeng Formation of middle Ordovician, and the secondary is the dolomite and dolomitic limestone of the lower Ordovician and Cambrian system.The recharge source of the spring karst groundwater system primarily originates from three sources: the first is the precipitation penetrating recharge in the carbonatite bare field to the north of the Fenhe River, which is the main recharge of the Jinci Spring; the second is the seepage recharge in the section of the carbonatite river that reaches Fenhe (the seepage reach of the Fenhe River in this spring mainly includes the Luojiaqu-Zhenchendi seepage section (the length is about 22 km), the Gujiao-Zhaishang seepage section (the length is about 5 km), and the Zhaishang-Shaoshi seepage section (the length is about 16 km)); and the third is the leaking recharge from the groundwater of the overlying formation.
The karst groundwater in the spring region migrates from northwest to southeast.The Malan syncline is located in the middle area of this spring region, which displays in the northwest~southeast direction and extends for 43 km.Because the Ordovician limestone is buried, the karst development is weak and has poor permeability in the synclinal axis.The region in the west Malan syncline and to the south of the Fenhe River is a stagnant zone, and the karst groundwater basin that is located on the east side of this syncline is in the strong runoff zone.In addition, a water-abundance zone with strong karstification and water conductivity runs along the Taiyuan-Jiaocheng fault zone.
Except for the natural discharge points at the Jinci Spring and the Ping Spring, the karst groundwater in this spring discharges a small amount of water into the Taiyuan Basin as subsurface flow.In addition, the exploitation of karst groundwater is another source of groundwater discharge [19].
The Jinci Spring flow has exhibited a reduction tendency year by year since the 1960s, especially in the late 1980s.It declined from 1.69 m 3 /s in the 1960s to 1.13 m 3 /s in the 1970s, to 0.46 m 3 /s in the 1980s, and to 0.15 m 3 /s in the 1990s.The Jinci Spring was cutoff on 30 April 1994.The groundwater level of the entire Jinci Spring sharply declined from 1994 to 2009.The water level at the spring vent of the Jinci Spring declined to its lowest point in history (774.94m), which is 27.65 m lower than the elevation at the spring vent of the Jinci Spring (802.59 m).However, it has exhibited a rising tendency from 2010 to 2012 and has recovered from 777.54 m in 2010 to 788.21 m in 2012 (Table A1 and Figure 2).The average annual precipitation of the Jinci Spring region in the 1960s, 70s, 80s and 90s are 507.7 mm, 456.1 mm, 424.1 mm and 396.7 mm respectively, which shows a descending trend.According to the precipitation and discharge data of the 1960s and 1990s, the precipitation decreased by 21.9%, and the average spring flow decreased by 91.8%, from 1.735 m 3 /s to 0.142 m 3 /s.According to this result, precipitation is not the main factor of the spring flow attenuation.The average annual The average annual precipitation of the Jinci Spring region in the 1960s, 70s, 80s and 90s are 507.7 mm, 456.1 mm, 424.1 mm and 396.7 mm respectively, which shows a descending trend.According to the precipitation and discharge data of the 1960s and 1990s, the precipitation decreased by 21.9%, and the average spring flow decreased by 91.8%, from 1.735 m 3 /s to 0.142 m 3 /s.According to this result, precipitation is not the main factor of the spring flow attenuation.The average annual precipitation from 2000 to 2012 is 493.9 mm, which is 7.9% higher than that from 1960 to 2012 (458 mm).Therefore, the change of precipitation has a small contribution to the change of the karst water level in the Jinci Spring region.The reason why the spring dried up is mainly attributed to groundwater over-abstraction.
The main reasons for the rise of the groundwater level in recent years are as follows: first, by means of utilizing the Yellow River surface water, the local government closed part of the groundwater exploitation wells to reduce the amount of groundwater exploitation; second, the high water consumption enterprises in this region have been shut down or moved; third, the implementation of the Fenhe River resuming flow engineering has achieved continuous flow of the Fenhe River for many years, which increases the recharge of groundwater; fourth, through the implementation of regional ecological environment comprehensive management, vegetation coverage has been raised, and water resources are conserved.After taking the above measures, the karst groundwater level in the spring area has risen remarkably in recent years.

Time-Series Test
To analyze climate change characteristics in this spring and the influence of climate on the karst groundwater recharge, the Mann-Kendall non-parametric test method was employed to conduct a trend analysis and abrupt change test [20].

Trend Test
The Mann-Kendall rank correlation test (referred to as the M-K trend test) is a non-parametric statistical test, which can accurately verify the sequence variation trend [21].
(2) The statistic is calculated, and the calculation formula is: (3) S is normal distribution, and its average value is 0. The variance is Var(S) = n(n − 1)(2n + 5)/18, and the calculation formula of the standard normal statistic is: then the original hypothesis is refused, which means the sequence has an obvious rising or declining trend.
(5) The tilt rate (quantization of monotonous trend) is calculated, and the calculation formula is: where 1 < j < i < n, and the median is the median function.When β > 0, the sequence exhibits a rising trend; when β < 0, the sequence exhibits a declining trend [22].

Abrupt Change Test
(1) The one time sequence is assumed to be X i = (x 1 , x 2 . . ., x n ), and a one-order sequence m i is constructed to show the accumulative quantity of a sample when x i > x j (1 ≤ j ≤ i).One statistic d k is defined as: (2) Under the assumption that the original sequence is stochastic independence, the average value and variance of d k are, respectively: (3) d k is standardized, and then: where u(d k ) is a standard distribution, and its probability can be achieved through calculation or looking it up in a table.Assuming one significance level to be a 0 , when a 1 > a 0 , the original hypothesis is accepted; however, when a 1 < a 0 , the original hypothesis is refused, which means that this sequence exhibits a strong rising or declining trend.All the u( The time sequence x k is ranked in reverse and is then calculated according to formula ( 7), and at the mean time: The variation trend of sequence x k can be further analyzed through analyzing the statistical series UF k and UB k , and the mutation time and area can be explicitly known.If UF k > 0, then the sequence exhibits a rising trend, and if UF k < 0, then the sequence exhibits a declining trend.If they are higher or lower than the critical lines, then the rising or declining trend is remarkable.If there is any point of intersection between the two curves of UF k and UB k , and the point of intersection is between critical lines, the time corresponding to the point of intersection is the mutation starting time [24].

Groundwater Recharge Amount Calculation
The main influencing factors of climate change on the karst groundwater recharge amount in the spring region are precipitation, land evaporation and river runoff [25].The change of the groundwater recharge amount due to precipitation, land evaporation, and river runoff can be analyzed based on a quantitative calculation of karst groundwater recharge amount in the spring region.The water balance method is a general calculation method used for the hydrological cycle elements.This research uses the Fenhe Reservoir hydrologic station as the basin import and the Zhaishang hydrologic station as the basin export to determine the scope of the study area, which is 1587 km 2 .
The length of the seepage zone in the Fenhe River from the Fenhe Reservoir to Zhaishang in the Jinci Spring is 27 km, and the length of the seepage zone in the Fenhe River from Zhaishang to Shaoshi is 16 km.The seepage zone from the Fenhe Reservoir to Zhaishang occupies 62.8% of the total seepage zone length.The exposed limestone area in this spring region is 391 km 2 , which is mainly located to the north of this spring and to the north of the Fenhe River channel.The exposed limestone area in the basin controlled by the Zhaishang hydrologic station is 317 km 2 , which occupies 81.1% of the total exposed area in this spring region.Thus, the research on the recharge of precipitation and the Fenhe River seepage of the zone from the Fenhe Reservoir to Zhaishang to the karst groundwater in the Jinci Spring region is highly important to the research on the karst groundwater recharge in the entire spring region.
The water balance in the basin from the Fenhe Reservoir to Zhaishang can be expressed as the following formula [26]: where Q 1 is the observed runoff amount in the Fenhe Reservoir hydrologic station; P is the basin precipitation; Q 2 is the observed runoff amount in the Zhaishang hydrologic station; E is the land evaporation amount; and R is the recharge amount of precipitation and river channel seepage on the groundwater in this spring region.

Results and Discussion
The data used in this research includes meteorological data and runoff data from two hydrologic stations located at the Fenhe Reservoir and Zhaishang from 1960 to 2012; the station location map is shown in Figure 1.

Climate Change Characteristic Analysis
The precipitation infiltration recharge is the most important recharge source in the Jinci Spring, and its main means of recharge include the direct precipitation infiltration recharge of the exposed carbonatite area, the indirect precipitation infiltration recharge of the overlay area, and the seepage recharge of the surface runoff yield formed by the precipitation in the carbonatite river channel.The precipitation will be transformed in three ways after landing on the ground: to form surface runoff, to infiltrate underground and then form groundwater, and to land-evaporate via vegetation and ground surface.In these three ways, the slope surface has a short lag time that allows the precipitation to quickly import into the channel and flow away, and the land evaporation and underground infiltration demonstrate the reciprocal relationship under the transduction function of all kinds of energy, including water and vapor in an unsaturated zone.Additionally, the temperature controls the land evaporation capacity to some extent.Thus, the precipitation amount and the temperature have important influences on the infiltration of karst groundwater in this area.

Precipitation Variation
The precipitation data in the research uses the annual precipitation data from ten precipitation stations in the Jinci Spring region, and the areal precipitation in the spring region can be obtained through the Thiessen polygons method.The annual precipitation in the spring varies from 217.87 to 676.63 mm, and the multi-year average value is 457.85 mm.The precipitation in the spring region exhibits greater fluctuation at different times.For example, the 1960s and 1970s have a relatively larger fluctuation, the 1980s and 1990s have relatively less fluctuation, and the precipitation amount since 2000 is a significant increase as compared to the 1990s (Table A2 and Figure 3a).The most obvious precipitation decline took place during the 1990s, and it is 111 mm less than that in the 1960s (as shown in Table 1).
Water 2017, 9, 267 8 of 16 676.63 mm, and the multi-year average value is 457.85 mm.The precipitation in the spring region exhibits greater fluctuation at different times.For example, the 1960s and 1970s have a relatively larger fluctuation, the 1980s and 1990s have relatively less fluctuation, and the precipitation amount since 2000 is a significant increase as compared to the 1990s (Table A2 and Figure 3a).The most obvious precipitation decline took place during the 1990s, and it is 111 mm less than that in the 1960s (as shown in Table 1).According to the sliding average method and the linear tendency estimate method, the number of years is 53, and by setting k = 5, the 5a sliding average value of the annual precipitation in the Jinci Spring region can be obtained, and the sliding curve can be charted (Figure 3a).The following conclusions were made from the analysis.The 5a sliding average precipitation from 2007 to 2011 is the largest (549.87 mm), whereas the value from 1997 to 2001 is the lowest (352.77mm).The annual precipitation in the Jinci Spring region generally exhibits an initial decrease and then subsequent increase trend.The variation amplitude of the precipitation in the Jinci Spring region shows volatility.The annual distribution is uneven.Additionally, the average precipitation is 457.85According to the sliding average method and the linear tendency estimate method, the number of years is 53, and by setting k = 5, the 5a sliding average value of the annual precipitation in the Jinci Spring region can be obtained, and the sliding curve can be charted (Figure 3a).
The following conclusions were made from the analysis.The 5a sliding average precipitation from 2007 to 2011 is the largest (549.87 mm), whereas the value from 1997 to 2001 is the lowest (352.77mm).The annual precipitation in the Jinci Spring region generally exhibits an initial decrease and then subsequent increase trend.The variation amplitude of the precipitation in the Jinci Spring region shows volatility.The annual distribution is uneven.Additionally, the average precipitation is 457.85 mm, the maximum precipitation occurred in 1964 (676.63 mm), and the minimum occurred in 1972 (217.87 mm); their quantitative difference is 458.76 mm.
The statistic Z, β, and significance level of the annual precipitation in the Jinci Spring region can be obtained by the Mann-Kendall method to conduct a trend test, and the test results are, respectively, −0.22, −0.41, and 0.05.These results reveal that the precipitation test value Z of the study area is less than 0, but it does not pass the significance test level of 0.05, which means the precipitation in the Jinci Spring region from 1960 to 2012 generally shows a declining trend, and the multi-year average declining degree of precipitation is 0.41 mm/a [27].
According to the M-K abrupt change test result of the annual precipitation time variation (Figure 3b), the precipitation in the spring region before 1971 exhibited a rising trend, whereas a declining trend has been evident since 1972.It passes the 0.05 significance test (U 0.05 = 1.96).According to the intersecting point of the UF k and UB k curves, it can be determined that the mutation points of the annual precipitation existed in 1967 and 2009 in the Jinci Spring region [28].The annual average temperature in the Jinci Spring region varies from 8.7 to 11.8 °C, and the multi-year average value is 10.09 °C over the last fifty years.As for the annual average temperature, the highest occurred in 2006 (11.8 °C), and the lowest occurred in 1967 (8.7 °C).The average temperature the Jinci Spring region increased over time (as shown according to decade in Table 1).The difference between the 1960s, 1970s, and 1980s is not very large, but the average temperature in the 1990s increased sharply to 0.88 °C (about 9.3%) higher than that of the 1960s, and in the 2000s it The annual average temperature in the Jinci Spring region varies from 8.7 to 11.8 • C, and the multi-year average value is 10.09 • C over the last fifty years.As for the annual average temperature, the highest occurred in 2006 (11.8 • C), and the lowest occurred in 1967 (8.7 • C).The average temperature the Jinci Spring region increased over time (as shown according to decade in Table 1).The difference between the 1960s, 1970s, and 1980s is not very large, but the average temperature in the 1990s increased sharply to 0.88 • C (about 9.3%) higher than that of the 1960s, and in the 2000s it was 1.48 • C (about 15.63%) higher than that of the 1960s.This conclusion is in accordance with the conclusion found in literature [29].According to the statistical temperature data, the 5a sliding average temperature from 1967 to 1971 in the Jinci Spring region is the lowest (9.06 • C), and the value from 2006 to 2010 is the highest (11.8 • C); the difference between them is 2.74 • C.
The multi-year average temperature trend test result, as tested by the Mann-Kendall method, is shown in Table 2.It can be seen that the statistic |Z| of the multi-year average temperature in the Jinci Spring region is larger than Z 1−α/2 , which means that the temperature in the entire spring region refused the null hypothesis, and it has a monotonous trend with statistical meaning.The β value of the multi-year average temperature in the table is 0.04, which is larger than 0; it also indicates that the temperature has exhibited a rising trend in the Jinci Spring over the last five decades [30].The M-K abrupt change test of the annual temperature time variation (Figure 4b) shows that the average temperature in the spring region from 1962 to 1978 exhibited a declining trend; however, the average temperature has been on a rising trend since 1979 (in addition to 1985 and 1986), and it is still rising today.This rising trend passes the 0.05 significance test (U 0.05 = 1.96), and the mutation occurred in 1994 (the intersection of the UF k -UB k curve in the figure).

Evaporation Variation
According to the annual precipitation and annual average temperature data from 1960 to 2012, the evaporation in the spring region is calculated based on the Turc land evaporation empirical formula [31]: where Q z is the land evaporation amount (mm), P is the precipitation (mm), and T is the temperature ( • C).Curves representing the annual evaporation and 5a sliding average evaporation time sequence in the spring region from 1960 to 2012 are shown in Figure 5a.The two curves generally show an initial decreasing and then subsequent increasing trend over time, which is similar to the precipitation variation trend (Figure 5a).
The annual land evaporation in the Jinci Spring region varies from 213.74 to 467.9 mm, and the multi-year average value is 369.29 mm.The lowest was 213.74 mm in 1972, and the highest was 467.9 mm in 2007; the difference between them is 254.15 mm.
where Qz is the land evaporation amount (mm), P is the precipitation (mm), and T is the temperature (°C).
Curves representing the annual evaporation and 5a sliding average evaporation time sequence in the spring region from 1960 to 2012 are shown in Figure 5a.The two curves generally show an initial decreasing and then subsequent increasing trend over time, which is similar to the precipitation variation trend (Figure 5a).The annual land evaporation in the Jinci Spring region varies from 213.74 to 467.9 mm, and the multi-year average value is 369.29 mm.The lowest was 213.74 mm in 1972, and the highest was 467.9 mm in 2007; the difference between them is 254.15 mm.
According to Table 1, the evaporation at each time period decreased first and then subsequently increased over time.As compared to the 1960s, the 1970s decreased 20.77 mm (about 5.4%), the 1980s decreased 33.6 mm (about 8.75%), the 1990s decreased 44.43 mm (about 11.57%), and 2000-2012 increased 16.43 mm (about 4.28%).As for the 5a sliding average evaporation, the value from 2007 to 2011 is the highest (428.01 mm), and the value from 1997 to 2001 is the lowest (317.85 mm); the difference between them is 110.15 mm.
The M-K mutation test of the annual land evaporation variation over time (Figure 5b) has shown that the land evaporation in the spring region from 1960 to 1969 exhibited a rising trend.From 1970 to 2010, it exhibited a declining trend, and 2011 and 2012 demonstrated a rising trend.According to the intersecting point of the UFk and UBk curves, it can be determined that the mutation existed in 1964 and 2007 [32].

Analysis of the Influence of Climate Change on Karst Groundwater Recharge in the Spring Region
The recharge amount of the Jinci Spring region from the precipitation and river channel seepage in the Fenhe Reservoir-Zhaishang zone can be obtained according to Formula (9).
The annual recharge amount in the Fenhe Reservoir-Zhaishang zone of the Jinci Spring region varies from 29.32 to 140.28 mm (about 11,500-54,900 thousand m 3 ), and the multi-year average value is 59.37 mm (about 23,200 thousand m 3 ).
The karst groundwater recharge amount in the study region has a larger fluctuation at different times.The 1960s, 1970s, and 1980s are relatively higher, whereas the 1990s is relatively lower, and since 2000 the recharge amount has slightly increased (Figure 6a).As compared to the 1960s, the

Analysis of the Influence of Climate Change on Karst Groundwater Recharge in the Spring Region
The recharge amount of the Jinci Spring region from the precipitation and river channel seepage in the Fenhe Reservoir-Zhaishang zone can be obtained according to Formula (9).
The annual recharge amount in the Fenhe Reservoir-Zhaishang zone of the Jinci Spring region varies from 29.32 to 140.28 mm (about 11,500-54,900 thousand m 3 ), and the multi-year average value is 59.37 mm (about 23,200 thousand m 3 ).
The karst groundwater recharge amount in the study region has a larger fluctuation at different times.The 1960s, 1970s, and 1980s are relatively higher, whereas the 1990s is relatively lower, and since 2000 the recharge amount has slightly increased (Figure 6a).As compared to the 1960s, the recharge amount of the 1980s decreased 1.42 mm (about 2.11%), that of the 1990s decreased 26.75 mm (about 39.68%), and that of 2000-2012 decreased 8.07 mm (about 11.97%).
The curves of the annual recharge amount and 5a sliding average recharge amount time sequence in the study area from 1960 to 2012 are shown in Figure 6a.The following conclusions were drawn from the analysis.The highest 5a sliding average recharge amount was from 1973 to 1977 at 84.98 mm (converted into 33,200 thousand m 3 ), and the lowest value was from 1997 to 2001 at 37.96 mm (about 14,850 thousand m 3 ).Similar to the precipitation variation trend, the karst groundwater recharge amount in the spring region generally exhibited an initial declining and then subsequent increasing trend [33].
The Mann-Kendall method was used to conduct the trend test, and the test result shows that the statistic Z, β, and the significance level of the annual recharge amount in the study region are, respectively, −1.04, −0.25, and 0.05.The result reveals that the karst groundwater recharge amount test value Z of the study region is less than 0, which means that the karst groundwater recharge amount in the study region from 1960 to 2012 generally exhibits a declining trend.The M-K mutation test shows that the karst groundwater recharge amount exhibited a rising trend before 1967, and then a declining trend predominated; however, it has been rising once again since 2007.According to the intersecting point of the UF k and UB k curves, it can be determined that the mutation existed in 1964 and 2007 (Figure 6b) [34].
The Mann-Kendall method was used to conduct the trend test, and the test result shows that the statistic Z, β, and the significance level of the annual recharge amount in the study region are, respectively, −1.04, −0.25, and 0.05.The result reveals that the karst groundwater recharge amount test value Z of the study region is less than 0, which means that the karst groundwater recharge amount in the study region from 1960 to 2012 generally exhibits a declining trend.The M-K mutation test shows that the karst groundwater recharge amount exhibited a rising trend before 1967, and then a declining trend predominated; however, it has been rising once again since 2007.According to the intersecting point of the UFk and UBk curves, it can be determined that the mutation existed in 1964 and 2007 (Figure 6b) [34].

Conclusions
This study employed the M-K non-parametric test method to analyze the variation trend and mutation characteristic of the annual precipitation, land evaporation, temperature, and karst groundwater recharge in the Jinci Spring region.The main conclusions are as follows: (1) By comparing the relationship between the change in precipitation and the discharge of the spring, the change in precipitation has a small contribution to the change in karst water level in the Jinci Spring region.The reason why the spring dried up is mainly attributed to groundwater over-abstraction.(2) A total of 53 years of statistics reveal that the annual precipitation of the Jinci Spring region generally follows an initial decreasing and then subsequent increasing trend, and the most obvious declining trend occurred in 1990s, which was 111 mm less than that of the 1960s.Likewise, the variation trend of land evaporation is similar to that of precipitation; however, the annual average temperature has maintained an increasing trend.(3) The variation process of the karst groundwater recharge in the Jinci Spring region is generally consistent with that of the precipitation.As compared to the climate of the 1960s, all the climate changes at the other time periods resulted in a reduction in the recharge amount.The climatic change had the strongest influence in the 1990s, which caused the karst groundwater recharge

Conclusions
This study employed the M-K non-parametric test method to analyze the variation trend and mutation characteristic of the annual precipitation, land evaporation, temperature, and karst groundwater recharge in the Jinci Spring region.The main conclusions are as follows: (1) By comparing the relationship between the change in precipitation and the discharge of the spring, the change in precipitation has a small contribution to the change in karst water level in the Jinci Spring region.The reason why the spring dried up is mainly attributed to groundwater over-abstraction.(2) A total of 53 years of statistics reveal that the annual precipitation of the Jinci Spring region generally follows an initial decreasing and then subsequent increasing trend, and the most obvious declining trend occurred in 1990s, which was 111 mm less than that of the 1960s.Likewise, the variation trend of land evaporation is similar to that of precipitation; however, the annual average temperature has maintained an increasing trend.(3) The variation process of the karst groundwater recharge in the Jinci Spring region is generally consistent with that of the precipitation.As compared to the climate of the 1960s, all the climate changes at the other time periods resulted in a reduction in the recharge amount.The climatic change had the strongest influence in the 1990s, which caused the karst groundwater recharge amount in the spring region to be reduced by 26.75 mm, about 39.68% lower than that of the 1960s.The reduction of precipitation is one of main reasons underlying the cutoff of the Jinci Spring.

Figure 2 .
Figure 2. The change trend chart of Jinci Spring discharge and the karst groundwater level of the spring vent.

Figure 2 .
Figure 2. The change trend chart of Jinci Spring discharge and the karst groundwater level of the spring vent.

Figure 3 .
Figure 3. Annual precipitation variation trend and mutation analysis in the Jinci Spring region: (a) annual precipitation variation trend; (b) annual precipitation mutation analysis.

Figure 3 .
Figure 3. Annual precipitation variation trend and mutation analysis in the Jinci Spring region: (a) annual precipitation variation trend; (b) annual precipitation mutation analysis.

4. 1
.2. Temperature Variation The annual average temperature and the 5a sliding average temperature time sequence from 1960 to 2012 in the Jinci Spring region are shown in Figure 4a.According to Figure 4a, the two curves generally exhibit a rising trend over time.Water 2017, 9, 267 9 of 16 4.1.2.Temperature Variation The annual average temperature and the 5a sliding average temperature time sequence from 1960 to 2012 in the Jinci Spring region are shown in Figure 4a.According to Figure 4a, the two curves generally exhibit a rising trend over time.(a) (b)

Figure 4 .
Figure 4. Annual average temperature variation trend and mutation analysis in the Jinci Spring region: (a) average temperature variation trend; (b) annual average temperature mutation analysis.

Figure 4 .
Figure 4. Annual average temperature variation trend and mutation analysis in the Jinci Spring region: (a) average temperature variation trend; (b) annual average temperature mutation analysis.

Figure 5 .
Figure 5. Annual evaporation variation trend and mutation analysis in the Jinci Spring region: (a) annual evaporation variation trend; (b) annual evaporation mutation analysis.

Figure 5 .
Figure 5. Annual evaporation variation trend and mutation analysis in the Jinci Spring region: (a) annual evaporation variation trend; (b) annual evaporation mutation analysis.

Figure 6 .
Figure 6.Annual karst groundwater recharge amount variation trend and mutation analysis in the Fenhe Reservoir-Zhaishang zone of the Jinci Spring region: (a) annual karst groundwater recharge amount variation trend; (b) annual karst groundwater recharge amount variation trend mutation analysis.

6 .
Annual karst groundwater recharge amount variation trend and mutation analysis in the Fenhe Reservoir-Zhaishang zone of the Jinci Spring region: (a) annual karst groundwater recharge amount variation trend; (b) annual karst groundwater recharge amount variation trend mutation analysis.

( 4 )
The results of the Mann-Kendall abrupt change test revealed that the karst groundwater recharge amount in the Fenhe Reservoir-Zhaishang zone of the Jinci Spring region generally exhibits a declining tendency.The karst groundwater recharge amount exhibited a rising trend before 1967.After that, a declining trend took place; however, a rising trend has been occurring since 2007.The mutations existed in 1964 and 2007.

Table 1 .
Annual average precipitation, temperature, and evaporation for each time period in the Jinci Spring region: Statistical Table.

Table 1 .
Annual average precipitation, temperature, and evaporation for each time period in the Jinci Spring region: Statistical Table.

Table 2 .
The trend and abrupt change test of the annual average precipitation, temperature, land evaporation, and karst groundwater recharge.

Table A2 .
The data of precipitation, temperature and observed runoff amount.