Open access peer-reviewed chapter

Hydrological Trend Analysis Integrated with Landscape Analysis at the Watershed Scale (Case Study: Langat Basin, Malaysia)

By Hadi Memarian and Siva K. Balasundram

Submitted: November 23rd 2015Reviewed: February 11th 2016Published: July 27th 2016

DOI: 10.5772/62463

Downloaded: 803

Abstract

In this study, the trends of water and sediment data collected from three hydrometer stations over the past 25 years of development in the state of Selangor, Peninsular Malaysia, were analyzed using the Mann–Kendall and Pettitt’s tests. Landscape metrics for establishing the relationship between land use changes and trends of hydrological time series were calculated. The hydrologic trends were also studied in terms of rainfall variations and man-made features. Results indicated upward trends in water discharge at the Hulu Langat sub-basin and sediment load at the Semenyih sub-basin. These increasing trends were mainly caused by rapid changes in land use. Upward trends of hydrological series at the Hulu Langat sub-basin matched its rainfall pattern. At the Lui sub-basin, however, trends of hydrological series and variations in rainfall and land use were not statistically significant.

Keywords

  • trend analysis
  • Mann–Kendall test
  • Pettitt’s test
  • landscape metrics
  • water discharge
  • sediment load

1. Introduction

Globally, increased sediment load (SL) and intense flooding due to land use changes at river basins are very challenging problems [14]. The impacts of human activities and climate change on hydrological processes occurring at river systems are well documented [510]. Understanding time series trends of water discharge (WD) and SL can be a key solution in determining how hydrological systems are affected by climate change and anthropogenic disturbances [4].

Zhang et al. [4] determined time series changes in water and sediment discharge at the Zhujiang (Pearl River) Basin in China. They applied Mann–Kendall (MK) as a gradual trend test and Pettitt’s as an abrupt change test on annual WD and SL from 1950 to 2004 at nine hydrometer stations. Their study showed that long-term changes in annual WD, which were originally controlled by variation in precipitation, were not significant. SL at all main hydrometer stations showed declining trends during the study period. Declining trends were principally influenced by the construction of reservoirs and dams. In a mountainous tributary of the lower Xinjiang in China, hydrological response to changes in precipitation and anthropogenic activities were tested using MK and Pettitt’s tests [11]. The power of MK and Spearman’s rho tests to assess the significance of hydrological trends has been studied by Yue et al. [12]. They demonstrated that the power of both tests is directly proportional to trend slope, sample size, and predetermined significance level and inversely proportional to time series variation.

Ouyang et al. [6] established a relationship between soil erosion and landscape metrics at the Logliu catchment in China. They showed that landscape pattern significantly impacted soil erosion and sediment transportation. In several other studies, landscape metrics were applied at the landscape and patch levels to determine how hydrological conditions of the basin are affected by human activities such as land use change [1317].

In recent decades, the Langat Basin has experienced rapid development towards urbanisation, industrialisation, and intense agriculture [18]. The Langat Basin is also a main source of drinking water for surrounding areas and a source of hydropower and has an important role in flood mitigation. Over the past four decades, the Langat Basin has served approximately 50% of the Selangor State population. However, the Selangor State is currently facing water shortage problems, especially in urban areas [19, 20].

This study was conducted to assess the impact of land use change, rainfall variation, and other anthropogenic manipulations on hydrological trends in selected upper catchments within the Langat Basin over a period of 25 years.

Advertisement

2. Methodology

2.1. Study area

The Langat Basin is located at the southern part of Klang Valley, which is the most urbanised river basin in Malaysia. It is believed that the Langat Basin is currently experiencing “spillover” effects due to the excessive development in the Klang Valley. Hydrometeorologically, the Langat Basin is affected by two types of monsoons, i.e., the Northeast (November–March) and the Southwest (May–September) [21]. The average annual rainfall is approximately 2400 mm. The wettest months are April and November, with an average monthly rainfall exceeding 250 mm, whereas the driest month is June, with an average monthly rainfall not exceeding 100 mm. Topographically, the Langat Basin can be divided into three distinct areas in reference to the Langat River, i.e., mountainous area in the upstream, undulating land in the centre, and flat flood plain in the downstream (Figure 1). The Langat Basin consists of a rich diversity of landform, surface feature, and land cover [21, 22].

