Seasonal and Inter-Annual Variability of Groundwater and Their Responses to Climate Change and Human Activities in Arid and Desert Areas: A Case Study in Yaoba Oasis, Northwest China

: Climate change and human activities have profound e ﬀ ects on the characteristics of groundwater in arid oases. Analyzing the change of groundwater level and quantifying the contributions of inﬂuencing factors are essential for mastering the groundwater dynamic variation and providing scientiﬁc guidance for the rational utilization and management of groundwater resources. In this study, the characteristics and causes of groundwater level in an arid oasis of Northwest China were explored using the Mann–Kendall trend test, Morlet wavelet analysis, and principal component analysis. Results showed that the groundwater level every year exhibited tremendous regular characteristics with the seasonal exploitation. Meanwhile, the inter-annual groundwater level dropped continuously from 1982 to 2018, with a cumulative decline depth that exceeded 12 m, thereby causing the cone of depression. In addition, the monthly groundwater level had an evident cyclical variation on the two time scales of 17–35 and 7–15 months, and the main periodicity of monthly level was 12 months. Analysis results of the climatic factors from 1954 to 2018 observed a signiﬁcant warming trend in temperature, an indistinctive increase in rainfall, an inconspicuous decrease in evaporation, and an insigniﬁcant reduction in relative humidity. The human factors such as exploitation amount, irrigated area, and population quantity rose substantially since the development of the oasis in the 1970s. In accordance with the quantitative calculation, human activities were decisive factors on groundwater level reduction, accounting for 87.79%. However, climate change, including rainfall and evaporation, which contributed to 12.21%, still had the driving force to change the groundwater level in the study area. The groundwater level of Yaoba Oasis has been greatly diminished and the ecological environment has deteriorated further due to the combined e ﬀ ect of climate change and human activities.


Introduction
Groundwater resources play a crucial role in social, economic, and agricultural development, particularly in arid and semi-arid regions where surface water is scarce [1][2][3]. It is not just the vital water

Research Methodology
The MK trend test was used to detect abrupt changes and the trends of climatic factors in this case study, which is a non-parametric statistical method based on rank correlation proposed by Mann in 1945 [28] and developed by Kendall in 1975 [29]. The test has a robust ability because it does not require the sample data to follow any specific distribution and is free from the disturbance of outliers [30]. Thereby, this trend test is widely used for the variation analysis of climatological variables over time. Polemio and Lonigro investigated the variations of monthly climatic time series (rainfall, wet days, rainfall intensity, and temperature) and the annual maximum of short-duration rainfall [31]. Chowdhury et al. [32] evaluated the tendency and step changes in annual and seasonal rainfalls over

Research Methodology
The MK trend test was used to detect abrupt changes and the trends of climatic factors in this case study, which is a non-parametric statistical method based on rank correlation proposed by Mann in 1945 [28] and developed by Kendall in 1975 [29]. The test has a robust ability because it does not require the sample data to follow any specific distribution and is free from the disturbance of outliers [30]. Thereby, this trend test is widely used for the variation analysis of climatological variables over time. Polemio and Lonigro investigated the variations of monthly climatic time series (rainfall, Water 2020, 12, 303 5 of 23 wet days, rainfall intensity, and temperature) and the annual maximum of short-duration rainfall [31]. Chowdhury et al. [32] evaluated the tendency and step changes in annual and seasonal rainfalls over more than 100 years in the South Australian region using MK test. Wu et al. [33] used the MK test to assess the trends in annual and seasonal rainfall since 1950s in Shaanxi Province of China. For details of the MK trend test we refer to Hamed [34] and Jia et al. [35].
This work next applied the wavelet analysis to identify the frequency characteristics and evolution trend of the monthly groundwater level. Wavelet transform is a more effective method than the Fourier transform to analyze multiple time scales of hydrological time series, which can truly reflect the periodicity of the given series and the distribution characteristics of the main cycles [36]. In particular, the discrete Morlet wavelet is simple to develop and requires less computation time for non-stationarity hydrological data such as rainfall [23], runoff [37] and groundwater [38,39]. The discrete form of the Morlet wavelet transform coefficient is defined as where ∆t represents the time interval; a represents the scale factor reflecting the period wavelet; b represents the time factor indicating the shift in time domain; and ψ(t) and f (k∆t) represent the mother wavelet function and the discrete time series, respectively [40]. W f (a, b) in the 2D contour map can describe the data information in different time scales, phase, location, and intensity. In addition, the wavelet variance is calculated by Equation (2). The coefficient squares indicate the distribution and intensity of fluctuating energy at different time scales, which can determine the primary periods of hydrological time series [37].
The quantitative analysis between the influencing factors and groundwater is conducted at the end of the research based on the PCA method. PCA is a multivariate statistical method that can aggregate the impacts of multiple variables into few factors [41]. The basic idea of this technology is that the original variables are transformed into fewer variables through dimension reduction under the premise of loss of little information [42]. In other words, the new variable (principal component) can reflect the original information and is not repeated. PCA can replace the inter-correlated variables by the independent uncorrelated principal components to eliminate the effect of multicollinearity; given these abilities, it is widely applied in the field of groundwater resources (i.e., groundwater level [43], quality [44], and geochemistry [45]). The detailed principle and calculation steps are given in previous study [43].

