Next Article in Journal
Shrinking Desert Channel Response to Increasing Human Interferences and Changing Natural Factors in the Upper Yellow River
Next Article in Special Issue
A CMIP6 Multi-Model Analysis of the Impact of Climate Change on Severe Meteorological Droughts through Multiple Drought Indices—Case Study of Iran’s Metropolises
Previous Article in Journal
Aquaculture Development in Nigeria: The Second Biggest Aquaculture Producer in Africa
Previous Article in Special Issue
Future Bioclimatic Change of Agricultural and Natural Areas in Central Europe: An Ultra-High Resolution Analysis of the De Martonne Index
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Climate Change Impacts on Agricultural and Industrial Water Demands in the Beijing–Tianjin–Hebei Region Using Statistical Downscaling Model (SDSM)

1
School of Economics and Management, North China Electric Power University, Beijing 100096, China
2
Hebei Institute of Water Science, Shijiazhuang 050057, China
3
State Key Laboratory of Desert and Oasis Ecology, Key Laboratory of Ecological Safety and Sustainable Development in Arid Lands, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi 830011, China
*
Author to whom correspondence should be addressed.
Water 2023, 15(24), 4225; https://doi.org/10.3390/w15244225
Submission received: 14 October 2023 / Revised: 27 November 2023 / Accepted: 5 December 2023 / Published: 8 December 2023
(This article belongs to the Special Issue Impacts of Climate Change on Water Resources: Assessment and Modeling)

Abstract

:
As a politically and culturally important city cluster, the Beijing–Tianjin–Hebei (BTH) region is the most prominent area in China where the imbalance between the supply and demand of water resources restricts the sustainable and healthy development of the regional social economy. In the context of global warming, research into water demand prediction that takes climate change into consideration would be more in line with the strategic goal of the low-carbon sustainable development of future cities. At the same time, the prediction of agricultural water demands against a background of climate change is urgently needed, while industrial water consumption is weakly correlated with climate change, an investigation of the statistical relationship between the two is needed. Thus, in this paper, future climate data from the BTH region under the scenarios RCP2.6, RCP4.5 and RCP8.5 were generated using a statistical downscaling model, and then coupled with agricultural and industrial water demand prediction models to simulate and analyze the impact of climate change on the agricultural and industrial water demands, respectively. The results show that during the forecast period (2020–2035), the reference crop evapotranspiration (ET0) growth rates in the Beijing, Tianjin and Hebei areas under the RCP2.6 scenario are 1.438 mm·a−1, 1.393 mm·a−1 and 2.059 mm·a−1, respectively. Under the RCP4.5 scenario, they are 2.252 mm·a−1, 2.310 mm·a−1 and 2.827 mm·a−1, respectively. Under the RCP8.5 scenario, they are 3.123 mm·a−1, 2.310 mm·a−1 and 2.141 mm·a−1, respectively. Furthermore, under each climate scenario, the increase in evapotranspiration in the Hebei area is the largest, followed by that in the Tianjin area, and that in the Beijing area is the smallest. For water consumption per CNY 10,000 of industrial added value during the forecast period, under the three different climate scenarios, a downward trend is seen in the Beijing area, with rates of 0.158, 0.153 and 0.110, respectively, but in the Tianjin area, there is an upward trend, with an upward tendency in rates of 0.170, 0.087 and 0.071, and an upward trend in the Hebei area, with an upward tendency in rates of 0.254, 0.071 and 0.036, respectively. This study will help the BTH region to rationally allocate agricultural and industrial water against the background of future climate change, and strengthen the coordination and cooperation between the different regions to promote the healthy and sustainable development of the cities.

1. Introduction

Climate change refers to changes in the average state of the climate system caused by natural factors and human activities. Climate change can have a significant impact on hydrological processes [1], and can also affect the demand for water resources. With successive increases in global warming, there are increasing regional changes in soil moisture, precipitation, and mean temperature, and changes in extreme weather are expected to be greater in intensity and frequency [2]. The uneven spatial and temporal distribution of water resources is a prominent feature of China’s water resources distribution. Under the dual impact of climate change and human activities, many city clusters in China are facing increasingly severe challenges to water security [3].
China is vigorously promoting the development strategy of the city cluster, and the development of city clusters is closely related to changes in water demand. As the economic core region in northern China, the Beijing–Tianjin–Hebei (BTH) city cluster has a poor natural endowment of water resources, but its population and economy are especially concentrated. Furthermore, it is affected by rapid industrialization and urbanization, and has become one of the regions most affected by the global water cycle disturbance [4]. On the whole, the contradiction between regional sustainable development and the serious shortage of water resources carrying capacity in the BTH region has existed for a long time and has been deepening.
Thus, this paper takes the BTH city cluster as the research area and tries to solve the contradiction.
Of all industries, agriculture is the most sensitive to the effects of climate change. Agricultural water includes water for crops and water for forestry, animal husbandry and freshwater aquaculture. The crop water requirement is the total amount of water consumed by transpiration by the plants and soil evaporation between the plants during the whole growth period and is an important factor that must be considered by irrigation decision-making and water resource planning [5]. Compared with that of crops, water consumption by forestry, animal husbandry and aquaculture industries is less affected by climate and is difficult to quantify. Therefore, in terms of agricultural water demand, this paper mainly considers the change in water demand related to crops.
Reference crop evapotranspiration (ET0) is the only factor that appears in both the water balance and the surface energy balance equations [6], and is an important part of the crop water requirement calculation performed in climate and hydrological research [7]. ET0 reflects the atmospheric evaporation capacity and the vegetation water demand. Only meteorological factors can affect ET0 without considering the influence of crop characteristics and soil conditions. Researchers have carried out in-depth studies of ET0 from general areas to regions with special climatic and geographical conditions, such as around the Mediterranean Sea [8,9]. In such studies, information describing crop water consumption is typically obtained by recording regional evapotranspiration and surface energy flux [10]. In the context of climate change, the study and analysis of temporal and spatial variations and complex periodic variations in evapotranspiration can provide a relevant basis for irrigation water management [11,12]. With progress in research and the continuous development of technology, regional evapotranspiration can now be calculated more accurately by integrating innovative technologies such as remote sensing, or by improving the precision of the crop coefficient, so as to ensure regional water security and crop yield through good water resource management [13,14]. Currently, the majority of studies on evapotranspiration primarily focus on investigating its sensitivity to various meteorological factors [15,16,17] and its variations across different geographical conditions and climates [18,19]. Specifically, research on evapotranspiration in the BTH area mainly centers around understanding its temporal and spatial changes as well as historical patterns [20]. This paper aims to predict future trends in evapotranspiration within the BTH area under diverse climatic conditions, while also providing practical recommendations for regional water resource management cooperation based on forecast results and the developmental status of the respective regions, thereby facilitating enhanced regional strategic coordination.
Industrial water management has an important impact on urban water consumption. Due to the continuous advancement of industrialization and urbanization, efficient urban water demand management is becoming more important, and the impact of industrial operations on the effective use and sustainability of water resources is attracting increasing attention [21]. At present, for industrial water demand, the main prediction methods are the quota method, regression analysis, grey theory prediction, the neural network method, and system dynamics [22,23,24]. Compared with socio-economic factors, meteorological factors have a relatively weak impact on industrial water demand. In view of this, most studies of industrial water demand focus on the impact of socio-economic factors [25,26] and ignore the impact of climate factors. Current relevant studies have included considerations such as industrial added value, industrial water repetition rate, an industrial technological progress index, industrial water price and policy-based water-saving rates, and various combinations of these have been incorporated into forecast models of industrial water demand [27,28]. In addition, a challenge arises in the process of establishing an industrial water demand prediction model in China due to excessively large forecast values. Therefore, this paper aims to focus on the impact of climate change on industrial water demand by incorporating meteorological factors into the existing model. This approach not only enables us to predict the changing trends in industrial water demand but also provides a method for quantifying the relationship between climate change and industrial water demand.
In a departure from the literature, our main contribution is that we apply the method of coupling an SDSM and a water demand model to China’s important region, BTH. First of all, in order to improve the spatial resolution of future climate information output by atmospheric circulation models, the statistical downscaling model (SDSM) was established for the BTH region, and a 95% confidence level was selected in the calibration process to reduce the uncertainty of simulation results. The calibration results show that the model can explain the daily temperature change in the study area well. Then, three climate scenarios RCP2.6, RCP4.5 and RCP8.5 were simulated based on the SDSM constructed by the CanESM2 atmospheric circulation model, and the daily maximum temperature, minimum daily temperature and daily precipitation data of meteorological stations in the BTH region from 2020 to 2050 were simulated. The data of Beijing, Tianjin and Hebei were obtained by Kriging interpolation of weather station data by ArcGIS. Secondly, the industrial water demand model established in this paper is different from other studies that mainly analyzed the impact of social and economic factors such as industrial water price, water use repetition rate and policy saving rate on industrial water demand. On that basis, this paper also fitted climate factors into the industrial water demand model through correlation analysis and multiple linear regression. Finally, meteorological data were fed into industrial and agricultural water demand models to simulate and predict the change in water demand against the background of climate change. The results can be compared across three climate scenarios to obtain the impact of future climate change on agricultural and industrial water demands. Additionally, inter-regional comparisons of forecast results can be conducted. This study focuses on the unique city cluster of BTH as the research region. After predicting the changing trend in water demand under various climate scenarios, this paper emphasizes the promotion of strategic cooperation in urban water resources and aims to provide relevant foundations and feasible recommendations for future water resource management collaboration by integrating the changing trends and developmental status of water demands across different regions.
On the whole, this study mainly used the SDSM to generate precipitation and temperature climate scenarios related to representative concentration path scenarios (RCP 2.6, RCP 4.5 and RCP 8.5). The simulation results were coupled with agricultural and industrial water demand models to compare the forecast results of water demand in different regions and under different climate scenarios, so as to provide a relevant basis for future agricultural and industrial water resources management in the BTH region, and at the same time make relevant policy recommendations to cope with climate change.

