Abstract

The changes of both climate and land use/cover have some impacts on water resources. In the Taoer River basin, these changes have directly influenced the land use pattern adjustment, wetland protection, connections between rivers and reservoirs, local social and economic development, and so forth. Therefore, studying the impacts of climate and land use/cover changes is of great practical significance. The Soil and Water Assessment Tool (SWAT) model is employed in this study. With historical measured runoff data and remote sensing maps of annual land use classifications, we analyzed the impacts of climate change on the runoff of the Taoer River. Based on the land use/cover classifications of 1990, 2000, and 2010, we analyzed the land use/cover change over the last 30 years and the contribution coefficient of farmland, woodland, grassland, and other major land use types to the runoff. This study can provide a reference for the rational allocation of water resources and the adjustment of land use structure for decision makers.

1. Introduction

Currently, climate change is a global problem. Climate change intensifies the global hydrological cycle [1]. The fifth climate change assessment report [2, 3] of the United Nations Intergovernmental Panel on Climate Change (IPCC, the Intergovernmental Panel on Climate Change) noted that, with the increasing cumulative greenhouse gas (CO2, etc.) emissions, the average global-surface temperature has been increased during the twenty-first century. Northeast China is located in the high latitudes of the Northern Hemisphere and the eastern end of Eurasia, and this region has shown high sensitivity to global climate change since the 1950s. Statistics from the meteorological department of China show that the phenomenon of rising temperature in this area is markedly greater than the national average [4]. Rising temperatures will accelerate the hydrological cycle through enhanced evaporation. Additionally, the intensity, frequency, and spatial distribution of precipitation will change [5, 6]. Rising temperature has a significant impact on water resources in a river basin [7]. Northeast China is an important commodity grain base, and it is a crucial component of China’s food safety system and agricultural production. Since agricultural water consumption covers approximately 70% of the total water consumption in this area, climate change will also indirectly affect agricultural production [8].

Land use/cover change (LUCC) seriously affects water resources [9, 10]. As the most direct expression of the interaction between human activities and the natural environment, land use/cover change is the interaction that links human activities and natural ecological processes [11]. Land use/cover changes directly affect the condition of water resources and agricultural economic growth. In recent years, frequent human activities, for example, deforestation, afforestation, and farmland reclamation, have led to significant agricultural-land-pattern changes (Henze, 1977; McCann, 1995). Land use/cover changes affect water resources mainly through vegetation interception, evapotranspiration, runoff, surface infiltration, soil moisture status, and so forth, thereby affecting the process of watershed hydrology and water resource cycles [1214] (Shi Xiaoliang et al., 2014). In addition, the rising temperature in Northeast China has led to an annual accumulated increase, which provides favorable climatic conditions for forests and grasslands reclaimed in the large arid regions [11]. In this way, climate change, water resource conditions, and land use/cover changes influence one another.

Both climate change and land use/cover change have an impact on water resources. Therefore, the study of the impacts of climate change and land use/cover change on stream flow within a river basin has become an important topic in hydrology and water resource fields. There are many scholars that have conducted studies on this topic [1518]. Currently, the most widely used method is the application of a hydrological model to quantify and assess hydrological conditions (Xu zongxue et al., 2010) [10]. The SWAT model has been demonstrated to be an effective tool for researching the impact of environmental changes on hydrology and water resources. Many scholars [1923] used the SWAT model to analyze the response of the hydrological cycle under the climate change and land use/cover change. Also some [2426] evaluated the impacts of increased atmospheric CO2 concentration and other potential climate change elements on the water cycle. Some [2628] also used SWAT to study the effects of climate change and land use/cover change on agriculture.

In this paper, the Taoer River basin was introduced as a case study. The data for the basin, including elevation, hydrology, slope, land use, soil type, and land coverage, were obtained by using RS and GIS technology. The SWAT model was employed to assess the effects of three main land use types (arable land, woodland, and grassland) on water resources considering climate change and land use/cover change processes in the Taoer River basin. The contribution coefficients [29] of the three land use types to the runoff depth were calculated. The results can help water conservancy administrative departments in land use planning and water resource management.

2. Study Area