Dynamic Characteristics of Annual Groundwater Level
With the monthly mean groundwater level from 1982 to 2018 of monitoring well ZB-02 in Yaoba Oasis taken as an example, the monthly average rainfall is found to be the most abundant from July to September, whereas the groundwater level is at the lowest phase of the whole year ( Figure 2a). Hence, the rainfall does not change the overall trend of groundwater level during the year. Meanwhile, the monthly groundwater level is uncorrelated with monthly mean evaporation, temperature, and relative humidity to an extent by the qualitative analysis. Thus, the influence of climatic factors on the monthly groundwater level is negligible (Figure 2b-d). By contrast, the monthly groundwater level obviously drops as the monthly average exploitation quantity increases (Figure 2e). In other words, the variation of monthly groundwater level in the study area is opposite the change in monthly exploitation. Therefore, the annual groundwater level in Yaoba Oasis belongs to the dynamic exploitation type.  As shown in Figure 3, the change fluctuation of monthly groundwater level in the ZB-01 well from 1982 to 2018 is low. The monitoring well lies in the living quarters where the domestic water for residents is supplied by the water supply company in the town. Thus, the exploitation quantity of groundwater is relatively few in the vicinity of the ZB-01 well. The variation curves of the other five monitoring wells are generally similar during the year and the monthly average groundwater level changes regularly as the seasonal exploitation amount varies in different seasons. As a result, the characteristic of annual groundwater level from 1982 to 2018 in the study area can be roughly summarized into four stages [27].
(1) Intermittent irrigation period from January to February: The intensity of groundwater exploitation remains low, and the groundwater is only supplied for residents' domestic water. The groundwater level begins to rise slowly after the end of winter irrigation in the previous As shown in Figure 3, the change fluctuation of monthly groundwater level in the ZB-01 well from 1982 to 2018 is low. The monitoring well lies in the living quarters where the domestic water for residents is supplied by the water supply company in the town. Thus, the exploitation quantity of groundwater is relatively few in the vicinity of the ZB-01 well. The variation curves of the other five monitoring wells are generally similar during the year and the monthly average groundwater level changes regularly as the seasonal exploitation amount varies in different seasons. As a result, the characteristic of annual groundwater level from 1982 to 2018 in the study area can be roughly summarized into four stages [27].
(1) Intermittent irrigation period from January to February: The intensity of groundwater exploitation remains low, and the groundwater is only supplied for residents' domestic water. The groundwater level begins to rise slowly after the end of winter irrigation in the previous year, reaching the highest level of the year in March and tending to be basically stable. The increased amplitude of groundwater level of the six wells is 0.25-0.79 m from January to February. (2) Spring irrigation stage from March to May: The intensity of groundwater exploitation increases.
Farmers begin large-scale cultivation and extract the groundwater for irrigating corn, sunflower, and millet. Groundwater level continues to drop to the lowest level in April, and the decline range is 0.06-1.35 m. The extraction intensity of groundwater is low due to the end of spring irrigation. The groundwater level gradually rises to 0.05-0.14 m during this stage. (3) Summer irrigation stage from June to August: This period lasts for the longest time and includes crop growth, which requires a large amount of groundwater to meet growth needs. The groundwater reaches the maximum intensity, and the exploitation quantity accounts for more than 65% of the whole year. The groundwater level drops rapidly to the lowest value of the year in August, forming a regional cone of depression, and the decline amplitude is 0.67-4.79 m at this period. (4) Intermittent irrigation stage from September to December: As the summer irrigation nears its end, exploitation intensity weakens. The groundwater exploitation remains at a low level during this period, and the groundwater level rises rapidly to 0.74-3.44 m. The increase in the rate of groundwater level fluctuates because of limited-scale winter irrigation. However, the overall trend continues to rise with the amplitude of 0.26-0.88 m from November to December.
The previous studies have indicated that the amount of groundwater recharge is maintained at approximately 31 million m 3 /year, while the total amount of exploitation is maintained at 40 million m 3 /year in Yaoba Oasis [46]. Although the groundwater level slowly recovers during the year, the recharge cannot still meet the exploitation, resulting in the gradual decline in groundwater level year after year. As a result, the groundwater level at the end of the year is always lower than that at the beginning of the year. Farmers begin large-scale cultivation and extract the groundwater for irrigating corn, sunflower, and millet. Groundwater level continues to drop to the lowest level in April, and the decline range is 0.06-1.35 m. The extraction intensity of groundwater is low due to the end of spring irrigation. The groundwater level gradually rises to 0.05-0.14 m during this stage. (3) Summer irrigation stage from June to August: This period lasts for the longest time and includes crop growth, which requires a large amount of groundwater to meet growth needs. The groundwater reaches the maximum intensity, and the exploitation quantity accounts for more than 65% of the whole year. The groundwater level drops rapidly to the lowest value of the year in August, forming a regional cone of depression, and the decline amplitude is 0.67-4.79 m at this period. (4) Intermittent irrigation stage from September to December: As the summer irrigation nears its end, exploitation intensity weakens. The groundwater exploitation remains at a low level during this period, and the groundwater level rises rapidly to 0.74-3.44 m. The increase in the rate of groundwater level fluctuates because of limited-scale winter irrigation. However, the overall trend continues to rise with the amplitude of 0.26-0.88 m from November to December.
The previous studies have indicated that the amount of groundwater recharge is maintained at approximately 31 million m 3 /year, while the total amount of exploitation is maintained at 40 million m 3 /year in Yaoba Oasis [46]. Although the groundwater level slowly recovers during the year, the recharge cannot still meet the exploitation, resulting in the gradual decline in groundwater level year after year. As a result, the groundwater level at the end of the year is always lower than that at the beginning of the year.