2. Study Area and Data Sources

2.1. Study Area

The BTH region is located between 113°04′~119°53′ east longitude and 36°01′~42°37′ north latitude, bordered by Bohai Bay to the east, the Taihang Mountain region to the west, the Yellow River area to the south, and the Mongolia Plateau to the north. The study area covers an area of 217,000 km2, accounting for 2.3% of the national land area. Within the BTH region, resource distribution, climate characteristics, and water resource systems are all similar and the region is a resource-based water-shortage area in China.
Total water resources in 2021 in the BTH region were 6.13 billion m3 (Beijing), 3.98 billion m3 (Tianjin), and 37.66 billion m3 (Hebei), accounting for about 1.6% of the national total water resources (2964 billion m3). Hebei is an important grain-producing area in China and agricultural water consumption in this province has reached 53.4% of the total regional water consumption. Tianjin has a relatively equal distribution of all kinds of water consumption, while Beijing has a relatively high proportion of domestic water consumption. The BTH region contributes about 10% of the national GDP, 6% of the food production and contains 8% of the population, and currently limited water resources restrict the sustainable development of the region [29].

2.2. Data Sources

The climate data used in this paper mainly include historical meteorological observational data and CMIP5 climate model data. The daily measured meteorological data of 20 meteorological stations in the BTH region were obtained from the Water Resources Bulletin, Water Statistical Yearbook and provincial and municipal statistical yearbooks of the BTH region, including daily minimum temperature, maximum temperature and precipitation for the period 1973 to 2005.
CMIP5 model data was obtained from the first set of NCEP reanalysis data on the website of https://psl.noaa.gov accessed on 12 March 2023, including 26 atmospheric circulation forecast factors. The atmospheric circulation model was HadCM3, the climate model was CanESM2. Based on China’s current carbon peaking and carbon neutrality policies in response to global warming, and considering the frequency of extreme weather in the future, this study selected representative concentration paths RCP2.6, RCP4.5 and RCP8.5 as proposed by the IPCC for future hypothetical predictions.
For the industrial water demand model, the data involved in the annual industrial water repetition rate and industrial water price in the BTH region from 2000 to 2019 were derived from the China Environmental Statistical Yearbook.

3. Methodology

3.1. Statistical Downscaling Model

The SDSM used in this study is a combination of multiple linear regressions and random weather generators proposed by Wilby and Dawson [30]. It is practical, simple and easy to operate [31]. The SDSM utilizes outputs from a global climate model (GCM). This is a mathematical model that describes the physical processes of the Earth’s atmosphere, ocean and land based on the Navier–Stokes equations of a rotating sphere. It is widely used in climate change research and weather forecasting [32]. However, current GCMs only provide future climate scenarios at large scales and cannot be directly applied at regional scales [33]. To overcome this limitation, the SDSM can use multiple regression technology to extract station-scale climate information from GCM output results. Using the coupling principle of random weather generators and multiple regressions, the SDSM can establish relationships between GCM variables and local scale variables. The SDSM is then used in two stages (Figure 1): (1) establishment of the statistical relationship between predictors from regional or station meteorological data and predictors from a GCM; (2) generation of future local-scale daily time series data using the GCM output data. The SDSM is widely used in climate change studies around the world [34,35].
Figure 1. Flowchart illustrating the operation of the SDSM. The predictands (i.e., daily observed maximum temperature, minimum temperature and rainfall) are first extracted for 20 meteorological stations in the BTH region. Time series of 26 predictor variables are then extracted from NCEP for the same time period (1973 to 2000). Correlation analysis is then performed between the predictands and the 26 predictor variables to obtain the statistical relationship between the predictors and the predictands. The results selection of the predictors for each station are presented in Table 1.
Figure 1. Flowchart illustrating the operation of the SDSM. The predictands (i.e., daily observed maximum temperature, minimum temperature and rainfall) are first extracted for 20 meteorological stations in the BTH region. Time series of 26 predictor variables are then extracted from NCEP for the same time period (1973 to 2000). Correlation analysis is then performed between the predictands and the 26 predictor variables to obtain the statistical relationship between the predictors and the predictands. The results selection of the predictors for each station are presented in Table 1.
Water 15 04225 g001
After the statistical relationships between predictors and predictands had been established, the SDSM was calibrated. Historical parameters were inputted to the model and the model was run. The simulation results were then obtained and compared to the actual observational data to verify the accuracy of the model. After model calibration, daily data for the period 2020 to 2050 were obtained from HadCM3 and used as input to the SDSM in order to predict daily maximum and minimum temperature time series at each of the 20 meteorological stations under the three future climate scenarios of RCP2.6, RCP4.5 and RCP8.5 (Table 2). ArcGIS 10.6 was then used to perform Kriging spatial interpolation using the meteorological station predictions.
The three climate scenarios selected in this paper correspond to different degrees of climate intervention policies and greenhouse gas emission levels. Through the comparative analysis of the forecast results under different climate scenarios, the influence of different degrees of climate policy intervention on them can be indirectly obtained. At the same time, the uncertainty of future forecast results would be evaluated better.

