Evapotranspiration is one of the hot issues of ecological hydrology. However, few studies have analysed the impact of vegetation changes on water consumption from the perspective of natural and artificial vegetation. Taking the Ziya River Basin as an example, the daily meteorological data from 2001 to 2015 were used to calculate the water consumption of vegetation based on the Penman–Monteith model. The results showed that the vegetation coverage increased. The total water consumption increased from 2001 (2.60 × 1,010 m3) to 2005 (2.65 × 1,010 m3) and decreased from 2005 (2.65 × 1,010 m3) to 2015 (2.40 × 1,010 m3). The water consumption per unit area in descending order was mixed forests (660 mm, annual average), croplands (640 mm), closed shrublands (581 mm), deciduous broadleaf forests (528 mm), grasslands (514 mm), savannas (459 mm), and woody savannas (454 mm). Finally, the regression equation between vegetation coverage change and water consumption was y = 0.377x + 84.516, which showed that there was a proportional relationship. Therefore, attention should be given to balancing local water allocation during vegetation restoration. The results can provide a reference for vegetation restoration policies

  • Few studies have analysed the impact of vegetation changes on water consumption to scientifically propose that there is a moderate threshold in vegetation restoration.

  • The quantitative relationship between vegetation change and water consumption was obtained by regression analysis. And we proposed that vegetation restoration has a moderate threshold.

Water consumption is the water that cannot be directly recycled in the basin. Vegetation water consumption in the basin is the most important form, and reducing water consumption in the basin is the most effective measure to achieve real water saving in the basin. In recent years, many watersheds have carried out afforestation to improve the environment, but it has also directly increased the water consumption of the basin, especially in water-deficient areas, often further exacerbating the tension of water resources. Water is the most important limiting factor for the growth and development of vegetation. What's more, water stress caused crop yield reduction, which can exceed the sum of other environmental stresses. In order to reduce the water pressure of plants, it is necessary to clarify the water consumption rules of different plants, identify the factors affecting the water consumption of plants, and establish the quantitative relationship between the water consumption of plants and the change rate of plants, so as to formulate a vegetation restoration plan with low water consumption and suitable for regional water resources, avoiding the phenomenon of planting high water-consuming plants in water-deficient areas, which not only reduces the pressure of regional water resources but also reduces the pressure of water shortage on plants. Therefore, in order to formulate a plant restoration policy suitable for water resources in the region, it is of great significance to study the impact of vegetation change on water consumption for regional sustainable development (Pereira et al., 2021).

Plant transpiration and soil evaporation, which are also known as evapotranspiration (ET), account for more than 90% of water consumption and are the main body of water consumption. The calculation of vegetation water consumption is the primary basis for formulating policies for the rational utilization of water resources (Lopez et al., 2019). Vegetation water consumption varies greatly in different phases and at different spatial scales in a watershed, which is affected by many factors. It was found that vegetation coverage and growth period were the most critical factors affecting vegetation water consumption in arid Arab. Wind speed was the main driving force in a 50-year (1961–2010) vegetation water consumption study in the Iranian region. The change of vegetation distribution in the Yanhe River Basin may be the main factor affecting the change of vegetation water consumption. Valipour and his team's monitoring of 18 sites found the complexity of climate change driving vegetation water consumption changes (Valipour et al., 2020). On the Loess Plateau from 2001 to 2020, the most important factors affecting vegetation water consumption were the normalized difference vegetation index (NDVI) and wind speed (Zhang et al., 2022). In the Pearl River Basin from 1980 to 2019, the main factor affecting vegetation water consumption was the average temperature (Zuo et al., 2022). The increase in vegetation water consumption in humid areas of China is caused by 91% of the climate and 9% by vegetation greening (Zhang et al., 2020).

Studies have shown that changes in vegetation (such as species, quantity, and distribution) affect changes in vegetation water consumption. The greater the vegetation density in arid Arab, the greater the water consumption (El-Shirbeny et al., 2022). The natural vegetation area increased at a rate of 12.37%/5a, and the water consumption increased at a rate of 12.82%/5a in the Tarim River Basin (Wang et al., 2022). In the lower reaches of the Yanhe River Basin located in the semi-arid area, the water consumption per unit area from high to low is woodland > shrubland > spare forest land > medium-coverage grassland > cultivated land > low-coverage grassland (Guo et al., 2022).

Accounting for vegetation water consumption is the core of water resource assessment. Therefore, the accurate calculation of vegetation water consumption needs to select the appropriate calculation method according to the actual situation of the study area. Although there are many calculation methods for vegetation water consumption, the FAO-56 Penman–Monteith (PM) equation is the most widely used. The FAO-56 PM method has accuracy, stability, and global acceptability (Pereira et al., 2021). The FAO-56 PM method has been widely used in Asia, America, and other regions (Wang, 2014; Sullivan et al., 2019; Zhang et al., 2021).

The above studies showed that there is a close relationship between vegetation change and water consumption, and there are regional factors affecting vegetation and water consumption in different regions. The PM model is an effective tool for studying regional evapotranspiration. However, few studies have comprehensively analysed the relationship between vegetation change and water consumption, and proposed that there is a moderate threshold in regional vegetation restoration; that is, vegetation restoration must be based on the regional water resource-carrying capacity to guide the formulation and implementation of water policy. This study focuses on the Ziya River Basin (ZYRB), which is primarily composed of croplands, grasslands, and so on. As one of the most important agricultural production areas in China, the ecological environment of the ZYRB is fragile and sensitive to climate change (Yu, 2018). Consequently, the calculation of water consumption is important to socioeconomic development and ecological environmental improvement in the region. In this paper, the land use of different phases is used to analyse interannual changes, and the NDVI is used to analyse annual changes. The objectives of this study are as follows: (1) to analyse vegetation changes and causes in different phases of the ZYRB; (2) to analyse the interannual variability of water consumption associated with dominant land-cover classes across the ZYRB; (3) to evaluate water consumption patterns in terms of crop type and water management strategies; and (4) to explore the influence of vegetation changes on water consumption according to the spatial and temporal distribution characteristics of water consumption in different phases of the basin and to propose that vegetation restoration has a moderate threshold.

2.1. Study area

The Ziya River Basin (112°20′–117°50′,36°03′–39°35′), which is located in the southern part of the Haihe River Basin in North China, is one of the seven major subbasins of the Haihe River Basin. Based on the administrative division, the basin is located in Hebei, Shanxi, and Tianjin (Xiong & Li, 2007), covering an area of 4.69 × 104 km2. The ZYRB has a semi-arid and semi-humid continental monsoon climate in the warm temperate zone (Li et al., 2021). Based on meteorological data from 2008 to 2012, the annual average temperature is 11.8–12.9 °C, and the annual average precipitation is 540 mm. Most precipitation is concentrated in June, July, and August. The precipitation in mountain areas is approximately 620 mm, while it is approximately 500 mm in plain areas. The potential evapotranspiration in the upstream region is approximately 700–1,200 mm, and the actual evapotranspiration is 300–500 mm (Li, 2012).

The vegetation in the ZYRB can be divided into eight categories: grasslands, closed shrublands, woody savannas, mixed forests, deciduous broadleaf forests, croplands, cropland/natural vegetation mosaics, and savannas. The total terrain of the basin is high in the northwest and low in the southeast. It is mainly located in the northwestern part of the basin in Shanxi Province, and this area is mountainous terrain. The eastern and central parts of the basin contain plains and are located in Hebei Province. Mountain areas account for approximately 67% of the total basin area, and plain areas account for approximately 33% (Figure 1). Crops in the ZYRB include cotton, corn, soybeans, peanuts, sorghum, sweet potato, etc., and crop growth characteristics are 1 year.
Fig. 1

Weather stations and topographic map of the ZYRB.

Fig. 1

Weather stations and topographic map of the ZYRB.

Close modal

Data

The data include geographic information data, remote sensing data, meteorological data, and hydrological data (Table 1). According to the availability of data, we choose 2000–2015 as the research period.

  • (1)

    Monthly NDVI data in 2001, 2005, 2010, and 2014 were obtained from the data cloud website. NDVI data include MODND1M China NDVI products calculated by MODND1D within 1 month and at a maximum daily temporal resolution, and spatial resolution is monthly and 500 m. Due to the incompleteness of the data in 2015, the NDVI data in 2014 were used as a replacement (Jiao, 2021).

  • (2)

    According to the FAO-56 Penman–Monteith formula, meteorological data, including daily average pressure, daily maximum temperature, daily minimum temperature, daily precipitation, daily evaporation, relative humidity, daily average wind speed, and sunshine hours, are needed.

Table 1

Data obtained in this article.

DataSourceYearSpatial resolution
The remote sensing data of vegetation and coverage https://lpdaac.usgs.gov/products(MCD12Q12001, 2005, 2010, 2015 463 m × 463 m 
DEM digital elevation www.gscloud.cn 30 m 
NDVI https://ladsweb.nascom.nasa.gov 2001, 2005, 2010, 2014 (monthly) 500 m 
Meteorological data http://data.cma.cn 2001, 2005, 2010, 2015 (daily) 
DataSourceYearSpatial resolution
The remote sensing data of vegetation and coverage https://lpdaac.usgs.gov/products(MCD12Q12001, 2005, 2010, 2015 463 m × 463 m 
DEM digital elevation www.gscloud.cn 30 m 
NDVI https://ladsweb.nascom.nasa.gov 2001, 2005, 2010, 2014 (monthly) 500 m 
Meteorological data http://data.cma.cn 2001, 2005, 2010, 2015 (daily) 

Penman–Monteith model

In theory, the crop water requirement is defined as the evapotranspiration of crops that grow in a large area without insect diseases and meet soil evaporation and plant transpiration under sufficient soil moisture and other conditions. The FAO-56 Penman–Monteith equation (hereafter referred to as the ‘P–M equation’) recommended by the Food and Agriculture Organization of the United Nations (hereafter referred to as ‘FAO’) is a relatively mature and widely used empirical formula method. The specific calculation formulas are as follows (Allen et al., 1998):
formula
(1)
where ET0 is the reference evapotranspiration, mm/d; Rn is the net radiation of the crop canopy surface, MJ/(m2·d); G is the soil heat flux, which is 0 for each day or every 10 days, MJ/(m2·d); T is the average daily air temperature at a 2 m height, °C; is the wind speed at a 2 m height, m/s; is the air saturated vapour pressure, kPa; is the actual vapour pressure of air, kPa; is the saturated vapour pressure difference, kPa; is the slope of the water vapour pressure curve, kPa/°C; and is the hygrometer constant, kPa/°C.
When the empirical formula method is used to calculate the crop water demand, the reference crop evapotranspiration needs to be multiplied by the crop coefficient. The crop coefficient, an important parameter for estimating crop water demand, reflects the crop species characteristics, cultivation, and soil conditions (Kang & Xiong, 1991; Zhang, 2006). Under suitable soil moisture conditions, the crop coefficient kc is the ratio of actual crop water demand to reference crop evapotranspiration, reflecting the influence of crop characteristics on water consumption. The calculation formula is as follows (Allen et al., 1998):
formula
(2)
where Kc is the crop coefficient (Liu, 2014); is the actual evapotranspiration, mm/d; and is the evapotranspiration, mm/d (Sheng, 2017).

Land-use change index model

The land-use change index can be expressed by the land-use dynamic degree.

  • (1)

    Single land-use type dynamic degree

The dynamic degree of a single land-use type is the quantitative change of a certain land-use type in a certain time range in a study area. The expression is:
formula
(3)
where K is the dynamic degree of a certain land-use type during the study period, Ua and Ub are the numbers of a certain land-use type at the beginning and end of the study period, respectively, and T is the study period. When the period of T is set to 1 year, the value of K is the annual changing rate of a certain land-use type in the study area.
  • (2)

    Comprehensive land-use dynamic degree

The dynamic degree of comprehensive land use in a research sample area can be expressed as follows:
formula
(4)
where LUi is the area of type i land use at the beginning of the monitoring period, is the absolute value of the area of type i land use converted to non-type i land use during the monitoring period, and T is the length of the monitoring period. When the period of T is set to be a year, the value of LC is the annual changing rate of land use in the study area.

The land-use dynamic degree quantitatively describes the speed of land-use change, which has a positive effect on predicting the future trend of land-use change (Wang et al., 1999).

Maximum value composite

The monthly NDVI value and the annual NDVI value of the study area were obtained by the maximum synthesis method to reflect the overall changing characteristics of vegetation in space. The maximum synthesis method can effectively reduce the influence of aerosols, cloud shadows, solar elevation angles, and other factors in the atmosphere.
formula
(5)
where MNDVIi is the maximum NDVI in month i; i = 1,2…, 12 is the monthly number; NDVI1 is the NDVI value in the first half of month i; and NDVI2 is the NDVI value in the second half of month i. Based on the monthly NDVI value, the annual NDVI was calculated.

Multiple linear regression model

Statistics are used to measure the fitting degree of the sample regression to the sample observations.
formula
(6)
where k is the number of explanatory variables, βj (j = 1, 2…, k) is the regression coefficient, and μ is the random error after removing the influence of k independent variables on Y.
Measuring the fitting degree of the sample regression to sample observations by using statistics:
formula
(7)
where R2 is the coefficient of determination. The closer this statistic is to 1, the higher the goodness of fit of the model.
The indigenous test of the overall linearity of the equation (F test) under the condition of the original hypothesis is H0.
formula
(8)

Following the F distribution of degrees of freedom (k, nk − 1), given the explicitness level α (generally 0.05 or 0.1, i.e., 95% confidence or 90% confidence, respectively), the critical value (k, nk − 1) is obtained by looking up the table. According to the sample, the value of the F statistic is calculated, and F > (k, nk − 1) is used to reject or accept the original hypothesis H0 to determine whether the linear relationship of the original equation is obvious in general.

Figure 2 shows the research ideas of this paper and the use of the above research methods in the paper.
Fig. 2

Technology roadmap of the paper.

Fig. 2

Technology roadmap of the paper.

Close modal

Distribution of land-use types

Through the analysis of the ZYRB land-use type distribution and area change (Figure 3), land-use types were mainly grasslands, croplands, and forests, accounting for more than 80% of the land-use area. Among them, grasslands accounted for approximately 30% and were scattered in the northwestern part of the ZYRB; croplands accounted for approximately 52% and were mainly distributed in the eastern and northwestern parts of the ZYRB; and forests were mainly scattered in the mountainous areas of the ZYRB, accounting for approximately 1% of the area. Water bodies and permanent wetlands accounted for approximately 0.1% of the total area of the basin. Wasteland, or barren land, accounted for 0.02% of the total area, and urban and built-up lands accounted for approximately 7–8%, were mainly distributed in the central and southern parts of the basin, and were mostly located in plain areas. In terms of change, the grassland, closed shrublands, and cropland areas decreased, and urban and construction lands, mixed forests, deciduous broadleaf forests, and woody savanna areas increased in land-use types. Other types, such as wastelands, cropland/natural vegetation mosaics, water bodies, and permanent wetlands, were changed. However, the land-use structure did not change in the basin and was still dominated by grassland, arable land, and woodland.
Fig. 3

Spatial distribution map of land-use types in the ZYRB from 2001 to 2015.

Fig. 3

Spatial distribution map of land-use types in the ZYRB from 2001 to 2015.

Close modal

Temporal and spatial evolution of land use

The main land-use types in the plain area were croplands, followed by some urban and built-up lands. Mountain areas had mainly grasslands and forests. The dynamic degree of land-use types in the ZYRB changed spatially from 2001 to 2015. Table 2 shows that the dynamic degree of different land-use types in the ZYRB was obviously different. The dynamic degree of cropland/natural vegetation mosaics was large, up to nearly 100%, increasing by 0.22 km2 from 2001 to 2015, and the area was doubled. Second, mixed forests, permanent wetlands, and deciduous broadleaf forests grew faster, with growth rates of 66.83, 66.64, and 64.27%, respectively, and their areas increased by 112.69, 13.34, and 247.75 km2, respectively. The growth rate of savannas was 46.64%, with a total increase of 920.94 km2. The total area of barren land, cropland/natural vegetation mosaics, and water bodies changed little to less than 3 km2, and the change in barren land showed that the development and utilization of barren land increased. The grassland area decreased the most, showing a continuous decreasing trend. The dynamic degrees in the three periods were −3.06, −4.26, and −3.41%, respectively. The dynamic degree in the whole period was −10.34%, and the total area decreased by 1,595.60 km2. With the development of the economy and urbanization, human demand for land increased gradually, and the area of urban and built-up lands increased continuously in the three periods, with an increase of 438.78 km2 in 2001–2015. The area of permanent wetlands continuously increased in the whole period from 2001 to 2015, with growth rates of 2.14, 41.05, and 15.67%, respectively. The ecological environment improved slightly with the increase in wetlands. Forests and savannas increased by 360.44 and 920.94 km2, respectively, with an average annual growth rate of approximately 66%. The area of croplands increased by 617.6 km2 from 2001 to 2010 and decreased by 698.53 km2 from 2010 to 2015. The area of deciduous broadleaf forests increased rapidly from 2001 to 2010, the area of mixed forests increased rapidly after 2005, and the area of savannas began to increase in large areas after 2010. This shows that in the early stage of returning croplands to forests and grasslands, the policy had just begun to advance, and the vegetation recovery was slow. After 2010, more remarkable results have been achieved. The dynamic degree of the water area was the smallest, and the area changed little (Zhang, 2005; Zhao et al., 2017; Li & Liu, 2022).

Table 2

Land-use changing range and dynamics in the ZYRB from 2001 to 2015.

 
 

Note: Red indicates a conspicuous decrease, and yellow indicates a conspicuous increase.

From the perspective of topography, the land-use area change in the plain area in the ZYRB was smaller than that in the mountainous area from 2001 to 2015. The main land-use types in the plain area were croplands, followed by urban and built-up lands. From 2001 to 2015, the urban and construction land area expanded from 3,603 to 3,942 km2, an increase of approximately 439 km2. The cropland area increased to 24,837 km2 in 2010 and decreased to 24,138 km2 in 2015. Mountain areas had mainly grasslands and forests. From 2001 to 2015, grassland areas decreased by 1,595.6 km2, mixed forests increased by 112.8 km2, and deciduous broadleaf forests increased by 247.75 km2.

To analyse the transformation of areas between land-use types, the land-transfer matrix of four periods, 2000–2005, 2005–2010, 2010–2015, and 2001–2015, was obtained by using ArcGIS technology. The area of grasslands, woody savannas, and croplands increased from 2001 to 2005, and the area of other land-use types remained basically the same. In terms of land inflow, the increase in the grassland area was mainly due to the conversion of croplands (365.88 km2) and savannas (232.08 km2). The increase in urban and built-up lands was mainly due to the conversion of agricultural land (82.46 km2). In terms of land outflow, grasslands were converted mainly to croplands (754.85 km2) and savannas (261.33 km2). Croplands were partially converted into grasslands (365.88 km2). Closed shrublands were partially transformed into the deciduous broadleaf forest (50.76 km2) and savannas (50.11 km2). From 2005 to 2010, grassland, urban and construction land, and farmland areas increased; savanna areas decreased; and the remaining land-type areas did not change much. From the inflow of land, the increase in the grassland area was transformed from croplands (301.76 km2). In terms of land outflow, grasslands were mainly converted into croplands (691.46 km2), savannas (432.95 km2), and closed shrublands (116.79 km2). From 2010 to 2015, the increase in the grassland area in terms of land inflow was mainly due to the conversion of croplands (676.10 km2) and savannas (99.16 km2). The increase in urban and built-up lands was mainly due to the conversion of croplands (213.58 km2). From the conversion out of land, grasslands were mainly converted to closed shrublands (83.66 km2), croplands (217.41 km2), and savannas (930.62 km2). Closed shrublands were mainly converted to deciduous broadleaf forests (66.03 km2) and savannas (200.01 km2). From 2001 to 2015, the types of land with greater internal conversions were grasslands, agricultural lands, dense scrubs, and savannas. The conversion types of grasslands were croplands (930.33 km2) and savannas (456.62 km2), and the conversion types of closed shrublands (159.59 km2) were croplands (1274.61 km2) and savannas (1487.25 km2). The decrease in the area of croplands was mainly due to the conversion into grasslands (930.33 km2) and urban and building lands (421.13 km2). It is inferred that the increase in urban and built-up lands was mainly due to urbanization and socioeconomic development, which made the original development space insufficient to meet human activities and thus occupied arable land, according to the land-use planning of the People's Government of Shanxi Province and Hebei Province (Shanxi, 2010; Hebei, 2015). The increase in the area of forest and grassland indicated that the policy of restoring agricultural land to forests and grasslands had been effective, with a significant increase in vegetation cover.

Analysis of the NDVI

Table 3 shows the annual and interannual NDVI values of the ZYRB in 2001, 2005, 2010, and 2014. The MODIS NDVI time series contains a large amount of information, which can reflect the changes in vegetation over time. In this paper, the traditional maximum value composite (MVC) method was used to smooth and reduce noise (Li et al., 2019).

Table 3

NDVI values for 2001–2014.

Month2001200520102014
0.20 0.23 0.21 0.34 
0.22 0.21 0.20 0.26 
0.31 0.27 0.23 0.38 
0.40 0.39 0.36 0.58 
0.51 0.55 0.57 0.62 
0.45 0.54 0.69 0.62 
0.70 0.74 0.78 0.79 
0.76 0.81 0.83 0.81 
0.72 0.74 0.81 0.76 
10 0.49 0.51 0.67 0.62 
11 0.35 0.36 0.50 0.51 
12 0.29 0.30 0.39 0.45 
Yearly mean values 0.45 0.47 0.52 0.56 
Yearly variance 0.0342 0.0406 0.0523 0.0293 
Month2001200520102014
0.20 0.23 0.21 0.34 
0.22 0.21 0.20 0.26 
0.31 0.27 0.23 0.38 
0.40 0.39 0.36 0.58 
0.51 0.55 0.57 0.62 
0.45 0.54 0.69 0.62 
0.70 0.74 0.78 0.79 
0.76 0.81 0.83 0.81 
0.72 0.74 0.81 0.76 
10 0.49 0.51 0.67 0.62 
11 0.35 0.36 0.50 0.51 
12 0.29 0.30 0.39 0.45 
Yearly mean values 0.45 0.47 0.52 0.56 
Yearly variance 0.0342 0.0406 0.0523 0.0293 

Note: The bold numbers in the table represent the maximum NDVI and the average NDVI of the corresponding year, so as to clearly see the intra-annual and inter-annual trends of NDVI.

The land-use data from different years only analysed its changes over years, while the use of the NDVI could particularly analyse monthly vegetation changes. The NDVI is represented by the monthly average value during the year. In 2001, 2005, 2010, and 2014 (Table 3 and Figures 4 and 5), the NDVI in the ZYRB increased from January to August and reached a maximum value of approximately 0.8, and the vegetation coverage was the highest at this time. The NDVI values in May, June, July, and August viewed from the spatial distribution, and those in September and October were higher than 0.5. Summer was the season with the highest NDVI value in a year.
Fig. 4

The annual distribution changes of the NDVI in the ZYRB in 2014.

Fig. 4

The annual distribution changes of the NDVI in the ZYRB in 2014.

Close modal
Fig. 5

The maximum value synthesis method was used to extract the annual maximum NDVI.

Fig. 5

The maximum value synthesis method was used to extract the annual maximum NDVI.

Close modal

The NDVI value in the plain area was higher than that in the mountainous area since the main vegetation in the plain area was croplands and the main vegetation in the mountain area was grasslands and forests, also indicating that the growth period of croplands was the main contributor to the NDVI value.

Water consumption calculation results and analysis

Figure 6 presents the results obtained from the calculation of the ZYRB water consumption. From the vegetation types (Figure 7), the highest actual water consumption was found in mixed forests, followed by croplands. The actual water consumption of woody savannas, savannas, grasslands, and deciduous broadleaf forests was medium, while the water consumption of closed shrublands was relatively low. In this paper, the calculation of vegetation water consumption in plain areas mainly considered croplands, and seven vegetation types in mountainous areas were mainly considered: deciduous broadleaf forests, mixed forests, closed shrublands, woody savannas, savannas, grasslands, and permanent wetlands. The results calculated by the area-weighted method show that the average annual water consumption in the plain area was larger than that in the mountainous area (Table 4). In 2005, the average annual water consumption in the plain area was 1.95 mm lower than that in the mountainous area. In 2010, the difference in the average annual water consumption between the plain area and the mountain area was 5.28 mm. The water consumption of cropland in the ZYRB was the highest in April, May, June, July, August, and September (Figure 8), indicating that the water consumption during this period in the plain area was higher than that in the mountainous area, while the vegetation water consumption in the mountainous area was higher than that in the plain area in January, February, March, October, November, and December. In addition, the variation in annual water consumption had the same trend as that of the NDVI. During the three periods from 2001 to 2015, vegetation water consumption first increased, then decreased, and finally increased, but overall water consumption throughout the period was reduced. The actual evapotranspiration of farmland increased slightly. The calculated results were consistent with the inference that the transpiration water consumption of general tree species was relatively large and that of shrubs was relatively small (Wang, 2014).
Table 4

Actual water consumption by terrain classification in the ZYRB in mm.

TopographyVegetation type2001200520102015
Mountainous area Deciduous broadleaf forests 367.62 370.07 416.02 367.55 
Mixed forests 
Closed shrublands 
Woody savannas 
Savannas 
Grasslands 
Plain area Croplands 365.68 375.35 426.04 390.06 
Difference value 1.95 − 5.28 − 10.02 − 22.50 
TopographyVegetation type2001200520102015
Mountainous area Deciduous broadleaf forests 367.62 370.07 416.02 367.55 
Mixed forests 
Closed shrublands 
Woody savannas 
Savannas 
Grasslands 
Plain area Croplands 365.68 375.35 426.04 390.06 
Difference value 1.95 − 5.28 − 10.02 − 22.50 

Note: The bold numbers represent the difference in annual water consumption between the mountainous area and the plain area. We can clearly see that the difference between the water consumption of the plain area and the water consumption of the mountainous area is gradually increasing.

Fig. 6

Vegetation water consumption results. (a) Vegetation water consumption in the ZYRB, 2001–2015 (b) Vegetation water consumption in the ZYRB, 2001–2015 (108m3).

Fig. 6

Vegetation water consumption results. (a) Vegetation water consumption in the ZYRB, 2001–2015 (b) Vegetation water consumption in the ZYRB, 2001–2015 (108m3).

Close modal
Fig. 7

Column chart of actual evapotranspiration (ETa) in the ZYRB from 2001 to 2015.

Fig. 7

Column chart of actual evapotranspiration (ETa) in the ZYRB from 2001 to 2015.

Close modal
Fig. 8

Column chart of monthly ETa in the ZYRB from 2001 to 2015.

Fig. 8

Column chart of monthly ETa in the ZYRB from 2001 to 2015.

Close modal

Effects of climatic factors on water consumption

