Abstract

Assessment of the effects of climate change and land use/cover change (LUCC) on the flow regimes in watershed regions is a fundamental research need in terms of the sustainable water resources management and ecosocial developments. In this study, a statistical and modeling integrated method utilizing the Soil and Water Assessment Tool (SWAT) has been adopted in two watersheds of northeastern Tibetan Plateau to separate the individual impacts of climate and LUCC on the flow regime metrics. The integrated effects of both LUCC and climate change have led to an increase in the annual streamflow in the Yingluoxia catchment (YLC) region and a decline in the Minxian catchment (MXC) region by 3.2% and 4.3% of their total streamflow, respectively. Climate change has shown an increase in streamflow in YLC and a decline in MXC region, occupying 107.3% and 93.75% of the total streamflow changes, respectively, a reflection of climatic latitude effect on streamflow. It is thus construed that the climatic factors contribute to more significant influence than LUCC on the magnitude, variability, duration, and component of the flow regimes, implying that the climate certainly dominates the flow regime changes in northeastern Tibetan Plateau.

1. Introduction

Climate change and land use/cover change (LUCC) are overlapped factors and function as two critical drivers that directly influence hydrological processes in various watersheds [1, 2]. This is especially true in arid and semiarid areas, where the ecosystems are fragile and more sensitive to climate change due to the restricted availability of water supply. Researchers have thus dedicated significant efforts towards broadening the understanding of the influence of these two variables on hydrological processes, as well as on the flow regimes [3]. As it is well-known, the increase in the concentration of the greenhouse gases emission is causing global climatic warming [4]. Consequently, the temporal spatial distribution and the magnitude of precipitation, frequency, and intensity of floods and recharge and discharge of streamflow are likely to be altered accordingly [5]. LUCC can alter the associations of streamflow generation, transpiration, interception, and ecosystem conservation [6, 7]. Deforestation, for instance, can produce higher streamflow generation and reduce lag time between precipitation and runoff leading to increase in peak flow, flood damage risks, and severe soil erosion [8, 9]. Conversely, reforestation and intensive agriculture rehabilitation can dramatically decrease the generation of streamflow and trigger an increase in the flow concentration time. Although several sources in literature have reported the interactions of forest change, climate, and hydrological activities [10], how these factors are incorporated into hydrological processes as well as the different flow regimes still remain unclear. It is thus important to gather further insights into the impacts of the climate and LUCC factors on the flow regimes in order to accurately quantify their respective roles in different watershed regions.

Available hydrological impact assessment methods commonly used in literature include the catchment paired experiments, statistics analysis, and measurements with hydrological models [11], but seldom have these been involved in streamflow regime studies [12, 13]. Moreover, the experimental and statistical methods treat the study basin as a black box and rarely examined the complexities of precipitation changes, underlying surface conditions, and the interactions between climate change and the respective hydrological processes [2]. Hydrological models, however, offer a viable framework for conceptualizing and investigating the relationships between climate, underlying surface, and hydrological processes in various categories of time and space [14, 15] and the approach of one-factor-at-a-time (OFAT) has been widely applied [16, 17]. For instance, Karlsson et al. [18] modeled the combined effects of land use and climate changes on the hydrology for a catchment located in Denmark. Zhao et al. [19] evaluated the climate variability and land use influences on green and blue water resource in Weihe River basin, NW China. These methods considered the hydrological processes and their interactions with the environment. It should be noted that hydrological models are more effective as they relate the model parameters directly to the physically observable land surface characteristics [20]. However, the assumption within OFAT method is that in the course of evaluating the influence of a given factor on hydrological processes, the other factors’ effects are not considered. In fact, the other factors not considered are changing concurrently over the entire period of observation which can contribute as an apparent bias in the separation results. Separation of the individual influence of the two factors on hydrological processes thus warrants the need for inputting the applicable status of two factors from the baseline period to the entire period. Yang et al. [21] recommended that a combined statistical and modeling method could be used to resolve this bias and further identified the climate and land cover change impacts on the hydrological processes in the Heihe River. In that study, the results were more reasonable and offered greater accuracy than the traditional OFAT and other statistical methods.