3.2. Agricultural Water Demand Model

In order to quantitatively evaluate the impact of climate change on crop water requirements, reference crop evapotranspiration was selected as the evaluation index. Crop water requirement can be estimated using the product of ET0 and a crop coefficient [36]. However, the crop coefficient varies in space and time, due to inherent variability in emergence date, land use pattern, antecedent precipitation, emissivity, vegetation amount, and atmospheric boundary conditions [37]. Hence, to avoid further error, we use ET0 alone to reflect the degree of impact of climate change on overall crop water requirements.
In this study, the Hargreaves formula was used to calculate ET0 at each of the 20 meteorological stations. The method estimates the reference crop evapotranspiration (ET0) value using temperature and solar radiation [38]. Compared with the traditional Penman–Monteith formula, this method does not require climate data that is unavailable from the observational meteorological stations [39,40]. Temperature and solar radiation are easy to measure quantitatively, and are used to estimate ET0 as follows:
E T 0 = 0.0023 R a λ × ( T m e a n + T o f f ) × T m a x T m i n   0.5
where ET0 is the reference crop evapotranspiration (mm·d−1); λ is the latent heat of water vapor (2.45 MJ·kg−1); Tmean is the daily average air temperature (°C); Toff is a temperature constant (recommended value is 17.8 °C [41]); Tmax is the daily maximum air temperature (°C); Tmin is the daily minimum air temperature (°C); and Ra is the radiation intensity of the upper atmosphere (MJ·m−2·d−1) and is given by:
R a = 118.08 π d r ω s sin ( φ ) sin ( δ ) + cos ( φ ) cos ( δ ) sin ( ω s )
ω s = arccos tan ( φ ) tan ( δ )
d r = 1 + 0.033 cos 2 π D J
δ = 0.409 sin 2 π D J 1.39
where dr is the reciprocal of the relative distance between the sun and the Earth; ωs is the solar hour angle (rad); φ is the geographical latitude (rad); δ is the solar magnetic decrement angle (rad); J is the ordinal number of a day in a year, ranging from 1 (1 January) to 365 (31 December) in common years and from 1 to 366 in leap years; and D represents the number of days in the current year, 366 in leap years and 365 in common years.
After calculating ET0 time series using the observed and predicted meteorological data for each of the observational stations located in the BTH region, average ET0 values for the Beijing, Tianjin and Hebei sub-areas were estimated using spatial interpolation. The data was then divided into the baseline period (2005–2019) and the forecast period (2020–2035) with a period of 15 years. Comparing the baseline period as a whole with the forecast period, the impact of future climate change on the overall evapotranspiration change trend in the forecast period can be more intuitive. And it would be easier to evaluate the ET0 trends for each sub-area under the three climate scenarios (i.e., RCP2.6, RCP4.5 and RCP8.5)

3.3. Industrial Water Demand Model

The quota method [42] was used to calculate the industrial water demand (Windustry). This method considers the influence of climate change on the water consumption per CNY 10,000 of industrial added value. Windustry (m3) is estimated as follows:
W i n d u s t r y = u i n d u s t r y × G D P i n d u s t r y
where GDPindustry is the industrial added value (CNY 10,000); and uindustry represents the water consumption per CNY 10,000 of industrial added value (m3·10,000 yuan−1).
The industrial water demand prediction model of Yingjie et al. [43] was used to estimate values of uindustry. This model utilizes an industrial water price, an industrial water reuse rate and the air temperature as follows:
ln u i n d u s r y = b 0 + b 1 ln x 1 + b 2 ln x 2 + b 3 ln x 3
where b0, b1, b2, b3 denote regression coefficients; x1 is the utilization rate of industrial water; x2 is the industrial water price (yuan·m−3); and x3 is the absolute change in temperature.
The statistical software package SPSS 26.0 was used to carry out correlation analysis and multiple linear regressions on the data from 2000 to 2019 in order to establish the regression coefficients b0, b1, b2, b3. Predictions of x1, x2 and x3 for the period 2020 to 2035 under the three climate scenarios (RCP2.6, RCP4.5 and RCP8.5) were obtained by SDSM and various planning data in relevant documents. They were inputted to the fitted formula to obtain the water consumption per CNY 10,000 of industrial added value and evaluate the simulation effect of the established model.

4. Results and Discussion

4.1. SDSM in the BTH Region

The SDSM was established using data for the period 1973 to 1990, and the model performance was then evaluated using data from the period 1991 to 2005. Two evaluation parameters were calculated for the validation period, the standard error (SE) and the explanatory variance (R-squared), and these are summarized in Table 3.
It can be seen from Table 3 that the explained variance of the model for daily maximum temperature and daily minimum temperature simulations varied between 0.559 and 0.705, with most values being greater than 0.6. This suggests that the SDSM explained the daily temperature variation well. Due to the complex factors affecting precipitation, the explained variance of daily precipitation was relatively lower, but to a certain extent the SDSM can be assumed to simulate precipitation to an acceptable degree. In general, the SDSM demonstrated reasonable performance and hence can be used to generate predictions of daily maximum air temperature, minimum air temperature and precipitation under RCP2.6, RCP4.5 and RCP8.5 climate scenarios.

4.2. Daily Maximum and Minimum Temperatures under Future Climate Scenarios

The simulated changes in daily maximum and minimum air temperatures in the BTH region under the three climate scenarios (RCP2.6, RCP4.5 and RCP8.5) for the period 2020 to 2050, are presented in Figure 2.
There is a difference of 0.14 to 1.55 °C between the forecast period (2020–2050) and the base period (2005–2019). The difference was greatest in the RCP2.6 scenario. The data for the forecast period are shown in the Table 4 below.
During the simulation period, the daily maximum and daily minimum temperatures in the BTH region fluctuate and show an overall upward trend. Under the RCP8.5 climate scenario, the daily maximum temperature and daily minimum temperature in each region increase the most, while under the RCP2.6 and RCP4.5 climate scenarios, the daily maximum temperature and daily minimum temperature increase less. This shows that governments need to implement policy interventions to effectively control carbon emissions, reduce the impact of greenhouse effect as much as possible, and control the temperature increase in future climate change.

