Abstract
During the summer in the western mountainous regions of China (WMR), the disasters such as mountain floods, landslides, and debris flows caused by heavy rain occur frequently, which often result in huge economic losses and many casualties. Therefore, it is of great significance to predict the precipitation accurately in these regions. In this paper, a statistical model is established to predict the precipitation in the WMR using the linear regression statistical method, in which the summer area-averaged precipitation anomaly in WMR is taken as the predictand and the prewinter Niño3 SST is taken as the predictor. The results of the return cross test for the historical years from 1979 to 2008 and independent sample return test from 2009 to 2018 show that this statistical model has a good performance in predicting the summer precipitation in the WMR, especially in the flood years. It has better skill in the prediction of WMR precipitation than the dynamical model SINTEX-F.
1. Introduction
Summer in China is a season when precipitation is most concentrated. More or less precipitation in summer often causes floods, droughts, and other natural disasters, which has large impacts on people in China and is a general concern of China’s meteorological departments and the whole society. Summer precipitation in China has significant characteristics of interannual and interdecadal variation, which is significantly related with ENSO [1–5]. Many studies have regarded ENSO as an important predictor for the prediction of precipitation in China [6–9]. However, some studies show that the correlation between ENSO and precipitation in most regions of China is weak and has interdecadal variation [10–12]. In addition to ENSO, the variation in precipitation is also affected by other factors, such as snow cover over the Tibetan plateau [13], monsoon [14, 15], soil moisture [16], and sea ice [17], which makes the precipitation prediction complicated and requires a comprehensive consideration of various external forcing factors. It is a pleasant surprise that the influence of ENSO on precipitation will be amplified due to the uplift of water vapor by mountain in the western mountainous regions (WMR) on the eastern edge of the Tibetan plateau of China, such as Qinling Mountains, Ta Pa Mountains, and Wushan Mountains [18]. Summer precipitation associated with ENSO is anchored by mountains firmly. It means that the predictability of summer precipitation in these regions has been strengthened greatly due to its close relationship with ENSO.
The Qinling Mountains, Ta Pa Mountains, and Wushan Mountains are large mountains located in the WMR. Among them, the Qinling Mountains are east-west mountains located in Shanxi Province, which are regarded as the geographical dividing line between the northeastern China and southeastern China. The altitudes of Qinling Mountains are mostly between 1000 meters and 3000 meters. The Ta Pa Mountains are the mountains in boundary region of three provinces (Shaanxi, Sichuan, and Hubei). They are in the south of the Qinling Mountains. Between the Ta Pa Mountains and Qinling Mountains lies the Han River Valley. Ta Pa Mountains are in the northwest-southeast direction, and the elevation of them is generally about 2000 meters. The Wushan Mountains are located at the border of three provinces (Hubei, Chongqing, and Hunan) and are continuous clusters of peaks in the northeast-southwest direction. The main peak of Wushan Mountains is located in the Chongqing Province, with an altitude of 2400 meters. Wushan Mountains are regarded as the dividing line between the second and third terrain steps of China.
Affected by the specific geological environment, the natural disaster in these mountainous regions is very active. In the rainy season, disasters such as mountain torrents, landslides, and mud-rock flows caused by rainstorm happen frequently, which always result in vast economic losses and causalities. For example, in the 10 years from 1981 to 1990, six large-scale natural disasters occurred successively in the Qinling and Ta Pa Mountains, causing more than 4,000 deaths and 1.5 million destroyed houses. The direct economic losses reach 56 million dollars [19]. The Wushan Mountains are also suffering from severe geological disasters. According to the statistics of nearly 10 years’ data, there are about 40 geological disasters occurring in the Wushan Mountains every year on average, and about 15,000 people are affected by the disasters. The direct economic losses amount to 2.8 million dollars every year. For example, due to the “8.31” continuous heavy rainfall in 2014, 68 geological disasters occurred. 6,548 households and 28,297 people in Wushan Mountains got affected seriously. The economic losses are about 5.6 million dollars [20]. Natural disasters have become a most important factor restraining the economic development and social progress in the WMR. Therefore, an accurate prediction of precipitation in WMR regions is necessary and of great significance for preventing disasters and reducing losses.
Prediction of summer rainfall in China is extremely challenging because of the atmospheric internal variability [21], diversity in ENSO evolution [22], and so on. The dynamical models usually have poor skills in predicting the summer rainfall in China due to its low resolution and incorrect representation of the complex terrain and some physical processes [23–26]. To address this issue, many studies use physically based statistical model to predict the summer precipitation in China and achieve good results [6];[8, 9, 24, 27]. However, most studies predict the precipitation in North China, or in Yangtze River valley, or in the whole China. No study focuses on the prediction of precipitation in WMR.
Therefore, it is necessary to develop a special statistical model which can focus on the prediction of summer precipitation in WMR. In order to do this, the characteristics of temporal and spatial variation in the summer precipitation in WMR and its relationship with ENSO are analyzed in this paper firstly. Then, a statistical prediction model is established and evaluated. The second part in this paper introduces the dataset and methods in details. The third part is the data analysis and prediction modeling. The fourth part is the comparison of results between the statistical model and a dynamical model. The last part is the conclusions and discussion.
2. Dataset and Methods
The dataset used in this paper includes monthly averaged precipitation data from 756 stations in China and Niño3 SST index. The precipitation dataset comes from the monthly dataset of Chinese surface climate data (V3.0) released by the National Meteorological Information Center, which contains the monthly average precipitation of 756 stations in China. It can be downloaded freely at http://www.nmic.cn/. The Niño3 SST dataset is from the climate prediction center (CPC) of National Oceanic and Atmospheric Administration (NOAA), which can be downloaded freely at https://www.esrl.noaa.gov/psd/data/climateindices/list/.
The statistical prediction model is established using the linear regression statistical method, which can be written as
Here, Y is the predictand (precipitation) and X is the predictor (SST and atmospheric circulation indices). The 30-year data from 1979 to 2008 are used to train the statistical model. Both the return cross test for the historical years from 1979 to 2008 and the independent sample return test of 10 years from 2009 to 2018 are conducted to evaluate the performance of the statistical model.
3. Results from the Statistical Prediction Model
Firstly, the relationship between the summer precipitation in the WMR and ENSO is analyzed. Here, the summer refers the months from June to August (JJA). The winter refers the months from December to February (DJF). Figure 1 is the distribution of correlation coefficients between the summer precipitation and prewinter Niño3 SST index from 1979 to 2008. From Figure 1, it can be found that the summer precipitation is significantly related to the prewinter Niño3 SST in the WMR and the upper reaches of Yangtze River. There are 15 stations with correlation coefficients above 95% confidence level (the black dots in Figure 1), which are mainly distributed in the Qinling Mountains, Ta Pa Mountains, and Wushan Mountains. The detailed information of the 15 stations is given in Table 1. The reason for the strong positive correlation relationship has been proposed by [18]. They showed that, in the post-El Niño years, the winter SST in the Niño3 region is higher than normal, an anticyclonic circulation anomaly exits over the tropical northwest Pacific, resulting in a west and strong western Pacific subtropical high, which further leads to strong warm wet southwesterlies to the western mountains. More precipitation is generated on the windward slope of the mountains due to the topographic uplift. In the post-La Niño years, when the Niño3 SST is lower than normal in winter, it is on the contrary. ENSO has a great impact on the precipitation in the WMR through the monsoon southwesterlies as the medium, which is amplified by the topographic uplift of water vapor in monsoon southwesterlies. Therefore, the correlation between summer precipitation and prewinter Niño3 SST in WMR is more obvious than other regions in China.