Impact of Climatic Factors on Groundwater Level
(1) Trend of temperature An analysis of the temperature data from 1955 to 2018 revealed that the average annual temperature in the study area was 8.33 • C, the lowest temperature was 6.25 • C in 1967, and the highest temperature was 10.02 • C in 1998. As shown in Figure 4a, the variability of inter-annual temperature exhibited a fluctuating rising trend, which was consistent with global warming. The trend equation of temperature was y = 0.0398x + 7.0372 (R = 0.812 **) using the linear trend analysis method, and the warming tendency indicated a decadal increasing rate of 0.398 • C ( Table 1). The MK trend test was used to test the mutation year of temperature for the 0.05 significance level and the ±1.96 critical lines [31]. Figure 5 clearly shows one intersection between the UF and UB curves, indicating that the mutation year of temperature occurred in 1986. Thus, the variability of annual temperature in the study area can be described as two stages: insignificant rise phase from 1955 to 1986 and significant rise phase from 1986 to 2018. The UF value was greater than 0 after 1955 and did not exceed the +1.96 critical line before 1986. That is, the temperature increased insignificantly during this period rising from 7.46 • C in 1955 to 7.99 • C in 1986. Next, the UF exceeded the +1.96 confidence interval after 1986. Thus, the temperature increased significantly rising from 7.99 • C in 1986 to 9.67 • C in 2018. The significant warming in temperature caused the irrigation water that was not absorbed by the crops to be consumed via evaporation, thereby causing drought in the farmland and increasing the irrigation times. This phenomenon further led to the increase in groundwater exploitation and the aggravation of soil salinization in the Yaoba Oasis.
Water 2020, 12, x FOR PEER REVIEW 8 of 24

