Next Article in Journal
Research on Characteristics of Groundwater Recharge in the Weishan Irrigated District Based on a Bromide Tracer
Previous Article in Journal
Disinfection Methods for Swimming Pool Water: Byproduct Formation and Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Assessment of Green Water Based on the SWAT Model: A Case Study in the Hai River Basin, China

1
School of Resources and Earth Science, China University of Mining and Technology, Xuzhou 221116, Jiangsu, China
2
State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin, China Institute of Water Resources and Hydropower Research, Beijing 100038, China
3
Uni Research Climate, Bjerknes Centre for Climate Research, Jahnebakken 5, 5007 Bergen, Norway
*
Authors to whom correspondence should be addressed.
Water 2018, 10(6), 798; https://doi.org/10.3390/w10060798
Submission received: 12 April 2018 / Revised: 20 May 2018 / Accepted: 14 June 2018 / Published: 16 June 2018
(This article belongs to the Section Hydrology)

Abstract

:
Green water accounts for two-thirds of precipitation, and the proportion could be even higher in dry years. Conflicts between water supply and demand have gradually become severe in the Hai River Basin (HRB) due to the socio-economic development. Thus, the exploitation and the utilization of green water have attracted increasing attention. By gathering the related hydrological, meteorological, and geographic data, the spatiotemporal distribution of green water in HRB and the impacts of land use types on green water are analyzed based on the SWAT (Soil and Water Assessment Tool) model in this study. Furthermore, three new indices are proposed for evaluation, including the maximum possible storage of green water (MSGW), the consumed green water (CGW), and the utilizable green water (UGW). The results show that (1) the MSGW is relatively low in plain areas and its spatial distribution is significantly associated with the soil type; (2) according to the evaluation results of CGW and UGW in HRB, a further improvement of utilization efficiency of green water could be achieved; (3) in general, the utilization efficiency of precipitation in farmlands is higher than other land use types, which means that the planting of appropriate plants could be helpful to enhance the utilization efficiency of green water. Our results summarize the spatiotemporal distribution of green water resource and provide a reference for water resources management in other water-short agricultural areas.

1. Introduction

Water is an indispensable natural resource for human beings, and it plays important roles in maintaining the sustainable development of the ecosystem and the stability of manufacturing activity [1,2]. Falkenmark proposed that water is separated in two main types: blue water and green water [3]. Blue water includes surface and groundwater runoff, and green water refers to precipitation that is stored in the root zone of the soil and that evaporates, transpires, or gets incorporated by plants. Traditional methods for evaluating water resources were mainly focused on blue water, whereas green water has been overlooked because its availability is limited (it cannot be directly exploited). However, green water is the principle water source that sustains the whole terrestrial ecosystem, including crops, grasslands, and forests, which guarantee the food security for people around the world [4,5,6,7]. About 60% of food production in the world depends on green water, and almost all meat products of animal husbandry and forest products rely on green water [8]. In the components of global annual average water footprint, green water footprint (consumption of green water) and blue water footprint (consumption of surface and ground water) account for 74% and 11%, respectively [9]. With such high contribution to food production, green water is a vital resource that requires more attention in terms of efficient utilization. Therefore, it is essential to study and evaluate green water.
Since green water was first proposed, it has been continually studied and its utility has improved. Falkenmark proposed that green water contains two parts: green water flow and the green water storage [10]. Rockström further divided green water flow into the productive part (the actual transpiration), which impacts plant biomass, and the non-productive part (the actual evaporation) [11]. In recent year, more research has been done on assessment of soil water resources and the corresponding evaluation indices system has been established [12,13]. The green water has gradually attracted attention as a significant resource in hydrology research and water resources assessment [5,14,15,16]. For example, some studies focused on the quantity and spatiotemporal distribution of green water [15,16,17,18,19], numerous studied have also been done on the impacts of land use types on green water [20,21,22]. Furthermore, the studies on water footprint provides a reference for the green water resource research [23,24]. These studies were particularly useful for developing an overview of the green water resources status and helped to find out areas where further analysis may be indispensable. Green water flow and green water storage were selected as the main evaluation indicators in those previous studies. However, these indicators may not be clear enough for the available water evaluation in practical applications. For instance, it is obvious that not all the soil water content could be absorbed by plants, therefore, the use of green water storage to estimate the quantity of green water resources is not reliable in a way and may even mislead the water resources management policy.
In this study, three new indicators (MSGW, CGW, and UGW) are proposed for the appraisal of the green water resource more holistically. The depth of plant roots was considered when analyzing the MSGW (the maximum possible storage of green water), which represented the theoretical utilization potential of green water in a watershed. Evaporation was regarded as the unproductive part of numerous studies [23,25,26]. However, evaporation and transpiration are closely interlinked hydrologic processes, and it is supposed that evaporated water has the potential to be partly used for food production [27]. Therefore, the total actual evapotranspiration in farmlands, grasslands, and woodlands was considered to be consumed green water (CGW) in this study. Evapotranspiration from the urban area and water body cannot be utilized and has been defined as non-utilized green water (NGW). CGW provide the consumption status of green water. To quantitatively evaluate the development and utilization value of green water storage, UGW (utilizable green water) was proposed to represent the part of green water, which can still be utilized by plants.
Agricultural droughts occur in almost all areas in the world. The problem is usually more severe because the potential evapotranspiration is considerably higher than precipitation in a dryland farming area [28,29]. The Hai River Basin (HRB) is one of the main agricultural areas in China. However, problems of agricultural drought existed, and the water consumption was mainly provided by groundwater in HRB [30]. With the development of the urban construction and agriculture, the basic demand for water kept increasing, which caused a serious overdrawing of the shallow groundwater. Many scholars have paid more attention to the precipitation, the surface water, and the groundwater in HRB than the green water [31,32,33,34]. Since green water is the main water resource consumed, this paper aimed at making a comprehensive assessment of the green water with the new indicators (MSGW, CGW, and UGW) in an agricultural region that is facing a water shortage problem. By using the SWAT (Soil and Water Assessment Tool) model, the hydrological cycle processes for the period of 1995–2004 in HRB were simulated. The evaluation of green water could provide practical and efficient references for water management.

2. Materials and Methods

2.1. Study Area

The study area is the Hai River basin with a drainage area of 3.18 × 105 km2, located in a semi-humid zone in northern China (35° N–43° N and 112° E–120° E). To make it convenient for management, the Hai River Water Conservancy Commission divided the whole basin into 15 hydrological regions according to the topography, distribution of river net and administrative region (Figure 1) [35]. The mean annual precipitation in HRB is 539 mm with a high monthly variability. It is rainless in autumn and winter and dry in spring with a high evaporation rate, which leads to frequent droughts. However, it is warm and wet in summer because of the abundant precipitation and rainstorms. The mean annual temperature is 7.8–13.5 °C and the relative humidity is 50%–70%. The mean annual land evaporation and water surface evaporation in the study area are 500 mm and 1100 mm, respectively.
Approximately 22% of total rainfall runoff is the water consumption in the watershed, and more than 60% of it is used in agricultural areas. HRB needs water from other river basins through the water diversion project to satisfy the water requirements during the past few years [30]. Maize and winter wheat are the main crops in HRB [30], the accurate quantification and spatiotemporal distribution of green water could potentially facilitate the selection of more appropriate copping patterns in HRB, to reduce evapotranspiration and utilize green water more efficiently. Furthermore, it could achieve the goal of reducing water consumption from surface water and groundwater, which is used for agricultural production, and alleviating the problem of the water shortage.

2.2. Data

The daily meteorological data from 43 meteorological stations (including air temperature, wind speed, solar radiation, and relative humidity) and precipitation data from 216 stations in the study region (1995–2004), were provided by the China Meteorological Administration (CMA) (http://www.cma.gov.cn/). The data is complete in most stations, dozens of precipitation stations, the data series were completed by using the data in adjacent station. The data quality check haven been made by CMA (http://data.cma.cn/).The hydrological data from 9 stations were collected from the Hai River Water Conservancy Commission (http://www.hwcc.gov.cn/). These data are the fundamental inputs for the SWAT model and served as the driving force for other simulation processes in this model. The distribution of these stations is shown in Figure 2. The soil data (1:1,000,000) were provided by the Institute of Soil Science of the Chinese Academy of Sciences. Figure 3a shows the distribution of different soil types, which can be divided into 27 categories. The soil classification is based on the Genetic Soil Classification of China (GSCC) [36]. Since the soil particle graded standard in GSCC differs from the standard in the SWAT model (the standard of United States Department of Agriculture), the standard conversion was completed by using the cubic spline interpolation [37]. Then, the soil parameters (such as the saturated hydraulic conductivity and bulk density) were estimated on the basis of the soil texture (the proportion of clay, silt, sand, and rock) using the Soil Water Characteristics Hydraulic Properties Calculator [38]. The DEM (Digital Elevation Model) data was collected from the Data Cloud of Chinese Academy of Sciences (http://www.csdb.cn/) and it was extracted from the SRTM DEM (Shuttle Radar Topography Mission Digital Elevation Model) with a resolution of 90 m. The land use data in 2004 with a spatial resolution of 250 m × 250 m were provided by Institute of Remote Sensing and Digital Earth of the Chinese Academy of Sciences. The source data is the Landsat TM/ETM (Thematic Mapper/Enhanced Thematic Mapper) remote sensing data and the land use data was generated through artificial visual interpretation. The land use classification is based on the relevant government standard in China [39]. The dominant types of land use are farmland, woodland, and grassland, which collectively account for 94% of the study region (Figure 3b).

2.3. Methodology

2.3.1. SWAT Model and Simulation Setup

The SWAT model has been widely applied in studies at different regions because of its high precision and wide applicability [40]. The SWAT simulator version 2000 and its ArcView interface [41,42,43] were used in this study. In SWAT modeling, the study basin is divided into many sub-basins based on the digital elevation model (DEM) data. Then, each sub-basin can be further subdivided into several hydrologic response units (HRU) which consist of similar land use types, soil, and slope characteristics [44]. There are 283 sub-basins in total and 2100 HRUs in the Hai river basin for SWAT modeling. Different hydrological components, for example, evapotranspiration, infiltration, and water surface evaporation are calculated in each HRU. Eventually, it obtains the hydrological variables, such as evapotranspiration and surface runoff, in each sub-basin. The modeling process adopts the water balance method, which is expressed as follows:
S W t = S W 0 + i = 1 t ( R d a y , i Q s e r f , i E i W s e e p , i Q g w , i )
where SWt is the final soil water content (mm), SW0 is the initial soil water content on day i (mm), t is the time (day), Rday,i is the amount of precipitation on day i (mm), Qsurf,i is the amount of surface runoff on day i (mm), Ei is the amount of evapotranspiration on day i (mm), Wseep,i is the amount of water entering the unsaturated zone from the soil profile on day i (mm), Qgw,i is the amount of return flow on day i (mm).
There are multiple options for runoff generation and evapotranspiration estimation, which are provided in the SWAT model. According to the previous research [31,45], the soil conservation service curve number method and the Priestley–Taylor method were applied to the runoff generation and evapotranspiration calculation in this study, respectively. The monthly assessment results of green water have practical significance to guide the agricultural production. Therefore, the time step of SWAT modeling is monthly in the study.
Because only about 6.3% areas were irrigated in the mountainous area [30], the irrigation has been ignored during model simulation. Due to the lack of exactly data about irrigation time and amount for irrigation zone in plain area, the irrigation in study area was triggered by a soil water deficit threshold [46]. In agricultural HRUs, the soil water deficit threshold was defined as 75% of the soil field capacity according to relevant studies [47,48].
A sensitivity analysis can dramatically reduce the number of parameters in model calibration. By using the LH-OAT (Latin hypercube one-factor-at-a-time) method, Sun selected the sensitive parameters for the mountains area and plain area in HRB, respectively [49]. Wang chose four main parameters to calibrate the SWAT model in HRB with observed runoff data [31]. Based on the former research, 12 sensitivity parameters were identified and selected for further model calibration of the sub-basins in this study (Table 1). Then, the final model calibration has been done by using SUFI-2 (sequential uncertainty fitting algorithm) method [50]. The Nash-Sutcliffe efficiency coefficients (NS) and the coefficient of determination (R2) are chosen as the indices for the model calibration and validation in the study [51]. Both are as much closer to be 1 as better, which means a good fitting result. The model performance is considered to be satisfactory if the NS > 0.5, according to the model evaluation guidelines for watershed simulations [52,53].

2.3.2. Evaluation Methods

Previous studies of evaluating green water were mainly focused on the green water flow (actual ET) and the green water storage (that is, soil water) [15,16,17,18,19]. To conduct a more comprehensive analysis of the quantity and distribution of green water, three evaluation indices are introduced in this study.
1. The maximum possible storage of green water (MSGW)
The soil moisture content could range from standard moisture content to the saturation moisture content. In general, it is widely recognized that the soil moisture could only be absorbed by plants when it ranges between field capacity (FC) and permanent wilting percentage (Wp). This kind of soil moisture is defined as the available water capacity (AWC). The maximum possible storage of green water can be calculated by using AWC and the corresponding depth of soil layer (D), where D is determined by the depth of the plant roots.
A W C i = F C i W P , i
W m s = i = 1 n A W C i × D i
where Wms is the maximum possible storage of green water, AWCi (mm) is the available water capacity in soil layer i, FCi (mm) is the field capacity in soil layer i, WP,i is the permanent wilting percentage in soil layer i (mm), Di is the depth of soil layer i (mm).
2. Consumed green water (CGW)
The land use types have been divided into farmland, grassland, woodland, urban area, and surface water body. In general, evapotranspiration from urban areas and water bodies cannot be utilized, which are usually regarded as nonproductive green water. This part of evapotranspiration has been defined as non-utilized green water (NGW) in this paper, and the rest of the evapotranspiration refers to the consumed green water (CGW). The formula is defined as the following:
W c = i = 1 n E i + i = 1 n T i
where Wc is the consumed green water (mm), Ti is the transpiration in different plants species (mm), and Ei is the evaporation between plants in the different plant species (mm).
3. Utilizable green water (UGW)
The utilizable green water is calculated according to the permanent wilting percentage (Wp) and soil water content (SW), and the formula is expressed as follows:
W u = i = 1 n S W i W p , i
where Wu is the utilizable green water, SWi is the soil water content in soil layer i (mm), WP,i is the permanent wilting percentage in soil layer i (mm).
Blue water can directly be obtained from the simulation results [40]. The three new indices represent the storage volume, utilization state, and the exploitation potential of green water. After the soil moisture parameters of the different soil types (including the available water capacity (AWC), field capacity (FC), permanent wilting percentage (Wp), and soil water content (SW)) have been calculated by the SWAT model, we can obtain the distribution of the MSGW and UGW in each sub-basin through area weighting. Using the same method, the CGW in a sub-basin equals the difference between the evapotranspiration (ET) and NGW.

3. Results

3.1. Analysis of the Simulation Results

3.1.1. Model Calibration and Validation

The model calibration was divided into two parts in the HRB according to the locations of sub-basin and its runoff data quality. For the sub-basins in the mountainous area with relatively low human impacts, the model calibrations are based on the monthly runoff. The model was calibrated from 1995 to 1999 and validated from 2000 to 2004 for the mountainous area based on the observed monthly runoff data. The final range of the parameters after model calibration are shown in Table 2. The NS > 0.65 and R2 > 0.8 for each discharge station during calibration are shown in Table 3. The model was validated using the calibrated parameters and the results are also satisfactory, as shown in Table 3. The monthly runoff at Zhangjiakou station and Dongyanghe station are further illustrated in Figure 4.
In this study, the yearly quantity of the river discharge into the sea in HRB was used for model calibration for the sub-basins at plain hydrological regions (Table 4 and Table 5). The monthly remote sensing ET and simulated ET in the plain hydrological region have been compared (Table 6). The ET data with a resolution of 1000 m in 2002–2004 were derived from the remote sensing data of MODIS (Moderate-resolution Imaging Spectroradiometer) by using the ETWatch model [54]. The relevant methodology and the results of the remote sensing data processing have already been published, which is based on a separate project by the Institute of Remote Sensing Applications, the Chinese Academy of Sciences [54,55]. The range of NS and R2 is 0.60–0.71 and 0.72–0.83, respectively, which show that the SWAT model performance works well in ET simulation. The final ranges of the calibrated parameters in the plain areas are shown in Table 4. Based on the calibration and validation, the water balance analysis and green water assessment are implemented in the following section.

3.1.2. Water Balance Analysis

In the study, the water balance is analyzed according to three components: precipitation, green water, and blue water. Figure 5 shows the annual precipitation, green water, and blue water in HRB from 1995 to 2004. The annual variation of soil water during the study period (1995–2004) is less than −2 mm, it indicates that the three components (precipitation, green water, and blue water), which are obtained from the SWAT model simulation, achieves a good balance (nearly all precipitation transforms into blue water and green water). The annual mean green water in HRB is about 365 mm, which accounts for 74% of the precipitation. In dry years (such as 1997), the proportion can reach up to 90%, it is because the precipitation is more likely to absorbed by soil and utilized by plants instead of formed the river discharge or infiltrated into the aquifer. Therefore, green water plays an important role in the water cycle. The improvement of the utilization efficiency of green water might be one of the keys to relieving the water shortages in HRB.

3.2. Assessment of Green Water

3.2.1. Distribution of the Maximum Possible Storage of Green Water

The average depth of the plant roots in HRB is obtained by calculating the area-weighted average of all the sub-basins. The value is 933 mm and the maximum possible storage of the green water (MSGW) is 106 mm (or about 33.2 km3). The spatial distribution of the MSGW is shown in Figure 6. According to Equation (3), the main influencing factors of MSGW are the depth of the plant roots and soil moisture capacity. From Figure 3 and Figure 6, we can see that the distribution of MSGW has a significant correlation with soil types. Loessial soil is the dominating soil type in the upper reaches of the Zhangwei River (CZWRB) and the Ziya River (CZYRB). The field capacity (FC) and permanent wilting percentage (Wp) are 157.6 mm and 40.9 mm, respectively. According to Equation (2), the available water capacity is 116.7 mm. Cinnamon soil is the dominating soil type in the piedmont plain areas except for the upper reaches of Luan River (CLRB), which is mixed with brown soil. The AWC in cinnamon soil and brown soil are 110.5 mm and 111.7 mm, respectively. Due to a low AWC of 95.1 mm and a shallow depth of plant roots, the MSGW in the plain area is lower than other areas. The upper reaches of Yongding River, including hydrological regions of the West of Yongding River basin (WDRB) and the East of Yongding River basin (EDRB), locates at the west part of HRB. The FC over there is 198.9 mm and WP is 98.7 mm, which is higher than that of its surrounding area. However, the AWC is comparatively low (100.1 mm), indicating that the MSGW is relatively poor in these areas.

3.2.2. Analysis of the Consumed Green Water and Utilizable Green Water

The consumed green water (CGW) and the utilizable green water (UGW) in HRB are calculated based on Equations (4) and (5). As is shown in Figure 7, the annual mean value of blue water in HRB is 129 mm, which accounts for 26.2% of the precipitation. The annual mean of CGW and UGW are 339 mm and 27 mm, respectively. The results show that green water accounts for 74.1% of the precipitation. The sum of blue and green water components is 100.3%. This is because the annual mean variation of the soil water content is −1.7 mm (accounts for the −0.3% of the precipitation). Therefore, it is unlikely to remit the water shortage in HRB by the integrated management of blue water. How to properly use the green water will become one of the keys to solving conflicts between the supply and demand of water in the HRB. According to the quantitative analysis, the annual mean value of UGW is 27 mm, which means that the green water resource in the whole HRB still has not been fully utilized. The UGW is still higher than 10 mm, even in dry years (that is, 1997, 1999, and 2002).
The distribution of green water is further analyzed and shown in Figure 8a,b, from which we can see that the CGW and UGW in the HRB have a similar spatial pattern, which indicates that the green water is abundant in CZWRB, CZYRB and CLRB. These areas have comparatively large values of CGW and UGW. Thus, there is a huge potential to boost the utilization efficiency of green water. The evaporation from urban areas cannot be utilized, and most of cities and towns are distributed in the plain areas. This helps to explain why the value of CGW is relatively low in those areas. A comparatively large amount of CGW (350 mm) can be observed in the East of the Yongding River Basin (EDRB), but the UGW is close to zero in this area. This indicates that the green water in the EDRB has been fully utilized in transpiration. The green water is relatively poor in the plain areas and particularly in the piedmont plain, where the CGW is less than 350 mm and the UGW is less than 40 mm. Considering that most of the cultivated lands are distributed in these areas, the green water cannot meet the requirements of crops, and irrigation is indispensable.

3.2.3. Impacts of Land Use Types on Green Water

The effective utilization of precipitation (the percentage of precipitation that transforms into green water) varied with land use/cover. The components of annual mean evapotranspiration are shown in Table 7. The efficient evapotranspiration, which equals to CGW, is 106.9 km3 (amounting to 340 mm), which accounts for about 70% of the precipitation or 93% of the green water. It is much more than the non-utilized green water (NGW) of 7.9 km3 (amount to 25 mm). Most of the CGW is in farmlands, woodlands, and grasslands (55%, 23% and 22%, respectively.)
Land use types in HRB mainly includes farmlands, woodlands, and grasslands. Table 8 shows the effective utilization of precipitation in those three land use types which demonstrates: (1) the annual averages of effective utilization of precipitation are similar in different land use types; (2) the effective utilization of precipitation in cultivated land is higher than that of grasslands and woodlands for most years; (3) the effective utilization of precipitation is comparatively larger in dry years, especially in 1997, 1999, and 2002, which exceed 80%, no matter what kind of land use type . In particular, it exceeds 90% in farmlands.

3.2.4. Monthly Variation of Green Water

The monthly variation of green water in the HRB was analyzed in this section. The annual precipitation in the HRB is 492 mm during 1995–2004, which is similar to the amount of precipitation (490 mm) in 2000. In this case, we selected 2000 as the representative year to analyses the intra-annual changes of green water. The variations of CGW, UGW, and the other main hydrological cycle elements are shown in Table 9. The green water, blue water, and the variation of soil water content in 2000 account for 66.8%, 28.0% and 4.8% of precipitation, respectively. The CGW in HRB kept rising in the first six months while the UGW fell. This is because the consumption of green water by plants was increasing during the transition from the cold season to the warm season. The CGW stays at a relatively high level in July and August, while the UGW is in a rapid growth because of the abundant precipitation of those two months, which accounts for 55% of the annual precipitation. From September to December, the CGW and UGW in the HRB decreased, while the rate of the latter was relatively low and remained above 35 mm. In this case, improvement of the utilization efficiency of green water can still be explored to increase the CGW and UGW in HRB.

4. Discussion and Conclusions

4.1. Discussion

Green water accounts for 74% of the mean annual precipitation in the HRB (it is higher in dry years), the proportion is similar in other study areas with a similar climate [17]. Numerous studies have demonstrated that there are positive correlations between the precipitation and green water at both spatial and temporal resolutions [56,57,58]. Other meteorological factors such as temperature can also have an impact on the green water flow and storage [40]. Besides, land use change can increase or decrease the green water storage [40]. According to the distributions of topography and land use types, the result indicates that green water has a clear spatial pattern. Based on the definition of green water, there must be positive correlation between soil moisture and green water, and this is supported by the spatial pattern of soil moisture in HRB [59]. Therefore, the soil moisture can be used for a quick preliminary estimation of green water with the use of remote sensing technology.
During the operational process of the SWAT model, a land use type will not be regarded as a hydrologic response unit if its area is less than 5% of the area of a sub-basin. Therefore, in this study, the area of farmland is 3.6% higher than the value directly obtained from land use data. According to the assessment results, the utilization rate of precipitation in farmlands (73.90%) is relatively higher than grasslands (73.38%) and woodlands (72.55%). As with the study at Weihe River Basin in northeast China [40], the land use change (an increase of woodland and grassland, a decrease of farmland) leads to the decrease of the green water flow, which supports our results.
Both CGW and UGW have high seasonal variability. The intra-annual distribution of precipitation and the process of plants growth are the dominant factors [60]. Maize and winter wheat are the main crops in the HRB [30], considering their relatively long growth period (2–5 months for maize and 8–10 months for winter wheat), the monthly assessment results of green water have practical significance to guide the agricultural production [61].
There are several limitations in this study, one of which is that the land use types in HRB have been divided into five categories broadly. Green water utilization varies among plant types [9,23]. Thus, to utilize green water more efficiently, it is necessary to analyze the consumption of green water in different plants (especially between various crops) in the future. The green water in HRB has been analyzed in a ten years period and the change trend of green water in a longer period can be further studied. Furthermore, identifying the impact of climate variability and land use change on green water can give us great insight into the effective utilization of green water. To provide more accurate guidance for the crops with a shorter growth period, the weekly (or even daily) assessment of green water could be further learned with corresponding initial data.
The SWAT model in HRB has been built at a relatively high spatial and temporal resolution with the limited freely available input data. The observed runoff and evapotranspiration data were used for testing the performance of the model; the water balance analysis provides further evidence for the rationality of the simulation results. With this model, several applications such as the variations of green water under different climate/land use scenarios in the past or future could be analyzed with befitting data. However, there are still some shortcuts in this model. For example, during model calibration, only nine hydrological stations were used in the mountain areas and the yearly data were used in plain areas, which may affect the simulation accuracy. Overall, the performance of the model (e.g., the accuracy of outcomes and the minimization of uncertainties) could be enhanced in future studies with the use of more refined datasets.

4.2. Conclusions

In this paper, three new indices, including the maximum possible storage of green water (MSGW), the consumed green water (CGW), and the utilizable green water (UGW), are proposed to quantitatively analyze green water use in the Hai River Basin. These indices were applied to the assessment of green water in HRB, their spatial and temporal distributions have been further investigated. The following conclusions are drawn from this study:
(1)
Green water is the dominant component of precipitation in the studied area. With the evaluation results of MSGW, CGW and UGW at a sub-basin level, it is feasible to abate the shortage of water by reducing the evapotranspiration in the sub-basins with high CGW and further developing the soil water content in the sub-basins with high UGW.
(2)
The distribution of CGW is significantly associated with the distribution of the land use type. Farmland is the main source of CGW in HRB (Table 7); growing low water consuming crops or returning farmland to forests and grasslands deliberately could help to reduce the CGW.
(3)
The monthly variate conditions of green water could be obtained from the model. With high accuracy land use data, the efficient utilization patterns of green water could be developed for the different vegetative stages of a specific plant. Model with daily inputs can provide more accurate results because it considers the growth period of a specific plant in detail.
Overall, this paper provided a feasible direction to complement the evaluation system of green water. For the water-deficient area (especially agricultural region), a full valuation of the green water resources could help in turning the proportion of blue and green water back to a reasonable range. Furthermore, our results can be potentially useful for alleviating the contradictions of water demand between the ecosystem and human beings. Additionally, the proper management mode could be prepared in advance for new challenges brought about by climate/land use change.

Author Contributions

K.Z., Y.Z. and F.L. came up with the idea and designed the study; Z.X. and X.S. (Xinyi Song) performed the analyses and wrote the paper; L.L. helped with language editing, and X.S. (Xiaomeng Song) helped in collecting data.

Funding

This research was supported by the National Key Research and Development Program (No. 2016YFC0401301), the Natural Science Foundation of China (No. 51409246 and 51679252), the Jiangsu Province Natural Science Foundation (No. BK20130178), and the Open Research Fund of State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin (China Institute of Water Resources and Hydropower Research) (No. IWHR-SKL-201603).

Acknowledgments

Thanks to editor and reviewers for their constructive comments and suggestions, which greatly help to improve the quality of our article. We would also like to thank Ravi Jain for his help to improve the English in our article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mooney, H.; Cropper, A.; Reid, W. Confronting the human dilemma. Nature 2005, 434, 561–562. [Google Scholar] [CrossRef] [PubMed]
  2. Oki, T.; Kanae, S. Global hydrological cycles and world water resources. Science 2006, 313, 1068–1072. [Google Scholar] [CrossRef] [PubMed]
  3. Falkenmark, M. Land-Water Linkages: A Synopsis: Land and Water Integration and River Basin Management; Land and Water Bulletin No. 1; FAO: Rome, Italy, 1995; pp. 15–17.
  4. Schyns, J.F.; Hoekstra, A.Y.; Booij, M.J. Review and classification of indicators of green water availability and scarcity. Hydrol. Earth Syst. Sci. 2015, 12, 4581–4608. [Google Scholar] [CrossRef]
  5. Rockström, J. Green water security for the food makers of tomorrow: Windows of opportunity in drought-prone savannahs. Water Sci. Technol. 2001, 43, 71–78. [Google Scholar] [CrossRef] [PubMed]
  6. Gordon, L.J.; Finlayson, C.M.; Falkenmark, M.; Molden, D.; De Fraiture, C. Managing water in agriculture for food production and other ecosystem services. Agric. Water Manag. 2010, 97, 512–519. [Google Scholar] [CrossRef]
  7. Vanham, D. A holistic water balance of Austria—How does the quantitative proportion of urban water requirements relate to other users? Water Sci. Technol. 2012, 66, 549–555. [Google Scholar] [CrossRef] [PubMed]
  8. Cosgrove, W.J.; Rijsberman, F.R. Challenge for the 21st century: Making water everybody’s business. Sustain. Dev. Int. 2000, 2, 149–156. [Google Scholar]
  9. Hoekstra, A.Y.; Mekonnen, M.M. The water footprint of humanity. Proc. Natl. Acad. Sci. USA 2012, 109, 3232–3237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Rockström, J. On-farm green water estimates as a tool for increased food production in water scarce regions. Phys. Chem. Earth 1999, 24, 375–383. [Google Scholar] [CrossRef]
  11. Falkenmark, M.; Rockström, J. The new blue and green water paradigm: Breaking new ground for water resources planning and management. J. Water Resour. Plan. Manag. 2006, 132, 129–132. [Google Scholar] [CrossRef]
  12. Xia, Z.Q.; Li, Q.F. Study on soil water resource and its evaluating methodology. Adv. Water Sci. 2001, 12, 535–540. (In Chinese) [Google Scholar]
  13. Wang, H.; Yang, G.Y.; Jia, Y.W.; Wang, J.H. Connotation and assessment index system of soil water resources. J. Hydraul. Eng. 2006, 37, 389–394. (In Chinese) [Google Scholar]
  14. Gerten, D.; Hoff, H.; Bondeau, A.; Lucht, W.; Smith, P.; Zaehle, S. Contemporary “green” water flows: Simulations with a dynamic global vegetation and water balance model. Phys. Chem. Earth 2005, 30, 334–338. [Google Scholar] [CrossRef]
  15. Zhao, A.Z.; Zhu, X.F.; Chen, S.C.; Li, M. Trend analysis for blue and green water resources in the Weihe River basin of Northwest China during the past thirty years and in the near future. In Proceedings of the International Conference on Agro-Geoinformatics, Beijing, China, 11–14 August 2014; pp. 1–5. [Google Scholar]
  16. Konar, M.; Dalin, C.; Hanasaki, N.; Rinaldo, A.; Rodriguez-Iturbe, I. Temporal dynamics of blue and green virtual water trade networks. Water Resour. Res. 2012, 48, 2360–2368. [Google Scholar] [CrossRef]
  17. Zhang, W.; Zha, X.; Li, J.; Liang, W.; Ma, Y.; Fan, D.; Li, S. Spatiotemporal change of blue water and green water resources in the headwater of yellow river basin, China. Water Resour. Manag. 2014, 28, 4715–4732. [Google Scholar] [CrossRef]
  18. Hanasaki, N.; Inuzuka, T.; Kanae, S.; Oki, T.; Hoff, H.; Falkenmark, M. An estimation of global virtual water flow and sources of water withdrawal for major crops and livestock products using a global hydrological model. J. Hydrol. 2010, 384, 232–244. [Google Scholar] [CrossRef]
  19. Hoff, H.; Falkenmark, M.; Gerten, D.; Gordon, L.; Karlberg, L.; Rockström, J. Greening the global water system. J. Hydrol. 2010, 384, 177–186. [Google Scholar] [CrossRef]
  20. Jewitt, G.P.W.; Garratt, J.A.; Calder, I.R.; Fuller, L. Water resources planning and modelling tools for the assessment of land use change in the Luvuvhu catchment, South Africa. Phys. Chem. Earth 2004, 29, 1233–1241. [Google Scholar] [CrossRef]
  21. Cao, X.C.; Wu, P.T.; Wang, Y.B.; Zhao, X.N. Assessing blue and green water utilisation in wheat production of China from the perspectives of water footprint and total water use. Hydrol. Earth Syst. Sci. 2014, 18, 3165–3178. [Google Scholar] [CrossRef] [Green Version]
  22. Quinteiro, P.; Dias, A.C.; Silva, M.; Ridoutt, B.G.; Arroja, L. A contribution to the environmental impact assessment of green water flows. J. Clean. Prod. 2015, 93, 318–329. [Google Scholar] [CrossRef]
  23. Zoumides, C.; Bruggeman, A.; Hadjikakou, M.; Zachariadis, T. Policy-relevant indicators for semi-arid nations: The water footprint of crop production and supply utilization of cyprus. Ecol. Indic. 2014, 43, 205–214. [Google Scholar] [CrossRef]
  24. Gobin, A.; Kersebaum, K.C.; Eitzinger, J.; Trnka, M.; Hlavinka, P.; Takáč, J.; Kroes, J.; Ventrella, D.; Marta, A.D.; Deelstra, J.; Lalić, B. Variability in the water footprint of arable crop production across European regions. Water 2017, 9, 93. [Google Scholar] [CrossRef]
  25. Bulcock, H.H.; Gush, M.B.; Jewitt, G.P.W. A comparison of productive and non-productive green water-use efficiency of Podocarpus henkelii and Pinus patula in the Kwazulu-natal midlands. South. For. J. For. Sci. 2014, 76, 75–84. [Google Scholar]
  26. Willaarts, B.A.; Volk, M.; Aguilera, P.A. Assessing the ecosystem services supplied by freshwater flows in mediterranean agroecosystems. Agric. Water Manag. 2012, 105, 21–31. [Google Scholar] [CrossRef]
  27. Schuol, J.; Abbaspour, K.C.; Yang, H.; Srinivasan, R.; Zehnder, A.J.B. Modeling blue and green water availability in Africa. Water Resour. Res. 2008, 44, 212–221. [Google Scholar] [CrossRef]
  28. Sposito, G. Green water and global food security. Vadose Zone J. 2013, 12, 1742–1751. [Google Scholar] [CrossRef]
  29. Stewart, B.A.; Peterson, G.A. Managing green water in dryland agriculture. Agron. J. 2015, 107, 1544–1553. [Google Scholar] [CrossRef]
  30. Water Resources Bulletin of the Hai River Basin. Available online: http://www.hwcc.gov.cn/hwcc/static/szygb/hhlyszygb.htm (accessed on 20 October 2017).
  31. Wang, Z.G.; Luo, Y.Z.; Zhu, X.J.; Rui, W.; Li, W.; Zhang, M.H. Watershed modeling of surface water-groundwater interaction under projected climate change and water management in the Haihe river basin, China. Br. J. Environ. Clim. Chang. 2013, 3, 421–443. [Google Scholar] [CrossRef]
  32. Xia, J.; Qiu, B.; Li, Y.Y. Water resources vulnerability and adaptive management in the Huang, Huai and Hai River basins of China. Water Int. 2012, 37, 523–536. [Google Scholar] [CrossRef]
  33. Zhang, D.D.; Yan, D.H.; Wang, Y.C.; Lu, F.; Wu, D. Changes in extreme precipitation in the huang-huai-hai river basin of China during 1960–2010. Theor. Appl. Climatol. 2015, 120, 195–209. [Google Scholar] [CrossRef]
  34. Yan, D.; Xu, T.; Girma, A.; Yuan, Z.; Weng, B.; Qin, T. Regional correlation between precipitation and vegetation in the Huang-Huai-Hai River Basin, China. Water 2017, 9, 557. [Google Scholar] [CrossRef]
  35. Ren, X.S.; Hu, Z.L.; Cao, Y.B.; He, S. Water Resources Assessment in the Haihe River Basin; Water Resources and Electricity Press: Beijing, China, 2007. (In Chinese) [Google Scholar]
  36. Shi, X.Z.; Yu, D.S.; Warner, E.D.; Pan, X.Z.; Petersen, G.W.; Gong, Z.G. Soil database of 1:1,000,000 digital soil survey and reference system of the Chinese genetic soil classification system. Soil Surv. Horiz. 2004, 45, 129. [Google Scholar] [CrossRef]
  37. Wei, H.B.; Zhang, Z.P.; Yang, J.P. Establishing method for soil database of SWAT model. Water Resour. Hydropower Eng. 2007, 38, 15–18. (In Chinese) [Google Scholar]
  38. Saxton, K.E. Soil Water Characteristics Hydraulic Properties Calculator. USDA-Agricultural Research Service in Cooperation with Department of Biological Systems Engineering, Washington State University, Pullman, WA. Available online: http://hydrolab.arsusda.gov/soilwater/Index.htm (accessed on 9 March 2018).
  39. Standardization Administration of China. Current Land Use Condition Classification. In China’s National Standard (GB/T21010-2007); Standards Press of China: Beijing, China, 2007. [Google Scholar]
  40. Zhao, A.; Zhu, X.; Liu, X.; Pan, Y.; Zuo, D. Impacts of land use change and climate variability on green and blue water resources in the Weihe river basin of northwest China. Catena 2016, 137, 318–327. [Google Scholar] [CrossRef]
  41. Neitsch, S.L.; Arnold, J.G.; Kiniry, J.R.; Srinivasan, R.; Williams, J.R. Soil and water assessment tool user’s manual version 2000. In GSWRL Report 202; Texas Water Resources Institute: College Station, TX, USA, 2002. [Google Scholar]
  42. Van Orshoven, J.; Coppens, J.; Terres, J.M. Possibilities and limitations of avswat2000 for the assessment of environmental impacts of farming practices. Int. J. Climatol. 2007, 32, 261–273. [Google Scholar]
  43. Zhang, Q.L.; Chen, Y.X.; Jilani, G.; Shamsi, I.H.; Yu, Q.G. Model avswat apropos of simulating non-point source pollution in Taihu lake basin. J. Hazard. Mater. 2010, 174, 824. [Google Scholar] [CrossRef] [PubMed]
  44. Arnold, J.G.; Moriasi, D.N.; Gassman, P.W.; Abbaspour, K.C.; White, M.J.; Srinivasan, R.; Santhi, C.; Harmel, R.D.; Van Griensven, A.; Van Liew, M.W.; et al. SWAT: Model use, calibration, and validation. Trans. ASABE 2012, 55, 1345–1352. [Google Scholar] [CrossRef]
  45. Luo, Y.; Zhang, X.; Liu, X.; Ficklin, D.; Zhang, M. Dynamic modeling of organophosphate pesticide load in surface water in the northern San Joaquin valley watershed of California. Environ. Pollut. 2008, 156, 1171–1181. [Google Scholar] [CrossRef] [PubMed]
  46. Neitsch, S.L.; Arnold, J.G.; Kiniry, J.R.; Williams, J.R. Soil & Water Assessment Tool; Theoretical Documentation Version 2009; Texas Water Resources Institute Technical Report No. 406; Texas A&M University System: College Station, TX, USA, 2011. [Google Scholar]
  47. Zhao, Y.; Liu, X.; Zhong, X.; Cao, J.; Ren, T.; Feng, D. Determination of soil depth and its water threshold for diagnosing water deficit of winter wheat based on grain yield. Trans. Chin. Soc. Agric. Eng. 2014, 30, 147–154. (In Chinese) [Google Scholar]
  48. Xue, C.; Liu, R.; Ma, Z. Drought grade classification of summer maize in Huang-Huai-Hai area. Trans. Chin. Soc. Agric. Eng. 2014, 30, 147–156. (In Chinese) [Google Scholar]
  49. Sun, C.; Ren, L. Assessment of surface water resources and evapotranspiration in the Haihe river basin of China using swat model. Hydrol. Process. 2013, 27, 1200–1222. [Google Scholar] [CrossRef]
  50. Abbaspour, K.C.; Johnson, C.A.; Genuchten, M.T.V. Estimating uncertain flow and transport parameters using a sequential uncertainty fitting procedure. Vadose Zone J. 2004, 3, 1340–1352. [Google Scholar] [CrossRef]
  51. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual models-part1—A discussion of principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  52. Moriasi, D.N.; Arnold, J.G.; Van Liew, M.W.; Bingner, R.L.; Harmel, R.D.; Veith, T.L. Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Trans. ASABE 2007, 50, 885–900. [Google Scholar] [CrossRef]
  53. Abbaspour, K.C.; Rouholahnejad, E.; Vaghefi, S.; Srinivasan, R.; Yang, H.; Kløve, B. A continental-scale hydrology and water quality model for Europe: Calibration and uncertainty of a high-resolution large-scale SWAT model. J. Hydrol. 2015, 524, 733–752. [Google Scholar] [CrossRef]
  54. Wu, B.; Yan, N.; Xiong, J.; Bastiaanssen, W.G.M.; Zhu, W.; Stein, A. Validation of ETWatch using field measurements at diverse landscapes: A case study in Hai Basin of China. J. Hydrol. 2012, 436–437, 67–80. [Google Scholar] [CrossRef]
  55. Xiong, J.; Wu, B.; Yan, N.; Zeng, Y.; Liu, S. Estimation and validation of land surface evaporation using remote sensing and meteorological data in north China. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2010, 3, 337–344. [Google Scholar] [CrossRef]
  56. Zang, C.F.; Liu, J.; Velde, M.V.D.; Kraxner, F. Assessment of spatial and temporal patterns of green and blue water flows in inland river basins in northwest China. Hydrol. Earth Syst. Sci. 2012, 16, 2859–2870. [Google Scholar] [CrossRef] [Green Version]
  57. Zang, C.F.; Liu, J.G. Trend analysis for the flows of green and blue water in the Heihe River basin, northwestern China. J. Hydrol. 2013, 502, 27–36. [Google Scholar] [CrossRef]
  58. Rouholahnejad, E.; Abbaspour, K.C.; Srinivasan, R.; Bacu, V.; Lehmann, A. Water resources of the Black Sea Basin at high spatial and temporal resolution. Water Resour. Res. 2014, 50, 5866–5885. [Google Scholar] [CrossRef] [Green Version]
  59. Moiwo, J.P.; Yang, Y.H.; Han, S.M.; Lu, W.X.; Yan, N.N.; Wu, B.F. A method for estimating soil moisture storage in regions under water stress and storage depletion: A case study of Hai River Basin, North China. Hydrol. Process. 2011, 25, 2275–2287. [Google Scholar] [CrossRef]
  60. Liu, J.G.; Yang, H. Spatially explicit assessment of global consumptive water uses in cropland: Green and blue water. J. Hydrol. 2010, 384, 187–197. [Google Scholar] [CrossRef]
  61. Liu, S.M.; Xu, Z.W.; Zhu, Z.L.; Jia, Z.Z.; Zhu, M.J. Measurements of evapotranspiration from eddy-covariance systems and large aperture scintillometers in the Hai river basin, China. J. Hydrol. 2013, 487, 24–38. [Google Scholar] [CrossRef]
Figure 1. The location of the Hai River basin (HRB) and its hydrological regions.
Figure 1. The location of the Hai River basin (HRB) and its hydrological regions.
Water 10 00798 g001
Figure 2. The distribution of the river net, sub-basins, and the locations of meteorological, precipitations and hydrologic stations which were used in the model simulations in HRB.
Figure 2. The distribution of the river net, sub-basins, and the locations of meteorological, precipitations and hydrologic stations which were used in the model simulations in HRB.
Water 10 00798 g002
Figure 3. The spatial distribution of soil types (a) and land use types (b) in HRB.
Figure 3. The spatial distribution of soil types (a) and land use types (b) in HRB.
Water 10 00798 g003
Figure 4. The comparison of the observed and simulated runoff at the Zhangjiakou station (a) and the Dongyanghe station (b) (1995–2004).
Figure 4. The comparison of the observed and simulated runoff at the Zhangjiakou station (a) and the Dongyanghe station (b) (1995–2004).
Water 10 00798 g004
Figure 5. The annual precipitation, blue water, green water, and soil water variation in HRB (1995–2004) (unit: mm).
Figure 5. The annual precipitation, blue water, green water, and soil water variation in HRB (1995–2004) (unit: mm).
Water 10 00798 g005
Figure 6. The spatial distribution of the maximum possible storage of green water (MSGW) in HRB.
Figure 6. The spatial distribution of the maximum possible storage of green water (MSGW) in HRB.
Water 10 00798 g006
Figure 7. The components of green water and blue water in HRB (1995–2004) (unit: mm).
Figure 7. The components of green water and blue water in HRB (1995–2004) (unit: mm).
Water 10 00798 g007
Figure 8. The spatial distribution of the consumed green water (a) and utilizable green water (b) in HRB.
Figure 8. The spatial distribution of the consumed green water (a) and utilizable green water (b) in HRB.
Water 10 00798 g008
Table 1. The parameters used in the model calibration and validation.
Table 1. The parameters used in the model calibration and validation.
ParameterDescriptionMaximum Theoretical Range
ALPHA_BFBase flow alpha factor (day)0–1
BLAIMaximum potential leaf area index0–12
CANMXMaximum canopy index (mm H2O)0–100
CH_K2Hydraulic conductivity in main channel alluvium (mm/h)0–150
CN2The SCS (soil conservation service) runoff curve number30–95
EPCOPlant evaporation compensation factor0.01–1
ESCOSoil evaporation compensation factor0.01–1
GW_DELAYGroundwater delay (day)0–500
GWQMNThreshold depth of water in the shallow aquifer required for return flow to occur (mm)0–5000
RCHRG_DPGroundwater recharge to the deep aquifer0–1.0
REVAPMNThreshold depth of water in the shallow aquifer for “revap” or percolation to the deep aquifer to occur (mm)0–500
SOL_AWCAvailable soil water content (mm/mm)0–0.5
Table 2. The value range of the sensitive parameters in the mountainous area after model calibration.
Table 2. The value range of the sensitive parameters in the mountainous area after model calibration.
ParameterFinal Range of Parameters After Model Calibration
CLRBCBRBEDRBWDRBCDRBCZYRBCZWRB
ALPHA_BF0.05–0.100.10–0.160.10–0.140.10–0.150.06–0.100.08–0.150.10–0.15
BLAI *−0.20 to −0.15−0.05 to 0.05−0.15 to −0.05−0.15 to −0.10−0.15 to −0.10−0.05 to 0.05−0.10 to −0.05
CANMX1.2–5.71.2–5.81.3–6.00.8–4.21.5–5.50.9–4.51.4–5.5
CH_K285–10070–9565–8570–9065–8570–8565–85
CN2 *−0.20 to −0.15−0.10 to 0.05−0.25 to −0.15−0.20 to −0.15−0.15 to −0.10−0.25 to −0.20−0.15 to −0.10
EPCO0.30–0.400.20–0.350.20–0.300.40–0.450.35–0.450.30–0.400.25–0.35
ESCO0.85–0.900.80–0.850.75–0.850.85–0.900.80–0.850.75–0.850.70–0.80
GW_DELAY30–4535–4530–4031–3535–3940–4542–44
GWQMN40–9050–10520–7560–9045–8555–7560–75
RCHRG_DP0.45–0.550.40–0.500.25–0.350.30–0.400.35–0.400.40–0.450.40–0.50
REVAPMN220–270350–380360–380290–320330–370360–390320–350
SOL_AWC *0.10–0.15−0.05 to 0.05−0.15 to −0.050.05–0.10−0.10 to −0.050.05 to 0.10−0.05 to 0.05
* Means a relative change to the initial parameter, the other parameters were replaced by a new value. The locations and full name of the hydrological regions in mountainous area (CLRB, CBRB, EDRB, WDRB, CDRB, CZYRB and CZWRB) have been shown in Figure 1.
Table 3. The comparison between the observed and simulated monthly discharge of the representative stations in the mountainous hydrological regions during 1995–2004.
Table 3. The comparison between the observed and simulated monthly discharge of the representative stations in the mountainous hydrological regions during 1995–2004.
Hydrological RegionStationCalibration (1995–1999)Validation (2000–2004)
NSR2NSR2
CLRBLuanxian0.710.880.610.82
CBRBZhangjiafen0.740.840.680.72
Xiahui0.680.830.630.78
EDRBDongyanghe0.700.920.650.89
Zhangjiakou0.690.960.650.90
WDRBCetian0.740.820.700.77
CDRBZhangfang0.700.860.640.72
CZYRBXiaojue0.790.890.700.82
CZWRBGuantai0.770.900.710.80
Table 4. The value range of the sensitivity parameters in the plain area after model calibration.
Table 4. The value range of the sensitivity parameters in the plain area after model calibration.
ParameterFinal Range of Parameters After Model Calibration
LRPTRPLRPBREDRP and WDRPHYP and ZYRP and ZWRP
ALPHA_BF0.20–0.240.22–0.250.25–0.300.25–0.300.20–0.25
BLAI *−0.05 to 0.15−0.05 to 0.050.05–0.10−0.05 to 0.100.05–0.15
CANMX0.8–6.02.0–5.61.2–6.30.5–3.52.5–4.6
CH_K2105–120100–120115–140100–115110–125
CN2 *−0.30 to −0.2−0.25 to −0.20−0.30 to −0.15−0.40 to −0.25−0.35 to −0.30
EPCO0.50–0.750.55–0.650.65–0.850.45–0.600.50–0.65
ESCO0.75–0.850.70–0.800.75–0.850.60–0.750.75–0.80
GW_DELAY25–3522–2525–3520–3020–25
GWQMN60–11070–9070–10045–8575–90
RCHRG_DP0.35–0.450.35–0.450.45–0.500.40–0.500.40–0.45
REVAPMN190–220160–180130–170160–190140–170
SOL_AWC *0.05–0.150.10–0.150.10–0.150.15–0.250.05–0.10
* Means a relative change to the initial parameter, the other parameters were replaced by a new value. The locations and full name of the hydrological regions in plain area (LRP, TRP, LRPBR, EDRP, WDRP, HYP, ZYRP and ZWRP) have been shown in Figure 1.
Table 5. The simulated and observed annual river discharge to the sea in the coastal hydrological regions (km3).
Table 5. The simulated and observed annual river discharge to the sea in the coastal hydrological regions (km3).
Hydrological RegionData TypeYear
19961997199819992000
LRPSimulated0.371.420.290.39
Observed a0.441.650.210.10
Relative error−0.15−0.140.412.95
TRPSimulated2.190.641.800.170.36
Observed a2.510.471.800.110.22
Relative error−0.130.370.000.590.65
LRPBRSimulated0.230.120.490.370.43
Observed a0.160.110.300.610.46
Relative error0.390.140.66−0.40−0.07
EDRP and WDRP *Simulated0.381.890.110.27
Observed a0.431.810.260.74
Relative error−0.120.04−0.55−0.64
HYP and ZYRP * and ZWRP *Simulated0.703.880.330.84
Observed a0.873.630.511.47
Relative error−0.190.07−0.36−0.43
A region name with a * is an inland hydrological region and it belongs to the same river system with the relative coastal region. The calibration results for those sub-basins in the same hydrological region group (such as EDRP and WDRP *) are consistent. For coastal hydrological regions, their discharges are summarized if there are multiple outlets to the sea. a The data was obtained from the Water Resources Bulletin of the Hai River Basin [30].
Table 6. The comparison between the monthly remote sensing ET and simulated ET in the plain hydrological regions during 2002–2004.
Table 6. The comparison between the monthly remote sensing ET and simulated ET in the plain hydrological regions during 2002–2004.
Hydrological RegionArea (km2)Average Observed ET (mm)Average Simulated ET (mm)NSR2
LRP10,84343.3338.480.710.74
TRP31,91945.1742.170.660.80
LRPBR15,46236.5934.760.670.78
EDRP and WDRP26,26355.2849.580.600.72
HYP and ZYRP and ZWRP46,86846.3646.040.660.82
Table 7. The evapotranspiration (ET) at different land use types in HRB (km3).
Table 7. The evapotranspiration (ET) at different land use types in HRB (km3).
TimeCGWNGWTotal ETRainfall
FarmlandGrasslandWoodlandUrban AreaWater Body
199558.5021.9424.316.392.09113.22188.70
199658.1122.6824.105.661.76112.30183.93
199753.7020.8522.654.912.13104.23115.31
199863.6025.6627.036.271.99124.54174.88
199954.7620.8023.714.982.05106.31119.72
200052.0121.1322.405.042.23102.81154.05
200156.4020.3922.865.442.32107.41131.61
200254.9523.4224.044.922.30109.64123.32
200364.4226.7326.736.102.13126.11184.73
200475.3827.6127.737.413.22141.33171.43
Average59.1823.1224.565.712.22114.79154.77
Percentage (%)51.5620.1421.394.971.94100
Table 8. The effective utilization rate of precipitation in different land use types.
Table 8. The effective utilization rate of precipitation in different land use types.
YearFarmlandGrasslandWoodland
RCGWPRCGWPRCGWP
1995585.17350.3859.88586.76341.9258.29606.92352.5558.08
1996556.54348.0562.54590.81353.4659.83606.63349.5057.62
1997346.97321.6392.70377.61324.9486.03387.07328.4884.85
1998543.30380.9370.11510.39399.9078.35587.78392.0066.70
1999363.50327.9890.23393.04324.1682.48390.98343.8587.96
2000489.34311.5163.66490.76329.3067.08452.91324.8571.73
2001413.75337.8181.64379.79317.7783.67437.53331.5275.77
2002362.36329.1290.83451.64364.9980.82393.30348.6388.62
2003598.65385.8464.45582.24416.5771.55522.08387.6574.27
2004536.90451.4984.09547.48430.2978.58523.53402.1576.82
Average479.65354.4773.90491.05360.3373.38490.87356.1272.55
R = Rainfall (km3), CGW = consumed green water (km3), and P = CGW/R (%).
Table 9. The monthly variation of green water resources in 2000 (unit: mm).
Table 9. The monthly variation of green water resources in 2000 (unit: mm).
MonthPrecipitationETBlue WaterSoil Water ContentCGWUGW
18.932.140.4194.041.8714.15
211.299.781.7092.368.9412.47
35.7314.410.4883.0513.533.14
411.9718.150.2377.2716.81−2.66
530.8133.921.0473.8031.89−6.12
640.4746.782.5165.8144.41−14.14
7148.6964.5249.3797.6359.3217.75
8121.9251.5737.05131.3948.2051.57
951.5437.3521.13128.9834.8449.16
1045.4026.7314.02134.1524.9054.34
1111.3813.856.23125.1312.7145.30
121.587.953.03117.297.2537.45
Total489.72327.15137.21 304.66

Share and Cite

MDPI and ACS Style

Zhu, K.; Xie, Z.; Zhao, Y.; Lu, F.; Song, X.; Li, L.; Song, X. The Assessment of Green Water Based on the SWAT Model: A Case Study in the Hai River Basin, China. Water 2018, 10, 798. https://doi.org/10.3390/w10060798

AMA Style

Zhu K, Xie Z, Zhao Y, Lu F, Song X, Li L, Song X. The Assessment of Green Water Based on the SWAT Model: A Case Study in the Hai River Basin, China. Water. 2018; 10(6):798. https://doi.org/10.3390/w10060798

Chicago/Turabian Style

Zhu, Kui, Zibo Xie, Yong Zhao, Fan Lu, Xinyi Song, Lu Li, and Xiaomeng Song. 2018. "The Assessment of Green Water Based on the SWAT Model: A Case Study in the Hai River Basin, China" Water 10, no. 6: 798. https://doi.org/10.3390/w10060798

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