The Taoer River is the largest tributary of the Nenjiang River (Figure 1). It originates in Horqin Right-Front-Country. Its stream length is 563 km, with a drainage area of 4.18 × 104 km2 [30, 31]. The Taoer River basin (TRB) terrain gradually rises from the southeast to the northwest; the upper reaches pass through the mountains of the Inner Mongolia Autonomous Region. The stream gradient is approximately 16.6%, while in the downstream area from Zhenxi to the plains, the stream gradient is approximately 0.2%. The TRB is located in a semiarid area with a temperate continental monsoon climate. The climate is cold, dry, and rainless in winter due to the Mongolia high-pressure system, whereas the climate is hot and rainy in summer due to the West Pacific subtropical high-pressure system. The precipitation, temperature, evaporation, and runoff vary greatly at interannual and intra-annual timescales. The mean air temperature is between −2.5 and 5°C (decreasing from north and south), with a mean annual precipitation of 430–900 mm (decreasing from southeast to northwest), potential evapotranspiration of 800–1,900 mm (decreasing from southeast to northwest), and annual average natural runoff in the TRB of 1600 km3 (i.e., from 1956 to 2000).

The upper mountainous area of the TRB, which is covered with forest and grassland, is the key area of runoff yield. The middle and downstream plain areas, which are covered with grassland and arable land, are an important part of the northeast commercial grain bases. The vegetation types and coverage of the entire basin in 1990, 2000, and 2010 are shown in Figure 2. Analysis of the main land cover type changes over the past 20 years shows that the arable land area increased by 13.94%, whereas woodland decreased by 13.14%, and grassland decreased by 10.23%. In addition, farmland, woodland, and grassland are the major land use types in the TRB, accounting for approximately 90% of the entire basin, followed by unused land and marshes, whereas the share of water and settlements account for a smaller proportion.

3. Applied Methodologies

3.1. Model Selection

To research on the impact of climate change and land use/cover change on hydrology and water resources, hydrological modeling techniques are currently the most widely used, and these methods can be mainly divided into three categories: statistical models, conceptual hydrological models, and distributed hydrological models. Experience statistical models are based on the relationship among the statistical data of runoff, rainfall, and air temperature. The disadvantage of these models is that, for long time series of data, the pure statistical model has some limitations in terms of projecting future water resources. Conceptual hydrological models, such as the Stanford model, the Tank model, and the Xin Anjiang model, are based on hydrological phenomena in the relationship of climate and runoff. The disadvantage of these models is that the river basin is assumed to be integral component, which neglects the spatial heterogeneity caused by the differences of topography, vegetation, and soil.

Distributed hydrological models, such as the SWAT model, TOPMODEL model, and VIC models, are large-scale basin models that have been widely used recently [3234]. As a typical distributed hydrological model, the SWAT model is globally used because the input variables can be easily obtained, it has high computational efficiency, it provides long-term watershed simulation, and it is open sourced. It can be modified based on the actual characteristics. By 2012, more than 1,000 peer-reviewed articles have been published that document its various applications [35, 36].

3.2. Method

In this paper, we use historical runoff data and previous studies on climate to analyze the impact of climate change on runoff in the Tao River. We also analyzed the cover change and its influence on runoff in Taoer River basin, on the basis of the land use classification data in different years from satellite remote sensing images. Then we analyze main types of land use change and runoff depth contribution coefficient, such as farmland, woodland, and grassland. This study provides a reference for the rational allocation of water resources and the adjustment of land use structure.

In order to calculate the contribution coefficient of runoff depth on the main land use types in the Taoer River basin, multiobjective decision-making method is used [29, 37, 38]. The multiobjective decision-making method is a method of decision analysis, developed in the mid-1970s. Multiobjective decision-making generally includes multiobjective decision making and multiattribute decision making. It has a wide range of applications. At present, it has achieved rapid development in theory, method, and application.

When the multiobjective decision-making method is used to calculate the contribution coefficient of runoff depth, the basic objective function can be expressed by the annual average runoff depth formula, as shown in formula (1):where represent the contribution coefficient of runoff depth in cultivated land (), woodland (), grassland (), and construction land () in each unit area (km2), respectively. The units are mm·km−2. is the land use scenario among different years, for example, 1990, 2000, and 2010.