The above analysis is based on the summer seasonal averaged data. If this relationship between the mountainous precipitation and ENSO is stable in an intraseasonal time scale, the further correlation analysis is needed. Figure 2 is the distribution of correlation coefficients between monthly precipitation in summer and prewinter Niño3 SST index from 1979 to 2008. It can be found from Figure 2 that the distribution of correlation coefficients changes from month to month. In June (Figure 2(a)), the significant correlation coefficients above 95% confidence level are mainly distributed in the Wushan Mountains and the upper reaches of Yangtze River, which is quite different from that in JJA (Figure 1). It indicates that the topographic uplift of the warm wet southwesterly in the WMR is not so strong in June. In July, the distribution of correlation coefficients becomes very close to that in JJA (Figure 2(b)), indicating strong orographic uplift of warm wet airflow in the WMR. In August (Figure 2(c)), the large correlation coefficients in Wushan Mountains have disappeared and move northward to the Qinling.

(a)

(b)

(c)
The month-to-month change in the correlation between precipitation and prewinter Niño3 SST is mainly related to the intraseasonal revolution of precipitation in the WMR (Figure 3). In June (Figure 3(a)), the area with rainfall amount greater than 200 mm is mainly located in the Wushan Mountains. The monsoon southwesterlies have not reached Ta Pa Mountains and Qinling Mountains. The main rain belt is located in the south China. It is the preflood season in south China in June. Thus, the correlation between June precipitation and prewinter Niño3 SST is not obvious in Qinling Mountains and Ta Pa Mountains. In July, the monsoon southwesterlies stretches northward to the Ta Pa Mountains and Qinling Mountains, and the precipitation increases greatly due to the orographic lift of warm wet southwesterlies in those mountains. Therefore, the correlation between July precipitation and prewinter Niño3 SST is also significantly enhanced in this month. In August, as the monsoon advances northward further, the precipitation in Wushan Mountains decreases greatly, whereas the precipitation in Qinling Mountains and Ta Pa Mountains increases (Figure 3(c)), leading to the disappearance of large correlation coefficients in Wushan Mountains (Figure 2(c)). The significant correlation coefficients above 95% confidence level only concentrate in Qinling Mountains, Ta Pa Mountains, parts of north China, and Taishan Mountain (Figure 2(c)). By comparing the distribution of JJA precipitation (Figure 3(d)) and the distribution of correlation coefficients between them and prewinter Niño3 SST in the WMR (Figure 1), it can be found that the distribution of large precipitation amount is basically consistent with the distribution of large correlation coefficients. The strong correlation between JJA precipitation and the ENSO is mainly caused by precipitation in July and August. ENSO has great influence on the precipitation in mountainous regions through the monsoon southwesterlies as the medium. The intraseasonal variation in correlation between rainfall and ENSO is due to the combined contribution of both the orographic uplift of warm wet southwesterlies and the northward advance of the East Asian monsoon.