The main factors affecting vegetation change in the ZYRB were human factors, climatic factors, and terrain factors. Vegetation changes in the basin were mainly represented by the vegetation cover area and NDVI changes in the vegetation index. The water consumption rate per unit area of vegetation was not only related to their own biological characteristics but also closely related to their ecological factors (Wang, 2014). From 2001 to 2015, the climatic factors in the ZYRB also fluctuated under the influence of global climate change (Figure 9).
Fig. 9

Variation trend of climatic factors of different land uses in the ZYRB.

Fig. 9

Variation trend of climatic factors of different land uses in the ZYRB.

Close modal
Stepwise multiple regression analysis was used in SPSS to analyse the correlation between vegetation water consumption changes and climatic factor changes. First, the correlation between water consumption per unit area and multiple climatic factors was analysed. From 2001 to 2005, the water consumption per unit area was positively correlated with sunshine (R2 = 0.585, DW = 2.404, and p = 0.045). From 2005 to 2010, the water consumption per unit area was positively correlated with humidity and pressure and negatively correlated with wind (R2 = 0.988, DW = 3.020, and p = 0.002). From 2010 to 2015, temperature (Pearson correlation: −0.565) and pressure (Pearson correlation: −0.533) had weak negative correlations with the water consumption per unit area. From 2001 to 2015, the water consumption per unit area was positively correlated with temperature (R2 = 0.697, DW = 1.643, and p = 0.019). A multiple linear regression was performed between the water consumption per unit area and meteorological elements for the 2001–2015 period, and the relationship between them was water consumption per unit area (mm) = 8.61, temperature (°C) +54.12, sunshine hours (h) +508.63, relative humidity (1) 66.28, wind velocity (m/s) 4.55 pressure (kPa) – 0.15 precipitation (mm) – 0.60 (Figure 10(a) and 10(b)).
Fig. 10

Relationship between the water consumption per unit area and meteorological elements. (a) Scatter plot. (b) Fitted line plot.

Fig. 10

Relationship between the water consumption per unit area and meteorological elements. (a) Scatter plot. (b) Fitted line plot.

Close modal

From 2001 to 2005, the total vegetation water consumption was positively correlated with precipitation (R2 = 0.887, DW = 1.818, and p = 0.001). From 2005 to 2010, the total vegetation water consumption was negatively correlated with precipitation and wind speed and negatively correlated with relative humidity (R2 = 0.990, DW = 2.288, and p = 0.002). From 2010 to 2015, the total water consumption of vegetation was negatively correlated with precipitation (R2 = 0.606, DW = 2.796, and p = 0.039). In the analysis of the correlation between the total amount of vegetation water consumption and precipitation in the ZYRB, because the water consumption of farmland was the highest, it played a leading role. The total amount of cropland water consumption was consistent with precipitation, and there was a positive correlation trend, and vice versa. The correlation between the total amount of water consumption and the overall trend of precipitation from 2001 to 2015 was negative (Figure 10(b)). From 2001 to 2015, the total water consumption of vegetation was positively correlated with sunshine hours (R2 = 0.871, DW = 1.995, and p = 0.002). Among them, R2 values indicated that the independent variable could explain the degree of the dependent variable. A DW close to 2 showed that the sample was independent, p < 0.05; that is, the level of the sample was less than 0.05, there were significant differences, and the sample was statistically significant.

In short, among the driving factors, the correlation between vegetation water consumption and meteorological factors was mainly annual precipitation and sunshine time, followed by relative humidity, wind speed, and pressure. The total water consumption was highly correlated with precipitation, and the impact of climate change on water consumption was mainly through precipitation and sunshine in regard to climatic factors.

Analysis of the influence of vegetation change on water consumption

In a calculation interval of 15 years, the water consumption of vegetation in the basin showed a decreasing trend from 2001 to 2015. The calculation results were the same as those of Feng (2020), who used the Fubaopu model based on the Budyko hypothesis to estimate the actual evapotranspiration in the upper reaches of the ZYRB from 1961 to 2013, which verified the correctness of the data calculations in this paper. Feng (2020) believed that climate was the main factor affecting the evapotranspiration of vegetation in the basin, and the underlying surface was the secondary factor affecting evapotranspiration. Croplands and urban and construction lands and even barren lands have a great impact on the actual water consumption of vegetation in the ZYRB. The increase in barren land would also reduce the water consumption of vegetation in the basin. For example, the barren area increased by 5.12% from 2001 to 2015, suggesting that the increase in the barren area in the ZYRB was an influencing factor in the reduction in vegetation water consumption in the basin.

The annual water consumption and NDVI values of vegetation in the ZYRB have the same trend. From March to August, there was an increasing trend; in the distribution of a single peak curve, the largest water consumption occurred in summer. The average annual water consumption in the mountainous and plain areas of the ZYRB was calculated (Table 4). The results indicated that the water consumption in plain areas of the ZYRB was greater than that in mountainous areas from 2001 to 2015. The water consumption in mountainous areas increased slightly by 1% from 2001 to 2005, increased by 12% from 2005 to 2010, and decreased by 12% from 2010 to 2015. The water consumption in plain areas increased by 3% from 2001 to 2005, increased by 14% from 2005 to 2010, and decreased by 8% from 2010 to 2015. From 2001 to 2015, the average annual water consumption in mountainous areas decreased by 0.07 mm and that in plain areas increased by 24.38 mm (7%). From the perspective of the changing range, the water consumption in plain areas increased more.