4.3. Evapotranspiration under Future Climate Scenarios

The predicted evapotranspiration results of the three sub-areas of the BTH region under the three climate scenarios are presented in Figure 3.
It can be seen from the data listed in the Table 5 that, under the three climate scenarios, the order of ET0 growth in the three regions of Beijing, Tianjin and Hebei is always RCP8.5 > RCP4.5 > RCP2.6 from 2005 to 2035.
A comparison of ET0 variation trends in the BTH region under the three climate scenarios shows that, although ET0 in the region is on the rise under all three climate scenarios, the increase in ET0 in each sub-area is the largest under the RCP8.5 climate scenario, followed by the RCP4.5 climate scenario, and under the RCP2.6 climate scenario, the increase in ET0 in each region is the smallest. The main reason for this situation is that in the RCP8.5 climate scenario, the future temperature increase is relatively large, and more water will be released by soil evaporation and crop transpiration. As a result, the increase in evapotranspiration will be large, and the agricultural water demand will also increase. At the same time, the differences between climate scenarios also indicate that for future climate change and even some extreme climates, the negative impact on agricultural water use will be effectively reduced under the condition of government policy intervention.
The regional comparison of ET0 in the three sub-areas shows that the overall increase in ET0 in Hebei is the largest, followed by that in Tianjin, with the smallest being in the Beijing area. The main reason for this result is that agricultural water accounts for the largest proportion of overall water consumption in the Hebei area, while agricultural water accounts for the smallest proportion in the Beijing area. Therefore, under the same climate scenario, the fluctuation in water demand caused by climate influence in the Hebei area is greater, while the fluctuation in the Beijing area is relatively small. The results can provide a reference to guide overall crop irrigation and water management in the BTH region in the future. According to the specific fluctuations in evapotranspiration in different regions under different climate scenarios in the future, combined with relevant factors such as regional water use structure and agricultural irrigation technology level, we can propose corresponding practical measures to deal with the negative impact of climate change on water demands in the future. In addition, the strategic cooperative relationships between city clusters can be fully utilized to coordinate the supply and demand of water resources between regions and formulate corresponding cooperation policies, so as to facilitate the rational distribution of water resources in order to guarantee the safety of agricultural water utilization.
Based on the above comparative analysis results of evapotranspiration changes in different climate scenarios and different regions, we propose the following suggestions. From the obvious gap between the increase in evapotranspiration under the three climate scenarios, it can be seen that under the trend of global warming in the future, if no active measures are taken to control carbon emissions and global temperature is allowed to rise, the increase in evapotranspiration will lead to an increase in agricultural water demand in many regions like the BTH region, especially in areas like Hebei where agricultural water consumption is dominant. However, effective human intervention in carbon emissions through the formulation of relevant policies, such as the two-carbon policy currently formulated in China, will achieve good results in the future when dealing with the negative effects of climate change on water resources. In view of the management of agricultural water use in the three regions against the background of future climate change, as the proportion of agricultural water consumption in Hebei and Tianjin is large, emphasis should be placed on improving irrigation technology to control agricultural water demand.

4.4. Industrial Water Consumption under Future Climate Scenarios

The water consumption per CNY 10,000 of industrial added value in the BTH region under the three different climate scenarios of RCP2.6, RCP4.5 and RCP8.5 is presented in Figure 4.
As we can see (Table 6), under the three climate scenarios of RCP2.6, RCP4.5 and RCP8.5, the overall change in water consumption per CNY 10,000 of industrial added value showed a downward trend in Beijing, but an upward trend in Tianjin and Hebei.
A comparison of the water consumption per CNY 10,000 of industrial added value in the BTH region under the three different climate scenarios shows that there are no significant differences in the change trend of industrial water demand in each sub-area. The reason for this is that the influence of temperature on industrial water demand is relatively weak. The change in industrial water demand is mainly controlled by regional water price, water repetition rate and water saving rate. By comparing the water consumption per CNY 10,000 of industrial added value between the three sub-areas, it can be concluded that the order of the water consumption per CNY 10,000 of industrial added value is Hebei > Tianjin > Beijing under all three climate scenarios. In addition, the future water consumption per CNY 10,000 of industrial added value in the Beijing area will decrease, while that in the Tianjin and Hebei areas will increase. The main reason for this is that there is a big difference between the water price and the industrial water repetition rate in the three sub-areas.
The above results show that the model of industrial water demand used in this study, which includes temperature as a factor, predicts a very limited impact from climate change on industrial water demand in the BTH region, but it also provides a reference calculation method for the change in industrial water demand in the BTH region under the extreme climate conditions that may occur in the future; that is, by using multiple linear regression and correlation analysis, climate factors such as temperature and rainfall are included in the calculation model of industrial water demand. In addition, the differences in the trends in industrial water demand among the three regions are mainly caused by social and economic factors such as the repetition rate of industrial water use and the price of industrial water. The industrial development status of the three regions is very different, which indicates that we should focus on combining the characteristics of industrial development and strategic cooperation in each region to provide a relevant basis for the overall industrial transfer in the BTH region in the future.
According to the industrial development of the different regions in recent years, the industrial water consumption in the Beijing area has been reduced due to the relocation of enterprises. The expansion of infrastructure investment in Tianjin has led to an increase in industrial water consumption. Hebei Province, as the sub-area with the highest water demand in the BTH region, has a much larger industrial scale than Beijing and Tianjin. Although its industrial water consumption has increased due to some relocating of manufacturing industry, improvements in both industrial development and water use efficiency have gradually reduced this consumption. According to the results of this study, industrial water usage in the three regions will continue to change with the same trend in each climate change context up to 2035. Therefore, for the industrial water allocation in the BTH region against the background of future climate change, targeted measures should be implemented according to the characteristics of industrial development in each region. For example, Tianjin should address the input–output imbalance of a number of resources, including water resources, while Hebei should focus on improving water efficiency. At the same time, with the continuous development of the coordinated development policy of the BTH region, the overall process of promoting industrial transfer in Beijing, Tianjin and Hebei should be “upgraded first and then transferred”, so as to achieve mutual benefit and a win-win situation.

4.5. Uncertainty Analysis

In this paper, a statistical downscaling method was used to simulate the variation trend of ET0 and the industrial water consumption in the BTH region under three typical concentration paths (RCP2.6, RCP4.5 and RCP8.5) using output data from the CanESM2 model in the CMIP5 model. It should be noted that, in this paper, only the CanESM2 model data was used to analyze the ET0 changes in the BTH region, which may have introduced some uncertainties. In subsequent studies, multiple model data, different downscaling methods and calculation models containing more climate impact factors should be used to evaluate the magnitude of the uncertainties and to increase the accuracy of the simulation results. In addition, some parameters involved in the ET0 model and in the industrial water consumption model (e.g., the latent heat value of water vapor and the constant value of temperature included in the Hargreaves formula; the industrial water repetition rate and industrial water price) have large associated uncertainties. Further improvement of the accuracy of these parameters would also make the simulation results of water demand more accurate.