The Yingluoxia catchment (YLC) and Minxian catchment (MXC), utilized in this study, are situated in the northeastern Tibetan Plateau, which are primarily influenced by the northern hemispheric middle latitude Westerlies and the Asia summer monsoon, respectively [22, 23]. The two watersheds are nearly located at the same longitude, but quietly different latitude. It means the further north the watersheds, the less the water vapor content received. Meanwhile, these regions are among the most sensitive zones in terms of climate change [24]. More than 1.3 million people and 266,000 ha of irrigated agricultural land rely on the streamflow from YLC region [25]. This flow is crucial for maintaining the oasis and agricultural ecosystems [26, 27]. The MXC is considered as an important water resource hub for the Yellow River, but it is also a source region addressing serious water shortage area in the Gansu Loess Plateau [28]. Therefore, the quantification of the impacts of climate change and LUCC on flow regimes in these regions is crucial in order to clarify the hydrological response to different climate types and understand the flow regime behaviors with different climate and land cover combinations. The objectives of this study are, thus, () to define a conceptual framework and to propose a modified OFAT method for separating the individual contribution of climate and LUCC to flow regimes, () to improve the separation accuracy and to compare the separation results with climatic conditions, and () to compare and gain further insights and understand the effects of the two major factors on streamflow. The results can provide a beneficial reference for an assessment and proper management of the water resources for the associated policymakers and stakeholders within the YLC and MXC regions.

2. Study Area and Data Collection

2.1. Study Area
2.1.1. The Yingluoxia Catchment (YLC)

YLC, the headwater region of the Heihe River basin, is located in the middle of the Hexi Corridor, China, with a total drainage area of 10,009 km2 lying between 99° and 101°E and 38° and 39°N. The elevation varies from 5058 m to 1668 m (Figure 1). YLC is the main region for the generation of streamflow in the entire basin. Approximately 90% of the water resources in the middle and lower reaches are therefore recharged by the streamflow from the YLC region. The climate, mainly controlled by the northern hemispheric middle latitude Westerlies, is characterized as cold and moist with large spatial and temporal heterogeneity [29]. The average annual precipitation is more than 400 mm and increases by 15.5–16.4 mm for every 100 m increase in the elevation. Nearly 70% of annual precipitation occurs from June to September and less than 10% from November to March. The annual average air temperature varies from −5°C to 4°C, and decreases 0.5°C along with the elevation increase for every 100 m [30]. The total annual runoff is 16.05 × 108 m3 with significant interannual variability [31]. Conditions within the Heihe River have attracted great attention in China due to increasing pressures in the nation, particularly in water resources and the degradation of natural environmental sectors [32, 33].

2.1.2. The Minxian Catchment (MXC)

MXC, the headwater of the Taohe River, is covered by a total drainage area of 14,912 km2 with an annual streamflow of 26.8 × 108 m3, lying between 101° and 105°E and 34° and 35°N, dominated by the Asia summer monsoon. The land cover and topography change from forests and mountains to grasslands and open valleys [28]. The annual average air temperature increases from 1 to 6°C at elevations ranging from 4562 to 2,411 m. The mean annual precipitation decreases from 640 to 560 mm from west to east of MXC region. The dominant climate varies from an alpine cold humid and subhumid climate to a temperate semiarid climate, and the terrestrial vegetation ranges from alpine grasslands and forest to arid grasslands and rain-fed cultivated lands [34].

2.2. Data Collection

The digital elevation model (DEM) (30 m × 30 m), map of soil types (1 : 1000000), and land use maps (100 m × 100 m) of 1980s and 2000s were collected for SWAT modeling. Soil properties were obtained from the Chinese Soil Database of the Institute of Soil Science, while land use properties were also sourced directly from the SWAT model database. Surface runoff data were received from the Scientific Database of the Yellow River Hydrology Service, which included the monthly time series from 1980 to 2010 and the meteorological data were collected from the China Meteorological Administration. The latter source provided crucial information on the daily precipitation, maximum and minimum temperature, sunshine hours, humidity, and wind speed. The long term measured gauge air temperature and precipitation were thus aggregated to determine their monthly values by using the daily datasets. In terms of data quality control, the homogenization of climate and streamflow data has been tested using the standard normal homogeneity test.

3. Methodologies

3.1. Flow Regime Metrics

Six representative streamflow regimes metrics, comprised of the flow magnitudes (MMF, Lw75, and Hg25), variability (CvMF), duration (DMF), and component, were selected to evaluate indices of the impacts of climate change and LUCC on streamflow (Table 1) [25, 35]. It is noteworthy that the base flow index (BFI) is the fraction of the base flow over the total streamflow [16], which is contributed by the groundwater. We have used BFI to reflect the low-flow characteristic of the streamflow component. Due to the restrictions by the simulation performance of the long term daily flow regimes in our study area, it was impossible to describe the fine-detail characteristics of daily flow regimes. Monthly flow metrics would thus be a valid alternative, and this has been widely used in the present water management strategies. However, the monthly flow regime could still capture several characteristics of daily flow regimes [36].

3.2. SWAT Model

SWAT is a physically and mathematically based hydrological model developed by the research centers at US Department of Agriculture [37]. The model requires significant quantities of specialized information, including weather, soil properties, topography, vegetation, and land management practices. The resulting data is used to determine parameters controlling hydrological processes in a given watershed. A detailed description of the model and its input/output variables can be found in the model’s public domain website at http://www.brc.tamus.edu/swat.

SWAT processes hydrological simulation by applying the principle of water balance in each of the hydrologic response units (HRUs) [38], which define areas with identical combinations of surface slope, land use, and soil type [18]. The water balance equation is given as where is the final soil moisture (mm), is the initial soil moisture at the time of (mm), is the precipitation at the time of (mm), is the streamflow at the time of (mm), is evapotranspiration at the time of (mm), is the water flow to the unsaturated zone from the soil profile at the time of (mm), and is the water flow from the watershed from underground at the time of (mm). In this study the SCS (Soil Conservation Service) curve number method is applied for computing streamflow volume. The Muskingum and Penman–Monteith methods are used for flow routing and estimating the potential evapotranspiration.

3.3. Framework for Separating the Effects of Climate Change and LUCC

In a given river basin, the changes in hydrological processes are caused by climate change and LUCC, both of which are assumed to be independent factors [39]. In this study, the hydrological regimes are treated as a function of these two factors. To separate the discrete effects of climate change on the flow regime, we have created a schematic in order to illustrate the variation in hydrological processes along with climate change under differing LUCC scenarios (Figure 2).

Due to the different response relationships between different land cover and hydrological processes, we show the schematic of flow regimes response to land covers and (Figure 2). Note that and are the responses of flow regimes to the climate change (Δ) under land cover conditions and . Normally, and are different, when the land cover change (Δ) is smaller, is closer to . Thus, we use the average of and to denote separate impacts of climate change on hydrological elements (): where , , , and donate the flow regimes value at points A, B, C, and D, respectively. Similarly, the responses of flow regime to LUCC are calculated by applying difference of hydrological component values of land use conditions and under the climate condition of and , denoted as and , respectively. The smaller the change in climate change (Δ), the closer is to . Thus, we accept the arithmetic mean of and as the separate impact of LUCC to hydrological elements () as

The total change in hydrological element equals the sum of the two factors contributions. The total change also can be estimated by the differences of the hydrological elements in the two periods, aswhere represents the total change of a given hydrological element; hydrological elements can serve as statistic mean values over annual, seasonal, or monthly time scales [21].

In this study, we select the period 1982–1995 as baseline situation and 1996–2010 as affection period, along with land use for each period reflected on a separate map. The land use map for the 1980s and 2000s is used to represent the land use patterns for the two periods. A calibrated SWAT model is then applied to each of the four permutations derived from the two time periods and the two land use maps; hereinafter these are referred to as four scenarios. The influences of the LUCC and climate change are quantified by comparing the SWAT outputs of the four scenarios by using the separation equations ((2) and (3)), as follows:()1980s land use and 1982–1995 climate2000s land use and 1982–1995 climate1980s land use and 1996–2010 climate2000s land use and 1996–2010 climate.

3.4. Model Calibration and Validation

The SWAT model is calibrated using historical data from January 1982 to December 1995. The validation period is based on the observed runoff from January 1996 to December 2010. And the warming up period is January 1980 to December 1981. Correlation coefficient (), Nash-Sutcliffe coefficient (NSE), and Relative Bias (BIAS) are used to evaluate the performance of SWAT during calibration and validation [40, 41]. The equations are defined aswhere is the number of data points, is the observed runoff, is the simulated runoff, and and are the mean values of the observed and simulated runoff, respectively. Model performance is considered satisfactory if NSE is greater than 0.5 and BIAS is less than 25% [42].

4. Results and Discussion

4.1. Calibration and Validation of SWAT Model

The simulated monthly streamflow has good correspondence to the observed series during calibration and validation (Figure 3). In accordance with Table 2 and Figure 4, there exist significant correlations between the observed and simulated streamflow at monthly and annual scales. That is, the BIAS is found to be approximately 1.14% and 6.80% for the calibration and −0.15% and 2.7% for the validation period in the YLC and MXC regions, respectively. The correlation coefficient in the YLC region is 0.95 and 0.96 for the monthly streamflow and approximately 0.94 and 0.96 for the yearly streamflow, in the calibration and the validation period, respectively. In the MXC region, the correlation coefficient is 0.90 and 0.94 for the monthly streamflow and approximately 0.91 and 0.96 for the yearly streamflow, in the calibration and the validation periods, respectively. The NSE in the YLC region is 0.91 and 0.90 for the monthly streamflow and approximately 0.88 and 0.76 for the yearly streamflow, in the calibration and the validation periods, respectively. In the MXC region, the NSE are 0.85 and 0.74 for the monthly streamflow and approximately 0.80 and 0.78 for the yearly streamflow during the calibration and the validation periods, respectively. Although for the both watersheds the performance of the SWAT model during the validation period is found to be less accurate compared to the calibration period and the monthly performance generally better than the yearly results, the results of the performance are still considered satisfactory, which indicates that the fundamental rainfall-runoff relationship and water balances including the intra-annual distribution are well captured. We thus aver that these results are within the range of “the very good performance” benchmarks [42].

4.2. Characteristics of Land Use/Cover Changes

In order to investigate the change in land use from 1980 to 2010, we integrated the different sources of land use conducts in 1980s and 2000s and compared the fraction of main land use in two time slices from 1980s to 2000s. Figure 5 shows the fraction of land use in YLC and MXC regions. The land use structure shows slight changes from 1980s to 2000s. Pasture is the dominate land use pattern, occupying about 52% and 61% of total area of YLC and MXC, respectively, followed by barren land and forest land, occupying 23% and 22% of total YLC, and forest land and farm land, occupying 27% and 7% of total MXC, respectively. The other land use types occupy less than 2.5%, which have limited effects on the runoff regimes.

Figure 6 shows the spatial distribution of land use conversion in the YLC and MXC regions. There are approximately 331 km2 in YLC and 176 km2 in MXC region that have changed their land use types occupying about 3.3% and 1.2% of their entire catchment that underwent type conversions, respectively. We mainly focus on the land use conversions that occupied large areas and may lead to significant change in runoff regimes. In the YLC region, the forest increased mainly due to the conversion from the pasture (63.65 km2), which is primarily attributed to the afforestation programs and ecological environment restoration programs in Qilian mountain launched in early 2000s [43]. Pasture area increased results from the conversion from barren land (226.03 km2), which is due to the global warming that led to pasture extending to higher altitude area [33]. The other two conversion types in YLC are pasture to farm land (13.3 km2) and glacier to barren land (28.37 km2), which occupies about 12% of total land use conversion area. The conversions of land use indicate the increase of vegetation dynamics in the YLC region, which may lead to increase in canopy interception and transpiration through vegetation, meanwhile decreasing streamflow. In the MXC region, the three main conversion types are forest to pasture (34.44 km2), pasture to farmland (92.37 km2), and farmland to pasture (37.23 km2) (Figure 6(b)). The forest conversion to pasture mainly occurred due to deforestation at three forest farms. The conversion of pasture to farmland is the largest land use conversion type in the MXC region, which, due to herdsmen, has settled down and grasslands are reclaimed [44]. The conversion of farmland to pasture mainly results from grain for green program implemented in early 2000s in loess area of low stream [34]. The other three conversion types occupy small area: forest to farmland (3.26 km2), pasture to residential (3.59 km2), and farmland to residential (5.89 km2). These conversions are relevant to urbanization and reclamation. From the analysis above, we can conclude that the climate change influences the land use conversion in the YLC region, and human activities influence the land use conversion in the MXC region over the 30 years.

4.3. Observed Changes in Precipitation, Air Temperature, and Streamflow

Figure 7 shows the yearly variation of precipitation, streamflow, and air temperature in the YLC region and MXC for the entire modeling period as well as a linear fitted to the observed values. Similar changes of these elements occurred in the two watersheds. The trends of precipitation and runoff decrease from 1980 to 1995 and increase from 1996 to 2010 in both watersheds. However these decrease or increase trends in YLC are slighter than that in MXC. This indicates that variation of water vapor sources due to climate change is more obvious in MXC than in YLC. This is mainly because the Asia summer monsoon tends to affect the MXC region more significantly than the YLC region resulting in more water vapor received in MXC than YLC region [22]. However, the results of precipitation and streamflow trend magnitudes did not pass the significance test of 0.05, except for the streamflow in the MXC region during the period 1980–1995. The increasing trend of air temperature in the YLC region is more obvious than in the MXC region with rate of 0.0381°C/a () in the YLC compared to 0.0265°C/a () in the MXC region from 1980 to 1995. It is important to note that the significant increasing of air temperature magnitude means significant effects on the watershed evapotranspiration. From the difference of the change rates between precipitation and runoff, we can find that evapotranspiration in YLC is not intensive compared to that in MXC. Based on the above analysis, it is evident that the precipitation and streamflow have changed to a notable extent around 1995 and 1996. In accordance with this, we divided the period 1980–1995 as baseline phase and the period 1996–2010 as the affection phase in our study.

4.4. Effects of Climate Change and LUCC on Streamflow Regimes