As shown in the table of water consumption calculation results (Table 5 and Figure 6), the water consumption of croplands was the largest by more than 50%, followed by grasslands, accounting for approximately 30%, followed by savannas, closed shrublands, and deciduous broadleaf forests. The woody savannas and mixed forests had the least water consumption, and their water consumption accounted for less than 0.1% of the total water consumption. The annual area of croplands changed greatly; therefore, it had the greatest impact on the total water consumption. The annual water consumption change accounted for 5.59% of the total water consumption from 2001 to 2015 (Table 2 and Table 5). The change in the area of grasslands, closed shrublands, woody savannas, mixed forests, and deciduous broadleaf forests had a positive effect on the change in water consumption, i.e., an increase in total water consumption due to an increase in the area and a decrease in total water consumption due to a decrease in the area. The change in water consumption caused by the cropland area had a negative effect, where the area of croplands decreased but the total water consumption increased because the water consumption per unit area increased by 24.38 mm from 2001 to 2015.

Table 5

Vegetation water consumption changes and dynamics in the ZYRB from 2001 to 2015.

 
 

Note: Red indicates a conspicuous decrease, and yellow indicates a conspicuous increase.

The bivariate correlation analysis between vegetation coverage change and vegetation water consumption change from 2001 to 2005 (Figure 11) showed that there was a positive linear correlation between them. Under the 95% confidence level, the correlation coefficient was 0.747, and the correlation was significant. From 2005 to 2010, the correlation coefficient was 0.171, and the correlation was so weak that the regression equation could not be obtained. This may be due to fewer data or because their relationship was not significantly correlated. In 2010–2015, the correlation coefficient was 0.708 under the 90% confidence level, and the correlation was significant. From 2001 to 2015, there was a positive correlation between the vegetation coverage change and water consumption. The specific regression equation from 2001 to 2015 was y = 0.377x + 84.516, where x is the change in the vegetation area, 106m2, and y is the change in vegetation water consumption, 106m3 (R2 = 0.622, DW = 1.198, and p = 0.035 < 0.05).
Fig. 11

Regression equation and scatter plot (2001–2015).

Fig. 11

Regression equation and scatter plot (2001–2015).

Close modal

For most vegetation, when the area decreased, the total water consumption tended to decrease. When the vegetation area increased, the total water consumption increased (Figure 11). Vegetation restoration increased the water consumption of the region. Therefore, vegetation restoration has a certain threshold and cannot be blindly expanded. It is necessary to combine the actual water resources of the basin to guide reasonable land-use planning. However, when the area of croplands decreased, the total water consumption increased, indicating that the water consumption of croplands per unit area was increasing. This may be due to the increase in high water-consuming crops. This was closely related to the change in cultivated land-planting structure and climate change. Therefore, cultivated crops and artificial plants have water consumption factors that are different from those of natural vegetation. They were the main water users in the basin and needed much irrigation water, which needs to be focused on in the water resource management of the basin.

The purpose of this paper is to quantitatively analyse the relationship between vegetation change and water consumption by analysing the water consumption of different types of vegetation under the ecological restoration scenario, so as to realize the scientific formulation of regional sustainable development water resource policy. The main conclusions are as follows:

  • (1)

    The water consumption per unit area in descending order was mixed forests (660 mm, annual average), croplands (640 mm), closed shrublands (581 mm), deciduous broadleaf forests (528 mm), grasslands (514 mm), savannas (459 mm), and woody savannas (454 mm).

  • (2)

    The main factors affecting vegetation water consumption in the ZYRB were vegetation and sunshine. What's more, there was a positive correlation between vegetation coverage changes and water consumption. The specific regression equation from 2001 to 2015 was y = 0.377x + 84.516, where x is the change in the vegetation area, 106 m2, and y is the change in vegetation water consumption, 106 m3. When vegetation restoration is increased, the water consumption will increase, which will increase the water pressure in the region.

  • (3)

    For water resource policymakers, it is necessary to rationally coordinate water resources and take into account regional ecological environment balance and socially sustainable development. Therefore, it is suggested that formulating water policies should combine the vegetation restoration policies of the ecological environment protection department to achieve the efficient use of water resources.

  • (4)

    At the same time, the deep-seated impact of vegetation restoration on water consumption needs to be further studied in the future. For example, what is the quantitative relationship between other factors and water consumption? How to quantitatively determine the reasonable threshold for regional vegetation restoration? These are the kinds of problems that should be solved.

The authors would like to acknowledge the financial support for this work provided by the National Natural Science Foundation of China (grant no. 52239004).

All relevant data are included in the paper or its Supplementary Information.

The authors declare there is no conflict.