5. Conclusions

In this paper, a climate downscaling model was used to simulate future regional climate change under the RCP2.6, RCP4.5 and RCP8.5 scenarios. The output from this model was used to drive models of agricultural and industrial water demand so as to simulate and predict future water demand in the BTH region. The following conclusions were reached:
(1)
During the forecast period (2020–2035), the ET0 growth rates in the Beijing, Tianjin and Hebei areas under the RCP2.6 scenario are 1.438 mm·a−1, 1.393 mm·a−1 and 2.059 mm·a−1, respectively. Under the RCP4.5 scenario, they are 2.252 mm·a−1, 2.310 mm·a−1 and 2.827 mm·a−1, respectively. Under the RCP8.5 scenario, they are 3.123 mm·a−1, 2.310 mm·a−1 and 2.141 mm·a−1, respectively. Under each climate scenario, the increase in evapotranspiration in the Hebei area is the largest, followed by that in the Tianjin area, and that in the Beijing area is the smallest. The order of increase in evapotranspiration in all three areas under the different climate scenarios is RCP8.5 > RCP4.5 > RCP2.6.
(2)
During the forecast period (2020–2035), under the three different climate scenarios, the water consumption per CNY 10,000 of industrial added value in the Beijing area shows a downward trend, with rates of 0.158, 0.153 and 0.110, respectively. The water consumption per CNY 10,000 of industrial added value in the Tianjin area shows an upward trend, with upward tendency rates of 0.170, 0.087 and 0.071 under the three climate scenarios, respectively. The water consumption per CNY 10,000 of industrial added value in the Hebei area also shows an upward trend, with upward tendency rates of 0.254, 0.071 and 0.036 under the three climate scenarios, respectively. The results of this study provide a statistical method of assessing the impact of climate change on industrial water demand in the BTH region.
(3)
In the context of future climate change, with the increase in temperature, the agricultural demand for water will increase significantly, and the industrial demand for water will be less affected, but will also rise. Therefore, in the future allocation process of agricultural and industrial water resources in the Beijing–Tianjin–Hebei region, on the one hand, targeted measures should be proposed to deal with the negative effects of climate change based on the development status and water use characteristics of each region, and on the other hand, the water use structure of each region should be fully considered to allow the strategic cooperative relationship between the city clusters to develop. Through the rational distribution of water resources in the different regions, the constraints of water supply and demand on social and economic development can be solved, so as to promote the healthy and sustainable development of the study region.
However, there are still some areas that need to be improved. First of all, compared with the traditional Penman–Monteith formula, the temperature-based Hargreaves formula adopted in this paper is relatively simple to calculate, and the meteorological parameters considered to change with climate change only include the minimum, maximum and average temperature. For the important parameter of solar radiation involved in the formula, it was calculated through parameters such as geographical latitude, which were easy to measure, and only its specific changes over time were taken into account. Therefore, the influencing factors of this method are relatively simple. As a future improvement, more climate-influencing factors could be added based on the Penman–Monteith method according to research needs. And it will make the prediction result more realistic. Secondly, there were relatively few variables in the industrial water demand model established in this paper. Future innovations and improvements to industrial models may use methods such as neural networks and system dynamics to include as many relevant variables as possible to improve the prediction results. Finally, the Coupled Model Intercomparison Project Phase 6 (CMIP6) provides a scenario framework for global socio-economic development and climate change by combining Shared Socioeconomic Pathways (SSP) and Representative Concentration Pathways (RCP). Unlike the RCP scenario in CMIP5, CMIP6 takes into account future socio-economic development, and the GCM data are more accurate and complex than the previous CMIP5 scenario. In future studies, CMIP6 data can be combined to make a more accurate prediction of water demand under different development scenarios in the future, and relevant studies will have more value to reflect future changes and hot spots.

Author Contributions

Conceptualization, Q.Z. and W.D.; methodology, Y.Z. and M.C.; software, Y.Z. and M.C.; validation, M.C. and Y.Z.; formal analysis, Y.Z., Q.Z. and M.C.; investigation, Y.Z. and M.C.; resources, M.C.; data curation, W.D. and M.C.; writing—original draft preparation, Y.Z.; writing—review and editing, Y.Z., Q.Z. and W.D.; visualization, Q.Z. and W.D.; supervision, Q.Z. and W.D.; project administration, Q.Z. and W.D.; funding acquisition, Q.Z. and W.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Beijing Social Science Foundation under grant number 19GLB014, the National Natural Science Foundation of China under grant number 42122004, and the Science Foundation of Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences under grant number E2500102.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the data sharing policy.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BTHBeijing–Tianjin–Hebei
SDSMStatistical downscaling model
CMIP5Coupled Model Intercomparison Project Phase 5
RCPRepresentative Concentration Pathway
SPSSStatistical Product and Service Solutions
NCEPNational Centers for Environmental Prediction
IPCCIntergovernmental Panel on Climate Change

