Assessment of Surface Water Availability under Climate Change Using Coupled SWAT-WEAP in Hongshui River Basin, China

: Climate change adversely affects the hydrological cycle at the basin level. This study integrated two models, the Soil and Water Assessment Tool (SWAT) for future climate prediction, and Water Evaluation and Planning (WEAP) for the simulation of water quantity in the Hongshui River Basin (HRB), to evaluate the impacts of climate change, which plays a signiﬁcant role in the lives of inhabitants downstream of the basin. Downscaled monthly rainfalls and temperatures under four Representative Concentration Pathways (RCPs) emission scenarios from ﬁve global circulation models (GCMs) were used to generate streamﬂow using the SWAT model. Streamﬂow data (1991–2001) were used to calibrate and validate, with the period of 1991–1997 used for calibration and that of 1998–2001 used for validation. Six scenarios were established to evaluate the response of the basin under socio-economic scenarios. The simulated results show that precipitation and streamﬂow would likely undergo a slight increase. The available water resources would be sufﬁcient to meet the existing needs until 2050. The results indicated that no water shortages exist under socio-economic, low, and medium climate change emission scenarios, however the basin will experience a water shortage under the high climate change emission scenario (RCP-8.5). The study proposed that, to ensure the sustainability of water resources, better long-term management policies are required to be implemented in the basin and to meet future downstream water needs.


