Spatiotemporal Evaluation of Blue and Green Water in Xinjiang River Basin Based on SWAT Model

: Poyang Lake is the largest freshwater lake in China. As an important tributary of Poyang Lake, Xinjiang River has an important inﬂuence on the water ecology and water resources of the Poyang Lake basin. Based on the hydrological simulation of the SWAT (Soil and Water Assessment Tool) model, the spatiotemporal distribution and evaluation of the blue and green water during the period (1982–2016) in the basin were explored by the Mann–Kendall test, precipitation anomaly percentage, and scenario simulation. It is found that the SWAT model presents a satisfactory perfor-mance in runoff simulation of the basin. The multi-year average blue water in the Xinjiang River basin is 1138 mm, and the green water is 829 mm, with a green water coefﬁcient of 0.42. The amount of blue water in wet years is about 1.5 times that in normal years and 2.4 times that in dry years. Compared with the green water, the blue water of the basin is more sensitive to the variations in precipitation. In spatial distribution, the blue and green water in the middle of the basin is obviously more than those in other parts of the basin. During the study period, the blue water in the basin shows a slight decreasing trend, and the green water shows a signiﬁcant decreasing trend. It is also found that climatic factors have a greater inﬂuence on the trend of blue and green water than land use, and the decrease in precipitation is the dominant cause for the trend of blue and green water. the in The evapotranspiration in the The results of the scenario simulation found that on the change trends in the green water compared to land use and are the cause of the temporal variation characteristics of blue and green water.


Introduction
Water is an indispensable resource to sustain human life and development. With the increasing population and rapid socio-economic development, human consumption of freshwater resources is increasing, and the imbalance between supply and demand of water resources not only restricts the balance of the ecosystem but also affects the sustainable socio-economic development and, therefore, is drawing increasing attention [1,2]. The current studies pay more attention to the surface streamflow, which can be directly used, but ignore the evaluation of ecological water [3,4]. In order to evaluate the role of water resources in the ecosystem more comprehensively, Falkenmark put forward the concept of blue water and green water and brought them into the consideration of water resources management [5]. The concept of blue and green water broadens the vision of water resources evaluation and research, especially the proposal of green water, which closely links the whole hydrological cycle system with the ecological cycle system [6,7]. In 2006, the concept of blue and green water was further improved by Falkenmark [8], which classified surface runoff, subsurface runoff, and interflow as blue water, classified the evapotranspiration of the land surface, vegetation, water area, and soil as green water flow, and classified the water content of the soil as green water storage. The sum of green water flow and green water storage is the so-called total green water. As the main part of water consumption, green water is a more active part of the land and water cycle system than blue water, and it plays an extremely important role in maintaining the balance, health, and stability of the global ecosystem [9,10].
As water consumption increases, it is projected that 1.8 billion people could face water shortages by 2025 [11,12]. Quantitative assessment of the availability and vulnerability of blue and green water plays a critical role in sustainable water resources management [13]. For example, Veettil et al. established the SWAT model in the Savannah River Basin of the United States and then quantitatively evaluated the water security of the basin based on the concept of blue and green water [12]. Rodrigues et al. proposed a watershed water security assessment method by quantifying the availability of blue and green water [14].
In the past few decades, human activities have had a strong impact on climate conditions, resulting in significant changes in precipitation, temperature, and other climate conditions, which would further affect the regional hydrological cycle and aggravate the shortage of blue and green water [15,16]. Meanwhile, the high intensity of human activities has changed the land use and land cover, resulting in significant changes in hydrological processes, such as infiltration and evaporation, which would also induce the imbalance between the supply and demand of the blue and green water [17,18]. A quantitative analysis of the impact of climate conditions and land use on blue and green water is of great significance for formulating scientific and reasonable water resources management methods [19,20]. Zuo et al. quantitatively evaluated the availability of blue and green water in different spatial scales of the Weihe River basin in China [21]. Lyu et al. conducted a detailed study on the spatial and temporal distribution characteristics of blue and green water in the Dongjiang River basin based on the SWAT model [3]. These studies on blue and green water can provide a valuable theoretical basis for regional water resources planning and management [22][23][24][25].
Poyang Lake, which is located on the southern bank of the middle and lower reaches of Yangtze River, is the biggest freshwater lake in China. It plays an irreplaceable role in regulating floods of Yangtze River, ensuring water supply for its surrounding areas, and providing a habitat for many endangered species [26]. There are many research studies on Poyang Lake, but they mainly focus on water quantity changes and the ecological environment [27][28][29]. Xinjiang River is an important tributary of Poyang Lake, and a change in its water resources can have a great influence not only on Poyang Lake but also on the regional ecological environment and economic development. However, the current research on water resources in the Xinjiang River basin is mainly focused on surface runoff and sediment [30], while the blue and green water in the basin is seldom addressed. Given what was mentioned above, we established the SWAT model in the Xinjiang River basin. Based on the simulation results of the model, we conducted a quantitative evaluation of the volumes and spatial distribution characteristics of the blue and green water on average in multiple years and in the typical years, and we also carried out an analysis of the variation trend of the blue and green water. Meanwhile, we set up a scenario simulation and studied the impacts of meteorological changes and land use changes on the blue and green water. It is hoped that the research results can provide a scientific basis for the watershed management and water resources planning in the Poyang Lake basin and Xinjiang River basin.