References

  1. Zhao, R.; Sun, H.; Xing, L.; Li, R.; Li, M. Effects of anthropogenic climate change on the drought characteristics in China: From frequency, duration, intensity, and affected area. J. Hydrol. 2023, 617, 129008. [Google Scholar] [CrossRef]
  2. Beans, E.W. A heat transfer model for global warming. Case Stud. Therm. Eng. 2023, 43, 102677. [Google Scholar] [CrossRef]
  3. Shao, W.; Liu, H.; Wang, H.-f.; Liu, J.; Yan, D.; Li, W.; Zhou, J.; Wang, H. Evaluation of regional water security in China and recommendations for counter measures. Arab. J. Geosci. 2020, 13, 107. [Google Scholar] [CrossRef]
  4. Xiang, M.; Zhang, S.-H.; Ruan, Q.; Tang, C.-H.; Zhao, Y. Definition and calculation of hierarchical ecological water requirement in areas with substantial human activity—A case study of the Beijing–Tianjin-Hebei region. Ecol. Indic. 2022, 138, 108740. [Google Scholar] [CrossRef]
  5. Zhao, L.; Zhao, X.; Pan, X.; Shi, Y.; Qiu, Z.; Li, X.; Xing, X.; Bai, J. Prediction of daily reference crop evapotranspiration in different Chinese climate zones: Combined application of key meteorological factors and Elman algorithm. J. Hydrol. 2022, 610, 127822. [Google Scholar] [CrossRef]
  6. Xu, C.Y.; Singh, V.P. Evaluation of three complementary relationship evapotranspiration models by water balance approach to estimate actual regional evapotranspiration in different climatic regions. J. Hydrol. 2005, 308, 105–121. [Google Scholar] [CrossRef]
  7. Wu, H.; Zhu, W.; Huang, B. Seasonal variation of evapotranspiration, Priestley-Taylor coefficient and crop coefficient in diverse landscapes. Geogr. Sustain. 2021, 2, 224–233. [Google Scholar] [CrossRef]
  8. de Moura Neto, A.; Tavares Fernandes, G.S.; Rafaelly de Almeida Lopes, J.; de Araujo Lima, E.; Marcos de Oliveira Lopes, P.; Rodrigues, L.S.; da Silva Gonçalves Juniora, A. Changes in land use and cover and their consequences on the spatial distribution of evapotranspiration and climatic elements in southwest Piauí, MATOPIBA. J. S. Am. Earth Sci. 2022, 120, 104084. [Google Scholar] [CrossRef]
  9. Ladlani, I.; Houichi, L.; Djemili, L.; Heddam, S.; Belouz, K. Estimation of Daily Reference Evapotranspiration (ET0) in the North of Algeria Using Adaptive Neuro-Fuzzy Inference System (ANFIS) and Multiple Linear Regression (MLR) Models: A Comparative Study. Arab. J. Sci. Eng. 2014, 39, 5959–5969. [Google Scholar] [CrossRef]
  10. Zitouna-Chebbi, R.; Jacob, F.; Prévot, L.; Voltz, M. Documenting evapotranspiration and surface energy fluxes over rainfed annual crops within a Mediterranean hilly agrosystem. Agric. Water Manag. 2023, 277, 108117. [Google Scholar] [CrossRef]
  11. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration: Guidelines for Computing Crop Water Requirements; FAO: Rome, Italy, 1998. [Google Scholar]
  12. Chen, H.; Sun, J. Changes in Drought Characteristics over China Using the Standardized Precipitation Evapotranspiration Index. J. Clim. 2015, 28, 5430–5447. [Google Scholar] [CrossRef]
  13. Gonçalves, I.Z.; Ruhoff, A.; Laipelt, L.; Bispo, R.C.; Hernandez, F.B.T.; Neale, C.M.U.; Teixeira, A.H.C.; Marin, F.R. Remote sensing-based evapotranspiration modeling using geeSEBAL for sugarcane irrigation management in Brazil. Agric. Water Manag. 2022, 274, 107965. [Google Scholar] [CrossRef]
  14. Feng, X.; Liu, H.; Feng, D.; Tang, X.; Li, L.; Chang, J.; Tanny, J.; Liu, R. Quantifying winter wheat evapotranspiration and crop coefficients under sprinkler irrigation using eddy covariance technology in the North China Plain. Agric. Water Manag. 2023, 277, 108131. [Google Scholar] [CrossRef]
  15. Li, C.; Wu, P.T.; Li, X.L.; Zhou, T.W.; Sun, S.K.; Wang, Y.B.; Luan, X.B.; Yu, X. Spatial and temporal evolution of climatic factors and its impacts on potential evapotranspiration in Loess Plateau of Northern Shaanxi, China. Sci. Total Environ. 2017, 589, 165–172. [Google Scholar] [CrossRef] [PubMed]
  16. Yang, C.; Lei, H. Climate and management impacts on crop growth and evapotranspiration in the North China Plain based on long-term eddy covariance observation. Agric. For. Meteorol. 2022, 325, 109147. [Google Scholar] [CrossRef]
  17. Han, J.; Wang, J.-h.; Zhao, Y.; Wang, Q.; Zhang, B.-C.; Li, H.; Zhai, J. Spatio-temporal variation of potential evapotranspiration and climatic drivers in the Jing-Jin-Ji region, North China. Agric. For. Meteorol. 2018, 256–257, 75–83. [Google Scholar] [CrossRef]
  18. Bakhshoodeh, R.; Ocampo, C.; Oldham, C. Evapotranspiration rates and evapotranspirative cooling of green façades under different irrigation scenarios. Energy Build. 2022, 270, 112223. [Google Scholar] [CrossRef]
  19. Zhang, K.-X.; Pan, S.-M.; Zhang, W.; Xu, Y.-H.; Cao, L.-G.; Hao, Y.-P.; Wang, Y. Influence of climate change on reference evapotranspiration and aridity index and their temporal-spatial variations in the Yellow River Basin, China, from 1961 to 2012. Quat. Int. 2015, 380–381, 75–82. [Google Scholar] [CrossRef]
  20. Shan, N.; Shi, Z.; Yang, X.; Gao, J.; Cai, D. Spatiotemporal trends of reference evapotranspiration and its driving factors in the Beijing–Tianjin Sand Source Control Project Region, China. Agric. For. Meteorol. 2015, 200, 322–333. [Google Scholar] [CrossRef]
  21. Wang, X.; Zhang, J.y.; Shahid, S.; Bi, S.h.; Elmahdi, A.; Liao, C.-h.; Li, Y. Forecasting industrial water demand in Huaihe River Basin due to environmental changes. Mitig. Adapt. Strat. Gl 2018, 23, 469–483. [Google Scholar] [CrossRef]
  22. Yan, B.; Xu, Y. Evaluation and prediction of water resources carrying capacity in Jiangsu Province, China. Water Policy 2022, 24, 324–344. [Google Scholar] [CrossRef]
  23. Gulfam-E-Jannat, S.; Golui, D.; Islam, S.; Saha, B.; Rahman, S.M.; Bezbaruah, A.N.; Iskander, S.M. Industrial Water Demand and Wastewater Generation: Challenges for Bangladesh’s Water Industry. ACS ES&T Water 2023, 3, 1515–1526. [Google Scholar]
  24. Li, K.; Ma, T.; Wei, G.; Zhang, Y.; Feng, X. Urban Industrial Water Supply and Demand: System Dynamic Model and Simulation Based on Cobb–Douglas Function. Sustainability 2019, 11, 5893. [Google Scholar] [CrossRef]
  25. Gumbo, B.; Mlilo, S.; Broome, J.H.; Lumbroso, D.M. Industrial water demand management and cleaner production potential: A case of three industries in Bulawayo, Zimbabwe. Phys. Chem. Earth 2003, 28, 797–804. [Google Scholar] [CrossRef]
  26. Vásquez-Lavín, F.; Vargas, O.L.; Hernández, J.I.; Ponce Oliva, R.D. Water demand in the Chilean manufacturing industry: Analysis of the economic value of water and demand elasticities. Water Resour. Econ. 2020, 32, 100159. [Google Scholar] [CrossRef]
  27. Gracia-de-Rentería, P.; Barberán, R. Economic Determinants of Industrial Water Demand: A Review of the Applied Research Literature. Water 2021, 13, 1684. [Google Scholar] [CrossRef]
  28. Lyu, F.; Zhang, H.; Dang, C.; Gong, X. A novel framework for water accounting and auditing for efficient management of industrial water use. J. Clean. Prod. 2023, 395, 136458. [Google Scholar] [CrossRef]
  29. Ministry of Water Resources of the People’s Republic of China. China Water Resources Bulletin 2021. (In Chinese). Available online: http://www.mwr.gov.cn/sj/tjgb/szygb/202206/t20220615_1579315.html (accessed on 21 April 2023).
  30. Wilby, R.L.; Dawson, C.W.; Barrow, E.M. SDSM—A decision support tool for the assessment of regional climate change impacts. Environ. Model. Softw. 2002, 17, 145–157. [Google Scholar] [CrossRef]
  31. Wilby, R.L.; Whitehead, P.G.; Wade, A.J.; Butterfield, D.; Davis, R.J.; Watts, G. Integrated modelling of climate change impacts on water resources and quality in a lowland catchment: River Kennet, UK. J. Hydrol. 2006, 330, 204–220. [Google Scholar] [CrossRef]
  32. Her, Y.; Yoo, S.H.; Cho, J.-P.; Hwang, S.; Jeong, J.; Seong, C. Uncertainty in hydrological analysis of climate change: Multi-parameter vs. multi-GCM ensemble predictions. Sci. Rep. 2019, 9, 4974. [Google Scholar] [CrossRef]
  33. Zhou, T.; Wu, P.; Sun, S.; Li, X.; Wang, Y.; Luan, X. Impact of Future Climate Change on Regional Crop Water Requirement—A Case Study of Hetao Irrigation District, China. Water 2017, 9, 429. [Google Scholar] [CrossRef]
  34. Eingrüber, N.; Korres, W. Climate change simulation and trend analysis of extreme precipitation and floods in the mesoscale Rur catchment in western Germany until 2099 using Statistical Downscaling Model (SDSM) and the Soil & Water Assessment Tool (SWAT model). Sci. Total. Environ 2022, 838, 155775. [Google Scholar] [PubMed]
  35. Mahmood, R.; Babel, M.S. Future changes in extreme temperature events using the statistical downscaling model (SDSM) in the trans-boundary region of the Jhelum river basin. Weather Clim. Extremes 2014, 5, 56–66. [Google Scholar] [CrossRef]
  36. Jia, K.; Yang, Y.; Dong, G.; Zhang, C.; Lang, T. Variation and determining factor of winter wheat water requirements under climate change. Agric. Water Manag. 2021, 254, 106967. [Google Scholar] [CrossRef]
  37. Ding, R.; Kang, S.; Zhang, Y.; Hao, X.; Tong, L.; Du, T. Partitioning evapotranspiration into soil evaporation and transpiration using a modified dual crop coefficient model in irrigated maize field with ground-mulching. Agric. Water Manag. 2013, 127, 85–96. [Google Scholar] [CrossRef]
  38. Gavilán, P.; Lorite, I.J.; Tornero, S.; Berengena, J. Regional calibration of Hargreaves equation for estimating reference ET in a semiarid environment. Agric. Water Manag. 2006, 81, 257–281. [Google Scholar] [CrossRef]
  39. Feng, Y.; Yue, J.; Cui, N.; Zhao, L.; Chen, L.; Gong, D. Calibration of Hargreaves model for reference evapotranspiration estimation in Sichuan basin of southwest China. Agric. Water Manag. 2017, 181, 1–9. [Google Scholar] [CrossRef]
  40. Gentilucci, M.; Bufalini, M.; Materazzi, M.; Barbieri, M.; Aringoli, D.; Farabollini, P.; Pambianchi, G. Calculation of Potential Evapotranspiration and Calibration of the Hargreaves Equation Using Geostatistical Methods over the Last 10 Years in Central Italy. Geosciences 2021, 11, 348. [Google Scholar] [CrossRef]
  41. George, H.H.; Zohrab, A.S. Reference Crop Evapotranspiration from Temperature. Appl. Eng. Agric. 1985, 1, 96–99. [Google Scholar] [CrossRef]
  42. Li, K.; Hajar, S.; Ding, Z.; Dooling, T.; Wei, G.; Hu, C.; Zhang, Y.; Zhang, K. Dynamic optimization of input production factors for urban industrial water supply and demand management. J. Environ. Manag. 2020, 270, 110807. [Google Scholar] [CrossRef]
  43. Chen, Y.; Jin, B.; Jin, J.; Wang, G.; Cao, M. System dynamics model-based water demand prediction under changing environment with consideration of physical mechanism. Hydro-Sci. Eng. 2021, 15, 84–95. (In Chinese) [Google Scholar]