Introduction
Climate change is a major challenge and a leading cause of global warming. Climate changes alter the average temperature of the Earth's climate system and induce related effects. These effects are likely to affect the hydrological cycle and, consequently, water availability for agriculture, industrial, and domestic sectors [1]. To tackle future water stress, the sustainable use of water resources, at local and global scales, is essential [2]. Sustainable water resources management requires analysis and quantification of several hydrological processes taking place within the watershed basin [3]. Anthropogenic effects have rapidly changed the global eco-system. These rapid changes can be seen as socio-economic development, population growth, and deforestation associated with the expansion of cities [4]. Climate change and urbanization are interconnected and have had an increasing environmental impact during the 21st century. The urban population, as a share of the total, recently reached 50%, and this proportion will be over 80% by 2050, according to Yan, et al. [5]. Global warming has altered precipitation and temperature patterns, which could have a significant impact on the local hydrological cycle. These changes in the hydrological cycle will increase or decrease stream water flows [6].
Climate change, population growth, and economic development will clearly affect the availability of potential water resources in numerous regions. Climate change impacts ISPRS Int. J. Geo-Inf. 2021, 10, 298 2 of 21 are a global issue, however, researchers are most interested in local effects at the basin level. One of the most influential findings of related research is that climate change will alter hydrological processes and water availability around the world [7]. Introducing proper water management strategies and understanding the interaction between climate change and water resources, will help researchers and policymakers reduce the harmful effects. Several studies of the impacts of climate change on hydrological processes indicate that streamflow variations are related to changes in temperature and rainfall [3,[8][9][10]. Researchers frequently analyze hydro-climatological processes using the global circulation model (GCM), and hydrological models have been demonstrated to support water resource management and prediction, especially at the basin scale [11][12][13][14]. In the Upper Colorado River Basin (UCRB), demand in most years exceeds supply, and is likely to continue to rise. Results of research using the GCM showed that the summer streamflow declines with a median decrease of 46%, and a total range of +22% until 2100 [15]. In a study of the uncertainty related to seven GCMs, the consequences for the Mekong River Basin of a 2 • C increase in the global average temperature were investigated. Result showed a relatively small, yet nonlinear, response of annual river discharge, from 4.5% to 5.4%, to a rise in the global average temperature [16]. Studies indicates that a temperature increase of 0.05 • C per decade has occurred as a consequence of urbanization in mainland China [17]. An investigation of climate change impacts on the Lanchang River (LR), in the Tibetan plateau in China, show that climate change contributed 57% to the increase in stream flows, and human activity contributed 95% to the post-impact period (2008-2014) [18]. Climate data were projected using a statistical downscaling model and three CMIP5 global circulation models under Representative Concentration Pathways for Xin River Basin, China. The average ensemble of the GCMs showed that streamflows under RCPs decrease significantly [12]. The characteristics of hydro-climatic changes in the Tarim River Basin were analyzed using data collected at 39 weather stations and 29 hydrological stations for the periods 1961-2008 and 1952-2008. Changes were explained by anthropogenic activities and climate change [19], and the hydrological response to climate change under the average of four RCPs showed a decline of 21.5-40% in annual average runoff [20]. During the past 50 years, runoff has decreased in most river basins in China. Examination of the contributions of related human activity and climate change to the hydrological response of the Weihe River Basin established that human activity has had a greater impact on basin runoff than climate change [21][22][23]. These studies have highlighted the strong impact on stream flows of climate change and illustrated the importance of incorporating such research into adaptive management.
Researchers have widely used hydrological models to assess climate change impacts on hydrological regimes. These tools can be used to calculate many hydrological parameters, including streamflow, runoff, evapotranspiration, and yield of water. Computing these hydrological parameters and selecting hydrological approaches are crucial to assessing the impacts of current and future climate change [2,9,24,25]. Integrated simulation models are needed to quantify water supply and demand balance [24]. Integrated water change management concepts have been well established as one of the best approaches to manage sustainability of water resources [25,26]. The Water Evaluation and Planning (WEAP) model provides a set of objects and procedures to address issues faced by water managers using a scenario-based approach that can be applied to natural watersheds, reservoirs, streams, and channels [27]. The WEAP tool has been used by several researchers to model water demand and supply, and to assess climate change impacts on hydrological processes, groundwater management, and water pricing [28][29][30][31][32]. These models are mainly restricted to water resource research in relation to runoff, however, and no analysis has been undertaken of other hydrological components for physical use, such as evapotranspiration (ET). The Soil and Water Assessment Tool (SWAT) model and other hydrological models can adequately represent hydrological processes and other physical hydrological elements. Nonetheless, traditional water-based simulation models are often unable to address contemporary water management problems [33,34]. Therefore, a robust water resource assessment and management tool is needed to consider hydrological responses to climate change, water resource allocation, the agricultural economy, and other factors in the design of water management options. In this study, we established a SWAT-WEAP combination approach that considered both hydrological processes and water management evaluations and analysis [35].
Historically, the Hongshui River Basin, China, one of the major streams of the Xijiang River Basin, did not undergo rapid land cover changes. However, government reforms in the early 1980s, undertaken to boost economic growth, directly improve living standards, and increase urbanization by 300 percent, had serious impacts on various hydro-climatic observations [36,37]. The overexploitation of water resources in the region has resulted in a series of environmental and ecological effects, including the drying up of rivers and lakes, groundwater depletion, grassland erosion, and ravaging dust storms. Additionally, water shortages have exacerbated disputes around provincial boundaries between upstream and downstream water users [38]. Guo, et al. [39] also stressed that climate change plays a key role in shifting the hydrology and flow of basins in the Xijiang River Basin. Studies have been conducted to examine climate change impacts on Hongshui River Basin water resources. Huang, Wang, Xiao, Chen, Zhou, Song and Wang [38] examined anthropogenic activities and climate change impacts on long-term streamflow patterns. The study reported negative effects, including reduced runoff, using a variable infiltration capability (VIC) model and hydrological sensitivity system in the upper reaches of Hongshui Basin, Southwest China, during the past 55 years. Touseef, Chen, Masud, Khan, Yang, Shahzad, Ijaz and Wang [9] applied the Soil and Water Assessment Tool (SWAT) to the upper reaches of the Xijiang Basin of Hongshui River and concluded that large negative shifts in low flows can be expected under future climate scenarios, showing a monthly decrease of 10% to 30% in water yield over the watershed. However, no study has been carried out under all four future scenarios to assess effects of climate change on streams and water demand using a combined SWAT-WEAP approach.
In this study, we introduced a SWAT-WEAP integrated approach, which considered both hydrological processes and water management assessments, and hydrological component analyses. Five global climate model (GCM) outputs under each of the four future scenarios were used to simulate the impact of climate change on streamflow and water demand across the Hongshui River Basin (HRB). The major objective of this research was to investigate the impacts of socio-economic and climate change scenarios on hydrological processes and water demand in the HRB. The study's sub-objectives were to calibrate and validate the hydrological module embedded in the WEAP model for streamflow simulation at Qianjiang and Tianer stations; to forecast and analyze changes in stream flows in the catchment under emission scenarios between 2010 and 2050; and forecast and analyze resulting water shortages in the HRB.

Study Area
The chosen study area is the Hongshui River Basin (HRB; shown in Figure 1), above the Qianjiang River Gauge Station. The Hongshui River Basin is the Pearl River Basin's main tributary, which lies in South China's subtropical region. In turn, the Pearl River Basin is China's third-largest river basin, with a population of more than 100 million. The Hongshui River Basin has a total drainage area of 98,500 km 2 [40]. The main channel of the Upper Xijiang River Basin has distinctive names for specific sections: Beipanjiang, Nanpanjiang, Hongshuihe, and Qianjiang [41]. The basin is dry and tropical with sufficient precipitation, and usually has a high air temperature. The mean air temperature is about 14-22 • C. The mean annual precipitation ranges from 760 to 1860 mm, with west-to-east variations. Precipitation occurs mainly from April to October, accounting for 72-86% of annual precipitation [42]. Significant changes have occurred in the natural cycle of the river, due to climate change and to anthropogenic activities, such as afforestation, reservoir operation and irrigation infrastructure [40]. Nanpanjiang, Hongshuihe, and Qianjiang [41]. The basin is dry and tropical with sufficient precipitation, and usually has a high air temperature. The mean air temperature is about 14-22 °C. The mean annual precipitation ranges from 760 to 1860 mm, with westto-east variations. Precipitation occurs mainly from April to October, accounting for 72-86% of annual precipitation [42]. Significant changes have occurred in the natural cycle of the river, due to climate change and to anthropogenic activities, such as afforestation, reservoir operation and irrigation infrastructure [40].

Data Availability
Hydrological modeling requires data on topography, soil characteristics, and land use/land cover (LULC) maps, as summarized in Table 1. The digital elevation model (DEM) is the most important input to provide topographical information in the WEAP system. This study used the Shuttle Radar Topography Mission Digital Elevation Model (SRTM-DEM) with a resolution of 90 m from the CGIAR Consortium for Spatial Information (CGIAR-CSI) (http://srtm.csi.cgiar.org/) [43]. The DEM is a simple input layer in the model configuration and was used to extract the slope and drainage network of the study area. The FAO-UNESCO global soil map, with a 5 km resolution, was obtained as the land map of the study area from the Food and Agriculture Organization (FAO) (http://www.fao.org/nr/land/soils/digital-soil-map-of-the-world/en/) [44].
LULC 2015, developed by the European Space Agency (ESA) as part of the Climate Change Initiative CCI-LC project, was analyzed to simulate stream flow discharge and water demand (http://maps.elie.ucl.ac.be/CCI/viewer/download.php) [45].
The National Meteorological Information Center (NMIC) of the China Meteorological Administration (CMA) provided daily precipitation data for 32 weather stations and two river discharge stations in the Hongshui River Basin [10].

Data Availability
Hydrological modeling requires data on topography, soil characteristics, and land use/land cover (LULC) maps, as summarized in Table 1. The digital elevation model (DEM) is the most important input to provide topographical information in the WEAP system. This study used the Shuttle Radar Topography Mission Digital Elevation Model (SRTM-DEM) with a resolution of 90 m from the CGIAR Consortium for Spatial Information (CGIAR-CSI) (http://srtm.csi.cgiar.org/ (accessed on 21 February 2021)) [43]. The DEM is a simple input layer in the model configuration and was used to extract the slope and drainage network of the study area. The FAO-UNESCO global soil map, with a 5 km resolution, was obtained as the land map of the study area from the Food and Agriculture Organization (FAO) (http://www.fao.org/nr/land/soils/digital-soil-map-of-the-world/en/ (accessed on 21 February 2021)) [44].
LULC 2015, developed by the European Space Agency (ESA) as part of the Climate Change Initiative CCI-LC project, was analyzed to simulate stream flow discharge and water demand (http://maps.elie.ucl.ac.be/CCI/viewer/download.php (accessed on 21 February 2021)) [45].
The National Meteorological Information Center (NMIC) of the China Meteorological Administration (CMA) provided daily precipitation data for 32 weather stations and two river discharge stations in the Hongshui River Basin [10].
GCMs are major tools used to test climate change scenarios, and to model the global climate system's reaction to changes in the atmosphere [46]. This study analyzed the climate data sets of five GCMs from the Inter-Sectorial Impact System Inter-Comparison Project (ISI-MIP)-GFDL-ESM2 M, HadGEM2-ES, IPSL-CM5A-LR, MIROC-ESM-CHEM, and NorESM1-M [47]-and three emission scenarios (RCP-2.6, RCP-4.5, and RCP-8.5). The raw GCM output was statistically downscaled, and bias-corrected using special disaggregation (BCSD) with the Climate Change Toolkit (CCT) [48][49][50]. The CCT package also includes historical climate data  from the Climate Research Unit (CRU TS 3.1) that can be used as an observed dataset. All climate datasets (precipitation, temperature, relative humidity, solar radiation, wind) are downscaled and available in text format. The CCT was used to undertake extraction, downscaling, bias correction, and interpolation of raw GCM outputs.

Integrated Hydrological Modeling Development
The combined SWAT-WEAP approach (shown in Figure 2) was used in this study to evaluate comprehensive water resource management options. The Soil and Water Assessment Tool (SWAT) is a physically based semi-distributed basin model [52]. The SWAT model uses the soil conservation services (SCS) approach and simulates hydrological processes, such as surface runoff, snowmelt runoff, lateral flow, and drainage into a river basin, in addition to water balance and streamflow [53][54][55]. The SWAT hydrological model is based on the basic water balance equation: where SW total is the total soil water content, SW 0 is the initial day soil water, t is the time in days, R d is the daily rainfall, Q s is the surface run-off, E a is the actual evapotranspiration, W seep is the quantity of water seep into the vadose zone on day t, and Q gw is the return flow on day t. The Water Evaluation and Planning (WEAP) software, developed by the Stockholm Environment Institute (SEI), is a practical semi-theoretical, semi-distributed, and deterministic water resource planning tool that incorporates both water supply and water demand issues, in addition to water quality and the conservation of ecosystems, as required by an integrated approach to basin management [56,57]. The WEAP model was used to simulate both the hydrological processes and the anthropogenic activities of water resources to analyze the availability of water in the basin. The WEAP model has been implemented at national and international levels because it offers a flexible and comprehensive policy analysis framework, in addition to a method for managing water supply and demand. The WEAP model uses five methods to simulate daily or monthly catchment processes. For this study, the soil moisture method was selected because it allows characterization of the impacts of land use and soil types on these processes [27,31]. WEAP simulates water supply and demand based on scenarios, such as demographic changes or water policies for a user-required period. These hypothetical predictions (or scenarios) are compared to the current year. Generally, the baseline is the year for which all the input data required by the WEAP model, including streamflow, climate, demographic, and land use/land cover, are available. It represents the actual condition in the catchment, also known as the current account. The baseline year was chosen because the model needs a current water availability value to be used to assess current and future water availability in different socio-economic, climate, or management scenarios [58]. The integrated approach was applied to simulate the future water resource status (2010-2050) by combining various water resource management options and climate change scenarios. catchment processes. For this study, the soil moisture method was selected because it allows characterization of the impacts of land use and soil types on these processes [27,31]. WEAP simulates water supply and demand based on scenarios, such as demographic changes or water policies for a user-required period. These hypothetical predictions (or scenarios) are compared to the current year. Generally, the baseline is the year for which all the input data required by the WEAP model, including streamflow, climate, demographic, and land use/land cover, are available. It represents the actual condition in the catchment, also known as the current account. The baseline year was chosen because the model needs a current water availability value to be used to assess current and future water availability in different socio-economic, climate, or management scenarios [58]. The integrated approach was applied to simulate the future water resource status (2010-2050) by combining various water resource management options and climate change scenarios.

Model Calibration and Validation
This study calibrated and validated the SWAT model using the SWAT-CUP (Soil and Water Assessment Tool-Calibration and Uncertainty Programme). The observed streamflow data were used for calibration and validation purposes at the Qianjiang and Tianer hydrological stations, calibrated for a period of 7 years (1991-1997) and validated for 1998-2001 from the available streamflow data presented in Figure 3 a,b. A future climate scenario was generated and embedded in the SWAT model and used to drive the SWAT model to simulate tributary streaming under various climate change scenarios. The simulated discharge output of SWAT was used to drive the WEAP model. The calibration and validation procedure of the WEAP model was carried out using the PEST (Parameter Estimation) routine. PEST is a nonlinear estimator of parameters and is considered a specific calibration method. Responsive parameter adjustment is undertaken by trial and error to determine the best value for a parameter. Hydrological output must be checked using objective functions.

Model Calibration and Validation
This study calibrated and validated the SWAT model using the SWAT-CUP (Soil and Water Assessment Tool-Calibration and Uncertainty Programme). The observed streamflow data were used for calibration and validation purposes at the Qianjiang and Tianer hydrological stations, calibrated for a period of 7 years (1991-1997) and validated for 1998-2001 from the available streamflow data presented in Figure 3a,b. A future climate scenario was generated and embedded in the SWAT model and used to drive the SWAT model to simulate tributary streaming under various climate change scenarios. The simulated discharge output of SWAT was used to drive the WEAP model. The calibration and validation procedure of the WEAP model was carried out using the PEST (Parameter Estimation) routine. PEST is a nonlinear estimator of parameters and is considered a specific calibration method. Responsive parameter adjustment is undertaken by trial and error to determine the best value for a parameter. Hydrological output must be checked using objective functions.

Statistical Performance Indices
Hydrological model performance must be validated using objective functions. This study used (R 2 , Nash-Sutcliffe efficiency (NSE), percent bias (PBIAS)) to evaluate the statistical performance of the SWAT and WEAP model, as presented in Tables 2 and 3. R 2 is the measurement of the co-relation between simulated and observed values, known as the coefficient of determination. The R 2 varies from 0 to 1, where approaching 1 means better performance. R 2 can be calculated as follows: where , are initially measured and simulated discharge, respectively. Nash-Sutcliffe efficiency (NSE) is a normalized dimensionless statistic that defines the relative magnitude of residual variance compared to the measured data variance [59]. NSE can be measured by the following equation:

Statistical Performance Indices
Hydrological model performance must be validated using objective functions. This study used (R 2 , Nash-Sutcliffe efficiency (NSE), percent bias (PBIAS)) to evaluate the statistical performance of the SWAT and WEAP model, as presented in Tables 2 and 3. R 2 is the measurement of the co-relation between simulated and observed values, known as the coefficient of determination. The R 2 varies from 0 to 1, where approaching 1 means better performance. R 2 can be calculated as follows: where Q m , Q s are initially measured and simulated discharge, respectively. Nash-Sutcliffe efficiency (NSE) is a normalized dimensionless statistic that defines the relative magnitude of residual variance compared to the measured data variance [59]. NSE can be measured by the following equation: where n is the total number of observations, Q o,i and Q s,i are the observed and simulated discharge at the ith observation, respectively, and Q mean is the mean observed data over the simulation period. NSE value ranges from −∞ to 1, with an optimal value of 1 [60]. PBIAS (percent bias) calculates the average tendency of simulated values to be larger or smaller than observed values [60]. PBIAS with a smaller magnitude shows better model performance. PBIAS can be represented as follows: where Q is discharge and the optimal value of PBIAS is zero. A positive PBIAS value indicates an underestimation of the model and negative values represent overestimation.

Scenario Development
In this study, the model was designed to simulate different socio-economic and climate change scenarios: Future climate change impacts were evaluated by comparing WEAP outputs between the reference and climate change extended scenarios. The baseline scenario considers the current account and sets the year as 2010, which provides actual water demand and supply for the system. The socio-economic reference-based scenario includes the current account data in the entire model and projected the outputs (2011-2050). The reference scenario represents the current actual situation that is modelled and projected under current year (2010) situations and conditions, considering assumptions of population growth rate, land use, and unit consumption.
The climate change extended scenario was designed by replacing climatic outputs by RCP outputs. The projections made by the climate models are not consistent with the "real" natural conditions, due to uncertainty and data errors in the models. CMIP5 results recently tried to fill this gap with finer model resolution and climate change scenarios. The projected precipitation and discharge generated by the SWAT model are shown in Figure 4 and indicate that there will be a slight increase in precipitation and discharge over HRB. The simulation results of SWAT model will be reflected in water demand for various sectors.

Baseline Scenario
Current water use was estimated based on current water supply and demand estimates for the year 2010. The reference scenario had four types of sectoral demand: agricultural, domestic, forests (including animal husbandry and livestock), and industries (including construction and fire water demand). Based on annual water use, the demand site data were established as 44,283 L per person annually for domestic purposes, including urban and rural. The total population according to the 2010 census was 6,350,000 in Hongshui River Basin (HRB). The total land area of agriculture was of 4,820,373 hectares, with 451.97 m 3 consumption per hectare annually. Forest area was 6,189,570 hectares, with 49.30 m 3 consumption per hectare annually. Figure 5 shows the monthly water demand for 2010 (current account) for the agriculture, domestic, forests, and industrial sectors with total demand of 4155 million cubic meters (Mm 3 ). Agriculture is the main sector of the basin with the highest annual water demand (2179 Mm 3

Reference Scenario
The reference scenario is the basic concept in the WEAP model, which can be used to manage various scenarios and reflects the actual role. In this model, the reference scenario was structured as a baseline scenario with population growth of 0.59% as per 2010 census data, and a 0.46% reduction in agriculture water consumption per year, due to advanced irrigation water techniques, 4.29% increase in forests, animal husbandry and livestock water consumption, and 1.24% increase in industrial water consumption based on previous years' data.

Reference Scenario
The reference scenario is the basic concept in the WEAP model, which can be used to manage various scenarios and reflects the actual role. In this model, the reference scenario was structured as a baseline scenario with population growth of 0.59% as per 2010 census data, and a 0.46% reduction in agriculture water consumption per year, due to advanced irrigation water techniques, 4.29% increase in forests, animal husbandry and livestock water consumption, and 1.24% increase in industrial water consumption based on previous years' data.

Reference Scenario
The reference scenario is the basic concept in the WEAP model, which can be used to manage various scenarios and reflects the actual role. In this model, the reference scenario was structured as a baseline scenario with population growth of 0.59% as per 2010 census data, and a 0.46% reduction in agriculture water consumption per year, due to advanced irrigation water techniques, 4.29% increase in forests, animal husbandry and livestock water consumption, and 1.24% increase in industrial water consumption based on previous years' data.
The overall outcome obtained from the reference scenario is shown in Figure 7. Based on the overall result, the value of unmet demand is zero. Therefore, water availability under the reference scenario (based on the previous population growth rate of 0.59 percent) is indicated to be reasonably sufficient. Figure 7 shows the projected water demand for the period 2010-2050. From the results, it was projected that agriculture water demand will likely decrease 49. Domestic demand will likely increase by 59 Mm 3 every five years and overall domestic demand will likely increased by 513.34 Mm 3 during 2010-2050. Industries and forestry, animal husbandry and livestock will likely have the highest water demand, with an average increase of 148 Mm 3 every five years, and livestock sector consumption will likely undergo an average increase of 105 Mm 3 every five years until 2050. The total water demand until 2050 will likely be 208 billion cubic meters (Bm 3 ), as shown in Table 4, with a breakdown of 39 percent for agriculture, 9.90 percent for forestry, animal husbandry, and livestock, 16.2 percent for domestic, and 34.87 percent for the industries sector. The overall outcome obtained from the reference scenario is shown in Figure 7. Based on the overall result, the value of unmet demand is zero. Therefore, water availability under the reference scenario (based on the previous population growth rate of 0.59 percent) is indicated to be reasonably sufficient. Figure 7 shows the projected water demand for the period 2010-2050. From the results, it was projected that agriculture water demand will likely decrease 49. Domestic demand will likely increase by 59 Mm 3 every five years and overall domestic demand will likely increased by 513.34 Mm 3 during 2010-2050. Industries and forestry, animal husbandry and livestock will likely have the highest water demand, with an average increase of 148 Mm 3 every five years, and livestock sector consumption will likely undergo an average increase of 105 Mm 3 every five years until 2050. The total water demand until 2050 will likely be 208 billion cubic meters (Bm 3 ), as shown in Table 4, with a breakdown of 39 percent for agriculture, 9.90 percent for forestry, animal husbandry, and livestock, 16.2 percent for domestic, and 34.87 percent for the industries sector.

Population Growth Scenario
A new scenario was introduced to predict the effect of potential situations on the model. The new scenario was generated to evaluate the impacts of positive and negative population growth rate for the Hongshui River Basin. This scenario also included a linear increase of 0.5% in agriculture area. The scenario was used to investigate the impact of a positive population growth rate ranging from 0.59% to 0.80% and a negative population growth rate ranging from 0.59% to 0.30%. Figure 8 shows the projections of water demand based on both positive and negative population growth scenarios. Water demand under positive population growth will likely be 218 Bm 3 , which is 10 Bm 3 more than under the reference scenario. Negative population growth results in a slight increase in projected water demand because a linear increase in agriculture area was added. The total water demand for this scenario was projected to be 215 Bm 3 , which is 7 Bm 3 more than in the reference scenario. It was noted that, in the scenario of a 0.50% increase in agriculture area and 0.80% increase in population, no unmet water demand was projected, implying that the water supplied each year will be sufficient.

Population Growth Scenario
A new scenario was introduced to predict the effect of potential situations on the model. The new scenario was generated to evaluate the impacts of positive and negative population growth rate for the Hongshui River Basin. This scenario also included a linear increase of 0.5% in agriculture area. The scenario was used to investigate the impact of a positive population growth rate ranging from 0.59% to 0.80% and a negative population growth rate ranging from 0.59% to 0.30%. Figure 8 shows the projections of water demand based on both positive and negative population growth scenarios. Water demand under positive population growth will likely be 218 Bm 3 , which is 10 Bm 3 more than under the reference scenario. Negative population growth results in a slight increase in projected water demand because a linear increase in agriculture area was added. The total water demand for this scenario was projected to be 215 Bm 3 , which is 7 Bm 3 more than in the reference scenario. It was noted that, in the scenario of a 0.50% increase in agriculture area and 0.80% increase in population, no unmet water demand was projected, implying that the water supplied each year will be sufficient.

Deforestation Scenario
Results of the deforestation scenario, shown in Figure 9, indicate a 4% decrease in total water demand. The forest area was linearly reduced by 1%, and this reduction in area will decrease the water demand stress. Future deforestation in the Hongshui River Basin is possible because of population expansion and rapid industrialization. (Table 5 shows the water demand for forests, animal husbandry, and livestock, for reference and deforestation scenarios.

Deforestation Scenario
Results of the deforestation scenario, shown in Figure 9, indicate a 4% decrease in total water demand. The forest area was linearly reduced by 1%, and this reduction in area will decrease the water demand stress. Future deforestation in the Hongshui River Basin is possible because of population expansion and rapid industrialization. (Table 5 shows the water demand for forests, animal husbandry, and livestock, for reference and deforestation scenarios.
The total water demand projected for forests, animal husbandry, and livestock were 25,796 Mm 3 for the deforestation scenario and 33,159 Mm 3 for the reference scenario, which represents a 22% reduction.  The total water demand projected for forests, animal husbandry, and livestock were 25,796 Mm 3 for the deforestation scenario and 33,159 Mm 3 for the reference scenario, which represents a 22% reduction.

Industrial Growth Scenario
Industrial water, which also includes construction and fire water, is the second-highest category of water consumption, following agriculture water consumption. This scenario provides important insight into the impact of industrialization on future water demand. The reference scenario included a 1.24% increase in water consumption until 2050. A further 2.5% increase in water consumption and a 1% linear increase in industrial sector area, due to rapid urbanization was added to this new scenario.
The total water demand projected for the industrial growth scenario was 258 Bm 3 , which is 19.6% more than in the reference scenario, as presented in Figure 10. The total water demand projected for industries will be 122 Bm 3 for the industrial growth scenario and 73 Bm 3 for the reference scenario, representing a 41% increase.

Industrial Growth Scenario
Industrial water, which also includes construction and fire water, is the second-highest category of water consumption, following agriculture water consumption. This scenario provides important insight into the impact of industrialization on future water demand. The reference scenario included a 1.24% increase in water consumption until 2050. A further 2.5% increase in water consumption and a 1% linear increase in industrial sector area, due to rapid urbanization was added to this new scenario.
The total water demand projected for the industrial growth scenario was 258 Bm 3 , which is 19.6% more than in the reference scenario, as presented in Figure 10. The total water demand projected for industries will be 122 Bm 3 for the industrial growth scenario and 73 Bm 3 for the reference scenario, representing a 41% increase.

Low Climate Change Extended Scenario (RCP-2.6)
Climate change scenarios are those scenarios developed by the Coupled Model Intercomparison Project Phase 5 (CMIP5) based on an unprecedented level of information on which to base projections, including new Earth system models with a more complete representation of force, new Representative Concentration Pathways (RCP) scenarios, and more output available for analysis. Global precipitation would certainly increase with increasing global mean surface temperature in the long term. Global mean precipitation will increase at a slower rate per degree Celsius than atmospheric water vapor. For scenarios other than RCP2.6, it will most likely rise by 1 to 3 percent • C −1 . The CMIP5 models' sensitivities for RCP2.6 range from 0.5 to 4 percent • C −1 at the end of the 21st century [61]. This scenario included statistically downscaled precipitation and temperature data of five GCMs, and the addition of streamflow data generated by the SWAT model was embedded in the WEAP model to project future water demand under the RCP-2.6 scenario.
The WEAP model predicted an increase in future water demand in the Hongshui River Basin according to the low emission scenario RCP-2.6. This scenario requires ambitious greenhouse gas emissions reductions and assumes a population of 9 billion at the end of 2100, an increase in croplands due to bio-energy production, and more intensive animal husbandry. In this scenario we considered these factors and results show an increase in water demand from domestic, agriculture, and forests, animal husbandry, and livestock sectors. The total water demand would likely increase under the low emission scenario RCP-2.6, as shown in Figure 11, to 231 Bm 3 , whereas in the reference scenario the total water demand was 208 Bm 3 with zero unmet water demand. Figure 11. The total water demand for the reference scenario and low climate change emission (RCP-2.6) scenario.

Medium Climate Change Extended Scenario (RCP-4.5)
RCP-4.5 is an intermediate emission scenario encompassing relatively ambitious emissions reductions, and the assumed outlook is consistent with strong reforestation, decreasing use of croplands and grasslands due to yield increases, and dietary changes and stringent climate policies. This study analyzed the ensemble output of five GCMs under RCP-4.5 and the WEAP model predicted an increase in all demand sectors.
This scenario included an increase in all demand sectors consistent with RCP-4.5 features. The total water demand was 258 Bm 3 , which is 24% more than in the reference scenario, as summarized in Figure 12 for each demand site. The maximum percentage projected changes were 54% in forests, animal, and husbandry, and 18% in the agriculture sector, to ensure the future food security. under RCP-4.5 and the WEAP model predicted an increase in all demand sectors.
This scenario included an increase in all demand sectors consistent with RCP-4.5 features. The total water demand was 258 Bm 3 , which is 24% more than in the reference scenario, as summarized in Figure 12 for each demand site. The maximum percentage projected changes were 54% in forests, animal, and husbandry, and 18% in the agriculture sector, to ensure the future food security.

High Climate Change Extended Scenario (RCP-8.5)
This scenario included the features of the RCP-8.5 high emission scenario. This RCP is consistent with a scenario of no future policy changes to reduce emissions. The scenario also assumes a global population of 12 billion by 2100, increased use of croplands and grassland, driven by an increase in population, a threefold increase in carbon dioxide emissions compared to current levels, and high energy intensity and high resilience on fossil fuels.
This scenario projected a significant amount of water demand as a result of the lack of implementation of new climate change policies, which leads to high methane and carbon emissions. The total water demand projected for RCP-8.5 will likely jump to 734 Bm 3 , which is approximately three times higher than in the reference scenario, as shown in Figure 13.

High Climate Change Extended Scenario (RCP-8.5)
This scenario included the features of the RCP-8.5 high emission scenario. This RCP is consistent with a scenario of no future policy changes to reduce emissions. The scenario also assumes a global population of 12 billion by 2100, increased use of croplands and grassland, driven by an increase in population, a threefold increase in carbon dioxide emissions compared to current levels, and high energy intensity and high resilience on fossil fuels.
This scenario projected a significant amount of water demand as a result of the lack of implementation of new climate change policies, which leads to high methane and carbon emissions. The total water demand projected for RCP-8.5 will likely jump to 734 Bm 3 , which is approximately three times higher than in the reference scenario, as shown in Figure 13. RCP-8.5 projections showed a likely significant increase in water demand in industries, and forests, animal, and husbandry, followed by agriculture, as shown in Figure 14, resulting in large quantities of carbon and methane emissions. The total water demand (2010-2050) projected by RCP-8.5 for the agriculture sector was 171 Bm 3 ; for forests, animal, and husbandry, projected demand was 100 Bm 3 ; and for industries, projected demand was 442 Bm 3 , which is five times higher than in the reference scenario. RCP-8.5 projections showed a likely significant increase in water demand in industries, and forests, animal, and husbandry, followed by agriculture, as shown in Figure 14, resulting in large quantities of carbon and methane emissions. The total water demand (2010-2050) projected by RCP-8.5 for the agriculture sector was 171 Bm 3 ; for forests, animal, and husbandry, projected demand was 100 Bm 3 ; and for industries, projected demand was 442 Bm 3 , which is five times higher than in the reference scenario. The total unmet water demand projected for RCP-8.5, as shown in Figure 15, was 5,428 Mm 3 in the period from 2039 to 2050. Forests, animal, and husbandry will likely undergo the highest unmet water demand, of 2,876 Mm 3 , followed by the unmet water demand of the agricultural sector of 2,372 Mm 3 .

Discussion
This study assessed the water demand focused on various socioeconomic, management, and climate change scenarios for sustainable availability of future water resources in Hongshui River Basin (HRB), China. To establish a sustainable water supply-demand system, we considered the current water supply-demand system in the study area and determined the potential impacts of climate change impacts of extreme drivers of demand for HRB water resources. Based on the current supply-demand situation, supply-and The total unmet water demand projected for RCP-8.5, as shown in Figure 15, was 5428 Mm 3 in the period from 2039 to 2050. Forests, animal, and husbandry will likely undergo the highest unmet water demand, of 2876 Mm 3 , followed by the unmet water demand of the agricultural sector of 2372 Mm 3 . The total unmet water demand projected for RCP-8.5, as shown in Figure 15, was 5,428 Mm 3 in the period from 2039 to 2050. Forests, animal, and husbandry will likely undergo the highest unmet water demand, of 2,876 Mm 3 , followed by the unmet water demand of the agricultural sector of 2,372 Mm 3 .

Discussion
This study assessed the water demand focused on various socioeconomic, management, and climate change scenarios for sustainable availability of future water resources in Hongshui River Basin (HRB), China. To establish a sustainable water supply-demand system, we considered the current water supply-demand system in the study area and determined the potential impacts of climate change impacts of extreme drivers of demand for HRB water resources. Based on the current supply-demand situation, supply-and

Discussion
This study assessed the water demand focused on various socioeconomic, management, and climate change scenarios for sustainable availability of future water resources in Hongshui River Basin (HRB), China. To establish a sustainable water supply-demand system, we considered the current water supply-demand system in the study area and determined the potential impacts of climate change impacts of extreme drivers of demand for HRB water resources. Based on the current supply-demand situation, supply-and demand-side management scenarios were applied in the model to optimize the water management system to achieve better results. In this study, a reference scenario was constructed based on the current water demand and supply in the Hongshui River Basin, which was estimated to be 4,155 Mm 3 . Similar studies were conducted by Matchett and Fleskes [62][63][64] using the WEAP model to analyze future water demand based on current and reference scenarios. In the current study, the total water demand for the reference scenario in 2050 was projected to be 55% more than current year water demand. Based on the reference scenario, water demand will increase by 38% for domestic, 60% for industries and 63% for forests, animal husbandry, and livestock. Similar results were reported by Asghar, et al. [65], who concluded that, based on the reference scenario, water demand in the central Indus Basin will increase by 11% for domestic and 55% for agriculture and livestock sectors. A similar study conducted on water demand in West Ethiopia's Didessa sub-basin found the existing water demand was 74 Mm 3 . Their study findings will help solve water management issues in the basin [66]. Another similar study conducted in the Mae Kong Basin in Thailand estimated an average unmet demands for agriculture of 62 and 17 Mm 3 under differing climate change scenarios [67]. In the current study, the WEAP model was also used to analyze external drivers (population growth, agricultural growth, deforestation, industrialization, and climate change), which place significant stress on the existing water supply system. These demands intensify water scarcity in the catchment, particularly downstream, and encourage water resource managers to implement water management policies in the Hongshui River Basin (HRB). An efficient management approach must be implemented based on the results of this report. A comparative analysis was designed to consider alternative efficient strategies for water security in the basin.
The integrated SWAT-WEAP methodology was applied in an efficient way to evaluate the technical management of water resources. This study shows that the HRB's high water management pressure stems predominantly from demand in agriculture and industry. The reference scenario has a negative effect on the HRB supply system, via rising water demand, which is expected to double by 2050. This rapid rise in demand is not only due to anthropogenic activities, but also to changes in rainfall and seasonal variation. Socio-economic scenarios, including population growth (positive and negative growth), industrial growth and deforestation scenarios, were analyzed in this study. Each of these socio-economic scenarios would likely result in an increase in water demand and zero unmet water demand. The available water will likely be sufficient for two cities (He chi and Laibin) in the basin for the period 2010-2050. Because downstream cities are more populated and support large industries, they may face water shortages, which will be covered in future studies. Touseef, Chen, Masud, Khan, Yang, Shahzad, Ijaz and Wang [9] concluded that the water yield is likely to decline by 30 percent for the future period (2020-2050) under all future scenarios, which supports the finding that water demand will likely increase. Sufficient water is likely to be available until 2050 in the socio-economic scenario and two climate change scenarios because major drought is not projected in the Xijiang River Basin. The Xijiang Basin likely underwent three significant cycles of hydrological drought (from the 1950s to the mid-1960s, during the late 1980s, and from the mid-1990s to the beginning of 2010) and two periods of drought relief (from the late 1960s to the mid-1980s, and the late 1990s to the early 2000s) [68]. Another study found that there will be growing pressure on the water supplies of Binhai New Area (BHNA), China, and numerous suggestions have been submitted to aid policymakers in water management to meet the potential water demand in the region [69]. A similar study conducted at Songhua River Basin, Northeastern China, concluded that the implementation of the WEAP model indicated a maximum decrease in water flows of 24% in potential climate change scenarios before the end of the 21st century. The drought indicators showed that serious to extreme drought events are likely to be most evident in 2059, 2060, and 2085, and mild to moderate drought events may occur in the study area in the remaining years [70].
The current study also modelled water demand using an ensemble of five GCMs under three RCPs: the low climate scenario (RCP-2.6), medium climate scenario (RCP-4.5) and high climate scenario (RCP-8.5). The results showed that water demand will likely increase under all climate change scenarios. The total water demand projected under low and medium climate scenarios was 231 and 258 Bm 3 , respectively, which represented increases in water demand of 11% and 24% compared to the reference scenario. The higher climate change scenario (RCP-8.5) will significantly increase water demand, by threefold relative to the reference water demand. Thus, it is important to implement policies to maintain the temperature rise and emissions, thereby lowering their effect and reducing water demand. Industrial policies will be important to ensure future water demands can be met. The water availability in Hongshui River Basin (HRB) will likely be sufficient for the two adjacent cities until 2050, but this finding will not uniformly apply to downstream cities. Further study will be required to analyze the water demand beyond 2050 and for downstream cities. The current results are supported by Wang, et al. [71], who summarized that the A1B scenario will increase water scarcity in river basins in the north of China, but the water situation will improve in river basins in Southern China. The results of another study differed slightly from our findings and showed that total water demand will increase in most basins of China, but with a declining share of agriculture, mainly due to commercial, domestic, and urban competition. In general, water availability for agriculture will decrease in southern China and remain stable in northern China [72].

Conclusions
This study integrated two models, the SWAT model for future climate prediction, and the WEAP model for the simulation of water supply and demand in the Hongshui River Basin (HRB), to evaluate the impacts of climate change and socio-economic scenarios. The results of the SWAT model showed that, in the HRB, temperature, precipitation, and water discharge will increase slightly until 2050. The WEAP model was developed and showed that the two major drivers of future water shortages in the basin will be climate variables and socio-economic factors. This study analyzed three socio-economic and three climate change scenarios. The results indicated that agriculture and industries are the major factors that will place pressure on the water resources in the future. Domestic water is sufficient for two cities located in the basin (He chi and Laibin) until 2050. A total unmet water demand of 5 Bm 3 would likely occur in the basin under RCP-8.5. All scenarios showed an increase in water demand. Therefore, we propose that better long-term management policies should be implemented in the watershed to ensure the sustainability of water resources. We also propose to use a range of GCMs with the latest RCPs, in addition to the five examined here, in future research on the Hongshui River Basin. These studies will help explain the uncertainties of modelling and simulation. Nonetheless, for research into future climate change scenarios and hydrological trials in this area, the high-resolution climate data used in this study performed well. Based on the results of the analysis of the socio-economic and climate changes scenarios, the study has high significance for water management and stakeholders of water supplies.
Author Contributions: Muhammad Touseef devised the idea, performed the analysis, designed the study, and wrote the paper. Lihua Chen provided technical support and supervision, and Wenzhe Yang helped in data collection. All authors have read and agreed to the published version of the manuscript.