Impact of Climatic Factors on Groundwater Level
(1) Trend of temperature An analysis of the temperature data from 1955 to 2018 revealed that the average annual temperature in the study area was 8.33 °C, the lowest temperature was 6.25 °C in 1967, and the highest temperature was 10.02 °C in 1998. As shown in Figure 4a, the variability of inter-annual temperature exhibited a fluctuating rising trend, which was consistent with global warming. The trend equation of temperature was y = 0.0398x + 7.0372 (R = 0.812 **) using the linear trend analysis method, and the warming tendency indicated a decadal increasing rate of 0.398 °C (Table 1). (a)   The MK trend test was used to test the mutation year of temperature for the 0.05 significance level and the ±1.96 critical lines [31]. Figure 5 clearly shows one intersection between the UF and UB curves, indicating that the mutation year of temperature occurred in 1986. Thus, the variability of annual temperature in the study area can be described as two stages: insignificant rise phase from 1955 to 1986 and significant rise phase from 1986 to 2018. The UF value was greater than 0 after 1955 and did not exceed the +1.96 critical line before 1986. That is, the temperature increased insignificantly during this period rising from 7.46 °C in 1955 to 7.99 °C in 1986. Next, the UF exceeded the +1.96 confidence interval after 1986. Thus, the temperature increased significantly rising from 7.99 °C in 1986 to 9.67 °C in 2018. The significant warming in temperature caused the irrigation water that was not absorbed by the crops to be consumed via evaporation, thereby causing drought in the farmland and increasing the irrigation times. This phenomenon further led to the increase in groundwater exploitation and the aggravation of soil salinization in the Yaoba Oasis.  (2) Trend of rainfall The average annual rainfall of the study area was 201.95 mm, minimum annual rainfall was 93.2 mm in 1963, and maximum annual rainfall was 330.1 mm in 2003 according to an analysis of rainfall data from 1954 to 2018. As shown in Figure 4b, the annual variation of rainfall greatly fluctuated with the trend equation of y = 0.8186x + 174.93. The increasing rate of rainfall was +0.8186 mm/year and R = 0.276 *. Three crossing points between UF and UB are listed in Table 1, and the mutation years of rainfall were 1989, 2005, and 2009, respectively. In accordance with the rainfall data from 1954 to 2018, the two points of 2005 and 2009 were mutation camouflage points, which were removed. Thus, the rainfall variation could be divided into two stages: insignificant reduction period from 1954 to 1989 and insignificant increase period from 1989 to 2018. The UF was less than 0, but it did not exceed the −1.96 critical line before 1989. Therefore, the rainfall trend had an insignificant decreasing fluctuation from 160.71 mm in 1954 to 144.32 mm in 1989. The UF was greater than 0, but it did not exceed the   Table 1, and the mutation years of rainfall were 1989, 2005, and 2009, respectively. In accordance with the rainfall data from 1954 to 2018, the two points of 2005 and 2009 were mutation camouflage points, which were removed. Thus, the rainfall variation could be divided into two stages: insignificant reduction period from 1954 to 1989 and insignificant increase period from 1989 to 2018. The UF was less than 0, but it did not exceed the −1.96 critical line before 1989. Therefore, the rainfall trend had an insignificant decreasing fluctuation from 160.71 mm in 1954 to 144.32 mm in 1989. The UF was greater than 0, but it did not exceed the +1.96 critical line after 1989. Hence, the increasing trend was insignificant, rising from 144.32 mm in 1989 to 237.6 mm in 2018. Rainfall is not only an essential element of the water cycle in the arid oasis but also an important recharge source for groundwater. It recharges groundwater by means of surface infiltration replenishment. Therefore, following the rule of rainfall change and using seasonal flood to regulate groundwater resources are feasible methods.  Table 1). The intersection between the UF and UB reflected the mutation year, that is, 2007. Thus, the evaporation change consisted of two phases: long-term insignificant reduction from 1954 to 2007 and short-term insignificant increased from 2007 to 2018. The UF was less than 0 after 1954, but it did not exceed the −1.96 critical line before 2007, indicating that the evaporation had an insignificant decline during this period, which decreased from 2400.7 mm in 1954 to 1987.09 mm in 2007. The UF then became greater than 0 but did not exceed the +1.96 critical line. Hence, the evaporation showed an insignificant increase afterwards from 1987.09 mm in 2007 to 2264.72 mm in 2018. According to previous research results, the intense evaporation in the arid region mostly exhausted the flood formed during the rainstorm season before it infiltrated into the groundwater [16,47]. In addition, water diversion and irrigation during the agricultural production increased the groundwater level and caused it to exceed the critical depth, resulting in continuous strong water evaporation in the shallow phreatic water area of the western part of study region, which worsened soil salinization [27,48,49].  Table 1). The intersection between UF and UB reflected that 1980 was the mutation year. The annual relative humidity could be divided into two stages: significant decreasing phase from 1954 to 1980, and insignificant decreasing phase from 1980 to 2018. The UF was less than 0, and it exceeded the −1.96 critical line after 1954. Hence, the relative humidity significantly decreased during this period from 43.27% in 1954 to 36.45% in 1980. The UF was less than 0, but it did not exceed the −1.96 critical line after 1980. Therefore, the relative humidity insignificantly declined from 36.45% in 1980 to 36.01% in 2018. The long-time reduction in relative humidity exacerbated the drought, resulting in increased groundwater extraction in the study area, which would cause a further decline in groundwater levels.

Impact of Human Activities on Groundwater Level
The impact of human activities on groundwater level dynamics in Yaoba Oasis is reflected by exploitation quantity, which is adjusted by the change in irrigated area, population, planting structure, and water-saving irrigation technology [26]. Taking the monitoring well ZB-02 as an example, the groundwater level in the non-irrigation (March) and irrigation (July) periods from 1982 to 2008 roughly declined as the exploitation quantity increased as evidenced in  Table 2. The other five wells showed the same variation law as the ZB-02 well. As a result, the inter-annual change in groundwater level in Yaoba Oasis affected by the exploitation amount can be roughly concluded as the following stages.
(1) Slowly dropping period from 1981 to 1997: In the 1960s, the irrigated area was only 2.40 km 2 , the amount of groundwater exploitation was only 1 million m 3 /year, and the depth of groundwater level was shallow. In 1973, agriculture began to rise rapidly with an irrigated area of 16.67 km 2 ( Figure 6). This rise caused the exploitation quantity to increase sharply to 2.10 million m 3 /year. The irrigated area increased to 25.33 km 2 in 1977 with an exploitation amount of 13.43 million m 3 /year due to the continuous influx of migrants and the expansion of land. However, the groundwater level was generally stable because the recharge of 22 million m 3 /year was greater than the exploitation quantity. The irrigated area then rose to 29.33 km 2 in 1979, and the exploitation increased to 24.21 million m 3 /year. The exploitation amount was slightly larger than the recharge. Therefore, the groundwater level slowly declined. In that year, the local government established the dynamic monitoring network because of the decline in groundwater (Figure 1c). Afterwards, the irrigated area was maintained at approximately 29.33 km 2 , and the exploitation ranged from 23 million m 3 Table 2. The other five wells showed the same variation law as the ZB-02 well. As a result, the inter-annual change in groundwater level in Yaoba Oasis affected by the exploitation amount can be roughly concluded as the following stages.
(1) Slowly dropping period from 1981 to 1997: In the 1960s, the irrigated area was only 2.40 km 2 , the amount of groundwater exploitation was only 1 million m 3 /year, and the depth of groundwater level was shallow. In 1973, agriculture began to rise rapidly with an irrigated area of 16.67 km 2 ( Figure 6). This rise caused the exploitation quantity to increase sharply to 2.10 million m 3 /year. The irrigated area increased to 25.33 km 2 in 1977 with an exploitation amount of 13.43 million m 3 /year due to the continuous influx of migrants and the expansion of land. However, the groundwater level was generally stable because the recharge of 22 million m 3 /year was greater than the exploitation quantity. The irrigated area then rose to 29.33 km 2 in 1979, and the exploitation increased to 24.21 million m 3 /year. The exploitation amount was slightly larger than the recharge. Therefore, the groundwater level slowly declined. In that year, the local government established the dynamic monitoring network because of the decline in groundwater (Figure 1c). Afterwards, the irrigated area was maintained at approximately 29.33 km 2 , and the exploitation ranged from 23 million m 3  The crops were mostly water-consuming crops, such as corn and wheat, due to the unreasonable planting structure [25]. Therefore, the reduction of exploitation quantity plays a positive role in controlling the decline in groundwater level. To reduce the exploitation amount of groundwater resources, the implementation of water-saving irrigation technology, adjustment of agricultural planting structure, contraction of irrigated area, and other methods should be considered.

Variation of Cone of Depression
The large-scale development of groundwater resources has been conducted in the study area since 1978. With the continuous overexploitation of groundwater, unbalanced extraction has destroyed the natural runoff field of groundwater, resulting in the formation of a regional cone of depression. Figure 7a-c clearly present that the cone of depression of groundwater level has changed from the emergence phase to the expansion phase to the recovery phase. (1) The emergence phase: As the decline rate of groundwater level in the east was greater than that in the west, the single local cone of depression had a prototype in the east in July 1985 with the groundwater level in the center of cone of depression at 1280.42 m (Table 3) Although the decline rate of groundwater level slowed significantly, and even the level of individual monitoring wells recovered from the historical minimum, the groundwater system still remained in over-exploitation status at this stage in the study area.
The difference in the cone of depression during the year was reflected in flow rate and area by comparing Figure 7c with Figure 7d. The shape of a cone of depression in the irrigation period was almost the same as that in the non-irrigation period, whilst the flow rate and direction of groundwater runoff had obviously changed. During the irrigation period in July, the variation in groundwater level was more complicated than that during the non-irrigation period in March. In particular, the center of the cone of depression was located in the strong extraction area with a central level of 1274.26 m, which mainly occurred in the middle and eastern parts of study area. The groundwater runoff flowed from the periphery to the center of cone of depression with a hydraulic gradient of approximately 1.91-8.30% . Subsequently, the flow rate of groundwater runoff decreased to 0.0245-0.168 m/day and the hydraulic gradient diminished to approximately 6.7% as the end of the irrigation period. The groundwater level was gradually recovered and stabilized in March of the following year with a central level of 1276.59 m. Meanwhile, the range of the cone of depression was significantly fewer than that during the irrigation period shrinking from 11.26 km 2 in July to 9.84 km 2 in March and the diminutive cone of depression in some areas began to disappear. Although the decline rate of groundwater level slowed significantly, and even the level of individual monitoring wells recovered from the historical minimum, the groundwater system still remained in over-exploitation status at this stage in the study area. The difference in the cone of depression during the year was reflected in flow rate and area by comparing Figure 7c with Figure 7d. The shape of a cone of depression in the irrigation period was almost the same as that in the non-irrigation period, whilst the flow rate and direction of groundwater runoff had obviously changed. During the irrigation period in July, the variation in groundwater level was more complicated than that during the non-irrigation period in March. In particular, the center of the cone of depression was located in the strong extraction area with a central level of 1274.26

Periodic Evolution Characteristics of Groundwater Level
In this work, Morlet wavelet transform was used to analyze the periodicity of the change in monthly groundwater level and the variation tendency of short-term seasonal groundwater level in the six monitoring wells from January 2010 to December 2018. The 2D isogram of the wavelet coefficient can visually present the signal strength of different time scales. The positive and negative changes in contour can further reflect the evolution and abrupt characteristics of the given sequence in near future. The positive value of the contour is represented by a solid line corresponding to the rise phase of the time series. By contrast, the negative value in the isogram is enclosed by a dashed line indicating the reduction phase [50]. As shown in Figure 8 The wavelet variance change in time scales was calculated using Equation (2) (Figure 9). As shown in the figure, the groundwater level of the ZB-04 well had two distinct peaks on the 12 months and 25 months scales. The 12 months corresponded to the maximum peak, indicating that the strongest oscillation occurred in this time scale, which was the first primary period of groundwater level. Therefore, the first and second main periods of groundwater level included 12 and 25 months, respectively. The above two main periods controlled the periodical evolution of groundwater level in the ZB-04 well from January 2010 to December 2018.
On the time scale of 20-35 months, the groundwater level was predicted to be in the rising period from September 2018 to October 2019 and from November 2020 to December 2021 and be in the dropping period from October 2019 to November 2020. On the time scale of 7-15 months, the groundwater level was predicted to be in the dropping period from August 2019 to March 2020 and in the rising period from March 2020 to August 2020. The oscillation intensity of the first main cycle was much higher than that of the second main cycle. Thus, the prediction of the groundwater level on the 7-15 months scale was more accurate. Thereby, the periodicity characteristics of the groundwater level is once a year, and the groundwater level will drop from August to November and rise from December to May every year. As a result, the periodicity characteristics of 12 months calculated by the wavelet variance coincide with the annual variation of groundwater level in Section The wavelet variance change in time scales was calculated using Equation (2) (Figure 9). As shown in the figure, the groundwater level of the ZB-04 well had two distinct peaks on the 12 months and 25 months scales. The 12 months corresponded to the maximum peak, indicating that the strongest oscillation occurred in this time scale, which was the first primary period of groundwater level. Therefore, the first and second main periods of groundwater level included 12 and 25 months, respectively. The above two main periods controlled the periodical evolution of groundwater level in the ZB-04 well from January 2010 to December 2018.
On the time scale of 20-35 months, the groundwater level was predicted to be in the rising period from September 2018 to October 2019 and from November 2020 to December 2021 and be in the dropping period from October 2019 to November 2020. On the time scale of 7-15 months, the groundwater level was predicted to be in the dropping period from August 2019 to March 2020 and in the rising period from March 2020 to August 2020. The oscillation intensity of the first main cycle was much higher than that of the second main cycle. Thus, the prediction of the groundwater level on the 7-15 months scale was more accurate. Thereby, the periodicity characteristics of the groundwater level is once a year, and the groundwater level will drop from August to November and rise from December to May every year. As a result, the periodicity characteristics of 12 months calculated by the wavelet variance coincide with the annual variation of groundwater level in Section 3.1. The wavelet analysis results of periodic change in groundwater level in the six monitoring wells are listed in Table 4. Differences in the periodic variation of groundwater level in the six monitoring wells from January 2010 to December 2018 were consistent and subtle. The groundwater level had almost the same multi-time scales in the phase structure, containing the time scale of 7-15 months and the first primary period of 12 months. The variation of groundwater level in the six wells affected by climatic and human factors had evident periodicity of one year and synchronization. The groundwater level in the study area also had the dominant seasonal characteristics during the year. Therefore, the prediction results showed that groundwater level will be in the declining stage from March to August and will enter the recovery stage from September to February every year.

Quantitative Analysis of Major Influencing Factors of Groundwater Level
The PCA method was used to study quantitatively the impact of climatic and human factors on the groundwater level in Yaoba Oasis in this work. The average rainfall (X1, mm), evaporation (X2, mm), temperature (X3, °C), relative humidity (X4, %), irrigated area (X5, km 2 ), population (X6, number), and exploitation amount (X7, million m 3 ) from 1973 to 2018 were used as the original variables, and the annual groundwater level (Y, m) was taken as the dependent variable. The data of The wavelet analysis results of periodic change in groundwater level in the six monitoring wells are listed in Table 4. Differences in the periodic variation of groundwater level in the six monitoring wells from January 2010 to December 2018 were consistent and subtle. The groundwater level had almost the same multi-time scales in the phase structure, containing the time scale of 7-15 months and the first primary period of 12 months. The variation of groundwater level in the six wells affected by climatic and human factors had evident periodicity of one year and synchronization. The groundwater level in the study area also had the dominant seasonal characteristics during the year. Therefore, the prediction results showed that groundwater level will be in the declining stage from March to August and will enter the recovery stage from September to February every year.

Quantitative Analysis of Major Influencing Factors of Groundwater Level
The PCA method was used to study quantitatively the impact of climatic and human factors on the groundwater level in Yaoba Oasis in this work. The average rainfall (X 1 , mm), evaporation (X 2 , mm), temperature (X 3 , • C), relative humidity (X 4 , %), irrigated area (X 5 , km 2 ), population (X 6 , number), and exploitation amount (X 7 , million m 3 ) from 1973 to 2018 were used as the original variables, and the annual groundwater level (Y, m) was taken as the dependent variable. The data of influencing factors were normalized, and their correlation coefficient matrix U was calculated as listed in Table 5. Correlation coefficients among the seven variables were high, especially between the exploitation and irrigated area. Hence, the seven variables showed multicollinearity which means the independent principal components need to be extracted in the next work to represent those original variables. The eigenvalue and contribution rate of U were listed in Table 6. In general, only the eigenvalue greater than 1 could be used as the principal component. As shown in Table 6, the maximum eigenvalue was 3.023, and its contribution rate was 43.189%. The second eigenvalue was 1.573, and its contribution rate was 22.465%. The third eigenvalue was 1.366, and its contribution rate was 19.511%. The cumulative contribution rate of the first three principal components reached 85.164% (>80%) [51]. In other words, three principal components can extract 85.146% of the original information. In addition, the slope of scree plots curve was noticeably slower reaching to the fourth principal component ( Figure 10). Thus, the original variable can be replaced by the first three principal components (i.e., Z 1 , Z 2 , and Z 3 ).
The eigenvectors corresponding to the eigenvalues are listed in Table 7. The principal component linear equation is given as follows:   As shown in Table 7, the X 2 (evaporation) had the largest coefficient of 0.534 in the first principal component Z 1 . Hence, the evaporation played a leading role in Z 1 . The first principal component Z 1 was defined as the evaporation climatic factor. The 0.765 coefficient of X 7 (exploitation amount) in the second principal component Z 2 was the largest value. Thus, the exploitation played a leading role in Z 2 . In summary, the second principal component Z 2 was defined as the human factor. The X 1 (rainfall) had the largest coefficient of 0.637 in the third principal component Z 3 . Hence, the rainfall played a leading role in Z 3 . Thus, the third principal component Z 3 was defined as the rainfall climatic factor. As listed in the linear representation of the three principal components (Equation (3)), Z 1 , Z 2 , and Z 3 , which are independent of one another and not collinear, can represent most original variable information. Therefore, multivariable linear regression analysis was performed in the next work.
In accordance with the principle of linear regression, the standardized annual groundwater level Y from 1982 to 2018 was taken as the dependent variable. The evaporation factor Z 1 , human factor Z 2 , and rainfall factor Z 3 were taken as the independent variables. The linear regression equation was obtained by SPSS 25 as follows: From Equation (4), the correlation coefficient R and the determination coefficient R 2 were 0.949 and 0.902, respectively. The test value F = 57.993, and the probability P = 0.000 < 0.001 significance level. On this basis, the regression equation performed well and could effectively quantify climate change and human activities on the impact of groundwater level in Yaoba Oasis. The contribution rate of human factor contributed 87.79%, evaporation factor accounted for 8.45%, and rainfall factor accounted for 3.76% according to the analysis on the relative influence percentage of various factors on the change in groundwater level. The influence of the human factor on the groundwater level was evidently much greater than that of the climatic factor. The reason was that the climatic factor affected the overall pattern of the ecological environment while the impact on a few areas was not as evident as the human factor. Therefore, the impact of the human factor had a clear indication of the change in groundwater level with the maximum coefficient of exploitation whereas rainfall and evaporation had a certain relatively weak influence on the groundwater level.
The positive and negative signs indicate the positive and negative effects of influencing factors on the groundwater level [44]. From Equation (4), human activities, such as unreasonable mass exploitation, rapid population growth, and gradual expansion of irrigated area, had a strong negative effect on the groundwater level in recent years, which caused a dramatic drop in groundwater level. Although the influence of rainfall on the groundwater level was much fewer than that of human activities, it had the driving force to increase the groundwater level. As a result, the groundwater level of the study area was mainly affected by human and climatic factors. The main influencing factors are exploitation, rainfall, and evaporation. In particular, the exploitation amount, as an important sensitive factor, aggravated the instability of groundwater system and significantly impacted groundwater circulation. In conclusion, affected by the combination of climate change and human activities, the groundwater level in Yaoba Oasis was greatly reduced, and the ecological environment tended to deteriorate.
In future development programs, the depressurization of exploitation amount in groundwater resources is vital for balancing the groundwater system. Feasible solutions, such as shrinkage of the irrigated area, the adjustment of the crop-planting structure, and the generalization of water-saving measures, can effectively reduce the exploitation of groundwater to achieve the sustainable development of the eco-environment, agriculture, and economy in Yaoba Oasis.