Allen
R. G.
,
Pereira
L. S.
,
Raes
D.
&
Smith
M.
(
1998
).
Crop Evapotranspiration: Guidelines for Computing Crop Water Requirements. FAO Irrigation and Drainage Paper 56. FAO, 56
.
El-Shirbeny
M. A.
,
Biradar
C.
,
Amer
K.
&
Paul
S.
(
2022
).
Evapotranspiration and vegetation cover classifications maps based on cloud computing at the Arab countries scale
.
Earth Systems and Environment
6
,
837
849
.
Feng
Y. R.
(
2020
).
Analysis of Actual Evapotranspiration Characteristics of Typical Watershed Based on Budyko Hypothesis
.
Masteral Dissertation
,
China University of Geosciences
,
Beijing
.
Jiao
Y. Z.
(
2021
).
Spatial and Temporal Distribution of Vegetation Drought and its Lag to Meteorological Factors in Ziya River Basin
.
Masteral Dissertation
,
Hebei University of Engineering
.
Kang
S. Z.
&
Xiong
Y. Z.
(
1991
).
Estimation of crop transpiration by Penman–Monteith model
.
Journal of Northwest A & F University (Natural Science Edition)
19
(
01
),
13
20
.
Li
L.
(
2012
).
Optimization Analysis on Rainwater Utilization Technology in Tianjin Urban District
.
Doctoral Dissertation
,
Tianjin University
.
Li
S. D.
&
Liu
M. C.
(
2022
).
The development process, current situation and prospects of the conversion of farmland to forests and grasses project in China
.
Journal of Resources and Ecology
13
(
01
),
120
128
.
Li
Y.
,
Zhang
C. C.
,
Luo
W. R.
&
Gao
W. J.
(
2019
).
Summer maize phenology monitoring based on normalized difference vegetation index reconstructed with improved maximum value composite
.
Transactions of the Chinese Society of Agricultural Engineering
35
(
14
),
159
165
.
Li
F.
,
Wang
Y.
&
Cao
R.
(
2021
).
Spatial pattern matching analysis of grain production and rainfall in the Ziya River Basin
.
Environmental Monitoring and Assessment
193
(
2
),
1
12
.
Liu
M. L.
(
2014
).
Study of Ecological Environment Evaluation and Ecological Water Requirement Based on Ziya River Basin
.
Doctoral Dissertation
,
Tianjin University
.
People ‘s Government of Hebei Province
(
2015
).
Notice of the People's Government of Hebei Province on Issuing the General Land Use Planning of Hebei Province (2006–2020)
.
People ‘s Government of Shanxi Province
(
2010
).
Notice of the People's Government of Shanxi Province on Issuing the General Land Use Planning of Shanxi Province (2006–2020)
.
Pereira
L. S.
,
Paredes
P.
,
López-Urrea
R.
,
Hunsaker
D. J.
&
Shad
Z. M.
(
2021
).
Standard single and basal crop coefficients for vegetable crops, an update of fao56 crop water requirements approach
.
Agricultural Water Management
243
,
106196
.
Sheng
C. H.
(
2017
).
Study on Water Requirement and Crop Coefficient of Typical Crops in Shule River Basin
.
Doctoral Dissertation
,
Tsinghua University
.
Sullivan
R. C.
,
Kotamarthi
V. R.
&
Feng
Y.
(
2019
).
Recovering evapotranspiration trends from biased CMIP5 simulations and sensitivity to changing climate over North America
.
Journal of Hydrometeorology
20
(
8
),
1619
1633
.
Valipour
M.
,
Bateni
S. M.
,
Sefidkouhi
M. G.
,
Raeinisarjaz
M.
&
Singh
V. P.
(
2020
).
Complexity of forces driving trend of reference evapotranspiration and signals of climate change
.
Atmosphere
11
(
10
),
1081
.
Wang
Y.
(
2014
).
The Forest Vegetation Change and Its Water Consumption Research in Chaobai River Under the Background of Climate Change
.
Doctoral Dissertation
,
Beijing Forestry University
.
Wang
Q.
,
Yang
Z. Y.
&
Sun
G. Y.
(
1999
).
Study on the methods of land use dynamic change research
.
Progress in Geography
18
(
1
),
81
87
.
Wang
Y. P.
,
Yang
P. N.
,
Zhou
L.
,
Xu
S. W.
,
Deng
X. Y.
&
Feng
S. Y.
(
2022
).
Spatio-temporal evolution of vegetation water consumption in lower Reaches of Tarim River
.
Bulletin of Soil and Water Conservation
42
, (
03
),
225
232
.
Xiong
Y.
&
Li
W. T.
(
2007
).
Thoughts on water resources protection of Ziya River
.
Haihe Water Resources
2007
, (
01
),
16
17
.
Yu
X. T.
(
2018
).
Characteristics of Meteorological Droughts in Ziya River Basin and its Impacts on Hydrological Droughts
.
China University of Geosceinces
,
Beijing
.
Zhang
Y.
(
2005
).
Study on Management of Conversion Cropland to Forest and Grassland in Shanxi
.
Doctoral Disseration
,
Shanxi Agricultural University
.
Zhang
D.
(
2006
).
Research on the equation of Penman-Monteith in reference crop evaporation
.
Journal of Anhui Agricultural Sciences
18
,
4513
4514
.
Zhang
D.
,
Liu
X.
,
Zhang
L.
,
Zhang
Q.
,
Gan
R.
&
Li
X.
(
2020
).
Attribution of evapotranspiration changes in humid regions of China from 1982 to 2016
.
Journal of Geophysical Research: Atmospheres
125
(
13
),
e2020JD032404
.
Zhang
Y. Q.
,
Kong
D. D.
,
Zhang
X. Z.
,
Tian
J.
&
Li
C. C.
(
2021
).
Impacts of vegetation changes on global evapotranspiration in the period 2003–2017
.
Acta Geographica Sinica
76
(
3
),
584
594
.
Zhang
R.
,
Song
X.
&
Zhang
J.
(
2022
).
Spatiotemporal responses of vegetation net primary productivity to drought/evapotranspiration in Yellow River Basin
.
Water Resources and Hydropower Engineering
53
(
9
),
57
69
.
Zhao
A. Z.
,
Zhang
A. B.
,
Liu
X. H.
,
Liu
Y. X.
,
Wang
H. F.
&
Wang
D. L.
(
2017
).
Spatiotemporal variation of vegetation coverage before and after implementation of grain for green project in the loess plateau
.
Journal of Natural Resources
32
(
03
),
449
460
.
Zuo
D.
,
Zang
C.
&
Wang
L.
(
2022
).
Temporal and spatial variation of potential evaporation and its influencing factors in the Pearl River Basin during 1980–2019
.
Pearl River
43
(
10
),
41
49
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).