This article was downloaded by: [NUS National University of Singapore]
On: 03 June 2014, At: 22:40
Publisher: Taylor & Francis
Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer House,
37-41 Mortimer Street, London W1T 3JH, UK

Hydrological Sciences Journal
Publication details, including instructions for authors and subscription information:
/>
Impact of climate and land-use changes on hydrological
processes and sediment yield—a case study of the Be
River catchment, Vietnam
a

Dao Nguyen Khoi & Tadashi Suetsugi

b

a

Faculty of Environmental Science, University of Science, Vietnam National University, Ho
Chi Minh City, Vietnam
b

Interdisciplinary Graduate School of Medicine and Engineering, University of Yamanashi,
Kofu, Yamanashi 400-8511, Japan
Accepted author version posted online: 04 Jul 2013.Published online: 29 Apr 2014.

To cite this article: Dao Nguyen Khoi & Tadashi Suetsugi (2014) Impact of climate and land-use changes on hydrological
processes and sediment yield—a case study of the Be River catchment, Vietnam, Hydrological Sciences Journal, 59:5,
1095-1108, DOI: 10.1080/02626667.2013.819433


Downloaded by [NUS National University of Singapore] at 22:40 03 June 2014

Interdisciplinary Graduate School of Medicine and Engineering, University of Yamanashi, Kofu, Yamanashi 400-8511, Japan

Received 9 June 2012; accepted 8 May 2013; open for discussion until 1 November 2014
Editor Z.W. Kundzewicz; Associate editor Q. Zhang
Citation Khoi, D.N. and Suetsugi, T., 2014. Impact of climate and land-use changes on hydrological processes and sediment yield—a
case study of the Be River catchment, Vietnam. Hydrological Sciences Journal, 59 (5), 1095–1108.

Abstract The impact of climate and land-use changes on hydrological processes and sediment yield is investigated in the Be River catchment, Vietnam, using the Soil and Water Assessment Tool (SWAT) hydrological
model. The sensitivity analysis, model calibration and validation indicated that the SWAT model could reasonably
simulate the hydrology and sediment yield in the catchment. From this, the responses of the hydrology and
sediment to climate change and land-use changes were considered. The results indicate that deforestation had
increased the annual flow (by 1.2%) and sediment load (by 11.3%), and that climate change had also significantly
increased the annual streamflow (by 26.3%) and sediment load (by 31.7%). Under the impact of coupled climate
and land-use changes, the annual streamflow and sediment load increased by 28.0% and 46.4%, respectively. In
general, during the 1978–2000 period, climate change influenced the hydrological processes in the Be River
catchment more strongly than the land-use change.
Key words climate change; hydrology; land-use change; sediment yield; SWAT model; Be River catchment, Vietnam

Impact des changements climatiques et de l’utilisation des terres sur les processus hydrologiques
et la production de sédiments—étude de cas du bassin versant de la rivière Be, Vietnam
Résumé L’impact des changements du climat et de l’utilisation des terres sur les processus hydrologiques et
l’apport de sédiments dans le bassin versant de la rivière Be (Vietnam) a été étudié en utilisant le modèle
hydrologique SWAT. L’analyse de sensibilité, l’étalonnage et la validation des modèles indique que le modèle
SWAT peut raisonnablement simuler l’hydrologie et la charge sédimentaire dans le bassin versant. C’est donc
avec cet outil que les réponses de l’hydrologie et des sédiments au changement climatique et au changement
d’utilisation des terres ont été étudiées. Les résultats indiquent que la déforestation a augmenté l’écoulement
annuel (1,2%) et la charge sédimentaire (11,3%), et que le changement climatique a également augmenté de

factors controlling the hydrological and sediment behaviours of catchments (Elfert and Bormann 2010). It is
important to understand the hydrological and sediment
responses to these changes in order to develop strategies for land-use planning and water resource management. Studies of the hydrological and water quality
impacts of climate change and land-use change are
desirable (Tong et al. 2012).
Many studies have considered the impact of climate change and land-use change on hydrology (Li
et al. 2009, 2012, Ma et al. 2009, 2010, Mango et al.
2011, Zhang et al. 2011). However, few studies have
investigated changes in hydrological processes and
water quality as well as sediment yield under the impact
of climate and land-use changes on a basin scale (Ward
et al. 2009, Tong et al. 2012). To assess the hydrological
and sediment impacts of environmental change, the
common methods used are the paired catchment
approach, statistical analysis and hydrological modelling (Li et al. 2009, 2012). Among these approaches,
the hydrological method is an appealing option, because
it is most suitable to be used as a part of scenario studies.
There are numerous hydrological models, such as the
Water Erosion Prediction Project (WEPP), Hydrologic
Simulation Program Fortran (HSPF), the Soil and Water
Assessment Tool (SWAT) and the physically-based
distributed hydrological model Système Hydrologique
Européen TRANsport (SHETRAN), that could be used
in simulating the runoff and transport of sediment and
pollutants in the catchment. The SWAT model has been
selected for the current study because it is widely used
to assess hydrology and water quality in agricultural
catchments around the world (see the SWAT literature
database: />Another reason for its selection is its availability and
user-friendliness in terms of handling input data

located in Dak Nong, Binh Phuoc, Binh Duong and
Dong Nai provinces, and has a catchment area of
about 7500 km2. The altitude varies from 1000 m a.m.
s.l. in the highland area to 100 m a.m.s.l. in the plain
area, in a northeast to southwest and south direction. The
origin of the branched-tree drainage system of the Be
River lies in Tuy Duc on the international border
between Vietnam and Cambodia, in Dak Nong
province. The study area is located in the steep area.
The degree of slope can be divided into three levels:
slopes of 0% to 7% account for 45% of the total area,
slopes of 8–15% account for 33% of the area, and slopes
greater than 15% account for 22% of the area. The
climate is tropical monsoon. The annual rainfall varies
between 1800 and 2800 mm, with an average of
2400 mm year-1. The area has two seasons: the rainy
season and the dry season. The rainy season lasts from
May to November and accounts for 85–90% of the total
annual precipitation. The average temperature is about
25.9°C, the maximum temperature is 36.6°C and the
minimum temperature is 17.3°C. The area has relatively
fertile land (75% basalt soil), consistent with agricultural
development. The main land-use types in this catchment
are forest and agricultural lands. The total population in
2010 was approximately one million inhabitants. The
mean annual flow of the catchment is about
7.51 × 109 m3. Similar to the distribution of rainfall,
the flow is distinguished by two distinct seasons: the
flood season (accounting for 67% of the total annual


S >0
S¼0
S

xj À xi
"i < j
β ¼ Median
jÀi

(5)


1098

Dao Nguyen Khoi and Tadashi Suetsugi

where 1 < i < j < n. The estimator β is calculated as
the median of all slopes between data pairs for the
entire data set.
The Pettitt test (Pettitt 1979) is a non-parametric
approach used for detecting the change point. There
are two samples (x1, x2, …, xt) and (xt+1, xt+2, …, xN)
that come from the same population (x1, x2, …, xN).
The test statistic Ut,N is given by:
Ut;N ¼

t X
N
X

À
Á
sgn xi À xj

Soil Loss Equation (MUSLE) to simulate the sediment yield for each HRU. The MUSLE (William
1995) is given as:


(8)

When p is smaller than the specific significance
level, the null hypothesis is not accepted. The time t
when Kt occurs is the change point time.
These methods have been commonly used to
detect changes in hydro-meteorological data (Ma
et al. 2008, Zhang et al. 2009, Zhang et al. 2011).
3.2 SWAT model
The SWAT model is a physically based, distributed,
continuous time model that is designed to predict the
effects of land management on the hydrology, sediment
and agricultural chemical yields in agricultural watersheds with varying soils, land-use and management
conditions (Arnold et al. 1998). In the SWAT model,
a catchment is divided into a number of sub-watersheds
or sub-basins. Sub-basins are further partitioned into
hydrological response units (HRUs) based on soil
types, land-use and slope classes that allow a high
level of spatial detail simulation. The model predicts
the hydrology at each HRU using the water balance
equation, comprising precipitation, surface runoff, evapotranspiration, infiltration and subsurface flow.
The SWAT model provides two methods for estimating surface runoff: the SCS curve number procedure
(USDA-SCS 1972) and the Green and Ampt infiltration
method (Green and Ampt 1911). SWAT calculates the
peak runoff rate using a modified rational method. The
potential evapotranspiration is estimated in the SWAT

image in 2001—obtained from the US Geological
Survey Earth Resources Observation and Science
Center. Land-use maps were generated using supervised
classification based on the maximum likelihood algorithm in the ENVI Version 4.4 image processing software. Overall accuracy and kappa statistic (κ) were used
to assess classification accuracy based on 256 ground
control points selected from the referenced land-use map


Impact of climate and land-use changes on hydrological processes

1099

Downloaded by [NUS National University of Singapore] at 22:40 03 June 2014

Table 1 Spatial model input data for the Be River catchment.
Data type

Description

Resolution Source

Topographic map
Land-use map
Soil map
Weather

Digital elevation map (DEM)
Land-use classification
Soil types
Daily precipitation, minimum and maximum temperature

root mean square error (RMSE) to the standard
deviation (STDEV) of measured data, RSR. The
RSR is calculated as (Moriasi et al. 2007):
sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
N
P
ðOi À Pi Þ2

RSR ¼

i¼1
RMSE
¼ sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
STDEVobs
N
P
 Þ2
ð Oi À O
i¼1

(10)

SRTM
Landsat TM, ETM+ (USGS/GLOVIS)
FAO
Hydro-Meteorological Data Center (HMDC)

where Oi is the observed value Pi is the simulated
 is the mean of the observed data, and N is
value, O

2002, 2008). The population of the Be River catchment was about 680 000 in 2000 compared to
400 000 in 1990 (SIWRP 2002). This represents a
population increase of about 170%.


Downloaded by [NUS National University of Singapore] at 22:40 03 June 2014

1100

Dao Nguyen Khoi and Tadashi Suetsugi

Fig. 2 Land-use maps of the Be River catchment.
Table 2 Statistics for land-use changes in the Be River catchment for the period of 1978–2007.
Land-use types

Agricultural land
Range land
Forest
Urban
Water
Total

1990

2001

2

2


91
7484

55.18
6.85
36.62
0.13
1.22
100

1114
–152
–1053
8
83

14.89
–2.03
–14.07
0.11
1.11

4.2 Change detection for hydro-meteorological
data
Annual temperature, precipitation and streamflow
were tested using the Mann-Kendall and Pettitt methods, as reported in Table 3 and illustrated in Fig. 3.
The results showed rises in annual temperature, precipitation and streamflow (by 0.035°C year-1,

Table 3 Summary of Mann-Kendall trend test and Pettitt
test statistics for annual rainfall, temperature and streamflow in the Be River catchment.

*

88
102
78

1989
1986
1989

*
*
*

*indicates significant at p < 0.05.

Change

20.613 mm year-1 and 3.142 m3 s-1 year-1, respectively) at the 5% significance level. In other words,
the null hypothesis H0 was not accepted for the
annual temperature, rainfall and streamflow time series. Change points in the annual rainfall and streamflow were detected as occurring around 1989, with a
significance level of 5%, while the change point in
annual temperature was statistically significant
in 1986.
The Mann-Kendall test was also applied to the
data series for monthly precipitation and temperature,
as summarized in Table 4. There are no significant
trends in most of the monthly precipitation time
series, except for October and December. The precipitation in October and December showed significant increasing trends of 2.24 and 1.73 mm year-1,
respectively. In the case of the monthly temperature,

December

0.66
1.58
0.63
0.58
1.85
–0.48
1.40
–1.11
–0.16
2.27
0.53
2.09

β
0.208
0.589
0.762
0.741
3.313
–2.740
2.783
–3.676
–0.233
2.244
1.325
1.733

Temperature

0.050
0.085

p
*

*
*
*
*
*
*
*

*indicates significant at p < 0.05.

catchment scale. The meteorological data were
divided into two periods, 1978–1989 and 1990–
2000, based on the change point analysis, and each
period included one land-use map. The land-use map
for 1990 was used to represent the 1978–1989 period, and that for 2001 was used to represent the
1990–2000 period. The following four scenarios
were investigated:
–
–
–
–

Fig. 3 Variations of mean values in (a) annual precipitation, (b) annual temperature and (c) annual discharge in the
Be River catchment (1978–2000).


1102

Dao Nguyen Khoi and Tadashi Suetsugi

Table 5 SWAT sensitivity parameters and calibrated values.
Simulation Parameter

Flow

Downloaded by [NUS National University of Singapore] at 22:40 03 June 2014

Sediment

CN2
ESCO
GQWMN
ALPHA_BF
SOL_Z
SOL_AWC
CH_K2
GW_REVAP
CH_N2
SOL_K
SPCON
SPEXP
USLE_P

Description of parameter


Calibrated value

–1
– 5000
–1
– 500
– 0.2
– 0.3
– 0.01

Phuoc Long

Phuoc Hoa

–0.29
0.95
456
0.11
0.03
0.23
184
0.17
0.04
–0.06
0.001

–0.36
0.42
2356
0.61

of observed sediment load, these data were only available from 07/1999 to 2004 at monthly levels. They were
divided into two periods for calibration (07/1999–2001)
and validation (2002–2004) using the land-use map for
2001. The flow calibration and validation was conducted
first, and then the sediment calibration and validation.
As a result of the calibration, the most sensitive flowrelated parameter of CN2 was adjusted to have values of
–0.29 for Phuoc Long and –0.36 for Phuoc Hoa, and the
most sensitive sediment-related parameter SPCON was
adjusted to have a value of 0.001. This value of SPCON
found here was similar to that in the study conducted by
Phan et al. (2011) in the Cau River watershed in northern
Vietnam. The details of the calibrated parameters are
presented in Table 5.
The SWAT flow simulations were calibrated
against the daily flow from 1981 to 1989 and validated
from 1990 to 1993 at the Phuoc Long gauging station,
as shown in Fig. 4. The simulated daily flow fit the

Fig. 4 Observed and simulated daily flow hydrograph at
the Phuoc Long station: (a) calibration and (b) validation.

observed data for the calibrated period well, with NSE,
PBIAS and RSR values of 0.77, 1.60% and 0.47,
respectively. For the validation period, the values of
NSE = 0.79, PBIAS = 3.30% and RSR = 0.45 suggest
that there was good agreement between the simulated
and observed streamflow during this period, based on


1103

PBIAS

RSR

Daily
Monthly
Daily
Monthly

0.77
0.87
0.79
0.91

1.60%
1.60%
3.30%
3.30%

0.48
0.36
0.45
0.30

Calibration (1981–1989)

Daily
Monthly
Daily
Monthly


Validation (1990–2000)

This match is shown in Table 6. Although the simulated
and observed streamflow followed the same trend, the
peak flow was overestimated for Phuoc Long station
and underestimated for Phuoc Hoa station. This may
have resulted from the uneven spatial distribution of the
rain gauges. In the study area, eight rain gauges are
located in the lower area of the catchment; however,
only one rain gauge located in the upper area of the
catchment has long-term records (Fig. 1). A further
reason can be attributed to the CN2, which is used to
simulate the surface runoff. The CN2 method assumes a
unique relationship between cumulative rainfall and
cumulative runoff for the same antecedent moisture
conditions (Betrie et al. 2011). Generally speaking,
these results reveal that the hydrological processes in
SWAT are modelled realistically for the Be River catchment, which is important for the simulation of sediment.
The simulated sediment load values were calibrated against monthly observed data from 07/1999 to
2001 and validated from 2002 to 2004 at the Phuoc Hoa
station, as presented in Fig. 6. The fit between the
simulated and observed sediment loads was acceptable,
according to Moriasi et al. (2007). The fit was indicated
by the values of NSE = 0.74, RSR = 0.51 and
PBIAS = –1.10% for the calibration period and
NSE = 0.55, RSR = 0.66 and PBIAS = 33.77% for
the validation period (Table 7). Although an underestimation of the monthly sediment yield by the model for
the validation period was within the satisfactory level of
acceptance, it can generally be said that the simulated

Time step NSE PBIAS

Calibration (07/1999–2001) Monthly
Validation (2002–2004)
Monthly

0.74
0.55

Fig. 8 Annual changes of hydrological components under
the impact of climate and land-use changes.

RSR

–1.10% 0.51
33.77% 0.66

conditions of the 1990s and the climate data of two
different periods: 1978–1989 and 1990–2000.
Figure 7 shows the absolute changes in the monthly
climate variables between the two periods. Compared
with the 1978–1989 period, the annual temperature
increased by 0.4°C and the annual precipitation
increased by 292.3 mm (12.8%). The increases in
temperature and precipitation were higher in the dry
season (0.5°C and 44.8%) than in the wet season
(0.4°C and 10.1%).
In the case of water balance components, climate
change caused increases in all water balance components, including a 7.9% increase in actual evapotranspiration, a 19.5% increase in groundwater discharge,
a 34.2% increase in surface runoff and a 56.1%

streamflow will increase the sediment load, while a
decrease in the streamflow will decrease the sediment
load.
4.6 Response to land-use change
The impact of land-use change on the water balance
components is illustrated in Fig. 8. Under the impact
of land-use change, surface runoff, soil water content
and sediment yield increased considerably, by 18.4%,
10.4% and 12.8%, respectively, while actual evapotranspiration and water yield increased slightly, by
approximately 1.3% and 1.1%, respectively. Aside
from this, the other water balance components
decreased, including a 5.8% decrease in groundwater
discharge, a 4.6% decrease in lateral flow, and a 4%
decrease in the amount of water percolating out of the
root zone. Deforestation and agricultural expansion
could be the cause of these changes. This is because
forest vegetation intercepts more water than other
land-use types (Ma et al. 2009), and the infiltration
rate of forest land is large compared with the other

1105

land-use types (Bruijnzeel 1990). Therefore, it is
likely that deforestation in the Be River catchment
caused an increase in runoff and decreases in groundwater discharge and lateral flow.
Under the impact of land-use change, deforestation and the increase in agricultural land resulted in
an increase in annual streamflow (1.2%) and sediment load (11.3%). Considering the seasonal change,
the streamflow decreased by 4.6% in the dry season
and increased by 1.8% in the wet season. In the case
of sediment load, it increased significantly in both the

land-use change occur simultaneously. In contrast,
when the directions of the changes affected by climate change alone and land-use change alone are
opposite, the change is reduced when climate change
and land-use change occur concurrently.


1106

Dao Nguyen Khoi and Tadashi Suetsugi

Table 8 Simulated streamflow at the Phuoc Hoa station under the impacts of climate and land-use changes.
Scenario

Downloaded by [NUS National University of Singapore] at 22:40 03 June 2014

1
2
3
4

Land-use

1990
1990
2001
2001

Climate

1978–1989

–
46

–
–
–
22.1

–
54
3
57

–
25.4
1.4
26.8

4.8 Limitations and recommendations
The SWAT hydrological model was successfully
applied in this study to the Be River catchment to
assess the impact of climate and land-use changes on
hydrology and sediment yield. However, there are
limitations in both the data and the model, which
are described as follows.
One of the limitations in this study comes from
the unavailability of data. Because of the lack of
sediment load data, sediment simulation is calibrated
and validated in monthly time steps for only a short
period of time, whereas hydrological modelling is

8 and 9, Table 8).