(a)

(b)

(c)

(d)
To establish the prediction model for the precipitation in WMR, the characteristics of the spatial and temporal variation in WMR precipitation need to be figured out. EOFs [28, 29] are convenient mathematical constructions to identify the basic characteristics of the precipitation in WMR. The EOF represents the spatial structure, and the principal component (PC) describes the temporal variation in the EOF. Figure 4 shows the first EOF and its corresponding PC of the summer precipitation anomaly based on the 15 stations data in the WMR. The variance contribution of the first EOF is 47%. It can be found from Figure 4 that all the eigenvector values of the first EOF are positive, indicating that the temporal variation in precipitation of the 15 stations in the WMR is basically synchronous. Therefore, it is feasible to take the area-averaged precipitation as the predictand in the prediction of WMR precipitation.

(a)

(b)
Figure 4 also shows that the evolution of the first leading PC is very similar with the evolution of DJF Niño3 SST. They both exhibit an obvious characteristic of quasibiennial oscillation. The correlation coefficient between the first leading PC of summer precipitation in the WMR and the DJF Niño3 SST is 0.64, which further illustrates the close relationship between the summer precipitation in WMR and ENSO.
The interdecadal variation in correlation between summer WMR precipitation and ENSO is also studied in this paper. Figure 5 shows the evolution of 21-year slide correlation coefficients between the area-averaged summer precipitation in the WMR and the prewinter Niño3 SST. As can be seen from Figure 5, the correlation between the summer precipitation in the WMR and ENSO has an abrupt interdecadal change in 1979. Before 1979, the correlation coefficients between them are not significant. However, they become significant after 1979. Many studies have demonstrated that an interdecadal variation in precipitation in China occurred in 1979 due to the interdecadal variation in monsoon and ocean conditions [30–33]. Thus, the interdecadal variation in correlation between the summer precipitation in the WMR and ENSO has close relationship to the interdecadal variation in 1979.

