Next Article in Journal
Propagation of Native Tree Species to Restore Subtropical Evergreen Broad-Leaved Forests in SW China
Next Article in Special Issue
Influences of Land Use Change on Baseflow in Mountainous Watersheds
Previous Article in Journal
Foundation Species Loss and Biodiversity of the Herbaceous Layer in New England Forests
Previous Article in Special Issue
Streamflow Regime Variations Following Ecological Management on the Loess Plateau, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling Ecohydrological Processes and Spatial Patterns in the Upper Heihe Basin in China

1
School of Water Resources and Environment, China University of Geosciences, Beijing 100083, China
2
State Key Laboratory of Hydroscience and Engineering, Department of Hydraulic Engineering, Tsinghua University, Beijing 100084, China
3
State Key Laboratory of Vegetation and Environmental Change, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China
*
Author to whom correspondence should be addressed.
Forests 2016, 7(1), 10; https://doi.org/10.3390/f7010010
Submission received: 29 September 2015 / Revised: 18 December 2015 / Accepted: 21 December 2015 / Published: 25 December 2015
(This article belongs to the Special Issue Forest Management and Water Resources in the Anthropocene)

Abstract

:
The Heihe River is the second largest inland basin in China; runoff in the upper reach greatly affects the socio-economic development in the downstream area. The relationship between spatial vegetation patterns and catchment hydrological processes in the upper Heihe basin has remained unclear to date. In this study, a distributed ecohydrological model is developed to simulate the hydrological processes with vegetation dynamics in the upper Heihe basin. The model is validated by hydrological observations at three locations and soil moisture observations at a watershed scale. Based on the simulated results, the basin water balance characteristics and their relationship with the vegetation patterns are analyzed. The mean annual precipitation and runoff increase with the elevation in a similar pattern. Spatial patterns of the actual evapotranspiration is mainly controlled by the precipitation and air temperature. At the same time, vegetation distribution enhances the spatial variability of the actual evapotranspiration. The highest actual evapotranspiration is around elevations of 3000–3600 m, where shrub and alpine meadow are the two dominant vegetation types. The results show the mutual interaction between vegetation dynamics and hydrological processes. Alpine sparse vegetation and alpine meadow dominate the high-altitude regions, which contribute most to the river runoff, and forests and shrub contribute relatively small amounts of water yield.

1. Introduction

Vegetation, topography and climate variables interactively influence multi-scale ecohydrological processes in a large catchment. Vegetation patterns are a key factor affecting the water balance and catchment water yield. Previous studies have demonstrated the impact of forest change on river discharge in China [1,2]. However, few studies have addressed the relationship between hydrological processes and vegetation patterns. In particular, the relationship between forest and water in the headwater catchments of arid inland basins remains unclear [3]. In arid inland basins, the river discharge generated from the mountain regions greatly affects the socio-economic development and ecosystem sustainability of the downstream regions. Therefore, understanding the complex relationship between spatial vegetation patterns and hydrological processes is important for integrated river basin management. Hydrological observations and small-scale experiments provide limited information for understanding the spatial patterns of ecohydrological processes in a large river basin [4]. Instead, a physically based distributed model that links ecohydrological processes across scales is needed to analyze the spatial variability of hydrology with vegetation pattern.
In the past few decades, with global climate change and population growth, the water shortage and ecosystem degradation in many river basins have gained increasing concern worldwide. Meanwhile, changes in natural river runoff in the headwater catchments worldwide have been reported in previous studies [5,6,7,8,9,10,11]. Because of the less direct influence of human activities in headwater areas, runoff changes are mainly caused by the mutual interactions among climate, vegetation and hydrology. Yang et al. [12] reported that climate contributions to river runoff showed a large spatial variation over China, and several studies attempted to attribute runoff change to the impacts of climate and vegetation changes [13,14,15]. However, vegetation is absent or simply parameterized in traditional hydrological models [16]. It is relatively difficult to evaluate the influence of vegetation change on runoff due to data availability and methodological limitations [17,18]. Therefore, it is important to develop the ecohydrological models for understanding and predicting changes in the regional water availability under the changing environment.
Modeling hydrological processes at the catchment scale requires a flexible distributed scheme to represent the catchment topography, river network and vegetation patterns [19]. Previous studies of distributed hydrological models focused on the representation of the heterogeneity of catchment landscape and physically based descriptions of hydrological processes, such as the Systeme Hydrologique Europeen (SHE) and the Variable Infiltration Capacity (VIC) model [20,21]. Yang et al. [22,23] proposed the geomorphology-based hydrological model (GBHM) considering sub-grid parameterization, which has been employed in macroscale studies [24,25] and mesoscale studies [13,14,26,27]. To better simulate the hydrological changes in a catchment, it is necessary to couple the vegetation pattern and vegetation dynamics in hydrological models [19]. Land surface models, such as the second version of the Simple Biosphere model (SiB2) [28], have mainly focused on the role of vegetation in the water-heat-carbon cycle. However, the catchment hydrology was poorly represented in SiB2. Further research embedded the SiB2 into the GBHM and developed a distributed biosphere hydrological model called the Water and Energy Budget-Based Distributed Hydrological Model (WEB-DHM) [29]. A comparative analysis of the WEB-DHM and GBHM has been carried out and indicated that the WEB-DHM could be used to predict streamflow and conduct the water and energy flux estimations [30]. Several land surface models can simulate the cryosphere processes, such as the Community Land Model (CLM) [31]. However, they typically simplify the catchment hydrological processes and cannot accurately simulate streamflow. Several studies attempt to develop physically based models to couple hydrological processes with ecological processes. Tague and Band [32] developed the Regional Hydro-Ecologic Simulation System (RHESSys) to simulate both hydrological processes and vegetation dynamics. However, in RHESSys, the frozen soil process is not considered and the soil temperature is estimated based on the average daily temperature. Maneta and Silverman [33] and Ivanov et al. [34] developed spatial distributed models to simulate the water, energy balance and vegetation dynamics. However, the frozen soil is not adequately considered, as there are only two soil thermal layers in this model. Therefore, it is important to develop a distributed model that couples hydrological processes, cryosphere processes and vegetation dynamics.
There are several previous studies focusing on the hydrological modeling in the Heihe River basin. Jia et al. [35] applied the water and energy transfer process (WEP) model to simulate and predict the annual runoff changes due to climate and land use changes. Wang et al. [36] and Zhang et al. [37] applied heat and water coupling models to a small experimental catchment in the upper Heihe basin to test the model applicability in the high, cold mountainous area. Zhou et al. [38] chose several modules for the hydrological processes in cold regions and linked them to a catchment hydrological model to improve the hydrological modeling capability in the cold mountainous area. Zang and Liu [39] applied the Soil and Water Assessment Tool (SWAT) in the Heihe basin to analyze the flow trends of green and blue water. Yang et al. [19] developed a distributed scheme for ecohydrological modeling in the upper Heihe River. However, this scheme needs further coupling of vegetation dynamics into the hydrological simulation and improvement of the cryosphere hydrological processes. Because of the complexity of ecohydrological processes in the upper Heihe basin, further research to develop a sophisticated ecohydrological model is still ongoing. A major research plan entitled “Integrated research on the ecohydrological processes of the Heihe basin” has been launched by the National Natural Science Foundation of China since 2010 [40]. One of the integrated research projects in this research plan aims to develop a distributed ecohydrological model for the cold mountainous regions of inland basins. With the help of the Heihe Watershed Allied Telemetry Experimental Research (HiWATER), which is a comprehensive ecohydrological experiment in this research plan [4], there are new opportunities for ecohydrological model development in the Heihe basin.
The major objectives of this study are to (1) develop a distributed ecohydrological model by coupling hydrological processes with vegetation dynamics and (2) analyze the relationship between the water balance and vegetation patterns in the upper Heihe basin.
In the following sections, the features of the study area are presented, followed by the data descriptions. The distributed ecohydrological model is then introduced. In the results and discussion section, the model validation is presented, followed by the analysis of water balance characteristics and the spatiotemporal variability of runoff. Finally, the catchment ecohydrological pattern is analyzed from the water balance and its relation to the vegetation distribution.

2. Study Area and Data Used

2.1. The Upper Heihe Basin

Heihe basin is the second largest inland basin in Northwest China. The upper reach of the Heihe River, which is gauged at the Yingluoxia hydrological station (see Figure 1), has a drainage area of 10,005 km2. The upper Heihe basin generates nearly 70% of the total river runoff, which supplies agricultural irrigation and benefits the social economy development in the middle and lower basin [19,41]. There are two tributaries in the upper Heihe basin, namely, the West Tributary and East Tributary, both of which originate in the Qilian Mountains. The annual precipitation is between 200 mm and 700 mm with highly seasonal variability, and nearly 60% of the total annual rainfall is concentrated in summer from June to September. The upper Heihe basin has an elevation of 1700–5200 m, low temperature and relatively abundant vegetation types [41,42,43,44]. The major vegetation types include coniferous forest (Picea crassifolia Kom.), shrub (Potentilla fruticosa Linn.), steppe (Stipa purpurea Griseb), alpine meadow (Kobresia pygmaea Clarke), alpine sparse vegetation (Saussurea medusa Maxim.), and desert (Sympegma regelii Bunge) (see Figure 2).

2.2. Data Used in the Study

Two types of data are used in this study: the first type is the data used to build and run the ecohydrological model, which include the geographic information and climatic forcing data, and the second type is the data used for model validation.
Meteorological observations are available at several stations within the study basin and its surrounding areas, as shown in Figure 1. The data are acquired from the National Meteorological Information Center affiliated with the China Meteorological Administration [45]. The observed meteorological data include daily precipitation, temperature, sunshine hour, wind speed and relative humidity. In this study, daily gridded precipitation is interpolated from the gauge data using the method developed by Shen and Xiong [46]. Other forcing data are interpolated using the inverse distance method. The hourly temperature is estimated from the daily maximum and minimum temperature using a sine curve. Hourly precipitation is estimated from the daily data according to the duration. Duration of precipitation within a day is determined by the precipitation amount estimated from the regional historical records. The starting time of the precipitation in a day is decided randomly and the hourly precipitation is specified using a normal distribution. The wind speed and relative humidity are assumed uniform in a day. The location of glaciers and their areas and volumes are obtained from the first and second glacier inventory datasets of China [47,48,49], which is downloaded from the Cold and Arid Regions Science Data Center at Lanzhou [50]. Digital elevation data with a resolution of 90 m resolution are downloaded from the Shuttle Radar Topography Mission (SRTM) database [51]. The vegetation map (see Figure 2) is obtained from the Institute of Botany, Chinese Academy of Sciences [52,53]. The leaf area index (LAI) is estimated from the remote sensing data by Fan [54]. The soil map of Heihe basin is produced by the second national soil survey [55]. We use the van Genuchten model to represent the retention curve of soil water. The soil water parameters used in this study, including the saturated hydraulic conductivity, residual soil moisture, saturated soil moisture, and parameters α and n in the van Genuchten model, are provided by the China soil dataset developed by Dai et al. [56]. The spatial resolution of this dataset is 1 km.
Figure 1. Location and topography of the upper Heihe basin.
Figure 1. Location and topography of the upper Heihe basin.
Forests 07 00010 g001
To validate the runoff simulation, the stream discharge data are obtained from the Hydrology and Water Resources Bureau of Gansu Province. Soil moisture measured by the wireless sensor network (WSN) in the East Tributary is obtained from the HiWATER [57,58]. The actual evapotranspiration dataset estimated by remote sensing is available in the study basin at a 1 km spatial resolution and monthly temporal scale [59,60], which is used for comparison with the simulated result.
Figure 2. Vegetation map of the upper Heihe basin.
Figure 2. Vegetation map of the upper Heihe basin.
Forests 07 00010 g002

3. Distributed Ecohydrological Model

3.1. Representation of the Landscape

The distributed scheme used to develop the ecohydrological model is originally from the geomorphology-based hydrological model (GBHM) [61,62,63,64] and was improved by Yang et al. [19] for ecohydrological modeling. The distributed ecohydrological model developed in this study is called the geomorphology-based ecohydrological model (hereafter referred to as the GBEHM). The major development in the GBEHM includes replacing the evapotranspiration estimation from the Penman-Monteith equation to a simple biosphere model used in SiB2 [65] and adding an energy balance based module for simulating the cryosphere hydrological processes.
A grid system with a resolution of 1 km × 1 km is used to discretize the study catchment, and the river network is extracted from the 1-km digital elevation model (DEM), by which the study catchment is divided into 461 sub-catchments. The main rivers of 461 sub-catchments are used to represent the streamflow pathway of the study catchment, and the Horton-Strahler ordering system is used to define the flow routing sequence. The mean terrain properties (slope length and gradient) and a particular soil type are assigned to each 1-km grid and the area ratios of vegetation types are determined for each grid.
Furthermore, each 1-km grid cell is represented by a number of topographically similar “hillslope-valley” systems. Length, gradient and aspect of the hillslope are estimated from the 90-m DEM and averaged on the 1-km grid [19]. The hillslopes within a 1-km grid are grouped according to vegetation type. The vertical structure of a hillslope is then subdivided into vegetation canopy, soil and bedrock.

3.2. Simulation of Ecohydrological Processes

Hillslope is the basic unit of ecohydrological simulation in the GBEHM, in which vegetation dynamics are coupled with the hillslope hydrological processes. The descriptions of hydrological processes, which mainly include the transfer of energy, water and carbon dioxide in the soil-plant-atmosphere continuum, are especially designed for better simulation in the cryosphere.
During the growing season, the photosynthesis process of vegetation is simulated together with canopy energy transfer and canopy evapotranspiration. The radiation transfer in the vegetation canopy layer includes interception, reflection, transmission and absorption and is described by the same scheme in SiB2 [28,66]. The energy balance equation in the canopy layer is expressed as [28]
C T c t = R n H λ E ξ
where C is the effective heat capacity of the canopy (J·m−2·K−1), Tc is the canopy temperature (K), t is the time (s), Rn is the absorbed net radiation (W·m−2), H is the sensible heat flux (W·m−2), λ is the latent heat of vaporization (J·kg−1), E is the evapotranspiration rate (kg·m−2·s−1), and ξ is the energy transfer due to water phase change (W·m−2), which is caused by the melting/freezing of the snow intercepted by the canopy. Thus, the energy budget is linked with water balance by the evapotranspiration component. The absorbed net radiation is calculated using the radiation transfer model described by Sellers et al. [65]. The sensible heat flux of the canopy H is calculated as
H = ρ c p T c T a r b
where ρ is the air density (kg·m−3), cp is the specific heat of air (J·K−1·kg−1), Tc is the canopy temperature (K), Ta is the air temperature (K), and rb is the bulk canopy boundary layer resistance (s·m1).
The actual evapotranspiration rate of the vegetation layer is influenced by the leaf stomatal conductance, which is related to environmental properties, such as the canopy temperature, carbon dioxide concentration and water vapor in the air. In this study, canopy conductance is estimated as [28]
g c = m A c C s h s p + b L T
where gc is the canopy conductance (m·s−1), m and b are the empirical coefficients of vegetation types, Cs is the carbon dioxide partial pressure at the leaf surface (Pa), hs is the relative humidity, p is the atmospheric pressure (Pa), LT is the total leaf area index (m2·m−2), and Ac is the canopy photosynthesis rate (mol·m−2·s−1) and is calculated as [67]
A c = A n 0 Π
Π = FPAR / k
where An0 is the net assimilation rate for the leaves at the top of the canopy (mol·m−2·s−1), FPAR is the fractional interception of photosynthetically active radiation, and k is the time-mean extinction coefficient for photosynthetically active radiation. The net assimilation rate of the leaf is calculated as
A n = A R d
where An is the leaf net assimilation rate (mol·m−2·s−1), A is the leaf photosynthetic rate (mol·m−2·s−1), and Rd is the leaf respiration rate (mol·m−2·s−1). A and Rd are estimated by the functions of the maximum catalytic capacity of the photosynthetic enzyme (Vmax), canopy temperature and other environmental factors [28,68,69]. Details about the calculation of gc, A and Rd were given by Sellers et al. [28]. Based on the canopy conductance gc, the canopy transpiration rate is expressed as [28]
λ E c = [ e * ( T c ) e a 1 / g c + 2 r b ] ρ c p γ ( 1 W c )
where Ec is the canopy transpiration rate (kg·m−2·s−1), e* is the saturation vapor pressure (Pa), ea is the vapor pressure in the canopy (Pa), rb is the canopy boundary layer resistance (s·m−1), ρ is the density of air (kg·m−3), cp is the specific heat of air (J·kg−1·K−1), and γ is the psychrometric constant (Pa·K−1), Wc is the canopy wetness-snow cover fraction, and the other parameters are the same as in Equations (1) and (3). The soil evaporation rate Eg of the surface soil layer and the evaporation rate of the canopy interception Ei are calculated using the same methods in SiB2 [28]. Eg and Ei are added to Ec as the total actual evapotranspiration over the hillslope unit.
In the non-growing season, ecohydrological simulation focuses on the cryosphere hydrological processes, which mainly include the snow melting process, soil freezing and thawing process, and glacier melting process. We use a similar scheme as in the Community Land Model version 4.0 (i.e., CLM4.0) [31] to represent the heat transfer in snow and frozen soil:
c T s t = z [ K T T s z ]
where c is the volumetric snow/soil heat capacity (J·m−3·K−1), Ts is the temperature (K) of the snow/soil layers, z is the vertical depth of snow/soil (m), and KT is the thermal conductivity (W·m−1·K−1). Equation (8) solves the snow/soil temperature with the boundary condition as the heat flux into the top surface layer of the snow/soil and zero heat flux at the bottom of the soil column. The surface layer heat flux from the atmosphere is expressed as
h = S g L g H g λ E g
where h is the upper boundary heat flux into the snow/soil layer (W·m−2), Sg is the solar radiation absorbed by the top layer (W·m−2), Lg is the long-wave radiation absorbed by the ground (W·m−2), Hg is the sensible heat flux from the ground (W·m−2), and λEg is the latent heat flux from the ground (W·m−2).
After solving Equation (8), the soil or snow layer temperature is evaluated to determine whether the phase change will take place [31]. For the soil or snow layers, melting takes place under the condition of
Ti > Tf and Wice > 0.
For the snow layers, freezing takes place under the condition of
Ti < Tf and Wliq > 0
and for the soil layers, freezing takes place under the condition of
Ti < Tf and Wliq > Wliq,max.
In Equations (10)–(12), Ti is the temperature of the soil or snow layers (K), Tf is the freezing temperature of water (K), Wice is the mass of ice in the soil or snow layers (kg·m−2), Wliq is the mass of liquid water in soil or snow layers (kg·m−2), and Wliq,max is the maximum mass of supercooled soil water, which is the liquid water that coexists with the ice (kg·m−2). The phase change rate is determined by the energy excess (or deficit) needed to change the temperature of soil or snow layers (Ti) to the freezing temperature (Tf). If the melting criteria in Equation (10) is met and energy excess is greater than zero, then the ice mass is calculated by [31]
W ice n + 1 = W ice n U i Δ t L f
where W ice n + 1 is the ice mass after the phase change (kg·m−2), W ice n is the ice mass before the phase change (kg·m−2), n + 1 and n refer to the time steps, Ui is the energy excess (W·m−2), Δt is the length of time step (s), and Lf is the latent heat of ice fusion (J·kg−1). If the freezing criteria in Equations (11) or (12) is met and energy excess is less than zero, then the ice mass in the snow layers is adjusted by
W ice n + 1 = min ( W ice n + W liq, n W ice n U i Δ t L f ) .
The ice mass in the soil layers is adjusted by
W ice n + 1 = { min ( W ice n + W liq n W liq,max , W ice n U i Δ t L f ) W ice n + W liq n W liq,max 0 W ice n + W liq n < W liq,max
and the mass of liquid water is adjusted by
W liq n + 1 = W liq n + W ice n W ice n + 1
where W liq n + 1 and W liq n are the mass of liquid water after and before the phase change (kg·m−2), respectively, in Equations (14)–(16), and the other parameters are the same as in Equation (13). Moreover, glacier melting is simulated using the degree-day model [70,71].
The surface runoff is from the infiltration excess and saturation excess calculated by solving Richards’ equation using an implicit finite difference method. The surface runoff flows through the hillslope into the stream via kinematic wave. The groundwater aquifer is treated as an individual storage corresponding to each grid. The exchange between the groundwater and river water is considered as steady flow and calculated using Darcy’s law [64]. The runoff generated from the grid is the lateral inflow into the river of the sub-catchment. The flow routing in the river network of the whole study catchment is solved using the kinematic wave approach:
{ A t + Q x = q Q = S 0 1 / 2 n r p 2 / 3 A 5 / 3
where Q is the discharge (m3·s−1), t is the time (s), x is distance along the river (m), A is the area of the cross-section (m2), q is the lateral inflow to the river from the hillslope (m3·s−1), S0 is the slope of the river bed, nr is the roughness of the river bed, and p is the wetting perimeter of cross-section (m). Equation (17) is solved using a nonlinear explicit finite difference method and Newton’s iteration scheme. The time step of the GBEHM model is one hour.

3.3. Model Calibration and Performance Evaluation Metrics

Most model parameters are estimated from field observations or remote sensing data, and some vegetation parameters are specified from previous studies [19]. The major parameters related to vegetation type are listed in Table 1. Considering the equilibration time for the hydrological state variables (e.g., soil moisture, soil temperature, groundwater table), a warm-up run of the model in the period of 1999–2000 is used to update the model state variables. We perform our analysis in the period from 2001 to 2012. The Nash-Sutcliffe efficiency (NSE) coefficient and relative error (RE) are used to evaluate the model performance.
Table 1. Major parameters related to vegetation and soil types in the geomorphology-based ecohydrological model (GBEHM).
Table 1. Major parameters related to vegetation and soil types in the geomorphology-based ecohydrological model (GBEHM).
Coniferous ForestShrubSteppeAlpine MeadowAlpine Sparse VegetationDesert
Dominant speciesPicea crassifolia Kom.Potentilla fruticosa Linn.Stipa purpurea GrisebKobresia pygmaea ClarkeSaussurea medusa Maxim.Sympegma regelii Bunge
Root depth (m)2.01.00.40.40.10.0
Surface storage (mm)30.025.010.015.015.05.0
Leaf reflectance to visible light0.1050.1050.1050.1050.105
Leaf reflectance to near-infrared radiation0.350.580.580.580.58
Leaf transmittance to visible light0.050.070.070.070.07
Leaf transmittance to near-infrared radiation0.100.250.250.250.25
Maximum Rubsico capacity of top leaf (10−5 mol·m−2·s−1)6.06.03.33.33.0
Intrinsic quantum efficiency (mol·mol−1)0.080.080.050.050.05
Mean (standard deviation) of saturated hydraulic conductivity (mm·hour−1)18.71 (7.86)20.32 (6.64)18.59 (8.89)24.64 (9.83)25.23 (6.48)11.18 (5.87)
Mean (standard deviation) of the saturated soil moisture content (cm3·cm−3)0.470 (0.010)0.465 (0.009)0.468 (0.011)0.465 (0.014)0.458 (0.011)0.475 (0.033)
Mean (standard deviation) of the residual soil moisture content (cm3·cm−3)0.100 (0.002)0.100 (0.001)0.100 (0.001)0.096 (0.019)0.099 (0.005)0.100 (0.000)
Mean (standard deviation) of n in van Genuchten Model1.195 (0.011)1.205 (0.019)1.197 (0.012)1.205 (0.023)1.215 (0.022)1.187 (0.005)
Mean (standard deviation) of α in van Genuchten Model0.020 (0.003)0.021 (0.003)0.021 (0.003)0.022 (0.004)0.022 (0.004)0.017 (0.003)

4. Results and Discussion

4.1. Model Validation

The model has been validated using the soil moisture measured in the East Tributary of Heihe River, streamflow discharge observed at three hydrological stations (Figure 1) and actual evapotranspiration estimated from the remote sensing observations.
At the catchment scale, the model is validated using the streamflow and spatial distribution of observed soil moisture. Figure 3 illustrates the monthly river discharge of the Yingluoxia, Zhamashike and Qilian hydrological stations, which are located at the outlet of upper Heihe River, the West Tributary and the East Tributary, respectively. The simulated streamflow generally shows a good agreement with the observed one. The NSE values for the three observation stations are 0.77, 0.80 and 0.67, respectively, and the absolute values of RE are smaller than 10%. The distributed soil moisture simulation of the top 5-cm layer is compared with the WSN observation at 4 cm from the surface of the HiWATER in the East Tributary of the upper Heihe basin. A comparison of the simulations and observations for the monthly averaged soil moisture over the East Tributary is shown in Figure 4. The areal average values of the simulated and observed soil moistures are highly similar at 0.41 and 0.40 in August, 0.40 and 0.39 in September, and 0.31 and 0.29 in October, respectively. The simulated soil moisture generally captures the spatial pattern of the observed soil moisture. However, the details need to be checked carefully in future research.
In addition, the simulated actual evapotranspiration (ET) is compared with the remote sensing-based estimation. Figure 5 illustrates the comparison of annual average ET between the model simulation and remote sensing-based estimation in the study catchment in the 2001–2012 period [59,60]. Both of the estimations have close long-term basin average ETs and similar spatial patterns over the study catchment. The areal average ET is 310.8 mm for the model simulation and 306.7 mm for the remote sensing-based estimation. However, the remote sensing-based estimation shows a higher spatial variability than the GBEHM simulation. The remote sensing-based estimation mainly considers the energy balance, and the land surface roughness plays an important role. In addition, the ecohydrological model considers both energy and water balances when estimating the actual evapotranspiration, in which vegetation plays an important role. Vegetation may re-distribute the precipitation and moderate the spatial variability of soil moisture and evapotranspiration.
Figure 3. Comparison of the monthly observed and simulated river discharge at the three hydrological stations: (a) Yingluoxia; (b) Zhamashike; and (c) Qilian, in the 2001–2012 period.
Figure 3. Comparison of the monthly observed and simulated river discharge at the three hydrological stations: (a) Yingluoxia; (b) Zhamashike; and (c) Qilian, in the 2001–2012 period.
Forests 07 00010 g003
Figure 4. Spatial distributions of monthly average soil moisture (m3·m−3) in the East Tributary: (a) observed soil moisture at 4 cm from the surface and (b) simulated soil moisture of the top 5-cm layer.
Figure 4. Spatial distributions of monthly average soil moisture (m3·m−3) in the East Tributary: (a) observed soil moisture at 4 cm from the surface and (b) simulated soil moisture of the top 5-cm layer.
Forests 07 00010 g004
Figure 5. Comparison of the annual average evapotranspiration between the (a) GBEHM simulation and (b) remote sensing-based estimation during the 2001–2012 period [59].
Figure 5. Comparison of the annual average evapotranspiration between the (a) GBEHM simulation and (b) remote sensing-based estimation during the 2001–2012 period [59].
Forests 07 00010 g005

4.2. Water Balance Characteristics and Spatio-Temporal Variability of Runoff

Based on the simulation, the annual average water balance during 2001–2012 is calculated for the East and West Tributaries and the entire catchment, and the results are given in Table 2. For the upper Heihe basin, the annual average precipitation, actual evapotranspiration (ET) and runoff are 479.9 mm, 310.8 mm and 169.0 mm, respectively. Comparing the water balance in the East and West Tributaries, the East Tributary has higher actual ET and runoff due to its higher precipitation. However, the West Tributary has a larger runoff ratio compared to the East Tributary due to the higher altitudes and relatively larger areas of glaciers.
Table 2. Water balance for the East and West Tributaries and the entire catchment during the 2001–2012 period.
Table 2. Water balance for the East and West Tributaries and the entire catchment during the 2001–2012 period.
CatchmentDrainage Area (km2)Precipitation (mm/a)Actual ET (mm/a)Runoff (mm/a)Runoff Ratio
East Tributary2457529.8344.9186.90.35
West Tributary4586485.3304.8178.30.37
Entire catchment10,005479.9310.8169.00.35
The seasonal characteristics of water balance for the East and West Tributaries and the entire catchment are also estimated and are shown in Table 3. In general, the East and West Tributaries and the entire catchment show similar seasonal patterns. For the entire catchment, the precipitation in winter (from December to February) is only 11.5 mm, and the actual evapotranspiration (ET) and runoff is also rather small (6.2 mm and 8.9 mm, respectively). The precipitation in spring (from March to May) is 81.4 mm, which is lower than the sum of the actual ET (66.0 mm) and runoff (19.3 mm). This implies that the runoff in spring is generated from not only precipitation but also snow and glacier melting. This characteristic is more obvious for the West Tributary due to the relatively larger area of glaciers (see Figure 2). As the glacier area is relatively small (less than 1% of the basin area), the contribution of glacier melting to river discharge in summer and autumn is quite low. In summer (from June to August), the precipitation is 281.6 mm, which is larger than the total of the actual ET (180.1 mm) and runoff (78.3 mm). This result implies that precipitation recharges the soil water and groundwater in summer. In autumn (from September to November), the precipitation is 105.4 mm for the entire basin, and it is lower than the total of the actual ET (58.5 mm) and runoff (62.5 mm). This result implies that the runoff is generated from the precipitation and also from the soil water and groundwater storage.
Table 3. Seasonal water balance for the East and West Tributaries and the entire catchment during the 2001–2012 period.
Table 3. Seasonal water balance for the East and West Tributaries and the entire catchment during the 2001–2012 period.
SeasonWater Balance ComponentsEast TributaryWest TributaryWhole Catchment
Spring (Mar.–May)Precipitation (mm/a)98.276.881.4
Actual ET (mm/a)71.164.266.0
Runoff (mm/a)19.018.219.3
Summer (Jun.–Aug.)Precipitation (mm/a)297.6294.3281.6
Actual ET (mm/a)201.8176.5180.1
Runoff (mm/a)83.880.278.3
Autumn (Sep.–Nov.)Precipitation (mm/a)124.8100.3105.4
Actual ET (mm/a)65.656.358.5
Runoff (mm/a)70.863.762.5
Winter (Dec.–Feb.)Precipitation (mm/a)9.213.911.5
Actual ET (mm/a)6.26.36.2
Runoff (mm/a)13.216.68.9

4.3. Spatial Pattern of Water Balance and Relation to Vegetation

To better understand the ecohydrological pattern in the upper Heihe basin, the spatial distributions of the water balance components are further analyzed, including the precipitation, evapotranspiration (ET), runoff and soil moisture in the top one meter during the vegetation growing season from May to October, based on a 30-year model simulation during 1981–2010. Figure 6 shows that annual precipitation and runoff have a similar spatial pattern, annual actual evapotranspiration has a similar pattern with the soil moisture of the top layer in the growing season. The mean soil moisture values are relatively high (0.22–0.41) because of the selected wet season and possible uncertainties of soil water parameters. As shown in Figure 6a, annual precipitation over the study catchment ranges from 220 mm to 630 mm, and the East Tributary has the highest annual precipitation. Comparing Figure 6b with Figure 2, the areas with relatively high ET correspond to the two major vegetation types, namely, steppe and alpine meadow. Figure 6b illustrates that ET has a similar pattern in the East and West Tributaries corresponding with the distribution of alpine meadow vegetation. However, the steep and narrow valley area along the main stream (i.e., downstream of the junction of East and West Tributaries) has relatively low ET and low soil moisture due to the lower precipitation, and this is also related to the steep hillslope and sparse vegetation. In general, the spatial distribution of runoff is mainly controlled by the precipitation and also affected by the topography and vegetation.
The area percentage of each vegetation type is calculated at a 200-m elevation interval in the basin. Accordingly, the water balance components (i.e., annual precipitation, ET and runoff) are also calculated in the same elevation intervals. Figure 7 shows the distributions of vegetation types, annual precipitation, ET and runoff along with the elevation. The vegetation areas of the four major vegetation types (namely, steppe, shrub, coniferous forest and alpine meadow) increase with altitude in the elevation range of 1700–3000 m. Vegetation coverage increases with elevation in the order of steppe, shrub, coniferous forest and alpine meadow. Desert dominates the lowest elevation area, and alpine sparse vegetation dominates the highest elevation area. This spatial vegetation pattern is closely related to the changes in precipitation and temperature along the elevation.
Figure 6. Spatial distributions of the mean values of (a) annual precipitation; (b) annual evapotranspiration; (c) annual runoff and (d) growing season soil moisture of the top 1 m during 1981–2010.
Figure 6. Spatial distributions of the mean values of (a) annual precipitation; (b) annual evapotranspiration; (c) annual runoff and (d) growing season soil moisture of the top 1 m during 1981–2010.
Forests 07 00010 g006
Figure 7. Changes in vegetation, annual average precipitation, evapotranspiration (ET), runoff, runoff/P, and ET/P ratios along with elevation during 1981–2010.
Figure 7. Changes in vegetation, annual average precipitation, evapotranspiration (ET), runoff, runoff/P, and ET/P ratios along with elevation during 1981–2010.
Forests 07 00010 g007
Precipitation and runoff have a similar pattern, both increasing with elevation. The actual evapotranspiration shares a similar pattern with the four major vegetation types (i.e., steppe, shrub, coniferous forest and alpine meadow) along the elevation. Actual evapotranspiration increases with elevation when the elevation is lower than 3000 m. The highest actual evapotranspiration is at the elevations of 3000–3600 m, where shrub and alpine meadow are the dominant vegetation types. Previous research has found that during the growing season, vegetation cover is the densest in the elevation range of 3200–3400 m in the Qilian Mountains [44]. This implies that vegetation dynamics are most intensive in the elevation range of 3200–3400 m. When elevation is higher than 3400 m, actual evapotranspiration gradually decreases and vegetation dynamics gradually weaken due to the decreasing temperature. The region with elevation above 4200 m has the highest runoff depth due to high precipitation and lower evapotranspiration. Glacier melting is also a reason for the high runoff depth. Regarding the annual water balance characteristics, the regions with elevations lower than 3200 m are water limited. In these regions, actual evapotranspiration and vegetation growth are controlled by water availability (precipitation). The regions with elevation higher than 3400 m are temperature/energy limited. In these regions, actual evapotranspiration and vegetation growth are mainly controlled by temperature, and the actual evapotranspiration decreases with increasing elevation. Vegetation grows best and the actual evapotranspiration has the highest value in the region with elevations ranging from 3200 to 3400 m. The above spatial patterns of vegetation and water balance components confirm that climate variability has a significant effect on ecohydrological patterns in the high mountainous region.
The water balance of the three different elevation intervals of the East Tributary and West Tributary are calculated and shown in Table 4. The water balance (actual ET and runoff) differs among vegetation types in each elevation interval where the precipitation is close. Both of the runoff depth and runoff ratio for the coniferous forest is less than other vegetation types within the same elevation interval. This implies that vegetation type enhanced the differences in water balance in addition to the climate variability.
Table 4. Water balance in different elevation intervals of the East Tributary and West Tributary during the 2001–2012 period.
Table 4. Water balance in different elevation intervals of the East Tributary and West Tributary during the 2001–2012 period.
East Tributary
Elevation Interval and the AreaMajor Vegetation TypesArea RatioPrecipitation (mm/a)Actual ET (mm/a)Runoff Depth (mm/a)Runoff Ratio
2800–2999 m (89 km2)Shrub16%413.2371.139.10.09
Steppe22%410.9358.050.90.12
Coniferous forest12%402.0376.725.30.06
Alpine meadow51%395.0359.537.80.10
3400–3599 m (458 km2)Shrub42%498.0428.475.90.15
Steppe2%460.1413.850.30.11
Coniferous forest1%468.9406.661.90.13
Alpine meadow52%513.4439.678.10.15
4200–4399 m (78 km2)Alpine meadow8%564.4400.9170.00.30
Alpine sparse vegetation87%606.0299.6307.80.51
West Tributary
Elevation interval and the areaMajor vegetation typesArea ratioPrecipitation (mm/a)Actual ET (mm/a)Runoff depth (mm/a)Runoff ratio
3000–3199 m (123 km2)Shrub24%436.5381.153.20.12
Steppe12%439.6398.637.70.09
Coniferous forest13%450.8421.934.10.08
Alpine meadow51%418.5383.836.00.09
3400–3599 m (590 km2)Shrub18%447.5384.468.80.15
Steppe2%466.6401.568.30.15
Coniferous forest1%458.1402.556.40.12
Alpine meadow78%437.8384.465.60.15
4200–4399 m (634 km2)Alpine meadow35%458.4337.3129.70.28
Alpine sparse vegetation60%526.8264.8259.00.49
The water balance of the entire catchment for each vegetation type is also calculated and shown in Table 5, which is the result of the combined effects of climate and vegetation. Because the annual precipitation increases with elevation and different vegetation types grow at different elevations, the annual precipitation of alpine meadow, alpine sparse vegetation and shrub are 488.5 mm, 547.3 mm and 495.9 mm, respectively, which are higher than those of coniferous forest (402.1 mm) and steppe (396.7 mm). The annual average actual evapotranspiration (ET) of coniferous forest, shrub, steppe and alpine meadow ranges from 331.5 mm to 355.0 mm, whereas the actual ET of alpine sparse vegetation is relatively lower (237.2 mm). Table 5 shows that the annual runoff depth varies for different vegetation types. Alpine meadow and alpine sparse vegetation have higher annual runoff (147.8 mm and 310.1 mm, respectively) than forest and steppe (70.5 mm and 65.2 mm, respectively). The runoff ratios of the four major vegetation types, namely, steppe, shrub, coniferous forest and alpine meadow, range from 0.16 to 0.30, whereas alpine sparse vegetation has higher runoff ratio of close to 0.5 due to the high altitude. The water yield per unit area from different vegetation type is in order of alpine sparse vegetation, alpine meadow, shrub, coniferous forest and steppe. The top three largest vegetation areas are covered by alpine meadow (with an area of 4549 km2), alpine sparse vegetation (with an area of 2009 km2) and shrub (with an area of 1652 km2), which are also located in relatively high-elevation regions. The total runoff amount produced by the top three largest vegetation areas are 6.72 × 108 m3/a (alpine meadow), 6.23 × 108 m3/a (alpine sparse vegetation) and 2.33 × 108 m3/a (shrub) (see Table 5). The runoff amount produced by forests (with area of 561 km2) is 0.40 × 108 m3/a. The coniferous forest has a small area (561 km2) and produces a small amount of runoff (0.40 × 108 m3/a).
Figure 8 shows the decadal changes in water balance components. Annual precipitation, actual ET and runoff show similar patterns along the elevation. However, decadal variability of the water balance is also observed. Compared with the 1980s, precipitation decreased and ET increased in the 1990s, which caused a reduction in runoff. However, in the 2000s, precipitation increased in the high-altitude region with elevations above 3600 m, which led to a significant increase in runoff compared to the 1990s.
Figure 8. Changes in annual average precipitation, evapotranspiration and runoff along with elevation during the 1980s, 1990s and 2000s.
Figure 8. Changes in annual average precipitation, evapotranspiration and runoff along with elevation during the 1980s, 1990s and 2000s.
Forests 07 00010 g008
Table 5. Water balance of the entire catchment for each vegetation type during the 2001–2012 period.
Table 5. Water balance of the entire catchment for each vegetation type during the 2001–2012 period.
Vegetation Type and Area Covered by Each Type (km2)Precipitation (mm/a)Actual ET (mm/a)Runoff Depth (mm/a)Runoff RatioRunoff Amount (108 m3/a)
Desert91253.1238.015.10.060.01
Shrub1652495.9355.0140.90.282.33
Steppe1063396.7331.565.20.160.69
Coniferous forest561402.1331.670.50.180.40
Alpine meadow4549488.5348.7147.80.306.72
Alpine sparse vegetation2009547.3237.2310.10.576.23
Snow or glaciers80586.782.7846.21.440.68

4.4. Comparison with Previous Studies in the Same and Similar Regions

In this study, the model simulation shows that the precipitation recharges soil water and groundwater in the summer. This result is consistent with a previous study by Yang et al. [72] on the Heihe River. Using the hydrochemistry approach, it is found that the precipitation contributes only slightly to the total surface runoff, instead mainly recharging the sub-surface soil and groundwater. The interflow and groundwater flow dominate the total runoff. Additionally, Wang et al. [73] inferred that as the air temperature has risen in the Heihe basin during the last 50 years, the active layer of frozen soil has increased in thickness. This may lead to the increase of both soil water storage and interflow and thus change the contribution of each runoff component in the future.
Our study shows that forests contribute only a small amount of the water yield. In a previous study, He et al. [3] analyzed the water balance of a small experimental catchment (2.91 km2) in the upper Heihe basin and found that the forests contribute only 3.5% to the total runoff, which is similar to this study. This result is also supported by model simulations. For example, using the Topography Driven Flux Exchange (FLEX-Topo) model, Gao et al. [74] report that forest hillslope generates only a small amount of runoff in the upper Heihe basin. Qin et al. [75] analyze the water balance components of different landscapes in the upper Heihe basin using the VIC model and find that glaciers contribute 3.57% of the total runoff and that the contribution of forests is also quite small (0.49%). This result also implies that the barren regions contribute most of the total runoff (52.46%) and steppe contributes 34.15% to the total runoff. This is in accordance with our results, as Qin et al. [75] consider the steppe and alpine meadow as the same vegetation type and consider regions in high elevation as barren.
Our study shows that the runoff ratio of coniferous forest has the lowest value (0.18) comparing with other vegetation types except the desert in the upper Heihe River basin. This is similar to the findings by Yaseef et al. [76], which shows that when precipitation larger than 300 mm ET accounts 85% of precipitation for Aleppo pine forest in Southern Israel. Wang et al. [77] analyzed the annual water balance of shrub at a station in the Inner Mongolian Highland Region with elevation of 1300 m and found that ET/P ratio is about 94%, which is higher than the results of the present study because of the differences of elevation and air temperature in the two study areas. Wang et al. [78] analyzed the water balance of different vegetation in a small watershed in the Liupan Mountains, Northwest China, based on field measurements. They found that the evaporation rate (ratio of ET to precipitation) was about 60% for grassland, 93% for shrubs, and >95% for forest. This also shows that the forest has the least water yield compared with other vegetation types in the Liupan Mountains.

5. Conclusions

A geomorphology-based ecohydrological model (GBEHM) is developed in the upper Heihe basin, and this model is validated using available observations, including soil moisture, streamflow discharge, and actual evapotranspiration estimated from remote sensing. The catchment water balance characteristics and their spatial-temporal variability are analyzed based on the ecohydrological simulation. The following conclusions can be drawn from the results of this study:
(1)
At the basin scale, the model provides a good simulation of streamflow discharge in the two tributaries and the entire catchment of the study area. It also captures the spatial pattern of soil moisture appropriately. In addition, the simulated actual evapotranspiration and remote sensing-based estimation have close long-term average values and similar spatial patterns over the entire study catchment. The GBEHM may be useful for ecohydrological simulation and prediction in cold high-altitude regions.
(2)
Analysis of the water balance characteristics shows that water balance characteristics are closely related to the altitude and vegetation patterns in the study catchment. Regarding the annual water balance characteristics, the low-altitude regions with elevations below 3200 m are water limited. The actual annual evapotranspiration and vegetation distribution and growth are controlled by water availability (precipitation). Seasonal analysis indicates that river runoffs are mainly in summer and autumn, and runoff in spring is generated from precipitation and snow melt.
(3)
In the upper Heihe basin, the precipitation and runoff share a similar pattern, increasing with elevation. Actual evapotranspiration has a similar pattern with the four major vegetation types (i.e., steppe, shrub, coniferous forest and alpine meadow) along the elevation. The highest actual evapotranspiration is at the elevations of 3000–3600 m, where shrub and alpine meadow are the two dominant vegetation types. Precipitation controls the spatial pattern of annual runoff and determines the spatial pattern of vegetation together with the air temperature. Climate variability in the high mountainous region has a significant effect on ecohydrological patterns.
(4)
At the same time, vegetation type enhanced the differences in annual runoff and actual evapotranspiration. In the same elevation interval with similar precipitation, differences in the runoff depth (and the actual evapotranspiration) were caused mainly by the vegetation types. For the whole study area, the water yield per unit area from different vegetation types is in order of alpine sparse vegetation, alpine meadow, shrub, coniferous forest and steppe. The three major vegetation types, namely, alpine meadow (with an area of 4549 km2), alpine sparse vegetation (with an area of 2009 km2) and shrub (with an area of 1652 km2), located in relatively higher elevation contribute most of the river runoff.
Several limitations remain in the current study. The GBEHM simulates the vegetation dynamics with the known leaf area index and other vegetation parameters. Further improvement of the model should include carbon partitioning to simulate the vegetation growth. The uncertainty of the model parameters should be assessed to apply this model to other catchments.

Acknowledgments

This research was supported by the major plan of “Integrated Research on the Ecohydrological Processes of the Heihe Basin” (Project Nos. 91225302, 91425303) funded by the National Natural Science Foundation of China (NSFC).

Author Contributions

Dawen Yang conceived and designed the research. Yuanrun Zheng performed the classification of vegetation types. Bing Gao performed the model programming and simulation. Yue Qin and Yuhan Wang analyzed the simulated data. Yue Qin and Bing Gao took the lead in writing the manuscript. All authors contributed to the revision of the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhou, G.; Wei, X.; Luo, Y.; Zhang, M.; Li, Y.; Qiao, Y.; Liu, H.; Wang, C. Forest recovery and river discharge at the regional scale of Guangdong Province, China. Water Resour. Res. 2010, 46, 5109–5115. [Google Scholar] [CrossRef]
  2. Liu, W.; Wei, X.; Liu, S.; Liu, Y.; Fan, H.; Zhang, M.; Yin, J.; Zhan, M. How do climate and forest changes affect long-term streamflow dynamics? A case study in the upper reach of Poyang River basin. Ecohydrology 2015, 8, 46–57. [Google Scholar] [CrossRef]
  3. He, Z.B.; Zhao, W.Z.; Liu, H.; Tang, Z.X. Effect of forest on annual water yield in the mountains of an arid inland river basin: A case study in the Pailugou catchment on northwestern China’s Qilian Mountains. Hydrol. Process. 2012, 26, 613–621. [Google Scholar] [CrossRef]
  4. Li, X.; Cheng, G.D.; Liu, S.M.; Xiao, Q.; Ma, M.G.; Jin, R.; Che, T.; Liu, Q.H.; Wang, W.Z.; Qi, Y.; et al. Heihe Watershed Allied Telemetry Experimental Research (HiWATER): Scientific Objectives and Experimental Design. Bull. Am. Meteorol. Soc. 2013, 94, 1145–1160. [Google Scholar] [CrossRef]
  5. Chang, H.J.; Jung, I.W. Spatial and temporal changes in runoff caused by climate change in a complex large river basin in Oregon. J. Hydrol. 2010, 388, 186–207. [Google Scholar] [CrossRef]
  6. Belmont, P.; Gran, K.B.; Schottler, S.P.; Wilcock, P.R.; Day, S.S.; Jennings, C.; Lauer, J.W.; Viparelli, E.; Willenbring, J.K.; Engstrom, D.R.; et al. Large Shift in Source of Fine Sediment in the Upper Mississippi River. Environ. Sci. Technol. 2011, 45, 8804–8810. [Google Scholar] [CrossRef] [PubMed]
  7. Taye, M.T.; Ntegeka, V.; Ogiramoi, N.P.; Willems, P. Assessment of climate change impact on hydrological extremes in two source regions of the Nile River Basin. Hydrol. Earth Syst. Sci. 2011, 15, 209–222. [Google Scholar] [CrossRef]
  8. Zhou, D.G.; Huang, R.H. Response of water budget to recent climatic changes in the source region of the Yellow River. Chin. Sci. Bull. 2012, 57, 2155–2162. [Google Scholar] [CrossRef]
  9. Wu, L.; Long, T.Y.; Liu, X.; Guo, J.S. Impacts of climate and land-use changes on the migration of non-point source nitrogen and phosphorus during rainfall-runoff in the Jialing River Watershed, China. J. Hydrol. 2012, 475, 26–41. [Google Scholar] [CrossRef]
  10. Ficklin, D.L.; Stewart, I.T.; Maurer, E.P. Climate Change Impacts on Streamflow and Subbasin-Scale Hydrology in the Upper Colorado River Basin. PLoS ONE 2013, 8, e71297. [Google Scholar] [CrossRef] [PubMed]
  11. Cuo, L.; Zhang, Y.X.; Gao, Y.H.; Hao, Z.C.; Cairang, L.S. The impacts of climate change and land cover/use transition on the hydrology in the upper Yellow River Basin, China. J. Hydrol. 2013, 502, 37–52. [Google Scholar] [CrossRef]
  12. Yang, H.B.; Qi, J.; Xu, X.Y.; Yang, D.W.; Lv, H.F. The regional variation in climate elasticity and climate contribution to runoff across China. J. Hydrol. 2014, 517, 607–616. [Google Scholar] [CrossRef]
  13. Ma, H.; Yang, D.W.; Tan, S.K.; Gao, B.; Hu, Q.F. Impact of climate variability and human activity on streamflow decrease in the Miyun Reservoir catchment. J. Hydrol. 2010, 389, 317–324. [Google Scholar] [CrossRef]
  14. Xu, X.Y.; Yang, H.B.; Yang, D.W.; Ma, H. Assessing the impacts of climate variability and human activities on annual runoff in the Luan River basin, China. Hydrol. Res. 2013, 44, 940–952. [Google Scholar] [CrossRef]
  15. Tang, Y.; Tang, Q.; Tian, F.; Zhang, Z.; Liu, G. Responses of natural runoff to recent climatic variations in the Yellow River basin, China. Hydrol. Earth Syst. Sci. 2013, 17, 4471–4480. [Google Scholar] [CrossRef]
  16. Peel, M. Hydrology: Catchment vegetation and runoff. Prog. Phys. Geogr. 2009, 33, 837–844. [Google Scholar] [CrossRef]
  17. Samaniego, L.; Bardossy, A. Simulation of the impacts of land use/cover and climatic changes on the runoff characteristics at the mesoscale. Ecol. Model. 2006, 196, 45–61. [Google Scholar] [CrossRef]
  18. Marshall, E.; Randhir, T.O. Spatial modeling of land cover change and watershed response using Markovian cellular automata and simulation. Water Resour. Res. 2008, 44, W044234. [Google Scholar] [CrossRef]
  19. Yang, D.W.; Gao, B.; Jiao, Y.; Lei, H.M.; Zhang, Y.L.; Yang, H.B.; Cong, Z.T. A distributed scheme developed for eco-hydrological modeling in the upper Heihe River. Sci. China Earth Sci. 2015, 58, 36–45. [Google Scholar] [CrossRef]
  20. Abbott, M.B.; Bathurst, J.C.; Cunge, J.A.; O’connell, P.E.; Rasmussen, J. An introduction to the European Hydrological System—Systeme Hydrologique Europeen, “SHE”, 2: Structure of a physically-based, distributed modelling system. J. Hydrol. 1986, 87, 61–77. [Google Scholar] [CrossRef]
  21. Liang, X.; Lettenmaier, D.P.; Wood, E.F. One-dimensional statistical dynamic representation of subgrid spatial variability of precipitation in the two-layer variable infiltration capacity model. J. Geophys. Res. 1996, 101, 21403–21422. [Google Scholar] [CrossRef]
  22. Yang, D.W.; Koike, T.; Tanizawa, H. Application of a distributed hydrological model and weather radar observations for flood management in the upper Tone River of Japan. Hydrol. Process. 2004, 18, 3119–3132. [Google Scholar] [CrossRef]
  23. Yang, D.W.; Li, C.; Ni, G.H.; Hu, H.P. Application of a distributed hydrological model to the Yellow River basin. Acta Geogr. Sin. 2004, 59, 143–154. (In Chinese) [Google Scholar]
  24. Xu, J.J.; Yang, D.W.; Yi, Y.H.; Lei, Z.D.; Chen, J.; Yang, W.J. Spatial and temporal variation of runoff in the Yangtze River basin during the past 40 years. Quat. Int. 2008, 186, 32–42. [Google Scholar] [CrossRef]
  25. Cong, Z.T.; Yang, D.W.; Gao, B.; Yang, H.B.; Hu, H.P. Hydrological trend analysis in the Yellow River basin using a distributed hydrological model. Water Resour. Res. 2009, 45, W00A13. [Google Scholar] [CrossRef]
  26. Valeriano, O.; Koike, T.; Yang, K.; Yang, D.W. Optimal Dam Operation during Flood Season Using a Distributed Hydrological Model and a Heuristic Algorithm. J. Hydrol. Eng. 2010, 15, 580–586. [Google Scholar] [CrossRef]
  27. Alam, Z.R.; Rahman, M.M.; Islam, A.S. Assessment of Climate Change Impact on the Meghna River Basin using geomorphology based hydrological model (GBHM). In Proceedings of the ICWFM 2011—The 3rd International Conference on Water and Flood Management, Dhaka, Bangladesh, 8–10 January 2011.
  28. Sellers, P.J.; Randall, D.A.; Collatz, G.J.; Berry, J.A.; Field, C.B.; Dazlich, D.A.; Zhang, C.; Collelo, G.D.; Bounoua, L. A Revised Land Surface Parameterization (SiB2) for Atmospheric GCMS. Part I: Model Formulation. J. Clim. 1996, 9, 676–705. [Google Scholar] [CrossRef]
  29. Wang, L.; Koike, T.; Yang, K.; Jackson, T.J.; Bindlish, R.; Yang, D.W. Development of a distributed biosphere hydrological model and its evaluation with the Southern Great Plains Experiments (SGP97 and SGP99). J. Geophys. Res. 2009, 114. [Google Scholar] [CrossRef]
  30. Wang, L.; Koike, T. Comparison of a distributed biosphere hydrological model with GBHM. Annu. J. Hydraul. Eng. 2009, 53, 103–108. [Google Scholar]
  31. Oleson, K.W.; Lawrence, D.M.; Bonan, G.B.; Flanner, M.G.; Kluzek, E.; Lawrence, P.J.; Levis, S.; Swenson, S.C.; Thornton, P.E.; Dai, A.G.; et al. Technical Description of Version 4.0 of the Community Land Model (CLM), NCAR Technical Note; National Center for Atmospheric Research: Boulder, CO, USA, 2010; p. 257. [Google Scholar]
  32. Tague, C.; Band, L. RHESSys: Regional Hydro-ecologic simulation system: An object-oriented approach to spatially distributed modeling of carbon, water and nutrient cycling. Earth Interact. 2004, 19, 1–42. [Google Scholar] [CrossRef]
  33. Maneta, M.; Silverman, N. A spatially-distributed model to simulate water, energy and vegetation dynamics using information from regional climate models. Earth Interact. 2013, 17, 1–44. [Google Scholar] [CrossRef]
  34. Ivanov, V.Y.; Bras, R.L.; Vivoni, E.R. Vegetation-Hydrology Dynamics in Complex Terrain of Semiarid Areas: I. A mechanistic Approach to Modeling Dynamic Feedbacks. Water Resour. Res. 2008, 44. [Google Scholar] [CrossRef]
  35. Jia, Y.W.; Wang, H.; Yan, D.H. Distributed model of hydrological cycle system in Heihe River basin I Model development and Verification. J. Hydraul. Eng. 2006, 37, 534–542. (In Chinese) [Google Scholar]
  36. Wang, L.; Koike, T.; Yang, K.; Jin, R.; Li, H. Frozen soil parameterization in a distributed biosphere hydrological model. Hydrol. Earth Syst. Sci. 2010, 14, 557–571. [Google Scholar] [CrossRef]
  37. Zhang, Y.L.; Cheng, G.D.; Li, X.; Han, X.J.; Wang, L.; Li, H.Y.; Chang, X.L.; Flerchinger, G.N. Coupling of a simultaneous heat and water model with a distributed hydrological model and evaluation of the combined model in a cold region watershed. Hydrol. Process. 2013, 27, 3762–3776. [Google Scholar] [CrossRef]
  38. Zhou, J.; Pomeroy, J.W.; Zhang, W.; Cheng, G.D.; Wang, G.X.; Chen, C. Simulating cold regions hydrological processes using a modular model in the west of China. J. Hydrol. 2014, 509, 13–24. [Google Scholar] [CrossRef]
  39. 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]
  40. Cheng, G.; Li, X.; Zhao, W.; Xu, Z.; Feng, Q.; Xiao, S.; Xiao, H. Integrated study of the water-ecosystem-economy in the Heihe River Basin. Natl. Sci. Rev. 2014, 1, 413–428. [Google Scholar] [CrossRef]
  41. Chen, Y.; Zhang, D.; Sun, Y.; Liu, X.; Wang, N.; Savenije, H. Water demand management: A case study of the Heihe River Basin in China. Phys. Chem. Earth 2005, 30, 408–419. [Google Scholar] [CrossRef]
  42. Herzschuh, U.; Kurschner, H.; Mischke, S. Temperature variability and vertical vegetation belt shifts during the last similar to ~50,000 yr in the Qilian Mountains (NE margin of the Tibetan Plateau, China). Quat. Res. 2006, 66, 133–146. [Google Scholar] [CrossRef]
  43. Wang, P.; Li, Z.Q.; Gao, W.Y. Rapid Shrinking of Glaciers in the Middle Qilian Mountain Region of Northwest China during the Last similar to 50 Years. J. Earth Sci. 2011, 22, 539–548. [Google Scholar] [CrossRef]
  44. Deng, S.F.; Yang, T.B.; Zeng, B.; Zhu, X.F.; Xu, H.J. Vegetation cover variation in the Qilian Mountains and its response to climate change in 2000–2011. J. Mt. Sci. 2013, 10, 1050–1062. [Google Scholar] [CrossRef]
  45. National Meteorological Information Center, the China Meteorological Administration. Available online: http://cdc.nmic.cn (accessed on 24 December 2015).
  46. Shen, Y.; Xiong, A. Validation and comparison of a new gauge-based precipitation analysis over mainland China. Int. J. Climatol. 2015. [Google Scholar] [CrossRef]
  47. Wu, L.; Li, X. Dataset of the First Glacier Inventory in China; Cold and Arid Regions Science Data Center: Lanzhou, China, 2004. [Google Scholar]
  48. Guo, W.; Liu, S.; Yao, X.; Xu, J.; Shangguan, D.; Wu, L.; Zhao, J.; Liu, Q.; Jiang, Z.; Wei, J.; et al. The Second Glacier Inventory Dataset of China (Version 1.0); Cold and Arid Regions Science Data Center: Lanzhou, China, 2014. [Google Scholar]
  49. Wei, J.F.; Liu, S.Y.; Guo, W.Q.; Yao, X.J.; Xu, J.L.; Bao, W.J.; Jiang, Z.L. Surface-area changes of glaciers in the Tibetan Plateau interior area since the 1970s using recent Landsat images and historical maps. Ann. Glaciol. 2014, 55, 213–222. [Google Scholar] [CrossRef]
  50. Cold and Arid Regions Science Data Center at Lanzhou. Available online: http://westdc.westgis.ac.cn (accessed on 24 December 2015).
  51. Jarvis, A.; Reuter, H.I.; Nelson, A.; Guevara, E. Hole-filled seamless SRTM data (Version 4). International Centre for Tropical Agriculture. 2008. Available online: http://srtm.csi.cgiar.org (accessed on 31 July 2015).
  52. Hou, X. China Vegetation Map (1:1,000,000); Science Press: Beijing, China, 2001. (In Chinese) [Google Scholar]
  53. Zhou, J.H.; Zheng, Y.R. Vegetation Map of the Upper Heihe Basin V2.0; Heihe Plan Science Data Center: Lanzhou, China, 2014. [Google Scholar]
  54. Fan, W. Heihe 1 km LAI Production; Heihe Plan Science Data Center: Lanzhou, China, 2014. [Google Scholar]
  55. Shi, X.; Yu, D.; Pan, X. A framework for the 1:1,000,000 soil database of China. In Proceedings of the 17th World Congress of Soil Science, Bangkok, Thailand, 14–21 August 2002; pp. 1–5.
  56. Dai, Y.J.; Shangguan, W.; Duan, Q.Y.; Liu, B.Y.; Fu, S.H.; Niu, G. Development of a China Dataset of Soil Hydraulic Parameters Using Pedotransfer Functions for Land Surface Modeling. J. Hydrometeorol. 2013, 14, 869–887. [Google Scholar] [CrossRef]
  57. Jin, R.; Li, X.; Yan, B.P.; Li, X.H.; Luo, W.M.; Ma, M.G.; Guo, J.W.; Kang, J.; Zhu, Z.L.; Zhao, S.J. A Nested Ecohydrological Wireless Sensor Network for Capturing the Surface Heterogeneity in the Midstream Areas of the Heihe River Basin, China. IEEE Geosci. Remote Sens. Lett. 2014, 11, 2015–2019. [Google Scholar] [CrossRef]
  58. Kang, J.; Jin, R.; Li, X.; Ma, M. HiWATER: WATERNET Observation Dataset in the Upper Reaches of the Heihe River Basin in 2013; Cold and Arid Regions Science Data Center: Lanzhou, China, 2014. [Google Scholar]
  59. Wu, B.F. Monthly Evapotranspiration Datasets (2000–2012) with 1 km Spatial Resolution over the Heihe River Basin; Heihe Plan Science Data Center: Lanzhou, China, 2013. [Google Scholar]
  60. Wu, B.F.; Yan, N.N.; Xiong, J.; Bastiaanssen, W.; Zhu, W.W.; Stein, A. Validation of ETWatch using field measurements at diverse landscapes: A case study in Hai Basin of China. J. Hydrol. 2012, 436, 67–80. [Google Scholar] [CrossRef]
  61. Yang, D.W.; Herath, S.; Musiake, K. Development of a geomorphology-based hydrological model for large catchments. Annu. J. Hydraul. Eng. 1998, 42, 169–174. [Google Scholar] [CrossRef]
  62. Yang, D.W.; Herath, S.; Musiake, K. Comparison of different distributed hydrological models for characterization of catchment spatial variability. Hydrol. Process. 2000, 14, 403–416. [Google Scholar] [CrossRef]
  63. Yang, D.W.; Herath, S.; Musiake, K. Spatial resolution sensitivity of catchment geomorphologic properties and the effect on hydrological simulation. Hydrol. Process. 2001, 15, 2085–2099. [Google Scholar] [CrossRef]
  64. Yang, D.W.; Herath, S.; Musiake, K. A hillslope-based hydrological model using catchment area and width functions. Hydrol. Sci. J. 2002, 47, 49–65. [Google Scholar] [CrossRef]
  65. Sellers, P.J. Canopy reflectance, photosynthesis, and transpiration. Int. J. Remote Sens. 1985, 8, 1335–1372. [Google Scholar] [CrossRef]
  66. Sellers, P.J.; Mintz, Y.; Sub, Y.C.; Dalcher, A. A Simple Biosphere Model (SiB) for Use within General Circulation Models. J. Atmos. Sci. 1986, 43, 505–531. [Google Scholar] [CrossRef]
  67. Sellers, P.J.; Berry, J.A.; Collatz, G.J.; Field, C.B.; Hall, F.G. Canopy reflectance, photosynthesis, and transpiration Part III: Are analysis using improved leaf models and a new canopy integration scheme. Remote Sens. Environ. 1992, 42, 187–216. [Google Scholar] [CrossRef]
  68. Collatz, G.J.; Ball, J.T.; Grivet, C.; Berry, J.A. Physiological and environmental regulation of stomatal conductance, photosynthesis and transpiration: A model that includes a laminar boundary layer. Agric. For. Meteorol. 1991, 54, 107–136. [Google Scholar] [CrossRef]
  69. Collatz, G.J.; Ribas-Carbo, M.; Berry, J.A. Coupled Photosynthesis-Stomatal Conductance Model for leaves of C4 plants. Aust. J. Plant Physiol. 1992, 19, 519–538. [Google Scholar] [CrossRef]
  70. Zhang, Y.; Liu, S.Y.; Ding, Y.J. Spatial Variation of Degree-day Factors on the Observed Glaciers in Western China. Acta Geogr. Sin. 2006, 61, 89–98. (In Chinese) [Google Scholar] [CrossRef]
  71. Gao, X.; Ye, B.S.; Zhang, S.Q.; Qiao, C.J.; Zhang, X.W. Glacier runoff variation and its influence on river runoff during 1961–2006 in the Tarim River Basin, China. Sci. China Earth Sci. 2010, 53, 880–891. (In Chinese) [Google Scholar] [CrossRef]
  72. Yang, Y.; Xiao, H.; Wei, Y.; Zhao, L.; Zou, S.; Qiu, Y.; Yin, Z. Hydrological processes in the different landscape zones of alpine cold regions in the wet season, combining isotopic and hydrochemical tracers. Hydrol. Process. 2012, 26, 1457–1466. [Google Scholar] [CrossRef]
  73. Wang, Y.; Yang, D.; Lei, H.; Yang, H. Impact of cryosphere hydrological processes on the river runoff in the upper reaches of Heihe River. J. Hydraul. Eng. 2015, 46, 1064–1071. (In Chinese) [Google Scholar]
  74. Gao, H.; Hrachowitz, M.; Fenicia, F.; Gharari, S.; Savenije, H. Testing the realism of a topography-driven model (FLEX-Topo) in the nested catchments of the Upper Heihe, China. Hydrol. Earth Syst. Sci. 2014, 18, 1895–1915. [Google Scholar] [CrossRef] [Green Version]
  75. Qin, J.; Ding, Y.; Wu, J.; Gao, M.; Yi, S.; Zhao, C.; Ye, B.; Li, M.; Wang, S. Understanding the impact of mountain landscapes on water balance in the upper Heihe River watershed in northwestern China. J. Arid Land 2013, 5, 366–384. [Google Scholar] [CrossRef]
  76. Yaseef, N.R.; Yakir, D.; Rotenberg, E.; Schiller, G.; Cohen, S. Ecohydrology of a semi-arid forest: Partitioning among water balance components and its implications for predicted precipitation changes. Ecohydrology 2009. [Google Scholar] [CrossRef]
  77. Wang, X.; Ronny, B.; Li, X.; Kang, E. Water balance change for a re-vegetated xerophyte shrub area. Hydrol. Sci. J. 2004, 49, 283–295. [Google Scholar] [CrossRef]
  78. Wang, Y.; Yu, P.; Xiong, W.; Shen, Z.; Guo, M.; Shi, Z.; Du, A.; Wang, L. Water-yield reduction after afforestation and related processes in the semiarid Liupan Mountains, Northwest China. J. Am. Water Resour. Assoc. 2008, 44, 1086–1097. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Gao, B.; Qin, Y.; Wang, Y.; Yang, D.; Zheng, Y. Modeling Ecohydrological Processes and Spatial Patterns in the Upper Heihe Basin in China. Forests 2016, 7, 10. https://doi.org/10.3390/f7010010

AMA Style

Gao B, Qin Y, Wang Y, Yang D, Zheng Y. Modeling Ecohydrological Processes and Spatial Patterns in the Upper Heihe Basin in China. Forests. 2016; 7(1):10. https://doi.org/10.3390/f7010010

Chicago/Turabian Style

Gao, Bing, Yue Qin, Yuhan Wang, Dawen Yang, and Yuanrun Zheng. 2016. "Modeling Ecohydrological Processes and Spatial Patterns in the Upper Heihe Basin in China" Forests 7, no. 1: 10. https://doi.org/10.3390/f7010010

APA Style

Gao, B., Qin, Y., Wang, Y., Yang, D., & Zheng, Y. (2016). Modeling Ecohydrological Processes and Spatial Patterns in the Upper Heihe Basin in China. Forests, 7(1), 10. https://doi.org/10.3390/f7010010

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