Sub-basinLuiHulu
Langat
Semenyih
Main riverLuiLangatSemenyih
Geographic
coordinate
3°07′–3°12′N and
101°52′–101°58′E
3°00′–3°17′N and
101°44′–101°58′E
2°55′–3°08′N and
101°49′–101°58′E
Drainage area
(km2)
68.25390.26235.62
Basin length
(km)
11.534.526.5
Average slope
(%)
3529.427.4
Max. altitude
(m)
120714791070
Min. altitude
(m)
612021
Ave. altitude
(m)
354277.4243.9
Ref. hydrometer
station
Sg. LuiSg. LangatSg. Semenyih
WD
(×106 m3/y)
55.05289.64146.11
SL
(×103 ton/y)
5.88146.636.81
Runoff
(mm/km2/y)
806.57742.16620.11
Sediment yield
(ton/km2/y)
86.22375.65156.22
Ref. rainfall
station
Kg. LuiUPM SerdangLdg. Dominion
Precipitation
(mm)
2188.324532548.8
Land covers*Forest 80.35%, cultivated
rubber 9.85%, orchards
2.6%, mixed horticulture
and crops,
urbanised area,
and mining land 7.2%
Forest 54.6%, cultivated
rubber 15.6%, orchards
2%, urbanised area 15%,
horticulture and crops,
oil palm, lake,
and mining land 12.8%
Forest 53.8%, cultivated rubber
17.4%, oil palm 6.3%, urban
area 5.6%, secondary forest 3.6%,
scrub land 2.4%, mining, other crops,
mixed horticulture, orchard, bare land,
marshland and aquaculture 10.9%

Table 1.

General information of the studied sub-basins.

*Based on the land use map dated 2006.


Based on the availability of hydrometric stations in the Langat Basin, three sub-basin (upstream of the Langat River) were investigated. The descriptions about these sub-basin are given in Table 1.

2.2. Data set

WD, SL, and precipitation data between 1984 and 2008 recorded at all three hydrometer and rain gauge stations under study (Table 1) were obtained from the Department of Irrigation and Drainage (DID) of Malaysia. The geographic location and general information of the hydrometer stations are presented in Figure 1 and Table 2. Land use maps dated 1984, 1988, 1990, 1995, 1997, 2001, 2002, 2005, and 2006 were obtained from the Soil Resource Management and Conservation Division, Department of Agriculture, Malaysia.

Figure 1.

Geographic locations of the three study sub-basins.

2.3. Trend analysis

In this study, non-parametric tests, such as MK and Pettitt’s, were used to detect gradual and abrupt changes in the hydrological data sets. According to Zhang et al. [4], non-parametric tests are preferred over parametric tests due to their strength in handling non-normally distributed data and missing data. The MK equation that is based on the Sstatistic is as follows:

S=i=1n-1j=i+1nsgnxj-xiE1

where xiand xjare sequential data values, nis the length of time series, and

sgn(θ)={1   if θ>00  if θ=01  if θ<1E2

Mann (1945) and Kendall (1975) (as cited by Yue et al. [12]) have posted that, when n≥ 8, Sis almost normally distributed with the following mean and variance:

ES=0E3
VS=nn-12n+5-i=1ntiii-12i+518E4

where tiis number of ties of the extent i.

The standard Zstatistic is calculated as follows:

ZMK=S-1Var(S)S>00S=0S+1Var(S)S<0E5

ZMK pursues the standard normal distribution with μ=0 and δ=1.

The probability (P) of the Sstatistic is estimated by normal cumulative distribution function as follows:

P=0.5Φ|Z|E6
ΦZ=12π0Ze-t22dtE6

Statistical significance of the data trend was based on 95% confidence level [12].

The MK test is not robust against autocorrelation [4, 11, 23]. As such, an autocorrelation test was performed on the data set to determine the degree of autocorrelation. The autocorrelation coefficient is estimated as follows [24]:

rk=1n-kt=1n-kXt-X¯Xt+k-X¯1n-1t-1nXt-X¯E7

where X=t=1nXtn, nis the sample size, and kis the lag.

For a completely random series, rk≈ 0 for all k≠ 0. If a series of rk(for k≠ 0) fall between the 95% confidence level estimated by ul=-1±Z1-/2n-2/n-1(where nis the length of tested time series, land uare the lower and upper limits, αis the significance level, and Zis the critical value of standard normal distribution for a given α), then the tested series will be independent at the 95% confidence level [23]. The WD and SL data showed significant autocorrelation. Therefore, Zhang’s method of data pre-whitening [25] was used to eliminate significant autocorrelation within the data.

Sen’s non-parametric method was used to estimate the change magnitude (i.e., slope of the linear trend). Sen’s method is robust against non-normally distributed data, missing values, and extreme outliers (Sen, 1968) (as cited by Zhang and Lu [11]).

Considering a sequence of random variables X1, X2, …, Xt, which have a change point at τ[Xtfor t= 1, 2, …, τhave a common distribution function F1(x) and Xtfor t=τ+1, …, Thave a common distribution function F2(x) and F1(x) ≠ F2(x)], Pettitt’s test (1979) (as cited by Zhang et al. [4] and Wolfe and Schechtman [26]) was used to detect one unknown change point in the pre-whitened WD and SL time series. In the Pettitt’s test, null hypothesis (H0): no change (τ=T) is tested against alternative hypothesis (Ha): change (1 ≤ τ<T)) by the non-parametric Kstatistic, as follows:

Kt=max1tTUt,T=maxKT+,KT-E8