To demonstrate this, the distributions of precipitation before 1979, after 1979, and difference between these two periods are analyzed (Figure 6). It can be found from Figure 6 that the precipitation in North China decreased after 1979. The precipitation in south China, the Yangtze River Basin, and the WMR increased after 1979. The precipitation in China changed from “north flood and south drought” to “north drought and south flood.” With the increase in precipitation in the WMR, the uplifting effect of topography to the warm moist flow on the windward slope also becomes stronger after 1979, which makes the relationship between WMR precipitation and ENSO closer. Therefore, it is better to use the data after 1979 to establish the prediction model rather than the data before 1979.

(a)

(b)

(c)
According to the above analysis, the 30-year data from 1979 to 2008 are used to establish the statistical model using the linear regression method. The summer area-average precipitation in the WMR is used as the predictand. The prewinter Niño3 SST index is used as predictors. The statistical model can be written as
In this equation, y is the summer area-averaged precipitation in the WMR. x is the prewinter Niño3 SST index.
In order to evaluate the prediction skill of the statistical model, both the return cross test for the historical years from 1979 to 2008 and the independent sample test for the years from 2009 to 2018 are conducted. Figure 7 shows the observed and predicted summer precipitation from 1979 to 2018 in the WMR. It can be found from Figure 7 that the observed precipitation shows an obvious characteristic of interannual oscillation. The prediction model can basically reproduce the interannual variation in summer precipitation in the WMR. The correlation coefficient between predicted precipitation and the observation is 0.64, indicating a satisfactory skill of the prediction model. However, it should also be noted that the prediction model has large bias in some years, such as 1980, 1992, 2015, and 2016, indicating that ENSO is not the only factor affecting summer precipitation in the WMR.