Figure 2. Predicted values of daily maximum and minimum temperatures under three climate scenarios in the BTH region. Left figures—daily maximum temperature. Right figures—daily minimum temperature.
Figure 2. Predicted values of daily maximum and minimum temperatures under three climate scenarios in the BTH region. Left figures—daily maximum temperature. Right figures—daily minimum temperature.
Water 15 04225 g002
Figure 3. Interannual evapotranspiration simulation results in the BTH region. The evapotranspiration fluctuations and trend lines of the base period (2002 to 2019) and forecast period (2020 to 2035) are also shown.
Figure 3. Interannual evapotranspiration simulation results in the BTH region. The evapotranspiration fluctuations and trend lines of the base period (2002 to 2019) and forecast period (2020 to 2035) are also shown.
Water 15 04225 g003
Figure 4. Water consumption per CNY 10,000 of industrial added value in BTH region. The fluctuation and trend lines for water consumption per CNY 10,000 of industrial added value in the base period (2005 to 2019) and forecast period (2020 to 2035) are also shown.
Figure 4. Water consumption per CNY 10,000 of industrial added value in BTH region. The fluctuation and trend lines for water consumption per CNY 10,000 of industrial added value in the base period (2005 to 2019) and forecast period (2020 to 2035) are also shown.
Water 15 04225 g004
Table 1. The predictors corresponding to the daily maximum temperature (Tmax), daily minimum temperature (Tmin) and daily precipitation (P) at 20 meteorological stations. Definitions of the predictor variables are given in the notes below the table.
Table 1. The predictors corresponding to the daily maximum temperature (Tmax), daily minimum temperature (Tmin) and daily precipitation (P) at 20 meteorological stations. Definitions of the predictor variables are given in the notes below the table.
Meteorological StationTmaxTminP
Beijing, Miyun, Baodingplth p1zh p5_z p500 p5zh tempp1zh p5_z p500 shum tempp1_f p1zh p500 p8_f prcp s850
Tianjin, Binhai, Tanggup1_u p1zh p5zh p500p8zh shum tempp1zh p5_z p500 shum tempp5_f p8_f p8_v prcp s500 s850
Chengde, Fengning, Weichangp1_v p1zh p5_z p5zh p8zh tempp1zh p5_z prcp s500 shum tempp1_f p1zh p500 p8_f prcp s500
Tangshan, Zunhua, Letingp1_u p1zh p5_z p500 p5zh tempp1zh p5_z p500 s500 shum tempp1_f p1zh p8_v prcp s500 s850
Qinglongp1_u p5zh p8th p8zh tempp1zh p5_z prcp s500 shum tempp1_f p1zh p8_v prcp s500 s850
Shijiazhuang, Nangong, Xingtai, Cangzhoup1_u p1th p500 p5zh p8zh shum tempp1th p1zh p5_z p500 shum tempp1zh p5_v p8_v prcp s500 s850
Huailai, Weixian, Zhangjiakoup1zh p5_z p500 p5zh p8_u p8zh tempp1zh p5_z prcp s500 shum tempp1_f p1zh p500 prcp s500 s850
Note: predictors are: p1_f—sea level wind speed; p1_u—sea level meridional temperature; p1_v—sea level zonal temperature; p1th—sea level wind direction; p1zh—sea level divergence; p5_f—500 hPa wind speed; p5_v—500 hPa zonal temperature; p5_z—500 hPa vorticity; p500—500 hPa geopotential height; p5zh—divergence of 500 hPa; p8_f—wind speed of 850 hPa; p8_u—meridional temperature of 850 hPa; p8_v—zonal temperature of 850 hPa; p8th—850 hPa wind direction; p8zh—850 hPa divergence; prcp—precipitation; s500—500 hPa relative humidity; s850—850 hPa relative humidity; shum—near surface humidity; temp—near surface temperature.
Table 2. Details of the three future climate scenarios (RCPs) used in this study. The paths of greenhouse gas (GHG) emissions differ in the three climate change scenarios.
Table 2. Details of the three future climate scenarios (RCPs) used in this study. The paths of greenhouse gas (GHG) emissions differ in the three climate change scenarios.
Climate
Scenario
Assumed Development Path
RCP2.6The lowest GHG emissions scenario. Emissions initially rise and then fall until reaching stability.
RCP4.5Intermediate GHG emissions scenario with moderate emissions. This scenario is consistent with the requirements of future Chinese economic development planning.
RCP8.5The highest GHG emissions scenario. Emissions continue to rise.
Table 3. Results of the SDSM performance evaluation at each of the 20 meteorological stations for the period 1991 to 2005. Values presented are the average percentage of explained variance (R2), and the standard error (SE) of daily predicted values of maximum temperature (Tmax), daily minimum temperature (Tmin) and daily precipitation (P).
Table 3. Results of the SDSM performance evaluation at each of the 20 meteorological stations for the period 1991 to 2005. Values presented are the average percentage of explained variance (R2), and the standard error (SE) of daily predicted values of maximum temperature (Tmax), daily minimum temperature (Tmin) and daily precipitation (P).
StationTmaxTminP
R2SE (°C)R2SE (°C)R2SE (°C)
Beijing0.5592.6400.7021.6700.4450.260
Miyun0.6102.4300.6991.6200.4100.290
Tianjin0.6222.4000.7091.5700.3650.260
Binhai0.5922.6000.5172.2800.3320.290
Tanggu0.5792.5000.6231.9100.3620.240
Baoding0.5932.6000.7071.5800.4830.320
Cangzhou0.6272.4800.7321.6000.5330.170
Chengde0.6952.3000.731.7800.5360.200
Fengning0.662.6000.6832.1400.4800.220
Huailai0.692.4000.7131.7900.5620.150
Leting0.642.2700.7111.8000.6220.210
Qinglong0.6552.3700.7361.8600.6270.220
Nangong0.6682.4300.6991.6100.4290.150
Shijiazhuang0.6252.6100.6881.6700.4110.230
Tangshan0.6532.2600.7411.6800.4090.220
Weichang0.7022.4400.7241.8600.2240.310
Weixian0.4563.4100.6512.4300.4150.140
Xingtai0.6612.5300.6911.6500.3290.180
Zhangjiakou0.7052.4200.7591.6800.6140.140
Zunhua0.6232.4700.7021.9600.3350.320
Table 4. The average daily maximum and minimum temperatures over the period 2020 to 2050 in the BTH under the three scenarios.
Table 4. The average daily maximum and minimum temperatures over the period 2020 to 2050 in the BTH under the three scenarios.
Climate ScenarioAverage Daily Maximum Temperature (°C)Average Daily Minimum Temperature (°C)
BeijingTianjinHebeiBeijingTianjinHebei
RCP2.618.9219.0118.786.949.356.81
RCP4.518.2918.7917.956.599.006.43
RCP8.519.1019.3119.057.339.697.20
Table 5. The mean ET0 of the BTH region in the base period (2005–2019) and forecast period (2020–2035), and the overall growth rate under each scenario in the forecast period.
Table 5. The mean ET0 of the BTH region in the base period (2005–2019) and forecast period (2020–2035), and the overall growth rate under each scenario in the forecast period.
Average ET0 (mm,
2005–2019)
Average ET0 (mm, 2020–2035)Growth Rate (mm·a−1, 2020–2035)
RCP2.6RCP4.5RCP8.5RCP2.6RCP4.5RCP8.5
Beijing1203.7471237.4471242.9861240.0961.4382.2523.123
Tianjin1151.6581189.9361206.2021213.1391.3932.3102.310
Hebei1095.8221224.1691224.1691240.6212.0592.8272.141
Table 6. The fluctuation range and growth rate of the water consumption per CNY 10,000 of industrial added value in the BTH region under three climate scenarios during the forecast period (2020–2035).
Table 6. The fluctuation range and growth rate of the water consumption per CNY 10,000 of industrial added value in the BTH region under three climate scenarios during the forecast period (2020–2035).
Fluctuation Range (m3)Growth Rate
RCP2.6RCP4.5RCP8.5
Beijing3.27–6.63−0.158−0.153−0.110
Tianjin13.82–28.280.1700.0870.071
Hebei12.45–17.890.2540.0710.036
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhou, Q.; Zhong, Y.; Chen, M.; Duan, W. Climate Change Impacts on Agricultural and Industrial Water Demands in the Beijing–Tianjin–Hebei Region Using Statistical Downscaling Model (SDSM). Water 2023, 15, 4225. https://doi.org/10.3390/w15244225

AMA Style

Zhou Q, Zhong Y, Chen M, Duan W. Climate Change Impacts on Agricultural and Industrial Water Demands in the Beijing–Tianjin–Hebei Region Using Statistical Downscaling Model (SDSM). Water. 2023; 15(24):4225. https://doi.org/10.3390/w15244225

Chicago/Turabian Style

Zhou, Qian, Yating Zhong, Meijing Chen, and Weili Duan. 2023. "Climate Change Impacts on Agricultural and Industrial Water Demands in the Beijing–Tianjin–Hebei Region Using Statistical Downscaling Model (SDSM)" Water 15, no. 24: 4225. https://doi.org/10.3390/w15244225

APA Style

Zhou, Q., Zhong, Y., Chen, M., & Duan, W. (2023). Climate Change Impacts on Agricultural and Industrial Water Demands in the Beijing–Tianjin–Hebei Region Using Statistical Downscaling Model (SDSM). Water, 15(24), 4225. https://doi.org/10.3390/w15244225

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop