Trend Analysis of Rainfall Time Series in Shanxi Province, Northern China (1957–2019)

Changes in rainfall play an important role in agricultural production, water supply and management, and social and economic development in arid and semi-arid regions. The objective of this study was to examine the trend of rainfall series from 18 meteorological stations for monthly, seasonal, and annual scales in Shanxi province over the period 1957–2019. The Mann–Kendall (MK) test, Spearman’s Rho (SR) test, and the Revised Mann–Kendall (RMK) test were used to identify the trends. Sen’s slope estimator (SSE) was used to estimate the magnitude of the rainfall trend. An autocorrelation function (ACF) plot was used to examine the autocorrelation coefficients at various lags in order to improve the trend analysis by the application of the RMK test. The results indicate remarkable differences with positive and negative trends (significant or non-significant) depending on stations. The largest number of stations showing decreasing trends occurred in March, with 10 out of 18 stations at the 10%, 5%, and 1% levels. Wutai Shan station has strong negative trends in January, March, April, November, and December at the level of 1%. In addition, Wutai Shan station also experienced a significant decreasing trend over four seasons at a significance level of 1% and 10%. On the annual scale, there was no significant trend detected by the three identification methods for most stations. MK and SR tests have similar power for detecting monotonic trends in rainfall time series data. Although similar results were obtained by the MK/SR and RMK tests in this study, in some cases, unreasonable trends may be provided by the RMK test. The findings of this study could benefit agricultural production activities, water supply and management, drought monitoring, and socioeconomic development in Shanxi province in the future.


Introduction
Rainfall is one of the most important climate elements that can directly affect the availability of water resources, can be used for diagnosing climate change under the background of global warming [1,2], and influence agricultural production activities and the ecohydrological environment [3][4][5]. Relevant studies have shown that water cycle modification, with rainfall being a key point in the process, is considered the most noticeable consequence of global atmospheric warming [6,7]. In addition, the spatial and temporal distribution of runoff, soil moisture, and groundwater reserves can be influenced by a change of rainfall distribution and the frequency of droughts and floods can also be altered [8]. Therefore, variations and trends around rainfall should be investigated and analyzed for planning and managing of water resources, and also for the purpose of designing suitable plans to cope with floods and drought conditions [3].
In recent years, many studies related to the trends and variability in precipitation under global warming have investigated different regions, such as India [4,5,[9][10][11][12][13], the Mediterranean [14][15][16][17][18], Ethiopia [19], China [1,2,[20][21][22][23][24][25], the Netherlands [26], South Korea [27], Brazil [28], Germany [29], and West Africa [30]. Although various studies on rainfall variation and trend analysis have been investigated in western [31,32], eastern [24], southern [33], north-west [34], south-east [2], and south-west China [23], little investigation has been carried out in the northern regions (e.g., Shanxi province, China). Meanwhile, relevant studies have suggested that various change patterns can be detected for different regions within a country [26,33,35,36]. Liu et al. [37] investigated the spatial and temporal patterns of trends of the precipitation in the Yellow River Basin (YRB) from 81 meteorological stations during the period 1961-2006. They found that most of the precipitation stations showed a downward trend, while only two meteorological stations displayed upward trends. Zhang et al. [33] analyzed the spatial and temporal changes of precipitation structure in the Pearl River Basin based on daily observed precipitation data at 42 rain gauging stations from 1960 to 2005. Decreasing trends in precipitation were found mainly in the middle and upper Pearl River Basin. However, a decreasing number of rainy days was identified almost over the entire basin. Wang et al. [24] analyzed the trend on annual and seasonal rainfall in the Yangtze River Delta, eastern China, using observed data from 14 rain gauging stations during 1961-2016, and found statistically significant upward trends in annual rainfall at all stations. Summer and winter presented increasing trends, while spring and autumn showed decreasing trends. In addition, precipitation extremes under global warming have also been investigated in various regions to better understand the rainfall regime changes due to climate change. Li et al. [38] applied an iterative-based Mann-Kendall trend test and examined the long-term trends in precipitation extremes . Frequency and intensity of precipitation in terms of inter-annual variability was also investigated. They found the frequency and intensity of precipitation extremes in Singapore showed a significant increasing trend. Three potential factors-El Niño-Southern Oscillation (ENSO), global mean temperature, and local temperature-are significantly correlated with precipitation extremes. Chow et al. [39] investigated the impacts of changing frequencies and intensities of extreme weather events on urban settlement, and concluded that the warm environment observed in Singapore and surface dryness have been more intense in recent periods. Change of occurrence of time and length of the extreme precipitation period, as well as influencing factors, were investigated by Guo et al. [40]. They found atmospheric circulation patterns have significant effects on the occurrence time and length of the extreme precipitation period. Changes of extreme precipitation induced by the global mean temperature and ENSO have been investigated by several researchers [41][42][43][44]. Meanwhile, the impacts of ENSO on regional rainfall have been examined around the world. Philippon et al. [45] analyzed the relationship between ENSO and South Africa austral winter rainfall using a long-term daily rainfall database, and found seasonal rainfall (May, June, and July) had a positive correlation with the Niño 3.4 index. Räsänen and Kummu [46] examined the impact of ENSO on large inter-annual variations between large floods and severe droughts. The results showed that both precipitation and discharge are influenced by ENSO. The precipitation and discharge decreased in El Niño years, while they increased in La Niña years. In order to understand the nature of the rainfall variations under climate change, more studies should be conducted for different geographical regions [11].
The non-parametric Mann-Kendall (MK) test and Spearman's Rho (SR) test are commonly employed to detect monotonic trends in time series of climatic and hydrologic variables. However, the results of the test may contain error if significant autocorrelation exists in the data series. To avoid this problem, a Revised Mann-Kendall test has been proposed by Hamed and Rao [47] using a variance correction method to eliminate the effect of autocorrelation in the time series data. The magnitudes of trend slope for monthly, seasonal, and annual rainfall series in many study areas have been estimated by using Sen's slope estimator [4,8,13,20].
Shanxi is a province with serious soil erosion concerns in northern China, and is also a fragile ecological zones which is sensitive to global climate change [48]. The per capita water resources in Shanxi province is only 456 m 3 , accounting for 1/5th of the national average and far below the severe water shortage limit of 1000 m 3 . The average amount of water used by cultivated land is 3255 m 3 /hm 2 , which is only 1/9th of the national average. Therefore, changes of rainfall in Shanxi province will have a negative effect on the water supply, agricultural irrigation, ecological systems, and socioeconomic development [24,25].
Most studies of rainfall changes in Shanxi province have been concerned with annual rainfall; little research has been undertaken on annual, seasonal, and monthly time series data using all available meteorological stations [48,49]. This study attempts to apply the Mann-Kendall (MK) test, Spearman's Rho (SR) test, Revised Mann-Kendall (RMK) test and Sen's slope estimator (SSE) to detect trends of rainfall on the annual, seasonal, and monthly scales in Shanxi province, based on the latest high-quality observed data from 1957 to 2019. Therefore, the main objectives of this study were to: (1) analyze temporal trend in annual, seasonal, and monthly rainfall time series data using the MK/SR and RMK tests; (2) investigate the magnitude of trend line in annual, seasonal, and monthly time series data using SSE; and (3) compare the difference between MK/SR and RMK in detecting the trends of long-term time series rainfall data.

Study Area and Data Availability
Shanxi province is located in the middle of China and lies in the east of the Loess Plateau, which is a typical semi-arid region in northern China. It lies between 34 • 34 N to 40 • 43 N latitude and 110 • 14 E to 114 • 33 E longitude ( Figure 1) and covers an area of 156,700 km 2 . Shanxi has a temperate continental monsoon climate. High temperatures and rainy days occur across the summer, while low temperatures and dry weather conditions characterize the winter. The average annual temperature ranges from 3 • C to 14 • C. The mean annual precipitation is between 350 mm and 700 mm, with 60% of it concentrated in June, July, and August, and the amount of precipitation is greatly affected by the topography. The average elevation of Shanxi province is 1000 m, with the highest elevation being 3059 m and the lowest elevation being 205 m ( Figure 1). There are four seasons according to the climatology, namely spring (March to May), summer (June to August), autumn (September to November), and winter (December to the following February). The location of Shanxi province in China, along with the topography and spatial distribution of meteorological stations, is shown in Figure 1.
Monthly rainfall data of 18 meteorological stations of Shanxi province for 63 years (1957-2019), except for Changzhi station (47 years), were collected from the China Meteorological Data Network. Quality checks and homogeneity were carried out and controlled before these data were released. The monthly rainfall data for the 18 meteorological stations were aggregated to total rainfall at seasonal and annual scale for trend analysis. The geographical location and elevation information of the meteorological stations and the corresponding length of the rainfall record is provided in Table 1. Quality checks and homogeneity were carried out and controlled before these data were released. The monthly rainfall data for the 18 meteorological stations were aggregated to total rainfall at seasonal and annual scale for trend analysis. The geographical location and elevation information of the meteorological stations and the corresponding length of the rainfall record is provided in Table  1.    [50] and Kendall [51], is a non-parametric test that is widely used to detect monotonic trends for time series of hydrology, meteorological, and environment data. A non-parametric test does not require data to be normally distributed and is flexible to outliers in the time series [52]. The test assumes a null hypothesis, H0, of no trend. The alternative hypothesis, H1, is that a monotonic trend is identified (increasing or decreasing). For a given time series X i = x 1 , x 2 , . . . , x n , the Mann-Kendall test statistic S is calculated as where n is the number of data points; x i and x j are the data values in time series i and j ( j > i), respectively; and sign x j − x i is the sign function as The variance is computed as where n is the number of data points; m is the number of tied groups; and t k denotes the number of extent k. Standardized test statistics, Z, are used to detect a significant trend, and expressed as: The monotonic trend can be identified by the values of Z. Negative values of Z indicate downward trends while positive values of Z indicate upward trends. The statistically significant trend test is undertaken at the specific α significance level. When |Z| > Z 1− α 2 , the null hypothesis is rejected, and a significant trend exists in the time series. Significance levels of α = 0.05 and α = 0.01 were used in this study. At the 5% significance level, the null hypothesis of no trend is rejected if |Z| > 1.96, and it is rejected if |Z| > 2.33 at the 1% significance level [3,4,53].

The Spearman's Rho Test
The Spearman's Rho (SR) test [54] is another widely used non-parametric test for detecting monotonic trends in time series data. Yue et al. [55] investigated the power of the MK and SR tests for detecting monotonic trends in hydrological series. The results indicated that these two non-parametric tests have similar power in detecting a trend. For a given time series X i = x 1 , x 2 , . . . , x n , the Spearman's correlation coefficient r SRC is given as . RX i is the rank of the variable X i ; RY is the chronological order of observations and i = 1, 2, . . . , n in series of size n. The test statistic t SRC is given by Equation (6).
The test statistic t SRC follows a t-distribution with degree of freedom ν and significance level α [43]. The null hypothesis of no trend is rejected when |t SRC | > t υ,1−(σ/2) .

Sen's Slope Estimator
The Sen's slope estimator [56], a robust non-parametric method, is used to estimate the magnitude of the trend slope. For a given time series X i = x 1 , x 2 , . . . , x n , with N pairs of data, the slope is calculated as Sen's [57] f (t) = Qt + B where Q indicates the slope, while B is a constant. The slope of Q i in the sample of N pairs of data is calculated: where x j and x k are the values of data pairs at times j and k.
The median of Sen's slope estimator is computed as The Q med sign reflects the data trend reflection, while its value represents the steepness of the trend [53]. A positive value of Q i means that there is an increasing trend, while a negative value of Q i indicates decreasing trend in the time series analysis.

Serial Correlation Effect
In a classical MK test, we assume that the observed time series data are serially independent. However, some significant autocorrelation coefficients may exist in the hydro-meteorological time series data [3,43]. The trend identified by the classical MK test would be increased (decreased) in a time series with positive (negative) autocorrelation coefficients. Generally, lag-1 autocorrelation is used to test for serial correlation in time series data [58]. The lag-1 autocorrelation coefficient is the simple correlation coefficient of the first observations N−1, X t , t = 1, 2, 3, . . . , N − 1 and the next observations, X t+1 , t = 2, 3, . . . , N [59]. The correlation between X t and X t+1 is given as where X = N t=1 X t is the overall mean. Coefficient r 1 is tested for its significance. The probability limits on the correlogram of an independent series of the two tailed test is given below [59].
where N is the sample size and k is the lag. If the value of r 1 is outside the confidence interval given above, the data are assumed to be serially correlated otherwise the sample data are considered to be serially independent [59]. In this study, the autocorrelation function (ACF) plot implemented in Python was used to explore the significant autocorrelation coefficients with varying lagged values (confidence level of 0.05).

The Revised Mann-Kendall Test
In this study, the revised Mann-Kendall test (RMK) proposed by Hamed and Rao [47] was employed to address serial autocorrelation issues in rainfall time series analysis. This modified MK test uses a variance correction approach to improve trend analysis. All significant lags in time series data are considered in the RMK test. Detailed information and calculation methods can be found in previous research [3,8].

Statistical Characteristics of the Rainfall Time Series
The basic monthly rainfall statistics of Shanxi province at 18 meteorological stations during the period from 1957 to 2019, such as the maximum, mean, standard deviation (SD), and coefficient of variation (CV), were analyzed and presented in Table 2. The maximum and mean of the long-term monthly rainfall at 18 meteorological stations in the study area varied from 264.70 mm to 468.60 mm and 31.93 mm to 64.27 mm, respectively. The SD ranged from 38.45 to 70.03 and CV ranged from 109% to 135%, which indicates that there is a great temporal variation for the long-term monthly rainfall for the whole period of 63 years (1957-2019) in Shanxi province. The highest mean value was found at Wutai Shan station, whereas the CV here was lowest among the meteorological stations. These values indicate that the regions covered by Wutai Shan station with greater rainfall had less variability than the regions with relatively lower rainfall. This finding is consistent with the results of Gajbhiye et al. [4]. The basic statistics of seasonal (spring, summer, autumn, and winter) and annual rainfall in the study period of 1957-2019, such as the mean, SD, and CV, were also analyzed, and are shown in Table 3. The mean and SD of annual rainfall at different stations varied from 427.82 mm to 771.28 mm and from 84.27 to 205.95, respectively. There was a low CV at the annual scale, which indicates a less inter-annual variability with varied meteorological stations. In terms of seasonal rainfall, the mean and SD of spring rainfall varied from 59.90 mm to 125.12 mm and from 28.08 to 49.11; summer rainfall varied from 230.87 mm to 458.38 mm and from 68.92 to 145.90; autumn rainfall varied from 84.68 mm to 158.28 mm and from 36.50 to 75.58; and winter rainfall varied from 7.42 mm to 34.69 mm and from 4.5 to 22.44 over the period from 1957 to 2019. Summer had the highest mean rainfall and lowest CV, while winter had the lowest mean rainfall and highest CV. These findings indicate that the regions with greater rainfall have less variability than the regions with relatively lower rainfall. This finding is similar to the results of Gajbhiye et al. [4,11].

Autocorrelation Analysis
In this study, a total of 306 time series (each station contains 12 monthly, four seasonal, and one annual series) generated from 18 meteorological stations were employed to investigate the serial correlation by using the ACF plot in Python. The results of the ACF plot demonstrated that most of the observed series were serially independent. For instance, the time series of rainfall in March at Changzhi and Yuanping stations, in June at Xingxian station, and in September at Yangquan station show slight autocorrelation (Figure 2), which were not expected to influence the results of the trend analysis when using the MK test. Detailed information can be found in the following sections. Statistically significant autocorrelation coefficients were identified in February (first and third lags) and winter (first and second lags) at Wutai Shan station (Figure 3), which resulted in a disproportionate result using a classical MK test [60].

Monthly, Seasonal, and Annual Rainfall Trend
The results of the Mann-Kendall (MK) test and Spearman's Rho (SR) test without considering

Monthly, Seasonal, and Annual Rainfall Trend
The results of the Mann-Kendall (MK) test and Spearman's Rho (SR) test without considering the effect of serial correlation for monthly rainfall series at each station in Shanxi province during the period from 1957 to 2019 are presented in Table 4. The significance level of 10%, 5%, and 1% were

Monthly, Seasonal, and Annual Rainfall Trend
The results of the Mann-Kendall (MK) test and Spearman's Rho (SR) test without considering the effect of serial correlation for monthly rainfall series at each station in Shanxi province during the period from 1957 to 2019 are presented in Table 4. The significance level of 10%, 5%, and 1% were considered at each station for Z-statistics of the MK and correlation coefficient of SR. Both positive and negative trends (either statistically significant or non-significant) depending on the stations, were identified using MK and SR tests for monthly rainfall series at 10%, 5%, and 1% levels of significance. However, negative trends dominated at the monthly time scale. As can be seen from Table 4, there was no significant trend detected in June or September for all meteorological stations. Neither significant nor non-significant trends were detected at the monthly scale using the MK and SR tests at 10%, 5%, and 1% significance levels at Youyu, Datong, Xingxian, Yuanping, Lishi, and Taiyuan stations. There was a significant increasing trend in May at Hequ and Wuzhai stations and a significant decreasing trend in March at Yuncheng and Yangcheng stations at the 10% level using MK and SR, and a significant decreasing trend in March at Yushe and Chagnzhi at the 10% level using MK and 5% level using SR. The Wutai Shan station showed significant decreasing trends in January, February, March, April, November, and December at the 1% level using MK and SR. For the remaining months at Wutai Shan station, non-significant decreasing trends were identified by MK and SR (Table 4). Two stations, Wuzhai and Changzhi, experienced significant increasing trends in February at the 5% level. Except for Yangquan and Yushe stations, neither significant positive nor negative trends were detected in August using the MK and SR tests. In July, only Yuncheng station experienced a significant decreasing trend at the 5% level, while the other stations showed non-significant decreasing trends. These finding are consistent with the findings of Gajbhiye et al. [4], which demonstrated that July in the Sindh Basin (1901-2002) experienced a decreasing trend (significant or non-significant at the 10% level). The results of precipitation trends in Jiangxi province during 1960-2008 indicated that June, July, and August showed upward trends as detected by the MK test, which present a contrasting result compared with our study in Shanxi province [2]. The largest number of stations showing decreasing trends occurred in March with 10 out of 18 stations at the 10%, 5%, and 1% levels. Overall, most stations experienced a decreasing trend (significant or non-significant) for the monthly rainfall series, as assessed using the MK and SR tests, during the period from 1957 to 2019, especially Wutai Shan station.
The results of the MK test and SR test, without considering the effect of serial correlation for seasonal and annual rainfall series at each station in Shanxi province over the period 1957-2019, are provided in Table 5. Across the four seasons, only Wutai Shan station showed significantly decreasing trends at the 10%, 5%, and 1% levels; the other stations did not have significant trends. In summer, all stations exhibited decreasing trends, except Lishi, based on the MK and SR tests. Ye [61] investigated the trend of China's summer precipitation; they showed a significant increasing trend in south China and a decreasing trend in north China. However, no significant trends were identified in the west and China overall. The trend in summer in Shanxi province is consistent with that of north China. In autumn, non-significant upward trends were observed at most stations (13 stations). Trends in spring and winter depended on stations. On the annual scale, only Wutai Shan and Yangcheng stations showed significant decreasing trends at the 1% and 10% levels, respectively. Similarly, Sanikhani et al. [3] reported that most stations in central India during 1901-2010 had decreasing trends on an annual time scale. Sharma and Singh [59] also found that the annual rainfall series over the state of Jharkhand, India, experienced a significant downward trend with a decrease of 14.11%. Wu and Qian [25] reported that annual rainfall shows a significant downward trend in the Wei River Basin and north of the Loess Plateau, where spring rainfall is the major contributor to the decline in annual rainfall.  Numbers in bold indicate significant values at the 10% level, * Significant at the 5% level, ** Significant at the 1% level.
From the trend analysis of the monthly, seasonal, and annual rainfall series above, we can conclude that MK and SR have similar power for detecting the trend of rainfall time series data. This result agrees with the finding of Yue et al. [55]. They demonstrated that these two tests have similar power in detecting monotonic trends in hydrological series.
The Revised Mann-Kendall (RMK) test was used to improve the trend analysis by considering the autocorrelation effect on rainfall time series data over the period of 1957-2019. The results of the RMK test for monthly rainfall series at each station are shown in Table 6. The number of significant trends for all stations have changed when using the RMK instead of the MK test. For instance, the Z-statistic value of Wutai Shan station in February and Datong station in June changed from −1.590 (MK) to −1.722 (RMK) and from 1.596 (MK) to 1.724 (RMK), respectively, which indicates that there is a significant trend at a significance level of 10% when using the RMK test. This is mainly contributed to the variance correction method for eliminating the effect of autocorrelation in rainfall time series. February at Wuzhai station showed a significant upward trend at level of 5% based on MK result, while the Z-statistic value of the RMK test indicated a significant increasing trend at the level of 1%. The Yuncheng station in July experienced a similar change; the Z-statistic value changed from −2.076 (MK test at significance level of 5%) to −2.974 (RMK test at a significance level of 1%). In addition, the values of the Z-statistic in August at Yushe station, in July at Xixian station, and in December at Hequ and Houma stations also changed when using the RMK instead of the MK test. This indicates that decreasing trends occurred more frequently in March, while increasing trends occurred more often in February (Tables 4 and 6).
The results of the RMK test, considering the effect of serial correlation for seasonal and annual rainfall series over the period of 1957-2019 are shown in Table 7. The trend in winter at Wutai Shan station was decreasing at a significance level of 5% when considering the serial correlation effect using the RMK test, while the MK test presented a downward trend at the level of 1%. The values of the Z-statistic based on MK test were −1.376, −1.56, and −1.75, for spring, summer, and annual at Yuncheng station, respectively, while the values of Z-statistic using RMK test were −2.301, −2.285, and −2.391, respectively (Tables 5 and 7). On the seasonal scale, summer for all stations showed a decreasing trend (significant or non-significant). There was no significant trend at the annual scale, except for Wutai Shan, Yushe, Yuncheng, and Yangcheng stations.
It should be noted that the Z-statistic value of Hequ station in December changed dramatically from −0.18 (MK) to −4.289 (RMK; Table 6). In winter, Jiexiu station presented an unreasonable Z-statistic value (10.179) when using the RMK test, which is not consistent with actual trend of the rainfall time series in winter ( Table 7). The Var(s) of MK in December at Hequn station and in winter at Jiexiu station was 27,813.667 and 28,423, respectively, while the RMK was only 48.92 and 116.79, respectively. Although similar trends and slopes obtained from monthly, seasonal, and annual rainfall series in Shanxi province over the period 1957-2019 were found using RMK test, mismatched results were also generated. Therefore, the RMK test may not be suitable for detecting monotonic trend in some cases. More trend identification methods should be considered, as suggested by Patakamuri et al. [52], to obtain accurate trend in hydro-meteorological variables.  Numbers in bold indicate significant values at the 10% level, * Significant at the 5% level, ** Significant at the 1% level.
As indicated by the above analysis, early spring (March) experienced a decreasing (significant or non-significant) trend at most stations. This indicated that the climate of Shanxi province will be drier in the future. March is the beginning of spring plowing in Shanxi province, so a negative influence on agricultural activities, the water supply, the natural environment, public health, and socioeconomic development can be expected. Stations with significant decreasing trends at the 10%, 5%, and 1% level are mainly distributed in the east, south, and south-east of Shanxi province. In the north-east of Shanxi province, only Wutai Shan station showed a significant decreasing trend at the 1% level in March. However, stations distributed in the west and north-west of Shanxi province presented non-significant decreasing trends. Therefore, adaptation pathways for changes of rainfall in Shanxi province should be formulated based on the sensitive regions and time period, such as optimizing the planting structure, building agricultural irrigation projects, and constructing water conservancy projects, and these should be flexible in order to respond to future knowledge of climate change to guards against the potential risks [62,63].

Magnitude of Trend
The magnitudes of the trends for the monthly, seasonal, and annual rainfall series over the period from 1957 to 2019 in the study area were calculated using Sen's slope estimator, and are provided in Table 8. In January, 6 out of 18 stations had no trend, while other stations had slight negative slopes (close to zero). Most of stations in February showed a positive slope or no trend, except Wutai Shan station (−0.104 mm/year), while all stations in March presented a negative slope, especially Wutai Shan station (−0.494 mm/year). There was an increasing slope of the rainfall trend in the monthly series (April to June, September) for most stations. The values of Sen's slope in Table 8 indicate that the falling slope of the rainfall trend existed in July, August, October, and November. There was no significant trend slope identified in December. A negative trend of rainfall slope was observed for all stations in summer. This finding is not consistent with the results of Gajbhiye et al. [4], who concluded that all stations in the Sindh River Basin in India had positive trend slope for the summer season. Huang et al. [2] applied the MK test to detect rainfall trends at the monthly, seasonal, and annual scale during 1960-2008 in Jiangxi province, south-east China. They found that significant increasing trends existed in January, August, winter, and summer.

Impacts of ENSO on Rainfall in Shanxi Province
Previous studies have demonstrated that atmospheric circulation patterns (e.g., ENSO) have an influence on global and regional precipitation [64][65][66][67][68][69][70]. In order to better understand the impacts of extreme climate events on rainfall in Shanxi province, seasonal and annual rainfall anomalies in the years of El Niño and La Niña were analyzed based on the average rainfall over the period of 1981-2010. As indicated by the National Oceanic and Atmospheric Administration (NOAA), El Niño and La Niña events can be identified by using Oceanic Niño Index (ONI), which is the running 3-month mean sea surface temperature (SST) anomaly for the Niño 3.4 region. According to the results of ONI, very strong El Niño patterns occurred in 1982-1983, 1997, strong El Niño patterns presented in 1957-1958, 1965-1966, 1972-1973, 1987-1988, and 1991-1992, and strong La Niña patterns were observed in 1973-1974, 1975-1976, 1988-1989, 1998-1999, 1999-2000, 2007-2008 and 2010-2011 (https://ggweather.com/enso/oni.htm). Rainfall anomalies in seasonal and annual rainfall in the specific El Niño and La Niña years were calculated and are presented in Table 9. As it can be seen from Table 9, both El Niño and La Niña have different degrees of impacts on rainfall. In general, negative annual rainfall anomalies occurred in El Niño years (e.g., 1997, 2015, 1957, 1965, 1972, 1987, and 1991), while positive annual rainfall anomalies presented in La Niña years, excluding 1974, 1999, 2000, 2008, and 2010. From the analysis of ENSO and annual rainfall anomalies, rainfall showed a reduction in El Niño years, and is mainly increased in La Niña years. This is consistent with the statistical analysis of Li et al. [71], which showed that there are many years with reduced annual precipitation in El Niño years, and the probability of negative anomalies is 58.8%, while many years with increased annual precipitation were La Niña years, with a probability of positive anomalies of 56.5%. However, annual rainfall increased in other El Niño years (Table 9), which is in line with results of Li et al. [71].  69.198 Author Contributions: F.G. conceived and designed the framework of this study; F.G. processed the data and analyzed the results; F.G. wrote the paper; Y.W., X.C., and W.Y. reviewed and edited the manuscript and resources for publishing the article. All authors have read and agreed to the published version of the manuscript.