Table 2 shows the predicted and observed summer rainfall anomaly and relative bias of independent sample return test from 2009 to 2018. The equation of relative bias can be written as follows:
In this equation, Rp is the predicted precipitation, Ro is the observed precipitation, and is the long-term mean precipitation from 1981 to 2010. It can be seen from Table 2 that the prediction model can reproduce the precipitation anomaly in most years. The signs of the Rp and Ro are consistent in 9 out of 10 years, indicating a good skill of the prediction model. The relative bias is generally less than 10% in most years. If the results with an absolute value of relative bias less than 10% are counted as correct, 8 out of 10 years are correct.
However, the bias of this statistical prediction model in some years is large, which are mostly concentrated on the drought years. For example, in the years of 1988, 1992, 2006, and 2015, the observed precipitation is far less than the average climatic value, and the relative bias is large as 20.9%, 31.7%, 27.0%, and 21.3%, respectively. On the contrary, the statistical model has better performance in predicting the summer precipitation during flood years than that during the drought years. The reason for this lies in that the precipitation is more than normal and the topographic uplift is stronger in the flood years, which strengthens the relationship between the summer precipitation and ENSO and the predictability of the WMR precipitation. This property of the statistical model is very helpful to reduce the losses and prevent the climatic disasters, such as mountain floods and mud-rock flows, in the flood years in western mountainous regions.
4. Results from a Dynamical Model
In addition to the statistical model, the ocean-atmosphere coupled general circulation models (OACGCMs) are wildly used to do the dynamical prediction of the future climate. Capturing the variation in the summer precipitation in the WMR is a notable problem. A dynamical model SINTEX-F is used to do the same prediction of WMR summer rainfall, and the results of SINTEX-F are compared with the results of the statistical model to further evaluate the performance of the statistical model.
The SINTEX-F model is developed from Japan Marine Science and Technology Center (JAMSTEC) [34, 35], which can make seasonal-interannual forecast for the global climate. The original version of the SINTEX-F is based on the European SINTEX model under the EU-Japan collaboration project [36, 37]. The ocean component is the reference version 8.2 of the Océan Parallélisé (OPA 8.2) [38]. Its longitude-latitude resolution is 2° × 2° cos (latitude) with increased meridional resolutions to 0.5° near the equation. There are 31 levels in the vertical direction. The atmosphere component is ECHAM4 [39]. It adopts a high horizontal resolution (T106) of about 1.1° × 1.1°. A hybrid sigma-pressure vertical coordinate is used with 19 levels. The atmosphere model and ocean model are coupled with the Océan Atmosphere Sea Ice Soil (OASIS 2.4) coupler. The initial condition of the atmosphere is provided by the 1 yr run force with the monthly climatological SST. The ocean is started from the Levitus annual mean climatologies with zero velocities. In this study, the SINTEX-F is started at 1 March every year from 1982 to 2018, and 24 months are forecast. The results in summer are used to compare with those of the statistical model.
Figure 8 shows the distribution of correlation coefficients between the predicted summer precipitation by SINTEX-F and prewinter Niño3 SST index from 1982 to 2018. Comparing Figure 8 with Figure 1, It can be found that there are obvious differences between them. The significant positive correlation coefficients between the predicted summer precipitation by SINTEX-F and prewinter Niño3 SST index are mainly distributed in the middle and lower reaches of Yangtze River Valley, whereas the correlation coefficients are not significant in the WMR. This kind of distribution of correlation coefficients predicted by SINTEX-F is not consistent with the observation. The impacts of ENSO on the summer precipitation in the middle and lower reaches of Yangtze River Valley are magnified by the SINTEX-F. However, the relationship between the ENSO and precipitation in WMR is not be represented correctly by the SINTEX-F.