where Ut,T=i=1tj=t+1Tsgn(XtXj),sgn(θ)={1(ifθ>0)0(ifθ=0)1(ifθ<0). KT+=max1tTUt,Tfor downward shift and KT-=-min1tTUt,Tfor upward shift. The significance level of KT+or KT-is estimated by P=exp-6KT2T3+T2. When KToccurs, the time twill be the point of change. When Pis smaller than the specific significance level, H0 is rejected.

The above procedures were performed using XLSTAT and R statistical packages.

2.4. Landscape analysis

To assess the changes in land use patterns over the period 1984–2006 (including nine records), Patch Analyst 3.0 (Grid) extension in ArcView 3.3 was applied to calculate landscape metrics [27], which are fundamental indices for the detection of trends in land use change [6]. A brief description of the six selected landscape metrics for this work is given below [28, 29]:

The number of patches (NUMP) is ≥1.

The patch size coefficient of variation (PSCOV) is the variability of the patch size relative to the mean patch size.

PSCOV=PSSDMPS100E9

where PSCOV 0, PSSD is the standard deviation in patch size and MPS is mean patch size of the corresponding patch type.

Edge density (ED) equals the sum of lengths (m) of all edge segments involving the corresponding patch type divided by the total landscape area (m2) in meters per hectare.

The Shannon’s diversity index (SDI; at the landscape level) is a measure between 0 and 1. SDI equals 0 when the landscape comprises only one patch (i.e., no diversity) and increases with increasing number of patch types. SDI equals 1 when the different patch types are distributed proportionally.

SDI=-i=1mPi.lnPiE10

where Piis proportion of the landscape occupied by the patch type (class) i.

The Shannon’s evenness index (SEI; at the landscape level) is a measure between 0 and 1. SEI equals 0 when the landscape comprises only one patch (i.e., no diversity) and approaches 0 as the areal distribution of patch types becomes uneven (i.e., dominated by one type). When the areal distribution of patch types is perfectly even, SEI equals 1.

SEI=-i=1mPi.lnPilnmE11

where mdefines the number of patch types (classes) present in the landscape including the landscape border.

The interspersion and juxtaposition index (IJI) is the observed interspersion over the maximum possible interspersion for a given number of patch types. IJI approaches 0 when the corresponding patch type is adjacent to only one other patch type and is 100 when the corresponding patch type is equally adjacent to all other patch types.

IJI=-k=1mk=1eikmeiklnk=1eikmeiklnm-110012E12

where eikis the total length (m) of edge in landscape between the patch types (classes) iand k, and m′ defines the number of patch types (classes) presented in the landscape.

3. Results

3.1. Hydrological trend analysis

Figure 2.

Autocorrelograms, resulted from autocorrelation test on WD and SL at the selected hydrometer stations.

The autocorrelation test reveals that the WD series (except that at Sg. Lui hydrometer station) and SL series have at least one autocorrelation coefficient that is significant at the 95% confidence level (Figure 2). The autocorrelation coefficients of the WD series at both Sg. Langat and Sg. Semenyih hydrometer stations are significant at the first lag. The SL series at Sg. Langat and Sg. Semenyih hydrometer stations have one significant autocorrelation coefficient at the second and first lags, respectively. The SL series at Sg. Lui hydrometer station is autocorrelated significantly only at the second lag.

The results of gradual trend analysis based on MK and Pre-Whitening MK (PWMK) tests for WD and SL data are shown in Table 2. At Sg. Langat, WD shows an increasing trend that is significant at the 95% confidence level; however, the ascendant trend of SL is not statistically significant. There are no significant trends in the hydrological time series of Sg. Lui and Sg. Semenyih.

Gradual changes in the hydrological time series at all three hydrometer stations are given in Table 2. At Sg. Lui, WD and SL are decreasing at a rate of 0.524×106 m3/y and 0.119×103 ton/y, respectively. At Sg. Langat, however, WD and SL are increasing at a rate of 9.899×106 m3/y and 1.415×103 ton/y. At Sg. Semenyih, WD and SL show declining tendencies at 3.686×106 m3/y and 0.316×103 ton/yr.

Station nameParameterMK and PWMK trend testsSen’s slope estimator
τPTrendTrendTrend_PLinearIntercept
Sg. LuiWD−0.1530.297Decreasing−0.524−13.105−0.52562.723
SL−0.0720.637Decreasing−0.119−2.986−0.4714.596
Sg. LangatWD0.3260.027Increasing9.899247.48511.531156.869
SL0.1300.385Increasing1.41535.38013.14246.278
Sg. SemenyihWD−0.1960.189Decreasing−3.686−92.145−3.912187.036
SL−0.0580.710Decreasing−0.316−7.909−3.61120.707

Table 2.

Results of MK and PWMK tests with the Sen’s slope estimator (at α=0.05) applied on WD and SL (data in bold are significant).

Trend: Sen’s slope (trend) per unit time; Trend_P: Sen’s slope (trend) over the time period; Linear: least-squares fit trend; Intercept: intercept of the Sen’s slope (trend).


The results of abrupt changes based on the Pettitt’s test for WD and SL are shown in Figures 3 and 4 and Table 3. The results show significant drastic changes in the hydrological time series at Sg. Langat and Sg. Semenyih hydrometer stations. At Sg. Langat, the mean level of WD (after 1998) shifted upward to 392.09×106 m3/y, which corresponds to a 77% increase, whereas the mean level of SL (after 1999) shifted upward to 297.27×103 ton/y, which corresponds to a 380% increase (Figure 3). At Sg. Semenyih, the mean level of WD (after 1993) shifted downward to 111.18×106 m3/y (44% decrease) and the mean level of SL shifted downward to 14.89×103 ton/y (78% decrease; Figure 4). At Sg. Lui, however, downward shifts in hydrological time series are not significant.

Figure 3.

Abrupt changes in the mean level of WD and SL for Sg. Langat at the significance level of 0.05.

Figure 4.

Abrupt changes in the mean level of WD and SL for Sg. Semenyih at the significance level of 0.05.

Station nameParameterKTPShiftT
Sg. LuiWater56.0000.430Downward1987
Sediment82.0000.084Downward1988
Sg. LangatWater104.0000.012Upward1998
Sediment108.0000.009Upward1999
Sg. SemenyihWater112.0000.010Downward1993
Sediment130.0000.005Downward1993

Table 3.

Results of Pettitt’s test applied on WD and SL (data in bold are significant at the level of 0.05).

Based on the PWMK test, the gradual change in the Lui and Semenyih sub-basin are showing declining tendencies. It is important to know whether these declining tendencies will prevail after the change point. Table 4 shows that all hydrological series have increasing tendencies after the change point, which are statistically significant only for SL series at Sg. Semenyih. Meanwhile, decreasing tendencies before the change point are statistically significant only for SL series recorded at Sg. Langat.

StationParameterTTime seriesMK trend test
τPTrend
Sg. LuiWD1987Pre-T*
Post-T0.0570.739Increasing
SL1988Pre-T−0.6660.308Decreasing
Post-T0.1840.269Increasing
Sg. LangatWD1998Pre-T−0.3310.112Decreasing
Post-T0.4050.127Increasing
SL1999Pre-T0.5820.004Decreasing
Post-T0.3880.175Increasing
Sg. SemenyihWD1993Pre-T0.0001.000
Post-T0.1430.511Increasing
SL1993Pre-T−0.3570.265Decreasing
Post-T0.4500.028Increasing

Table 4.

Results of the MK and PWMK tests on WD and SL before and after the change points (data in bold are significant at the level of 0.05).

*Limitation in number of records.


3.2. Landscape analysis

The relationships between landscape metrics and hydrological variables are given in Table 5. At the Lui sub-basin, all metrics with the exception of IJI and SEI are negatively correlated with WD and SL. SL correlates significantly with PSCOV and SEI, whereas WD correlates significantly with ED and NUMP.

At the Semenyih sub-basin, all metrics are negatively correlated with WD and SL. Correlations between these metrics and WD are statistically significant. Correlations between these metrics with the exception of IJI and PSCOV and SL are statistically significant. At the Hulu Langat sub-basin, on the contrary, all metrics are positively correlated with WD and SL but are not statistically significant.

In comparison, correlations between hydrological series and landscape metrics are more pronounced after the change point (Table 5). For example, at the Semenyih sub-basin, correlations between hydrological series and all metrics change from negative to positive and are only significant for PSCOV. Also, at the Hulu Langat sub-basin, correlations between WD and the metrics ED, SDI, and NUMP are statistically significant.

Table 6 shows the trend analysis of landscape metrics during the period 1984 to 2006. At the Lui sub-basin, only NUMP shows a significant increasing trend. At the Hulu Langat sub-basin, all metrics show increasing tendencies but are statistically significant only for PSCOV, ED, SDI, and NUMP. At the Semenyih sub-basin, all metrics with the exception of IJI show significant increasing trends.

Landscape
metric
All recordsRecords after the change points
LuiHulu LangatSemenyihLuiHulu LangatSemenyih
WDSLWDSLWDSLWDSLWDSLWDSL
PSCOV−0.625−0.788*0.3580.2430.921**−0.584−0.4080.0810.4600.6740.887*0.932**
ED0.680*−0.4710.2890.2440.957**0.790*−0.549−0.0260.941*0.7840.5520.548
IJI0.0910.2380.2060.1910.973**−0.6600.0480.1370.6740.6590.3800.575
SDI−0.556−0.0990.2570.2280.885**0.865**−0.532−0.0260.981*0.8460.7890.670
SEI0.4740.951**0.3400.3790.940**0.873**−0.1500.1680.4810.6460.3890.398
NUMP0.703*−0.4810.3620.5230.954**0.704*−0.584−0.0250.955*0.8700.7620.759

Table 5.

Correlations between the different landscape metrics and hydrological series using the Pearson correlation method.

*Correlation is significant at the 0.05 level.


**Correlation is significant at the 0.01 level.


Landscape metricLuiHulu LangatSemenyih
PSCOVNS ↑* ↑* ↑
EDNS ↑* ↑* ↑
IJINS ↓NS ↑NS ↑
SDINS ↑* ↑* ↑
SEINS ↓NS ↑* ↑
NUMP* ↑* ↑* ↑

Table 6.

Trend analysis of the landscape metrics during 1984 to 2006 at the studied sub-basins.

*Trend is significant at the 0.05 level.

NS, not significant. ↑, increasing; ↓, decreasing.

The categorisation of the landscape metrics using the clustering technique for all three sub-basins is given in Figures 57.

At the Lui sub-basin, PSCOV, ED, and NUMP corresponding to years 1984, 1988, and 1990 are categorised in the first cluster, whereas landscape metric values of the remaining years are classified in the second cluster. Meanwhile, SEI (1984 and 1988) and SDI (1988 and 1990) are categorised in the first cluster, whereas the values of the remaining years are classified in the second cluster (Figure 5).

Figure 5.

Change trends and classification of the landscape metrics at the Lui sub-basin during 1984 to 2006.

All landscape metric values at the Hulu Langat (except for IJI) and Semenyih sub-basins corresponding to years 1984, 1988, and 1990 are grouped in the first cluster and the values of the remaining years are grouped in the second cluster (Figure 6).

Figure 6.

Change trends and classification of the landscape metrics at the Hulu Langat Sub Basin during 1984-2006.

Figure 7.

Change trends and classification of the landscape metrics at the Semenyih sub-basin during 1984 to 2006.

3.3. Land use change detection

Change detection was based on the land use maps dated 1984 and 2006. The results show noticeable gains in agriculture, bare land, mining, oil palm, and urban acreage and a remarkable loss in rubber acreage (Table 7).

Land useDifference (2006–1984) in hectares
LuiHulu LangatSemenyih
Agriculture436.391545.32451.57
Bare land6.30370.60175.50
Forest359.75112.04-237.72
Grassland5.00−471.70382.30
Marshland/swamp46.42148.93-39.73
Mining16.67395.27147.57
Oil palm3.4421.811051.89
Rubber−809.43−7817.72-3860.35
Urban/built-up area111.095649.901501.61
Water body0.0045.55427.36

Table 7.

Land use change detection between 1984 and 2006 at the studied sub-basins.

Increasing gradual trends (as determined by the MK test using the 1984–2006 time series) in agriculture and mining acreages are significant only at the Hulu Langat sub-basin, whereas increasing gradual trend in oil palm acreage is significant only at the Semenyih sub-basin. Decreasing gradual trend in rubber acreage is significant at both Hulu Langat and Semenyih sub-basins. At all three sub-basins, a change trend in forest area is not significant, whereas increasing trend in the urbanised area is significant (Table 8).

Land useLuiHulu LangatSemenyih
AgricultureNS ↑* ↑NS ↑
ForestNS ↓NS ↓NS ↓
Urban* ↑* ↑* ↑
RubberNS ↓* ↓* ↓
Oil palmNS ↑NS ↑* ↑
MiningNS ↑* ↑NS ↑

Table 8.

Trend analysis of land use change during 1984–2006 at the studied sub-basins.

*Trend is significant at the 0.05 level.

NS, not significant. ↑, increasing; ↓, decreasing.

4. Discussion

Based on the hydrological trend analysis, the Hulu Langat sub-basin showed a significant increasing trend of WD. However, the Semenyih and Lui sub-basins showed decreasing tendencies of WD and SL. Gradual increase in hydrological series after the change points is significant only in the case of SL at Sg. Semenyih. In the following sections, hydrological alterations are discussed in relation to land use change, rainfall fluctuations, and other anthropogenic manipulations.

4.1. Effect of land use/cover change (LUCC)

Based on the data from the Department of Statistics, Malaysia (2001) and the National Urbanisation Policy of Malaysia (1981), rapid development in the state of Selangor started in 1981. The rapid development was aimed at attracting approximately 18% of Malaysia’s population to be settled in the state of Selangor by the year 2000 [30, 31]. The Langat Basin appears as a suitable barometer to measure urbanisation and agricultural/industrial development in the state of Selangor.

Based on Tables 7 and 9, rubber acreage at the Hulu Langat and Semenyih sub-basins decreased significantly between 1984 and 2006. During this period, at the Hulu Langat sub-basin, 34% of rubber acreage was transformed into urban areas, whereas another 11% was used for other agricultural production. Similarly, at the Semenyih sub-basin, 21% of rubber acreage was transformed into urban areas, whereas 15% and 6% were used for oil palm and other agricultural productions, respectively (Table 7). Based on these results and the work of Noorazuan et al. [22] and Juahir et al. [19], it is expected that these changes will affect the stream flow behaviour and characteristics.

19842006
Lui
ForestRubberAgricultureScrub/idle
grassland
Oil
palm
Urbanised/
industrial area
Water
bodies
Forest111
Rubber215
Scrub/idle
grassland
4711334
Hulu Langat
Forest43.5
Rubber1134
Agriculture1726
Scrub/idle
grassland
38171031
Semenyih
Forest2.51.42.7
Rubber61521
Agriculture7105615
Scrub/idle
grassland
701019
Oil palm8834
Swamp/
marshland
78013
Bare land67

Table 9.

Land use change matrix for important transitions (frequencies in %) between the years 1984 and 2006 at the studied sub-basins.

Tables 6 shows the increasing trend in landscape change during 1984–2006, which is confirmed by correlation analysis (Table 5), especially after the change points [6, 32].

At the Hulu Langat sub-basin, cluster analysis shows that discriminant points between the clusters of landscape metrics (except for IJI) are within the period 1990–1997 (Figure 6). These points correspond to the points of change in WD and SL. At the Semenyih sub-basin, the change point in WD and SL (i.e., 1993) matches the discriminant point in all landscape metric clusters, with the exception of SDI (Figure 7). At the Lui sub-basin, points of change in WD and SL (1987 and 1988) match only the discriminant point in SEI. As indicated in Table 6, at the Lui sub-basin, change trends in landscape metrics are not statistically significant. This could have contributed to the insignificant impact of LUCC on the basin hydrological conditions.

4.2. Effect of rainfall variations

The rainfall stations (i.e., Kg. Lui, UPM Serdang, and Ldg. Dominion) were analysed for rainfall change trend. These three rainfall stations represent rainfall events at the Lui, Hulu Langat, and Semenyih sub-basins, respectively (Table 10). Increasing trend in the rainfall time series is only significant at UPM Serdang, which corresponds to increasing trends in WD and SL at Sg. Langat. Although the change point in rainfall series at UPM Serdang (i.e., 1998) is not statistically significant, it matches the change point in WD and SL at Sg. Langat. This point matches the critical water level at Langat Reservoir [33], which has been reported by Shaaban and Low [34]. At Kg. Lui, the change tendency of rainfall series after the hydrological change point matches the trend of WD and SL at Sg. Lui (Table 11). At Ldg. Dominion, the change tendency in rainfall does not match the available tendency in hydrological series at Sg. Semenyih, especially for SL after the change point.

Station nameMK trend testPettitt’s test
τPTrendKTPShiftT
Kg. Lui0.1440.333Increasing70.0000.188Upward2003
UPM Serdang0.3410.020Increasing78.0000.079Upward1998
Ldg. Dominion0.1010.503Increasing68.0000.231Upward1990

Table 10.

Results of the PWMK and Pettitt’s tests applied on the rainfall time series at the representative stations (data in bold are significant at α=0.05).

x
Rainfall stationT_Hydro seriesTime seriesPWMK trend test
τPTrend
Kg. Lui1988Pre-T−0.3330.734Decreasing
Post-T0.1110.528Increasing
UPM Serdang1998Pre-T0.2560.246Increasing
Post-T0.0660.858Increasing
Ldg. Dominion1993Pre-T0.4440.117Increasing
Post-T0.0001.000

Table 11.

Results of PWMK test applied on the rainfall time series before and after the hydrological change points.

Thus far, results reveal the significant impacts of land use and rainfall variations on WD at the Hulu Langat sub-basin. However, the impacts on the sub-basin SL are not statistically significant. At the Lui sub-basin, the change trends in rainfall and landscape variables are not statistically significant. Hence, the insignificance of hydrological series trend is expected. At the Semenyih sub-basin, the impact of land use change on hydrological series is driven by the significant increasing trend in SL after 1993. However, rainfall variations do not impact the trend of hydrological series.

From the preceding discussion, two questions are important. First, have the Semenyih Reservoir and its connected water treatment facilities at the Semenyih sub-basin impacted the basin WD significantly? Secondly, despite the significant impact of land use change on the change trend in WD at the Hulu Langat sub-basin, why is the change trend in SL not statistically significant? In the following discussion, these questions are addressed.

4.3. Effect of man-made structures

There are two strategic dams in the Langat Basin. The Langat Dam, constructed in 1979, has a drainage catchment area of 41.5 km2 and a reservoir capacity of 37.5 Mm3. The Semenyih Dam, built in 1985, has a drainage catchment area of 56.7 km2 and a reservoir capacity of 62.6 Mm3. Both these dams supply domestic and industrial water. The Langat Dam is also used to generate power supply at moderate capacity for consumption within the Langat Valley. Currently, there are three major water treatment plants (WTP; operating 24 hours a day) within the study area. The Sg. Langat and Cheras WTPs along the Langat River produce 386.4 and 27 million litres per day (MLD) of clean water, respectively. The Semenyih WTP along the Semenyih River produces 545 MLD of clean water [30, 33].

Figure 8.

Cumulative double mass plot at Sg. Semenyih.

To evaluate the impact of Semenyih Dam construction on trends of hydrological time series, a double mass curve was plotted. As illustrated in Figure 8, at Sg. Semenyih, the hydrological time series trend after 1985 is seriously affected by the dam and WTP construction. The mean WD level changes from 231.8×106 m3/y before 1985 to 141×106 m3/y after 1985. The results of the Pettitt’s test for WD and SL during the period 1975–1993 were statistically significant (P=0.018 and P<0.001, respectively). This confirms 1985 as the point of change during the period 1975 to 1993.

Shaaban and Low [34] showed that drought events reduced WD at the Semenyih sub-basin, particularly in the period 1993–1998. As such, WTP and dam together with the effect of drought have been able to reduce the increasing trend of WD, especially after 1993.

At the Hulu Langat sub-basin, due to the significant trend in urbanisation and agricultural activities (Tables 7 and 9), the number and size of natural or artificial ponds are expected to increase dramatically. Field observation from this study confirms that the quantity of natural and artificial ponds is higher at Hulu Langat compared to that at Semenyih. The ponds are believed to affect the sedimentation process by increasing the deposition rate, hence resulting in the reduction in SL of the basin (Figure 9). It is clear that the Langat Dam and other sediment trapping features (i.e., natural and artificial ponds) contributed to the insignificant trend of SL at the Hulu Langat sub-basin.

Figure 9.

Ponds arisen from urban and agricultural development at the Hulu Langat sub-basin (extracted from SPOT 5 satellite images, dated 2006).

5. Conclusion

Increasing trend in WD at the Hulu Langat sub-basin was originally controlled by significant variations in land use and rainfall. However, increasing trend in SL was not significant due to dam construction and increase in the number and size of sediment trapping features, which is due to urbanisation and agricultural activities. At the Semenyih sub-basin, increasing trend in SL after 1993 was closely related to significant trends in landscape metrics and land use changes. However, WD did not increase significantly after 1993, primarily due to the impact of dam and water treatment facilities and continuous drought until 1998. At the Lui sub-basin, trends in land use and rainfall variations were mostly insignificant, thus causing an insignificant change in hydrological series.

This study demonstrates the power of the PWMK and Pettitt’s tests for trend evaluation of hydrological time series. The results obtained in this work are consistent with studies done by other researchers [3545]. Also, integrating landscape analysis with statistical analyses as emphasised in the work of several others [6, 1317, 32, 46, 47] could increase the depth of interpretation with regard to the complex hydrological conditions of developed basins.

Advertisement

Acknowledgments

The authors acknowledge Universiti Putra Malaysia for procuring land use maps and satellite imagery and the Department of Irrigation and Drainage, Malaysia, for supplying hydrological data.

© 2016 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Hadi Memarian and Siva K. Balasundram (July 27th 2016). Hydrological Trend Analysis Integrated with Landscape Analysis at the Watershed Scale (Case Study: Langat Basin, Malaysia), Landscape Ecology - The Influences of Land Use and Anthropogenic Impacts of Landscape Creation, Amjad Almusaed, IntechOpen, DOI: 10.5772/62463. Available from:

chapter statistics

803total chapter downloads

1Crossref citations

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

The Anthropic Pressure on the Landscapes of the Subcarpathian and Piedmont Basin of Dâmboviţa River (Romania)

By Mihaela Sencovici and Gica Pehoiu

Related Book

First chapter

Introductory Chapter: A General Reading Process on Landscape Architecture

By Amjad Almusaed

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us