3

the streamflow beyond the time period of observed
sediment data using the calibrated SWAT model for
streamflow simulation in order to calibrate and validate the sediment simulation. These decrease the
accuracy of the model performance in sediment simulation. Therefore, to improve the simulation results,
collecting additional data on sediment load should be
considered to improve the model performance in
streamflow and sediment yield simulations.
Another limitation comes from the SWAT model.
The SWAT model uses a number of empirical and quasiphysical equations that were developed based on the
climate conditions in the United States, and those equations may not be appropriate for the tropical climate in
Vietnam. For example, the CN2 equation was a product
of more than 20 years of studies involving rainfall–runoff relationships in small rural watersheds across the
United States (Neitsch et al. 2011). In addition, the
MUSLE was also developed based on the hydrological
conditions throughout the United States. In the tropical
area, the heavy rainfall that may accompany a storm has
the potential to erode as much surface soil in the catchment as the subsequent runoff, but the MUSLE does not
account for such factors (Phomcha et al. 2011). It is
suggested that some parameters in the empirical equation should be modified to suit the tropical climate area
in order to improve the simulation results. The use of
oversimplified sediment routing algorithms to simulate
both landscape and in-stream erosion is a further limitation of the SWAT model. Aside from this, the SWAT
model allows all the soil eroded by runoff to reach the
channel directly, without considering sediment deposition remaining on surface catchment areas (Oeurng et al.
2011). Even though the SWAT model has some limitations, the simulation results were within the performance
criteria provided by Moriasi et al. (2007).

than the land-use change in the catchment during the
1978–2000 period. Therefore, when planning and
managing for water resources, the importance of
increasing adaptation to climate change is emphasized. However, with the considerable changes in
the surface runoff and sediment load under the
impact of land-use change, the effect of land-use
change should be accounted for in water resource
management in the Be River catchment.
Investigating not only the separate but also the
combined impacts of climate and land-use changes
helps to enhance our understanding of these impacts
on hydrological processes and soil erosion in the
catchment. The results obtained from this study
could be of value to managers/decision-makers in
integrated river basin management, as well as to the
development of adaptation and mitigation strategies
regarding climate and land-use changes.

Acknowledgement The authors thank their colleagues in Vietnam for assisting with the data. They
are grateful for the comments of two anonymous
reviewers, which greatly enhanced the quality of the
manuscript.

1107

Funding The authors acknowledge the Global Center
of Excellence (GCOE) program of the University of
Yamanashi, which funded this study.

REFERENCES

Ma, H., et al., 2010. Impact of climate variability and human activity
on streamflow decrease in the Miyun Reservoir catchment.
Journal of Hydrology, 389, 317–324. doi:10.1016/j.
jhydrol.2010.06.010.
Ma, X., et al., 2009. Response of hydrological processes to landcover and climate changes in Kejie watershed, South-West
China. Hydrological Processes, 23, 1179–1191. doi:10.1002/
hyp.7233.
Ma, Z., et al., 2008. Analysis of impacts of climate variability and
human activity on streamflow for a river basin in arid region of
Northwest China. Journal of Hydrology, 352, 239–249.
doi:10.1016/j.jhydrol.2007.12.022.
Mango, L.M., et al., 2011. Land use and climate change impacts on
the hydrology of the upper Mara River Basin, Kenya: results of
a modeling study to support better resource management.
Hydrology and Earth System Sciences, 15, 2245–2258.
doi:10.5194/hess-15-2245-2011.
Mann, H.B., 1945. Non-parametric tests against trend. Econometrica,
13, 245–259. doi:10.2307/1907187.
MONRE (Ministry of Natural Resources and Environment), 2009.
Climate change, sea level rise scenarios for Vietnam. Hanoi:
Ministry of Natural Resources and Environment.
Monteith, J.L., 1965. Evaporation and the environment. In: G.E. Fogg,
ed. The state and movement of water in living organisms. XIXth
symposium on the society of experimental biology (Swansea,
UK). Cambridge: Cambridge University Press, 205–234.


Downloaded by [NUS National University of Singapore] at 22:40 03 June 2014

1108

Monthly Weather Review, 100, 81–92. doi:10.1175/1520-0493
(1972)100<0081:OTAOSH>2.3.CO;2.
Ranzi, R., Le, T.H., and Rulli, M.C., 2012. A RUSLE approach to
model suspended sediment load in the Lo River (Vietnam): effects
of reservoirs and land-use changes. Journal of Hydrology, 422–
423, 17–2929. doi:10.1016/j.jhydrol.2011.12.009.
Rossi, C.G., et al., 2008. Hydrologic calibration and validation of the
soil and water assessment tool for the Leon River watershed.
Journal of Soil and Water Conservation, 63, 533–541.
doi:10.2489/jswc.63.6.533.
SIWRP (Southern Institute for Water Resources Planning), 2002.
Integrated water resources planning for Be River catchment.
Ho Chi Minh city: SIWRP (in Vietnamese).

SIWRP (Southern Institute for Water Resources Planning), 2008.
Integrated water resources planning for Dong Nai River
basin. Ho Chi Minh city: SIWRP.
Tong, S.T.Y., et al., 2012. Predicting plausible impacts of sets of climate
and land use change scenarios on water resources. Applied
Geography, 32, 477–489. doi:10.1016/j.apgeog.2011.06.014.
Trinh, M.V., 2007. Soil erosion and nitrogen leaching in northern
Vietnam—experimentation and modeling. Thesis (PhD).
Wageningen University.
Tu, J., 2009. Combined impact of climate and land use changes on
streamflow and water quality in eastern Massachusetts, USA.
Journal of Hydrology, 379, 268–283. doi:10.1016/j.
jhydrol.2009.10.009.
USDA-SCS (US Department of Agriculture—Soil Conservation
Service), 1972. National engineering handbook. Section 4: hydrology. Washington, DC: US Department of Agriculture, United State.
Van Griensven, A., et al., 2006. A global sensitivity analysis tool

China. Journal of Hydrology, 410, 239–247. doi:10.1016/j.
jhydrol.2011.09.023.