Comparing the predicted precipitation in WMR from SINTEX-F with the observation (Figure 7), it can be found that the SINTEX-F underestimates the WMR precipitation in most years. The average simulated rainfall amount from 1982 to 2018 is 431 mm, whereas the observation is 518 mm. The reason for this bias of SINTEX-F lies in that the resolution of the GCM is low, and it cannot represent the detailed topography in the WMR and the topographic uplift of the water vapor well. Therefore, the performance of the SINTEX-F in predicting the interannual variation in summer precipitation in the WMR is not as good as the statistical model. The correlation coefficient between the predicted precipitation from SINTEX-F and the observation from 1982 to 2018 is 0.21, which is much lower than that of the statistical model (0.62) and cannot pass the significance test. The root mean square error of the predicted precipitation by SINTEX-F is 155 mm, which is larger than that of the statistical model (72 mm).
However, it can be noted from Figure 7 that the predictions of SINTEX-F are better than those of the statistical model in some years. The statistical model is not absolutely better than the dynamical model. Therefore, using the ensemble mean results of dynamical models and statistical models to predict the precipitation in the WMR is the recommended method.
5. Conclusions
Based on the observation data of 756 stations in China and Niño3 SST index, the relationship between the summer precipitation in the WMR with ENSO has been analyzed. The basic characteristics of the spatial structure and temporal variation in the summer precipitation in the WMR are analyzed using the EOFs. A statistical model is established to predict the summer precipitation in the WMR. The conclusions are as follows.
5.1. Correlation between the Summer WMR Precipitation and ENSO
It is found that the summer precipitation in the WMR has significant positive correlation with the prewinter Niño3 SST. The significant positive correlation coefficients are mainly distributed in the Qinling Mountains, Ta Pa Mountains, and Wushan Mountains.
The correlation between summer precipitation and ENSO in the WMR has a characteristic of intraseasonal variation, which is mainly caused by the northward advance of East Asian summer monsoon. The ENSO imposes its influence on the summer precipitation in the WMR through the summer monsoon southwesterlies as the medium. Therefore, both the northward advance of summer monsoon and the orographic uplift of the warm moist air flows result in the intraseasonal variation in the correlation between precipitation and ENSO.
The correlation of summer precipitation in the WMR and ENSO also has an obvious characteristic of interdecadal variation. The correlation between the summer WMR precipitation and ENSO is not significant before 1979, whereas it becomes significant after 1979. This interdecadal variation in correlation between summer precipitation in WMR and ENSO is related to the well-known interdecadal change of climate occur in 1979. The precipitation in the WMR increased after 1979. With the increase in precipitation in the WMR, the uplifting effect of topography to the warm moist wind on the windward slope also becomes stronger after 1979, which makes the relationship between WMR precipitation and ENSO closer.
5.2. Statistical Model
The area-averaged summer precipitation in the WMR is taken as the predictand. The prewinter Niño3 SST is used as the predictor. Then, the statistical model is established using the linear regression method. The results of the return cross test for the historical years from 1979 to 2008 and the independent sample test for the years from 2009 to 2018 show that this statistical model has a good performance in predicting the summer precipitation in the WMR, with low relative bias less than 10% in most years.
The good performance of this statistical model in predicting the summer precipitation in the WMR is mainly due to the stronger correlation between the WMR precipitation and ENSO than that in other regions in China. It is well known that ENSO has great impacts on the precipitation in China. However, the correlation between ENSO and precipitation in most regions of China is weak because the summer precipitation is also affected by other factors, such as snow cover over the Tibetan plateau, soil moisture, and sea ice, which makes the prediction complicated and requires a comprehensive consideration of various external forcing factors. However, in the WMR, the influence of ENSO on the WMR precipitation is amplified due to the topographic uplift of water vapor by the mountains. The statistical model can achieve good results using only prewinter Niño3 SST as the predictor. Because the ENSO often reaches its mature phase in winter, we can predict the precipitation in WMR several months in advance.
In addition, another advantage of this statistical model is that the prediction results in the rainy years are more accurate than that in the drought years. This property of the statistical model can be very helpful for reducing the losses of mountain floods, mud-rock flows, and other climatic disasters in the flood years of western mountainous region.
Compared the statistical model with the dynamic model SINTEX-F, it is found that the statistical model in this study has better performance in predicting the precipitation in the WMR. The SINTEX-F magnifies the impacts of ENSO on the precipitation in the middle and lower reaches of Yangtze River Valley and cannot represent the positive correlation between the WMR precipitation and ENSO. This defect leads to a lower skill of the SINTEX-F in predicting the WMR precipitation during summer. However, the statistical model is not absolutely superior to the dynamical model. In some years, the predicted rainfall amounts of SINTEX-F are more accurate than the statistical model in some years. Therefore, the ensemble method which can consider the advantages of both dynamical models and statistical models is recommended.
Data Availability
The precipitation dataset comes from the monthly dataset of Chinese surface climate data (V3.0) released by the National Meteorological Information Center, which contains the monthly average precipitation of 756 stations in China. It can be downloaded freely from http://www.nmic.cn/. The Niño3 SST dataset is from the climate prediction center (CPC) of National Oceanic and Atmospheric Administration (NOAA), which can be downloaded freely from https://www.esrl.noaa.gov/psd/data/climateindices/list/.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was sponsored jointly by the National Natural Science Foundation of China (41675077, 41975090, and 41605072), National Key Research and Development Project (2018YFC1505803 and 2018YFC1505900), Scientific Research Program of National University of Defense Technology (ZK17-03-22), and Jiangsu Collaborative Innovation Center for Climate Change.