Study Area
Xinjiang River is located in the northeast of the Poyang Lake basin. The longitude and latitude are between 116 • 38 E-118 • 36 E and 27 • 33 N-28 • 59 N, respectively ( Figure 1). It has a total length of 313 km and a catchment area of about 17,600 km 2 . Meigang hydrological station controls a catchment area of 15,319 km 2 , accounting for about 87% of the total area of the Xinjiang River basin. The upstream part of the Xinjiang River basin is dominated by hills, the midstream part of the basin is a bason, downstream of the basin are alluvial plains, and the elevation of the whole basin is about 16-2136 m. The basin is located in the subtropical monsoon humid climate zone, with a mild climate and abundant rainfall. The annual average temperature is 17.6 • C, and the annual average precipitation is 1876 mm. are alluvial plains, and the elevation of the whole basin is about 16-2136 m. The basin is located in the subtropical monsoon humid climate zone, with a mild climate and abundant rainfall. The annual average temperature is 17.6 °C, and the annual average precipitation is 1876 mm.

Data
The construction of the SWAT model requires topographic data, land use data, soi data, meteorological data, and hydrological data. The topographic data used in this study are from the Geospatial Data Cloud, with a resolution of 90 m (https://www.gscloud.cn/ accessed on 23 March 2022). Land use data of 1990 and 2010 were downloaded from the Resource and Environmental Science and Data Center platform of the Chinese Academy of Sciences (https://www.resdc.cn/, accessed on 24 March 2022). Soil data were obtained from the Harmonized World Soil Database v1.2 (HWSD) (https://www.fao.org/, accessed on 19 March 2022). The observed daily runoff data of the Meigang hydrological station from 1979 to 2016 were obtained from the Jiangxi Provincial Bureau of Hydrology.
The meteorological data used are CMADS (the China Meteorological Assimilation Driving Datasets for the SWAT model) long time-series data (http://www.cmads.org/, accessed on 20 March 2022). The data were created using various data processing tools, such as LAPS\STMAS assimilation algorithm, resampling, and bilinear interpolation, which now contains CMADS V1.0, CMADS-L V1.0, CMADS-ST V1.0, and other datasets with different spatial resolution and time series length [31]. The CMADS meteorological data used in this study are derived from the CMADS-L V1.0 dataset, which contains 23 CMADS stations covering the Xinjiang River basin and its surroundings, with a spatial resolution of 1/3° and a period of 1979 to 2018. Compared with traditional weather station data, CMADS weather data have the advantages of more uniform station distribution and easier data processing, and some studies have also shown that the CMADS-weather-data-

Data
The construction of the SWAT model requires topographic data, land use data, soil data, meteorological data, and hydrological data. The topographic data used in this study are from the Geospatial Data Cloud, with a resolution of 90 m (https://www.gscloud.cn/, accessed on 23 March 2022). Land use data of 1990 and 2010 were downloaded from the Resource and Environmental Science and Data Center platform of the Chinese Academy of Sciences (https://www.resdc.cn/, accessed on 24 March 2022). Soil data were obtained from the Harmonized World Soil Database v1.2 (HWSD) (https://www.fao.org/, accessed on 19 March 2022). The observed daily runoff data of the Meigang hydrological station from 1979 to 2016 were obtained from the Jiangxi Provincial Bureau of Hydrology.
The meteorological data used are CMADS (the China Meteorological Assimilation Driving Datasets for the SWAT model) long time-series data (http://www.cmads.org/, accessed on 20 March 2022). The data were created using various data processing tools, such as LAPS\STMAS assimilation algorithm, resampling, and bilinear interpolation, which now contains CMADS V1.0, CMADS-L V1.0, CMADS-ST V1.0, and other datasets with different spatial resolution and time series length [31]. The CMADS meteorological data used in this study are derived from the CMADS-L V1.0 dataset, which contains 23 CMADS stations covering the Xinjiang River basin and its surroundings, with a spatial resolution of 1/3 • and a period of 1979 to 2018. Compared with traditional weather station data, CMADS weather data have the advantages of more uniform station distribution and easier data processing, and some studies have also shown that the CMADS-weather-data-driven SWAT model has a better simulation effect than the traditional SWAT model driven by weather station data [32,33].

Methodology
The research framework of this paper is shown in Figure 2. Firstly, the CMADS meteorological data were used to drive the SWAT model. Then, the blue and green water of the Xinjiang River basin can be separated based on the simulation results of the SWAT. Then, a statistical analysis of the multi-year average blue and green water and their spatial distribution characteristics in the basin was performed. The precipitation anomaly percentage was used to discriminate typical years, and the amount of blue and green water and their spatial distribution characteristics in different typical years were analyzed. The Mann-Kendall test and linear regression were employed to test the changing trend of the blue and green water. Finally, the influence of climatic factors and land use factors on the changing trends of the blue and green water were explored based on scenario simulation. driven SWAT model has a better simulation effect than the traditional SWAT model driven by weather station data [32,33].

Methodology
The research framework of this paper is shown in Figure 2. Firstly, the CMADS meteorological data were used to drive the SWAT model. Then, the blue and green water of the Xinjiang River basin can be separated based on the simulation results of the SWAT. Then, a statistical analysis of the multi-year average blue and green water and their spatial distribution characteristics in the basin was performed. The precipitation anomaly percentage was used to discriminate typical years, and the amount of blue and green water and their spatial distribution characteristics in different typical years were analyzed. The Mann-Kendall test and linear regression were employed to test the changing trend of the blue and green water. Finally, the influence of climatic factors and land use factors on the changing trends of the blue and green water were explored based on scenario simulation.

SWAT Model
The SWAT model is a distributed hydrological model with a strong physical mechanism. The model consists of several sub-modules and can be used for various simulations, such as runoff simulation, sediment change simulation, and non-point source pollution simulation [34,35]. The simulation of the hydrological cycle is mainly composed of land surface hydrological cycle calculation and river confluence calculation. The calculation of the land surface hydrological cycle includes the calculation of surface runoff, interflow, and deep groundwater recharge to the main river, and the calculation methods include the SCS (soilconservation service) curve method and the Green and Ampt method. The river confluence evolution of runoff from the river network to the basin outlet is calculated by methods such as the Muskingum method. The SWAT model can fully consider the heterogeneity in soil properties, land use, and topographic slope by dividing the whole basin into many hydrological response units (HRUs). The model uses HRU as the basic unit for simulation and aggregates the simulation results of HRU to sub-basins, and then

SWAT Model
The SWAT model is a distributed hydrological model with a strong physical mechanism. The model consists of several sub-modules and can be used for various simulations, such as runoff simulation, sediment change simulation, and non-point source pollution simulation [34,35]. The simulation of the hydrological cycle is mainly composed of land surface hydrological cycle calculation and river confluence calculation. The calculation of the land surface hydrological cycle includes the calculation of surface runoff, interflow, and deep groundwater recharge to the main river, and the calculation methods include the SCS (soilconservation service) curve method and the Green and Ampt method. The river confluence evolution of runoff from the river network to the basin outlet is calculated by methods such as the Muskingum method. The SWAT model can fully consider the heterogeneity in soil properties, land use, and topographic slope by dividing the whole basin into many hydrological response units (HRUs). The model uses HRU as the basic unit for simulation and aggregates the simulation results of HRU to sub-basins, and then aggregates the simulation results of each sub-basin to the total basin outlet [36,37]. The simulation of land surface hydrological cycle processes by the SWAT model is based on the water balance equation of the soil layer as follows: where SW t is the final water content of the soil, mm; SW 0 is the initial water content of the soil, mm; t is time, d; R day is the precipitation of day i, mm; Q sur f is the surface runoff of day i, mm; E a is the evaporation of day i, mm; W seep is the infiltration and lateral flow at the bottom of soil profile of day i, mm; Q gw is the subsurface water content of day i, mm. More information about the principles and structure of the SWAT model can be found in the literature [38,39].

Model Assessment
There are many calibration methods for watershed hydrological models, including generalized likelihood uncertainty estimation (GLUE), Sequential Uncertainty Fitting program (SUFI-2), Bayesian parameter estimation methods, and particle swarm optimization algorithms [17,40]. The latest more advanced calibration methods include deep learning (DL). Traditional calibration methods are stable and reliable but time-consuming (slow convergence of parameters). DL can narrow down the optimal range of model parameters in a short time, which greatly improves the efficiency of model calibration and also obtains good calibration results [40,41].
In this work, the model calibration and validation were performed in SWAT-CUP using the SUFI-2 [42]. The SUFI-2 can provide the optimal parameter values and optimal parameter intervals at 95% confidence intervals. Nash efficiency coefficient (NSE), correlation coefficient (R 2 ), percentage bias (PBIAS), and ratio of root mean square error to standard deviation (RSR) were used to evaluate the calibration effect and simulation effect of the SWAT model. NSE indicates the degree of similarity between the simulated and measured values of the model, and a larger value indicates a better simulation of the model. R 2 indicates the degree of consistency between the simulated value and the measured value in the changing trend: the closer its value is to 1, the closer the changing trend is to the same. PBIAS indicates the average relative size of the measured and simulated values, PBIAS > 0 means the measured value is large compared to the simulated value, and PBIAS < 0 means the measured value is small compared to the simulated value. RSR reflects the magnitude of the model calibration error, and a smaller value indicates a smaller calibration error. It is generally considered that the model simulation is good when NSE > 0.75, R 2 > 0.8, |PBIAS| < 25%, and RSR < 0.7. The calculation equations of each indicator are as follows [43,44]: where Q 1 is the monthly measured runoff, m 3 /s; Q 2 is the monthly simulated runoff, m 3 /s; Q 1 is the monthly measured average runoff, m 3 /s; Q 2 is the monthly simulated average runoff, m 3 /s; n is the number of observations.

Blue Water, Green Water, and Green Water Coefficient
The SWAT model is able to separate the components of blue and green water at the HRU or sub-basin scale. Based on the output results of the SWAT model, blue water, green water, and green water coefficient can be directly calculated as follows: where BW is blue water, mm; GW is green water, mm; GWC is the green water coefficient; WYLD is the amount of water produced, mm; DA_RCHG is the amount of recharge of the deep aquifer, mm; ET is the actual evapotranspiration within the time step, which is also known as the green water flow, mm; SW is the soil water content, which is also known as the green water storage, mm.

Mann-Kendall Test
Mann-Kendall test is a non-parametric statistical test used for trend determination. The analysis results of this method are intuitive and reliable, and it has been widely used in hydrology, meteorology, and other fields. Supposing the time series X = {x 1 , x 2 , . . . , x n }, the order column S k can be calculated by: where the order column S k is the cumulative number of times that the value of the statistic r i at the ith moment is greater than the value at the jth moment. Then, the statistic UF k can be calculated by: where Calculate the statistic UB k : Calculate Z: A positive value of Z indicates an increasing trend of the time series and vice versa. When |Z| > 1.64, it means that it has passed the 5% significance test [45]. At the same time, the trend of the series can also be judged by UB k and UF k . Set the significance level (e.g., α = 0.05) and plot UB k and UF k curves; when the UB k or UF k curve is above the zero horizontal line, it indicates that the series has an increasing trend. On the contrary, when the UB k or UF k curve is below the zero horizontal line, it indicates that the series has a decreasing trend. If the two curves cross the critical line of significance (U α=0.05 = ±1.96), it indicates a significant increasing or decreasing trend of the series [46].

Precipitation Anomaly Percentage
The precipitation anomaly percentage is a relatively simple method to determine the typical year of hydrology, and the method is more accurate for long time series with annual and monthly scales. The calculation formula is as follows: where P a is the precipitation anomaly percentage; P i is the precipitation of the ith year; P is the average annual precipitation during the study period. were the calibration and validation periods, respectively. In this study, 23 runoff-related parameters were calibrated, and three indicators, including NSE, R 2 , PBIAS, and RSR, were used to evaluate the simulation effect of the SWAT. The values of each indicator during model calibration and validation periods are shown in Figure 3. It can be seen that the SWAT model exhibits a rather good performance; hence, the blue and green water in the Xinjiang River basin will be separated based on the simulation results of the model.

Precipitation Anomaly Percentage
The precipitation anomaly percentage is a relatively simple method to determine the typical year of hydrology, and the method is more accurate for long time series with annual and monthly scales. The calculation formula is as follows: where is the precipitation anomaly percentage; is the precipitation of the ith year; is the average annual precipitation during the study period.

Model Evaluation
The simulation period for the SWAT model is from 1979 to 2016, among which the period 1979-1981 was used as the warm-up period and the periods 1982-2000 and 2001-2016 were the calibration and validation periods, respectively. In this study, 23 runoffrelated parameters were calibrated, and three indicators, including NSE, , PBIAS, and RSR, were used to evaluate the simulation effect of the SWAT. The values of each indicator during model calibration and validation periods are shown in Figure 3. It can be seen that the SWAT model exhibits a rather good performance; hence, the blue and green water in the Xinjiang River basin will be separated based on the simulation results of the model. SWAT-CUP can also apply T-stat and p-value to analyze the sensitivity of each parameter. The value of T-stat represents the degree of sensitivity of the parameter, and a larger absolute value means that the parameter is more sensitive. The p-value represents the significance of the parameter sensitivity, and a value closer to zero means that the parameter sensitivity is more significant. The twelve most sensitive parameters in this study are listed in Table 1. SWAT-CUP can also apply T-stat and p-value to analyze the sensitivity of each parameter. The value of T-stat represents the degree of sensitivity of the parameter, and a larger absolute value means that the parameter is more sensitive. The p-value represents the significance of the parameter sensitivity, and a value closer to zero means that the parameter sensitivity is more significant. The twelve most sensitive parameters in this study are listed in Table 1.

Multi-Year Average Blue and Green Water and Their Spatial Distribution Characteristics
The annual average precipitation, blue water, green water, green water flow, green water storage, and green water coefficient of the Xinjiang River basin from 1982 to 2016 were 1867 mm, 1138 mm, 829 mm, 729 mm, 100 mm, and 0.42, respectively. The annual average blue water is about 1.38 times green water; hence, the water resources in the basin are dominated by blue water. As can be seen from Figure 4, the interannual variability in blue water is roughly consistent with that of precipitation. Green water is dominated by green water flow, which accounts for about 88% of the total green water, and green water storage accounts for about 12%. The interannual variation in green water, green water flow, and green water storage is relatively small, and green water and green water flow present similar variations. The interannual variation in the green water coefficient is large, with the smallest coefficient of 0.28 in 1998 when the blue water was about 2.61 times the green water and the largest coefficient of 0.57 in 2011 when the blue water was about 0.75 times the green water.  The spatial distribution characteristics of multi-year average precipitation, blu ter, and green water in the Xinjiang River basin are shown in Figure 5. It can be seen the spatial distribution characteristics of blue water are similar to precipitation, an blue water in the south-central part of the basin is obviously more than in other parts spatial distribution characteristics of green water and green water flow are similar the central bason covers both the most green water and green water flow, followed b northeast hills, southeast hills, northwest plains, and southwest hills. Green water st The spatial distribution characteristics of multi-year average precipitation, blue water, and green water in the Xinjiang River basin are shown in Figure 5. It can be seen that the spatial distribution characteristics of blue water are similar to precipitation, and the blue water in the south-central part of the basin is obviously more than in other parts. The spatial distribution characteristics of green water and green water flow are similar, and the central bason covers both the most green water and green water flow, followed by the northeast hills, southeast hills, northwest plains, and southwest hills. Green water storage is relatively evenly distributed throughout the basin, with most areas in the basin having green water storage between 95 and 104 mm, and very few areas being below 85 mm. The spatial distribution characteristics of the green water coefficient are opposite to that of blue water; the green water coefficient is smaller in areas with abundant blue water and larger in areas with little blue water. Green water is particularly important for maintaining the ecological balance in areas with large green water coefficients [10].

The Amount of Blue and Green Water in Typical Years
In this study, the precipitation anomaly percentage ( ) was used to discriminate the typical years, and the years within the study period were classified into five typical years based on values, including wet year ( > 20%), sub-humid year (10% < ≤ 20%), normal year (−10% ≤ ≤ 10%), sub-arid year (−20% ≤ < −10%), and dry year ( < −20%). As can be seen from Figure 6, wet years, sub-humid years, and sub-arid years each include 5 years, normal years include 14 years, and dry years include 6 years. The sum of the proportion of wet years and dry years is nearly 1/3, and the proportion of extreme precipitation years is high. At the same time, it can be found that there are frequent rapid transitions between dry and wet between 1995 and 1998 and between 2009 and 2013.

The Amount of Blue and Green Water in Typical Years
In this study, the precipitation anomaly percentage (P a ) was used to discriminate the typical years, and the years within the study period were classified into five typical years based on P a values, including wet year (P a > 20%), sub-humid year (10% < P a ≤ 20%), normal year (−10% ≤ P a ≤ 10%), sub-arid year (−20% ≤ P a < −10%), and dry year (P a < −20%). As can be seen from Figure 6, wet years, sub-humid years, and sub-arid years each include 5 years, normal years include 14 years, and dry years include 6 years. The sum of the proportion of wet years and dry years is nearly 1/3, and the proportion of extreme precipitation years is high. At the same time, it can be found that there are frequent rapid transitions between dry and wet between 1995 and 1998 and between 2009 and 2013. In order to avoid the one-sidedness and errors brought by the selection of a single typical year to the analysis results, we improved the previous research method by counting the multi-year average data for each type of typical year. The multi-year average blue and green water for typical years in the Xinjiang River basin are shown in Figure 7. The amount of blue water in the wet years is 2412 mm, which is about 1.5 times that in normal years and 2.4 times that in dry years. This analysis result is similar to the Ganjiang River basin, which is also a sub-basin of the Poyang Lake basin (the blue water in Ganjiang River basin is about 1.54 times normal years and 2.01 times dry years [19]). The variation in green water in different typical years is smaller compared with blue water, and the green water in sub-humid years is about 1.08 times that in dry years. Green water flow occupies the vast majority of green water (green water flow accounts for about 88% of green water), with an average of about 720 ± 30 mm per year. The green water coefficient is 0.33 in wet years, 0.43 in normal years, and 0.53 in dry years. The green water coefficients of the Xinjiang River basin were 0.33, 0.43, and 0.53 in wet, normal, and dry years, respectively, which were significantly lower than those of the Ganjiang River basin (0.40 in wet years, 0.51 in normal years, and 0.57 in dry years). This is related to the more abundant precipitation in the Xinjiang basin because it is closer to the ocean. It can be seen that the proportion of green water in dry years is significantly larger than that in wet years; hence, green water is particularly important for maintaining ecosystem balance in dry years. In order to avoid the one-sidedness and errors brought by the selection of a single typical year to the analysis results, we improved the previous research method by counting the multi-year average data for each type of typical year. The multi-year average blue and green water for typical years in the Xinjiang River basin are shown in Figure 7. The amount of blue water in the wet years is 2412 mm, which is about 1.5 times that in normal years and 2.4 times that in dry years. This analysis result is similar to the Ganjiang River basin, which is also a sub-basin of the Poyang Lake basin (the blue water in Ganjiang River basin is about 1.54 times normal years and 2.01 times dry years [19]). The variation in green water in different typical years is smaller compared with blue water, and the green water in sub-humid years is about 1.08 times that in dry years. Green water flow occupies the vast majority of green water (green water flow accounts for about 88% of green water), with an average of about 720 ± 30 mm per year. The green water coefficient is 0.33 in wet years, 0.43 in normal years, and 0.53 in dry years. The green water coefficients of the Xinjiang River basin were 0.33, 0.43, and 0.53 in wet, normal, and dry years, respectively, which were significantly lower than those of the Ganjiang River basin (0.40 in wet years, 0.51 in normal years, and 0.57 in dry years). This is related to the more abundant precipitation in the Xinjiang basin because it is closer to the ocean. It can be seen that the proportion of green water in dry years is significantly larger than that in wet years; hence, green water is particularly important for maintaining ecosystem balance in dry years.

Spatial Distribution Characteristics of Blue and Green Water in Typical Years
Three types of typical years were selected to analyze the spatial distribution characteristics of blue and green water in the Xinjiang River basin. It can be seen from Figure 8 that the spatial distribution characteristics of blue water in wet years, normal years, and dry years are similar to those of precipitation. The central bason has the most abundant blue water, and the southeast mountains and northwest plains have the least blue water. In wet years, the blue water in the basin presents an uneven distribution such that the central bason has the most blue water of about 1772 mm while the northwest plain has the least blue water of only 1558 mm. In normal years, the central bason has the most blue water of about 1221 mm, and the northwest plain has the least blue water of about 994 mm. In dry years, the blue water of the basin is more evenly distributed than in wet and normal years.

Spatial Distribution Characteristics of Blue and Green Water in Typical Years
Three types of typical years were selected to analyze the spatial distribution charac teristics of blue and green water in the Xinjiang River basin. It can be seen from Figure that the spatial distribution characteristics of blue water in wet years, normal years, an dry years are similar to those of precipitation. The central bason has the most abundant blu water, and the southeast mountains and northwest plains have the least blue water. In we years, the blue water in the basin presents an uneven distribution such that the central baso has the most blue water of about 1772 mm while the northwest plain has the least blue wate of only 1558 mm. In normal years, the central bason has the most blue water of about 122 mm, and the northwest plain has the least blue water of about 994 mm. In dry years, th blue water of the basin is more evenly distributed than in wet and normal years.
The spatial distribution characteristics of green water are similar in the three types o typical years. In wet years, the central bason has the most green water of about 840 mm and the southwest hills have the least green water of about 821 mm. In normal years, th green water in most parts of the basin is more than 800 mm, and the uneven spatial dis tribution of green water in the whole basin is the most obvious in three types of typica years. In dry years, the green water is significantly reduced, and the green water in th whole basin is less than 820 mm.
The spatial distribution characteristics of the green water coefficient are exactly op posite to those of the blue water. It is shown that the green water coefficients in the south east and northwest are obviously larger than that in the central bason. In wet years, th green water coefficient is between 0.32 and 0.35 in all parts of the basin, and there is littl difference in green water coefficient across the basin. In normal years, the green wate coefficient is in the range of 0.41 to 0.46 throughout the basin, with small differences from place to place. In dry years, the green water coefficient in the whole basin varies from 0.5 to 0.58, and green water is the dominant water resource in the whole basin. The decreas in precipitation and the increase in evapotranspiration lead to a significant increase in th green water coefficient in dry years.
The Xinjiang River basin has abundant topographic features, consisting of fores The spatial distribution characteristics of green water are similar in the three types of typical years. In wet years, the central bason has the most green water of about 840 mm, and the southwest hills have the least green water of about 821 mm. In normal years, the green water in most parts of the basin is more than 800 mm, and the uneven spatial distribution of green water in the whole basin is the most obvious in three types of typical years. In dry years, the green water is significantly reduced, and the green water in the whole basin is less than 820 mm.
The spatial distribution characteristics of the green water coefficient are exactly opposite to those of the blue water. It is shown that the green water coefficients in the southeast and northwest are obviously larger than that in the central bason. In wet years, the green water coefficient is between 0.32 and 0.35 in all parts of the basin, and there is little difference in green water coefficient across the basin. In normal years, the green water coefficient is in the range of 0.41 to 0.46 throughout the basin, with small differences from place to place. In dry years, the green water coefficient in the whole basin varies from 0.51 to 0.58, and green water is the dominant water resource in the whole basin. The decrease in precipitation and the increase in evapotranspiration lead to a significant increase in the green water coefficient in dry years.
The Xinjiang River basin has abundant topographic features, consisting of forest land, basin, and plain. Its spatial distribution of precipitation is affected jointly by geographical location, topography, and climatic conditions. The spatial distribution characteristic of blue water is dominated by precipitation, while that of green water is associated with soil types and land cover in addition to precipitation. The blue water in the middle of the basin has always been prominent in typical years for the precipitation in most of the middle reach. The green water in the middle reach is possibly related to the regional evapotranspiration increased by flat terrain, extensive plowland, and grassland in the middle reach apart from the distribution characteristic of precipitation.
with soil types and land cover in addition to precipitation. The blue water in the middle of the basin has always been prominent in typical years for the precipitation in most of the middle reach. The green water in the middle reach is possibly related to the regional evapotranspiration increased by flat terrain, extensive plowland, and grassland in the middle reach apart from the distribution characteristic of precipitation.

Temporal Variation Characteristics of Blue and Green Water
Meteorological factors and land cover are two important factors that may affect the temporal change characteristics of blue and green water. Before studying the temporal

Temporal Variation Characteristics of Blue and Green Water
Meteorological factors and land cover are two important factors that may affect the temporal change characteristics of blue and green water. Before studying the temporal change characteristics of blue and green water, it is necessary to analyze the changes in meteorological factors and land use in the whole basin.

Changes in Precipitation and Temperature
This study applied the Mann-Kendall (MK) test and linear regression to detect the changes in the precipitation and temperature of the basin. It can be seen from Figure 9 that the Z value of the MK test for precipitation is −0.072, which fails to pass the 5% significance test and shows a slight decrease. The Z value of the MK test for average temperature is 0.405, which passed the 5% significance test and shows a relatively significant increase. The analysis result of linear regression shows that precipitation decreased at a rate of 3.51 mm/year, and average temperature increased at a rate of 0.026 • C/year. change characteristics of blue and green water, it is necessary to analyze the changes in meteorological factors and land use in the whole basin.

Changes in Precipitation and Temperature
This study applied the Mann-Kendall (MK) test and linear regression to detect the changes in the precipitation and temperature of the basin. It can be seen from Figure 9 that the Z value of the MK test for precipitation is −0.072, which fails to pass the 5% significance test and shows a slight decrease. The Z value of the MK test for average temperature is 0.405, which passed the 5% significance test and shows a relatively significant increase. The analysis result of linear regression shows that precipitation decreased at a rate of 3.51 mm/year, and average temperature increased at a rate of 0.026 °C/year.

Land Use Change
The land use of the Xinjiang River basin in 1990 and 2010 is shown in Table 2. It can be seen that the land use of the Xinjiang River basin is dominated by woodland and plowland, accounting for about 67% and 26% of the total area of the basin, respectively. The grassland, water area, construction land, and bare land are all relatively small, which are all no more than 4% of the total area of the basin. The area of construction land in the whole basin increased rapidly between 1990 and 2010, with an increase of more than 75 and an area change rate of 0.5%. From 1990 to 2010, the area of plowland reduced by more than 65 and the area of grassland reduced by more than 25 km , part of which were transformed into construction land, woodland, and water area. The proportion of bare land area in the basin is negligible and has not demonstrated a significant change.

Temporal Variation Characteristics of Blue and Green Water
The MK test was used to analyze the trend of blue water, green water, green water flow, and green water coefficient in the Xinjiang River basin, and the results are shown in

Land Use Change
The land use of the Xinjiang River basin in 1990 and 2010 is shown in Table 2. It can be seen that the land use of the Xinjiang River basin is dominated by woodland and plowland, accounting for about 67% and 26% of the total area of the basin, respectively. The grassland, water area, construction land, and bare land are all relatively small, which are all no more than 4% of the total area of the basin. The area of construction land in the whole basin increased rapidly between 1990 and 2010, with an increase of more than 75 km 2 and an area change rate of 0.5%. From 1990 to 2010, the area of plowland reduced by more than 65 km 2 and the area of grassland reduced by more than 25 km 2 , part of which were transformed into construction land, woodland, and water area. The proportion of bare land area in the basin is negligible and has not demonstrated a significant change.

Temporal Variation Characteristics of Blue and Green Water
The MK test was used to analyze the trend of blue water, green water, green water flow, and green water coefficient in the Xinjiang River basin, and the results are shown in Figure 10. The Z value of the MK test for blue water is −0.008, which fails to pass the 5% significance test and shows a very slight decrease. The Z values of the MK test for green water and green water flow are both less than −1.64, which suggests significant decreasing trends at a 5% significance level. Meanwhile, it can be seen from the UF curve that the UF values of green water and green water flow are all less than zero most of the time and are positive only in some years from 1983 to 1990, which indicates that green water and green water flow had a decreasing trend most of the time. The UF and UB curves of green water and green water flow have an intersection point in 1992, and the intersection point is within the confidence interval, so it indicates that green water and green water flow have been decreasing since 1992. The Z value of the MK test for the green water coefficient is −0.096, which indicates a slight decrease and failed to pass the 5% significance test. It can also be found that the results of linear trends of blue water, green water, and green water coefficients are consistent with the results of the MK test. Figure 10. The Z value of the MK test for blue water is −0.008, which fails to pass the 5% significance test and shows a very slight decrease. The Z values of the MK test for green water and green water flow are both less than −1.64, which suggests significant decreasing trends at a 5% significance level. Meanwhile, it can be seen from the UF curve that the UF values of green water and green water flow are all less than zero most of the time and are positive only in some years from 1983 to 1990, which indicates that green water and green water flow had a decreasing trend most of the time. The UF and UB curves of green water and green water flow have an intersection point in 1992, and the intersection point is within the confidence interval, so it indicates that green water and green water flow have been decreasing since 1992. The Z value of the MK test for the green water coefficient is −0.096, which indicates a slight decrease and failed to pass the 5% significance test. It can also be found that the results of linear trends of blue water, green water, and green water coefficients are consistent with the results of the MK test.

Analysis of the Reasons for the Changing Trend of Blue and Green Water
In order to further explore the factors that induce the trend of blue and green water in the Xinjiang River basin, this study sets three different scenarios of land use and climate conditions (scenario settings are shown in Table 3). By comparing the simulation results of scenario 1 and scenario 2, which have the same climatic conditions but different land use, we can explore the influence of land use on the changing trend of blue and green water. By comparing the simulation results of scenario 2 and scenario 3, which have the same land use but different climatic conditions, we can explore the influence of climatic factors on the changing trend of blue and green water. The simulation results are shown in Figure 11. It can be found that the simulation results of scenario 1 and scenario 2 are relatively close, but the simulation results of scenario 2 and scenario 3 are quite different.

Analysis of the Reasons for the Changing Trend of Blue and Green Water
In order to further explore the factors that induce the trend of blue and green water in the Xinjiang River basin, this study sets three different scenarios of land use and climate conditions (scenario settings are shown in Table 3). By comparing the simulation results of scenario 1 and scenario 2, which have the same climatic conditions but different land use, we can explore the influence of land use on the changing trend of blue and green water. By comparing the simulation results of scenario 2 and scenario 3, which have the same land use but different climatic conditions, we can explore the influence of climatic factors on the changing trend of blue and green water. The simulation results are shown in Figure 11. It can be found that the simulation results of scenario 1 and scenario 2 are relatively close, but the simulation results of scenario 2 and scenario 3 are quite different. Therefore, we can conclude that climatic factors have a greater impact on blue and green water than land use and are the dominant reason for the temporal variation characteristics of blue and green water. According to the above analysis, it is apparent that the decrease in blue water in the basin is related to the decrease in precipitation, while the decrease in green water is related to the decrease in precipitation and evapotranspiration area, and the decrease in precipi tation is the main reason for the decrease in green water.   At the same time, it can be seen from Figure 11 that blue water in scenario 2 is 1.5 mm more than that in scenario 1, which may be related to the change in land use. The increase in construction land leads to the increase in impervious surface area, and the rain falling on the impervious surface cannot penetrate the soil and directly forms surface runoff, which leads to a slight increase in the blue water. It can also be found that the green water in scenario 2 is 1.6 mm less than that in scenario 1, which may be related to the decrease in the area of evapotranspiration in the basin. From 1990 to 2010, the increase in construction land and the decrease in plowland and grassland in the basin led to partial precipitation that could not be converted into ecosystem water through infiltration and vegetation interception, which led to a slight decrease in the green water. At the same time, the small difference between the simulation results of scenario 1 and scenario 2 is related to the small change in land use in the Xinjiang River basin. Compared to the simulation results of scenario 2 and scenario 3, it can be found that the blue water and green water in scenario 3 are 154.8 mm and 48.5 mm less than those in scenario 2, respectively, which is related to the decrease in precipitation in scenario 3. The decrease in precipitation directly leads to the decrease in surface runoff and evapotranspiration, which leads to a significant decrease in both blue and green water.
According to the above analysis, it is apparent that the decrease in blue water in the basin is related to the decrease in precipitation, while the decrease in green water is related to the decrease in precipitation and evapotranspiration area, and the decrease in precipitation is the main reason for the decrease in green water.

Conclusions
In this study, we applied CMADS long time series ) meteorological data to the study of blue and green water and conducted a comprehensive analysis of the spatiotemporal distribution characteristics of blue and green water in the Xinjiang River basin based on the simulation results of the SWAT model. Firstly, the annual average blue and green water and their spatial distribution characteristics in the basin were quantitatively evaluated. Then, the precipitation anomaly percentage was used to discriminate the typical years, and the blue and green water in typical years were assessed. The variation trend of blue and green water was analyzed by using the MK test and linear regression. Finally, three different scenarios were set up to explore the effects of climatic factors and land use on the trend of the blue and green water. The following conclusions are drawn from this study: (1) The SWAT model driven by CMADS long-term meteorological data has achieved satisfactory simulation results, with an NSE of 0.97 in the calibration period and an NSE of 0.93 in the validation period.
(2) The annual average blue water, green water, green water flow, green water storage, and green water coefficient of the basin were 1138 mm, 829 mm, 729 mm, 100 mm, and 0.42, respectively. The water resources in the basin are dominated by the blue water. However, the annual average green water coefficient in dry years is 0.53, and green water is the main water resource in dry years. Green water flow accounts for about 88% of the total green water and is the main part of the green water.
(3) The spatial distribution of blue and green water in the basin is uneven, with the central bason having significantly more blue and green water than other parts of the basin. The spatial distribution of the blue water is consistent with that of precipitation, but the spatial distribution of the green water coefficient is the opposite to that of precipitation. Compared with wet years and normal years, the distribution of the blue and green water in the basin is more uniform in dry years.
(4) The blue water in the Xinjiang River basin shows a slight decreasing trend, and the green water shows a significant decreasing trend. The decrease in blue water is mainly related to the decrease in precipitation. The decrease in evaporable water and evapotranspiration area resulted in a decrease in the green water. The results of the scenario simulation found that climatic factors have a greater influence on the change trends in the blue and green water compared to land use and are the dominant cause of the temporal variation characteristics of blue and green water.

Conflicts of Interest:
The authors declare no conflict of interest.