Conclusions
This study analyzed the annual and inter-annual variation characteristics of groundwater levels from 1982 to 2018 in YaoBa Oasis via Morlet wavelet analysis. The results are as follows: (1) The annual variations of groundwater level contained a spring irrigation decline period from March to May, a summer irrigation decline period from June to August, and an intermittent irrigation recovery period from September to February in the following year. (2) The inter-annual variations of groundwater level were divided into four stages: a slowly dropping period from 1981 to 1997, a rapidly declining period from 1997 to 2004, a slowly declining period from 2004 to 2008, and a stable period from 2008 to 2018. The MK test and linear trend method were also used to identify the variation trend of influencing factors in groundwater level. In accordance with the change in climatic factors from 1954 to 2018 in the study area, the temperature showed a significant upward trend, which was consistent with the global warming trend. Furthermore, the rainfall increased slightly, and the evaporation and relative humidity showed an insignificant downward trend. In the past 60 years, human activities, such as exploitation and irrigated areas, had intensified the utilization of groundwater resources in the oasis. Therefore, as a qualitative result, the combination of climatic factors and human activities caused a significant drop in groundwater level.
Moreover, PCA was further proposed to quantify the contribution percentage of climate change and human activities to the groundwater level changes. The calculation showed that human factors accounted for 87.79%, whereas climatic factors contributed 12.21%. Accordingly, the overall variation in groundwater level could be attributed to human activities, which outpaced the impacts of climate change. Thus, exploitation amount could be considered a major factor influencing the characteristics of groundwater level. In summary, the groundwater level dropped as exploitation increased in Yaoba Oasis.
Not only were the temporal and spatial variation characteristics of groundwater levels in a typical oasis in the arid northwest region of China realized, but the influencing degree of climate change and human activity response to groundwater was also quantified in this study. The research results can guide scientific groundwater exploitation planning and management for sustainable development of the oasis and provide novel ideas for similar research in other similar arid oases.
Author Contributions: H.L. and Y.L. were responsible for the original concept and writing the paper. C.Z. and X.Z. processed the data and conducted the program design. B.Z. and J.W. revised the manuscript and shared numerous comments and suggestions to improve the study quality. All authors have read and agreed to the published version of the manuscript.