Table 3 shows the results of streamflow simulated by the SWAT model under the four hypothetical scenarios as presented in Section 3.3. Compared to , the simulated average annual runoff in appeared to increase by approximately 5.6 mm in the YLC region and appeared to decrease by approximately 8.0 mm in the MXC region. These changes consequently result in, respectively, 3.1% and 4.3% of their total streamflow changes. It represents the combination of climate and land use change effects during the period 1980 to 2000 in the two watersheds. The difference between and and between ) and ) shows the influence of land use change during the two subject periods. The separation equation (3) is used to assess the impact of land use change. The results also suggest that the land use change leads to a decrease in average streamflow by approximately 0.4 mm and 0.5 mm, which then accounts for approximately 7.30% and 6.25% of the annual average streamflow changes in the YLC and MXC regions. The influence of climate change in the two periods can be calculated by the difference between and and between and (. The climate variation increases streamflow by 6 mm, which accounts for about 107.30% of total streamflow change in the YLC region and a decrease in streamflow by approximately 7.5 mm, which accounts for approximately 93.75% of total streamflow change in the MXC region. This indicates that the climate change shows opposite effects on streamflow and the impact of climate exceeds by far that of LUCC in two studied watersheds.

4.4.1. Effects of LUCC on Flow Regimes

In this paper, the simulated results instead of the measured data are used to compare hydrological model outcomes. The flow duration curve (FDC) is applied to the monthly stream flow to further quantify the temporal variation of the hydrological regimes. Flow duration curves are considered as a simple yet comprehensive approach to stream flow analysis that gives a graphical view of the overall variability associated with stream flow over a period of time [16]. Flow duration curves for the YLC and the MXC region are shown in Figure 8 with different time slices. The contrast between and and between and shows the influence of land use change in the two study periods. From Figures 8(a) and 8(b), we can see the differences between the two land use situations in the YLC region. The curves for the land use situation within 2000s are almost under the below of the curves for the land use situation of 1980s. This indicates that the LUCC leads to a decrease in monthly stream flow in the YLC region although the decreasing effects are very small. Figures 8(c) and 8(d) present the differences of LUCC impacts on the flow duration curve in the MXC region. Evidently, the two curves for the land use situation of 1980s and 2000s are almost overlapped. This means the LUCC impacts on the stream flow in MXC are almost negligible. Based on the above analysis, it is evident that the LUCC has limited impacts on flow regimes in both watersheds, primarily attributable to the area of land use changes that are small compared to entire watershed area.

4.4.2. Effects of Climate Change on Flow Regimes

Figure 9 presents the climate change impacts on stream flow under different land use situation in two catchments. The curves for two slices are intertwined in each land use situations and respective catchments. This reflects the uncertainty and complexity of the climate variation within these regions. Compared to the impacts of LUCC, climate change impacts on the stream flow are more clear, indicating that the influences of climate factor on streamflow are predominant in these two study catchments. However, the impacts of climate change and LUCC on stream flow are various in different phases of flow duration curve. Thus, we selected six representative streamflow regime metrics to identify the separate impacts of climate change and LUCC on stream flow in two watersheds.

Due to the combined climate and land use change influences, the six flow regimes appear to have changed accordingly. From Table 4, we can see that climate change leads to all streamflow metrics increasing in the YLC region, except for the decreasing of CvMF. This implies that climate change increases both the high-flow and low-flow, while the variation of streamflow is decreasing. In contrast, climate factors have an opposite influence on streamflow metrics in the MXC region. That is, climate has led to a decrease in all streamflow metrics, except the BFI, indicating that the climate change was able to induce a decrease in the high-flow, low-flow, and variation of streamflow. LUCC cause the MMF, Lw75, Hg25, and CvMF decrease and DMF and BFI increase in both two watersheds. These changes imply that LUCC can cause a decrease in the magnitudes and variations of streamflow while increasing the base flow. Meanwhile, climate change factor has larger contributions to all streamflow metrics compared to LUCC in both watersheds. In accordance with this, we can conclude that all streamflow metrics have a larger contribution from climate change than LUCC, indicating that climate change is more relevant to a shift in the flow regimes. It can thus be averred that the contribution of climate change acts to dominate the change in flow regimes and the contribution of land use deduction is relatively small.

This study attempted to disentangle the effects of climate change and LUCC on streamflow in two selected watersheds and to compare their different effects towards the dominant climate types and land cover conditions. In published literature, there have been some previous research works that have focused on the investigation of climate and land cover changes impacts [21, 26, 4548]. The results of these studies agree with our study with the primary finding that climate change has influenced the streamflow regime more significantly than land use/cover change in the study regions [21, 48, 49]. It is worthwhile to note that our study also emphasize the impacts of climate and land use/cover change on flow regimes and the latitudinal effect of climate on streamflow regimes. The results of this study further confirmed that the climate change offers a dominative effect on the magnitude, variability, duration, and component of streamflow regimes.

5. Conclusions

Deficit of water resources in China has raised the importance of investigating the causality of hydrological processes changes. It is widely accepted that these processes are controlled mainly by climate conditions. Notwithstanding this, it is also believed that they are strongly influenced by LUCC. The latter effect is produced by anthropogenic activities, which, in part, are the causes of climate change. In the present study, the SWAT method is successfully applied to separate the climatic and land cover impacts on the flow regime changes in two watersheds of northeastern Tibetan Plateau. This study has confirmed the usefulness of SWAT model for separating the impacts of climate change compared to those from LUCC on flow regimes.

The integration of the effects of LUCC and climate change was seen to increase annual streamflow by 5.6 mm in the YLC region and a subsequent decrease of annual streamflow by 8.0 mm in the MXC region, of approximately 3.2% and 3.4% of the total streamflow, respectively. It was evident that climate change shows an opposite effect in two watersheds. The results showed that climate change is confirmed to outweigh the impacts of LUCC on streamflow, leading to an increase in streamflow in the YLC and a decline in the MXC region, measured as approximately 107.3% and 93.75% of the total streamflow changes, respectively. The LUCC influence showed changes of approximately −7.3% and 6.25% of total runoff. These results outcomes were based on analysis of data collected over the last 30 years, showing that the climate change has significant impact compared to LUCC on flow regimes in the study region. Climate change is thus more relevant to shifting the flow regimes and can be seen to dominate the hydrological processes changes in northeast Tibetan Plateau.

Evident from study, we conclude that it is critical to apply an integrated approach, which combines the effects of climate and LUCC, for an accurate assessment of flow regimes. Future investigations should thus address the projection of different environmental scenarios for the study area and also take into account the uncertainty inherent in predicting both climate change and LUCC.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This study was supported by National Key R&D Program of China (2017YFC0404302, 2016YFC0400908), the Key Research Program of Frontier Sciences, CAS (QYZDJ-SSW-DQC031), the National Natural Science Foundation of China Project (41601038), and the China Postdoctoral Science Foundation (2015M572620). The authors thank the Scientific Data Center in Cold and Arid Regions and the China Meteorological Data Sharing Service System for providing the data.