Transforming formula (1) can obtain formula (2):

Different subscripts represent different years.

After the establishment of the multiobjective decision-making function, the function is solved by MATLAB, and we can obtain the runoff depth contribution coefficient of different land use types.

3.3. Data Sources

To set up the SWAT model, the input data required include DEM (Digital Elevation), land use/cover, soil type, and hydrometeorology.

The DEM data (SRTM90m) were obtained from the Chinese Academy of Sciences Computer Network Information Center and the international scientific data mirror sites (http://www.gscloud.cn/). Because the middle and lower reaches of the basin are located in the plains, there will be a larger error when using DEM to extract the river catchment. “Burn in” was applied as a corrective technique by importing a drainage map to extract the river catchment. The 1 : 4,000,000 electronic water system diagram was obtained from the National Fundamental Geographic Center (http://ngcc.sbsm.gov.cn).

Land use data were obtained from the “Western China Environment and Ecological Science Data Center” (http://westdc.westgis.ac.cn) of the National Natural Science Foundation. The original land use/cover types of the data classification method adopted six categories of 25 subcategories. To make the land use types of the study area the same as the land use parameter library, the original ones were reclassified; the code of each type was the same as the standard code of the SWAT model. The classification scheme is shown in Figure 3(a).

Soil data were obtained from the “Black River Planning Data Management Center” (http://westdc.westgis.ac.cn). The original soil data were reclassified by ArcGIS. The various types of soil were converted into those codes that can be identified by the SWAT model. The soil classification scheme is shown in Figure 3(b).

Hydrological and meteorological data were obtained from the China Basin Meteorological Data Sharing Service Network (http://data.cma.cn) and hydrological stations in the Taoer river basin, including the Aershan, Ulanhot, Solon, Baicheng, and Tailai weather stations. The meteorological parameters included daily precipitation, maximum and minimum temperature, relative humidity, wind speed, and solar radiation. Data were also obtained from the Solon, Chersen, Taonan, and Zhenxi hydrological stations. The geographical locations are shown in Figure 1.

3.4. Simulation

According to the Taoer River Drainage Network, the Solon, Charlson, Zhenxi, and Taonan hydrological stations were selected as subbasin outlets, and the entire river basin was divided into 27 subbasins (Figure 4). On the basis of the subbasin, land use data and a soil type map were input into the SWAT model, with a set basin slope rating. Combined with actual observations of the Taoer River basin, land use, soil type, and slope of the threshold values were set at 10%, and the Taoer River was divided into 1,663 hydrological response units (HRUs). After meteorological data and other database files were input to the model, the watershed hydrological simulation process was achieved.

3.5. Parameter Sensitivity Analysis

SWAT-CUP was developed by the Swiss Federal Institute of Water Science and Technology (EAWAG), which specializes in SWAT calibration, validation, and uncertainty analysis. They integrated the continuous uncertainty matching algorithm (SUFI-2), particle swarm optimization (PSO), generalized likelihood uncertainty algorithm (GULE), parameter Solution (ParaSol), Markov Chain Monte Carlo (MCMC), and many other programs. The parameter sensitivity analysis, calibration, verification, and uncertainty analysis can be performed by comparing the results of the SWAT model to those of the SWAT-CUP [39]. Based on the SWAT-CUP program, the SUFI-2 algorithm (Sequential Uncertainty Fitting algorithm) [40] was chosen for the sensitivity analysis of the model parameters.

Snowmelt, runoff confluence, evapotranspiration, and subsurface flow are the main hydrological processes involved in SWAT model simulation (Xia Jun et al., 2003; Liu Changming et al., 2004; Wang Zhonggen et al., 2004). Many parameters may impact the simulation; some of the less sensitive parameters can be removed through the model parameter sensitivity analysis to improve the efficiency of the model calibration. Based on the actual situation of the Taoer River basin, 19 parameters were chosen, and the calibrated values of parameters are shown in Table 1. The sensitivity analysis results are shown in Tables 25 (top ten sensitivity parameters).

3.6. Model Simulation and Verification

Hydrometeorological data from 1988-1989 were used to pretest the model, data from 1990–1993 years were used for model calibration, and data from 1994–1996 were used for model validation. In this paper, the Nash-Sutcliffe efficiency coefficient , the correlation coefficient , and the relative error RE were used to evaluate the adaptation degree between the SWAT model simulation results and the measured values. The computational formulas of , and RE are as follows:where is the measured flow rate, m3/s; is the analog flow rate, m3/s; is the average measured flow rate, m3/s; is the average analog flow rate, m3/s; and is the number of samples.

indicates that the analog value and measured value are equal. is usually between 0 and 1; a higher value indicates a higher coincidence of the simulated and measured values. If is negative, then the credibility of the analog value is less than the measured average. A correlation coefficient value closer to 1 indicates that the analog value and measured value trends are more synchronized; that is, the analog is better. The relative error RE reflects the deviation degree between the overall simulated and measured values. An RE closer to 0 indicates a greater match between the simulated and measured values. For runoff simulations of large basins over long time periods, a combination of these three evaluation criteria is used; when , , and RE < 25%, the simulation result is excellent.

3.7. Analysis of the Simulation Results

Comparison of simulated and measured monthly runoff is shown in Figure 5; the results of the simulation evaluation are shown in Table 6.

As seen from Table 3, the simulation results are fairly accurate at the Solun and Chaersen stations, which are in the upstream region of the river. In the calibration and verification periods, the efficiency coefficients were greater than 0.70; the correlation coefficients were also greater than 0.70; and the relative errors RE were less than 20%. In general, the simulation accuracy is high. At the Zhenxi and Taonan stations, which are also in the upstream region, the simulation results of the calibration period are acceptable. However, in the verification period, the correlation coefficient of the Zhenxi station was 29%, and the efficiency coefficient of the Taonan station was 0.42. These values suggest relatively low accuracy. This result is mainly because the description of the surface runoff yield of the SWAT model is in relatively great detail, but the groundwater process is simply generalized, and the manner of simulating groundwater is one-dimensional and conceptual. Groundwater flow between the different subbasins is thus not considered. Exchange between shallow groundwater and deep groundwater is also lacking in terms of a dynamic description. Therefore, it is difficult to accurately describe the transformation process between the complex surface water and groundwater in the SWAT model [41]. The Taonan and Zhenxi stations are located in the middle and lower plains of the basins. The conversion mode between surface water and groundwater is very complicated. The river recharges groundwater, while the SWAT model only considers the shallow groundwater undirectionally recharging the river. The model therefore has inaccurate assumptions that are inconsistent with observations. These inherent defects can lead to a certain degree of error in the SWAT model simulation of the Zhenxi and Taonan stations.

4. Hydrological Response in the Assumed Future Climate Scenarios

4.1. Future Climate Change

Climate change scenarios mainly consider temperature and precipitation changes. Previous studies have debates in the mechanism of its influence on the river basin. Wang Ni et al.   [42] analyzed historical hydrological and meteorological data of the Taoer River basin between 1958 and 2006. The results show that the average temperature increases, and the general trend of precipitation in the basin is not significant. Jiang et al. [43] noted that the precipitation and temperature of the Taoer River basin showed a rising trend. Xiaofan et al. [44] predicted the future temperature and precipitation in the Songhua River basin (the Taoer River basin is in the Songhua River basin range) with ECHAM5/MPI-OM, which is provided by the Climate Change Center of the China Meteorological Administration and the German Max Planck Meteorological Research Institute, using IPCC SRES 3 emission scenarios (A2, A1B, and B1). The results show that, regardless of the scenario, there is a significant increase in the temperature, but the annual precipitation shows no obvious change. Lu et al. [45] used the Co-Kriging interpolation method to study the spatial distribution of temperature in the Songhua River basin and the linear trend estimation method and Mann Kendall rank correlation method to analyze the temporal variation of the regional temperature. The study showed that annual average temperature in the Songhua River basin significantly increased at a rate of 0.039 degrees per year, which is above the global and national rates.

Under climate change scenarios, to investigate the changes of Taoer River basin runoff, this paper imported the results of the IPCC fifth climate change assessment report for future climate prediction. Considering 1990–1996 as the base period, the change in temperature and precipitation over the next 10 years (2020s), 20 years (2030s), and 30 years (2040s) was modeled (Table 7). Climate change mainly considers two factors, that is, temperature and precipitation, and it is assumed that the monthly distribution of the input factors remains unchanged.

4.2. Hydrological Response Analyses of Future Climate Change

Based on the assumption of future climate change scenarios, the SWAT model was used to estimate runoff in the Taoer River basin on a monthly scale for the future periods. Table 8 lists the runoff in the reference and future periods. The results show that the average annual runoff in the basin is slightly increased over the baseline period (1990–1996) during the next thirty years. However, in the 2040s, the average annual flow is reduced compared with the 2030s. Overall, the average runoff of each station has a tendency of increase and then decrease. In summary, there is a close relationship between runoff and temperature, rainfall, and other climate factors in the Taoer River basin.

5. Analysis of the Impact of Land Use Types on Water Resources

5.1. Runoff Contribution Coefficient Analysis of Different Land Use Types
5.1.1. Land Use of the Basin Controlled by the Hydrological Stations

For each subbasin of a given hydrological station, the vegetation type is shown in Figure 6. As shown in the figure, regardless of each subbasin, cultivated land (AGRL), woodland (FRST), and grassland (PAST) are the main vegetation types. Suolun station and Chaerson station are located in the upper reaches of the Taoer River basin and in mountainous areas, which are at high altitudes. The areas controlled by stations were 5273 km2 and 7182 km2. The area covered by cultivated land, woodland, and grassland accounts for more than 98% of the basin controlled by the Chaersen station. The Zhenxi and Taonan stations are located in the middle and lower reaches of the river basin, which is flat terrain at low altitude. The areas controlled by stations were 18702 km2 and 29154 km2. The area covered by cultivated land, woodland, and grassland accounts for more than 95% of the basin controlled by Taonan station.

5.1.2. Simulation Results of the SWAT Model for Different Land Use Types

Using the calibrated SWAT model, all parameters are kept constant, respectively, and the land use data (1990, 2000, and 2010) are incorporated. According to the model output, the Suolun station, the Chaersen station, the Zhenxi station, the Taonan station, and some other output results were analyzed. The runoff of each station is shown in Table 9.

5.1.3. Analysis of the Simulation Results

Different land use types have different effects on watershed runoff. To assess the three major land use types (AGRL, FRST, and PAST), this study calculates the runoff depth contribution coefficients using the multiobjective decision-making method in MATLAB. The results are shown in Table 10.

Analysis of the results shows that the contribution coefficients of the same land use type are different in the mountains area compared to the plains area.

The contribution coefficients of arable land are all negative; that is, arable land has an inhibitory effect on runoff generation. The runoff decreased with increasing of cultivated area. However, in decreasing cultivated area, the runoff increased. In arable land, plowing makes the upper soil loose, which is conducive to water infiltration, thus inhibiting runoff.

The contribution coefficients of woodland are all positive; that is, woodland enhances runoff generation. As woodland is normally located in mountains or hilly areas, it is difficult for the rainfall to remain in these steep slope areas and the runoff can rapidly form. Furthermore, forested areas are generally unsuitable for farming, the land is barren, the vegetation coverage is low, and the litter and humus layers are small; these characteristics are not conducive to rainfall infiltration. Therefore, the formation of woodland played a facilitating role for runoff generation. The runoff increased with increasing woodland.

For grassland, the contribution coefficients are positive in basins controlled by the Soren and Chersen hydrological stations, while those in basins controlled by the Taonan and Zhenxi hydrological stations are negative. The main reasons are as follows: grasslands in the area controlled by Soren and Chersen station are located in mountainous regions where the rainfall infiltration is not significant; in the downstream plains area, lush grass is not conducive to water flow; and the grassland inhibits runoff generation, in which case the runoff decreased with increasing grassland area.

5.2. Impacts of Land Use Change on Water Resource Analysis
5.2.1. Land Use of the Area between Chaersen and Taonan

The Taoer River basin is located in arid and semiarid area, where water resource shortage is a pressing concern. Based on the impact analysis of land use type changes on runoff, the impact of land use type changes on water demand in the region between the Chaersen hydrological station and the Taonan hydrological station (Figure 7) was also studied. The area of this region is 10,400 km2, and the land use types are shown in Figure 8. Over the past 20 years, the area of arable land increased by 1103 km2, and grassland decreased by 1206 km2; areal changes of other land use types are less obvious.

5.2.2. Impacts of Land Use on Agricultural Water Demand

The water demand of different land use types (such as arable land, grassland, and woodland) is different. In the study area the areal changes of land use types were not obvious, except for arable land and grassland. Therefore, arable land and grassland were the two types that have been studied in terms of the effects of land use change on agricultural water demand.

As the water demands of different farming types are not the same, the water quotas of agricultural water demand of the respective arable land areas were calculated. According to the results of actual investigation and a remote sensing data analysis, paddy acreage accounts for 8.1% of the total arable land, and dry land accounts for the remaining 91.9%. The irrigation quota of the average paddy is 15000 m3/hm2, and the average irrigation quota for dry land is 4125 m3/hm2. Amount of agricultural water demand increased by 552 million m3 over the last 20 years. The growth of grass does not generally require irrigation. Because grassland does not need irrigation, the water demand change of grassland was not calculated in this study. Therefore, with the increase of the cultivated area, the agricultural water demand also showed an increasing trend.

Furthermore, not only the agricultural water demand but also the industrial water demand and domestic water demand have increased, with the gradual increase of socioeconomic development and living standards. Thus, the gap between supply and demand of water resources will be increasingly prominent.

6. Discussion

The simulation, calibration, and validation results of the SWAT model applied to the Taoer River basin show that the SWAT model has strong applicability; this model can be used to simulate the monthly runoff process in this basin. The simulation accuracy of the Solon and Chaerson stations in the upstream region is much higher than that at other stations. Due to the limitations of rainfall data and generalized simple simulation results of groundwater, there were some errors in the simulation. These errors are within the acceptable range for the study purpose of assessing the long-term effects [24].

At present, emissions scenarios (Special Report on Emissions Scenarios, SRES) and a number of climate models are used in future climate predictions. However, the simulation results of different climate models are not the same, and the uncertainty of future climate forecasting is increased [46]. In this study, the future climate change assumption is based on the study results of local climate predictions. As the future climate is uncertain, many other factors will affect hydrology and water resources, in addition to rainfall and temperature. Thus, the predictable results are based on the assumption of the ideal situation.

Research on the relationship between land use and runoff has been extensively conducted. As the relationship between land use and water resources is very complex, topography, soil type, and crop management all impact runoff, and the impact of land use types on runoff differs from region to region [47]. The results of this study can be applied only when the geographical and climatic conditions are similar to this area.

This paper studied the impacts of climate change and land use changes on runoff, and a combined comprehensive analysis will be performed in future studies.

7. Conclusions

Based on the assumption of future climate scenarios, the change on the predicted runoff of the Taoer River compared with that in the base period is not obvious. It increases slightly in the 2020s and 2030s but decreases in the 2040s. The relevant authorities should consider that climate change will affect water resource changes, and the use of local engineering measures and nonengineering measures should be employed to minimize the adverse impacts of water resources on social and economic development.

The land use types in the TRB are mainly farmland, woodland, and grassland. The SWAT model and multiobjective decision-making method were used to quantitatively analyze the contribution coefficients of these three major land use types to the runoff. The results showed that, in the TRB, the contribution coefficient of farmland is negative; that is, farmland inhibits runoff; the woodland contribution coefficient is positive; that is, woodland enhances runoff, while, for grassland, in the upper reaches of the mountainous area, grassland plays a catalytic role in runoff, whereas grassland inhibits runoff in the downstream plains. In the past 20 years, the land use type changes caused an increase in the agricultural water demand, which has intensified the disparity between the water supply and demand.

To fill the gap between the supply and demand of water resources and to ensure sustainable and stable development of the local economy, authorities should address this issue through multiple measures. Construction of rivers and lakes has been proposed, where the Yueliangpao reservoir at the end of Taoer River will be used to connect the Nenjiang River to the Taoer River.

Conflicts of Interest

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

Acknowledgments

This research is supported by the Natural Sciences Foundation of China (51679026, 51209030) and the National Key Research and Development Program of China (2016YFC0400903).