Abstract

Variability and availability of water resources under changing environment in a regional scale have been hot topics in recent years, due to the vulnerability of water resources associated with social and economic development. In this paper, four subbasins in the Dongjiang basin with a significant land use change were selected as case study. Runoffs of the four subbasins were simulated using the SCS monthly model to identify the quantitative impacts of land use and climate change. The results showed that (1), in the Dongjiang basin, temperature increased significantly, evaporation and sunlight decreased strongly, while precipitation showed a nonsignificant increase; (2) since the 1980s, land uses in the Dongjiang basin have experienced a significant change with a prominent increase in urban areas, a moderate increase in farmlands, and a great decrease in forest areas; (3) the SCS monthly model performed well in the four subbasins giving that the more significant land use change in each subbasin, the more runoff change correspondingly; (4) overall, runoff change was contributed half and half by climate change and human activities, respectively, in all the subbasins, in which about 20%~30% change was contributed by land use change.

1. Introduction

Environment Change (including land use change and climate change) and its impacts on water resources have always been the hot issues in recent years. The direct or indirect impacts on the hydrological regime brought by land use and climate change both have contributed to some water problems, such as water shortage, flooding, and water logging to different extent. Some researches have been conducted to study the impacts of land use and climate change on water resources in different basins [15]. Particularly in the humid region of south China, six hydrological models were used to simulate the hydrological impact of some climate change scenarios, and response strategies for water supply and flood control due to climate change were analyzed in the Dongjiang basin by Jiang et al. [6]. South China has identified an increasing trend of extreme rainfall events, and the correlation between such events and flood events was studied by Fu et al. [7]. Climate change makes flooding more frequent in some regions [8, 9]. Additionally, land use change also has notable impact on the hydrological cycle [1012]. Overall, hydrological process and the variability and availability of water resources would change a lot due to land use change and climate change.

As the biggest developing country in the world, China has experienced an explosive economic growth over a couple of decades, which results in the prominent land use change. This change can alter evaporation patterns and potentially affect water resources in a region. The impacts of climate change on water resources and agriculture in China were analyzed by Piao et al. [13], and the results indicated that agriculture relying on the resultant increase in glacier runoff, especially in western China, would face a challenge, although its climate trends remained moderate compared to natural variability [13]. Even worse was the more frequent occurrence of drought and flood extremes largely due to climate change in different parts of China [14]. Some researches have focused on the variability of runoff in the Dongjiang basin under the impacts of climate change and human activities, especially land use change in recent decades [15]. Definitely, there is an increasingly serious challenge in the availability of water resources, and its exploitation and management will be more or less affected.

As to a basin, climate change and human activities both contribute to the hydrological cycle, and this finding has been supported by many studies [1619]. Many studies have been carried out to distinguish the roles of land use and climate change on water resources. The approach to the calculation for the runoff change due to precipitation and potential evaporation was proposed [2022]. Based on which, the separated impacts of human activities and climate change on natural runoff change in the Poyang basin from 1992 to 2000 were analyzed quantitatively by Ye et al. [23]. In recent years, much attention has been paid to analyze the separated quantitative effects of land use and climate change, and some progress has been made, although further research still needs to be carried out. A simple approach to distinguish land-use and climate-change effects on water resources was developed with a coupled water-energy budget analysis in US Midwest area by Tomer and Schilling [11]. Combined impacts of land cover and climate changes on hydrological processes of the Kejie watershed in the eastern Himalayas were assessed with a SWAT model, which revealed that land-cover change had more effects than those of climate change in the short and middle terms [24]. Furthermore, in many river basins of China, many studies have been conducted to assess the impacts of climate change and human activities on the runoff variation. As to Tarim River, the impact of human activities had increased the runoff with a ratio of 41%~75% for different ages [25], and similar research can be found in the Wuding River [26]. In addition, some relative studies also took place in the Heihe catchment, Chaobai River, and Mian River basin [2729], which came to a conclusion that climate change and human activities have the separated impacts on the runoff change varying from place to place. There is no exception for the Dongjiang basin [30]. However, few researches have been conducted to identify the quantitative effects of human activities and climate change on runoff in the study area. Moreover, most of the previous studies focused on a large scale basin and had no cases for comparison, which tends to limit our understanding of the impacts of climate change and human activities on water resources in the study area. The quantitative effect of land use change on runoff of the basin has not been revealed yet. Therefore, it is desirable to separate the impacts of climate change and human activities, specifically the land use change, on water resources under a changing environment in the study area.

Hydrological models have been widely used to study many practical and pressing issues that arise during planning, design, operation, and management of water resources systems [31, 32] and also to quantify the impacts of land use and climate change on the hydrological cycle. For example, a conceptual rainfall-runoff model was applied to the Rhine basin for the purpose of modeling the effect of land use change on the runoff. The results suggested that increased urbanization led to an increase in the lower peak runoff, compared to a considerable reduction of both the peak runoff and the total runoff volume resulted from intensified forestation [33]. Besides, the hydrological model considering other elements such as biogeochemistry was applied to reveal how land use change affects hydrological regimes at the watershed scale [34]. According to the available climate and hydrological data and the hydrological characteristics of the basin studied, many conceptual or distributed hydrological models were introduced to analyze the hydrological response under land use and climate change [3538]. Trend analysis of climate variables is necessary to detect the variability of climate variables, such as temperature and precipitation [39, 40], which provides supports for study of the impacts on the hydrological cycle of land use change and climate change. Overall, hydrological model combined with statistic methods has been a prevalent and useful tool to clarify such an interesting phenomenon for a long time.

The Dongjiang River, which lies in south China, is an important fresh water source of the Pearl River Delta (one of the most developed areas in China). Specifically, it supplies 70% fresh water for Hong Kong. However, water shortage and water pollution in the Dongjiang basin appear to be more and more serious in recent years, due to the fast and persistent economic development and urbanization. Land use in this region has changed prominently since 1980s and contributed to hydrological response in the Dongjiang basin. Meanwhile, climate change under the global warming also plays an important role in the variation of hydrological processes. The combined effects of land use and climate changes lead to a series of conflicts between water use and water supply. In this paper, four sub-basins, the Shuntian, Lantang, Jiuzhou, and Yuecheng, within the Dongjiang basin are selected for study (Figure 1). Major climate variables were taken to trend analysis by using the Mann-Kendall test method, and runoffs in two different periods (natural and human activity periods) were simulated with the SCS monthly model for each basin. The aims of this study are to explore the temporal and spatial characteristics of land use change and climate variables in each sub-basin, to study the possible causes of the change, more specifically, to identify the roles played by land use change and climate change, as well as other human activities, on the runoff change, respectively. A comparison between the independent effect of land use and climate change in the four sub-basins will be conducted.

2. Basic Knowledge and Data

2.1. Study Area

The Dongjiang River (Figure 1) springs from Jiangxi Province and flows into the Pearl River estuary from northeast to southwest, which forms the Dongjiang basin. The Dongjiang basin is located between 113°52′ and 115°52′ E in longitude and 22°38′and 25°14′ N in latitude. This basin has a drainage area of 35,340 km2, about 90% of which is in Guangdong Province. The mainstream of the basin is about 562 km long with an average slope of about 0.039% [30].

The Dongjiang basin has a subtropical climate with a mean annual temperature of 21°C. The annual rainfall over the basin varies between 1500 mm in the dry season from October to March and 2400 mm in the wet or monsoon season from April to September. The basin has some different soil types and the dominated type is alluvial, which is centered largely in the central and southern part of the basin. The major land use type of this area is forest, although the urban area has extended year by year with the fast pace of urbanization since the 1980s. Locations of the four sub-basins are shown in Figure 1.

2.2. Development of Economy and Society

The Dongjiang basin has experienced a fast development of economy and society over the recent 30 years. Population and GDP (gross domestic product) of the three regions (Huizhou, Heyuan, and Dongguan) within the basin from 1980 were analyzed. Detailed information can be seen in Figure 4.

Population and economy in the three regions have kept a rapid and persistent development due to China’s reform and opening policy. Refering to Figure 4, population of each region has been expanding rapidly since 1980. In which, population of Huizhou had increased from 1.92 million in 1980 to 3.6 million in 2006, with the average growth rate of 1.08%. The GDP of each region also increased rapidly, especially Dongguan, whose GDP increased from 8 billion Chinese Yuan in 1980 to 262.7 billion Chinese Yuan in 2006 with the average growth rate of 14.34%. There is no doubt that the growth of population and the development of economy depend largely on land resources. Therefore, we can infer that land use has changed greatly for decades in accordance with the social and economic development in the study area.

2.3. Data

The climatic data used in the study, including time series of annual average precipitation, evaporation, and temperature from 1956 to 2008 of the 21 meteorological stations in the Dongjiang basin, was provided by the Guangdong Meteorological Bureau. The monthly runoff time series for 4 hydrological stations (Shuntian, Lantang, Jiuzhou, and Yuecheng) from 1970 to 2008 were acquired from Hydrological Bureau of Guangdong Province. The daily precipitation, evaporation, and runoff time series of the main rainfall and hydrological stations in the Dongjiang basin from 1970 to 1978 were extracted from Water Conservancy and Electric Power Bureau of Guangdong Province.

Two multitemporal satellite sensor images, Landsat Thematic Mapper (TM) imagery of the 1980s and 2000s, were downloaded from Global Land Cover Facility. Based on the images, two periods (1980s and 2000s) of land use and vegetation cover maps of the Dongjiang basin (Figure 2) were gained with the ArcGIS spatial analysis technique. A digital elevation model (DEM) with a spatial resolution of 30 m and soil type data with a spatial resolution of 90 m (Figure 3) of the Dongjiang basin were downloaded from China Soil Scientific Database (CSSD) and CGIAR-CSI, respectively.

3. Methodologies

3.1. Mann-Kendall Test for Time-Series Trend

Developed by Mann and later improved by Kendall [41, 42], highly recommended for general use by the World Meteorological Organization [43], the Mann-Kendall test was widely used to detect time-series trends in hydrological and climatic data in many researches. The Mann-Kendall test has some advantages including its simplicity, ability to deal with nonnormal and missing data distributions, and robustness to the effects of outliers and gross data errors [4446].

Firstly, the Mann-Kendall test defined a variable as:

Then the presence of a statistically significant trend is evaluated using the value:

Besides, the standard normal cumulative distribution function under the significance level is given by . If , the hypothesis of for a two-sided test can be accepted. Positive value indicates an upward trend while negative value a downward trend. More descriptions of this method can be found in many researches [4749]. In this paper, trend analysis of climatic variables and runoff in the Dongjiang basin were conducted with the Mann-Kendall test based on over 50 years of climatic and hydrological data.

In addition, mutation detections of the climatic data were also conducted by the Mann-Kendall mutation analysis test method; a statistic variable of this method is where is the counts for the series when is greater than .

is the standard normal distribution and can be calculated based on the series: where =0, and are the expected value and variance of which can be calculated on condition that the series is mutually independent with the same continuous distribution:

3.2. SCS Rainfall-Runoff Model

Developed by the Soil Conservation Service of U.S. Department of Agriculture (USDA) early in 1954, the SCS model is widely used in the USA and many other countries [50, 51]. Several advantages of this model are as follows: the characteristics of the land cover such as soil, slope, vegetation and land use in a basin can be considered with the SCS model; the possible change on the rainfall-runoff relationship can be preliminarily estimated according to the change of land use, and it has the advantages of simplicity in structure and convenience in use with very low dependence on data.

The runoff yield for the SCS model is where is the runoff, is the precipitation, is the initial loss, and is the present probable maximum retention in the basin which is the top limit of the later loss.

Initial loss can be calculated by an empirical correlation with  : where is the coefficient of the initial losses which is related to the initial soil water content of the basin.

To eliminate the large variation scope of the value, there is an empirical relation for (a non dimensional parameter) and  :

The curve number CN is a key and comprehensive parameter within SCS model. CN describes the watershed features before a rain and is affected by AMC (antecedent moisture condition), slope, vegetation, soil type, and land use condition with a value of 0~100. AMC can be divided into 3 classes: AMC I-arid condition, AMC II-normal condition, and AMC III-moist condition. The specific basis for the classification and the estimation of the CN value can be found in related references [52].

Input by the observed monthly precipitation and evaporation, the SCS model yields the monthly runoff. A key issue here is to calculate the actual evaporation by using the SCS monthly model. Based on the research results from the two-parameter monthly water balance model [53], the monthly actual evaporation can be calculated by where is the actual monthly evaporation, is the observed evaporation from the evaporating dish, is the monthly precipitation, and is one of the parameters in the model (nondimensional).

The total discharge consists of two parts: the surface flow () and the baseflow (). can be gained by (6), while can be calculated with where is the coefficient of the groundwater flow and is the soil moisture content which can be calculated by

A schematic diagram for the water balance in a basin is shown in Figure 5. The SCS model is used to simulate the monthly runoff process in this study. There are two parameters, and , in the model.

3.3. Approach to Distinguish the Respective Impact of Land Use and Climate Change to Water Resources

The total runoff change can be obtained by the difference between the observed runoffs in two different periods, respectively, that is, the intensive human activities period and the low human activities period (natural condition). It is assumed that the runoff change in the low human activities period is affected only by climate change.

Therefore, the impacts of climate change, land use change, and other human activities on runoff change can be separated according to the following calculations: where is the runoff change due to land use change, from time period of phase I to phase II; and are the simulated runoffs corresponding to the land uses of phase II and phase I, respectively, and both of them can be calculated with the SCS monthly model; is the runoff change caused by climate change; is the natural runoff of phase I; is the total runoff change; is the runoff change due to other human activities except for land use change; , , and are the percentages for the roles played by land use change, climate change, and other human activities, respectively.

4. Results and Discussion

4.1. Time-Series Analysis

The variation trends of the climatic variables over 50 years in the study basin were analyzed with the linear regression and Mann-Kendall trend test method. Table 1 listed the detailed results obtained by the Mann-Kendall trend test method. Figure 6 demonstrated the results for mutation analysis with the Mann-Kendall trend test and the linear regression method, respectively.

It can be seen from Table 1 that each climatic variable showed different degree of change in the Dongjiang basin during the past 50 years (1959–2008), according to its statistics by the Mann-Kendall trend test method. Precipitation showed a nonsignificant increasing trend at 99% confidence level, which is largely due to local atmospheric circulation and topography. Temperature in the same period showed a significant increasing trend . While sunlight showed a significant decreasing trend. Meanwhile, very fast urbanization caused quick increase of impervious area in the study basin in the last 50 years. These all contributed to a decreasing trend of evaporation. Temperature showed an increasing trend with a distinct mutation taking place around 1997, and sunlight and evaporation both had an abrupt change in 1982. Overall, the change of climatic variables in the Dongjiang basin during the past 50 years was significant.

It can be found from the value of each climatic variable that there was a significant internal variation for precipitation but not for temperature. It can be illustrated by both the rainfall pattern and the atmospheric circulation in the study area: Dongjiang basin is located in the southern humid region in China and deeply affected by the monsoon circulation with the rainfall patterns of frontal rain and convectional rain. Therefore, the internal variation of precipitation is significant.

To clarify the relations among each climatic variable and the runoff, the correlation coefficients and relevancies were analyzed. According to the spatial differences of land cover and climatic variables, the Dongjiang basin was divided into three parts: the upper part centered with Longchuan, marked D1; the middle part centered with Heyuan, marked D2; and the whole basin centered with Boluo, marked D3. Tables 2 and 3 listed the correlation coefficients and relevancies, respectively. It can be seen from Table 2 that precipitation was positively related to runoff, and their correlation coefficient was the largest, which revealed that runoff in the Dongjiang basin depended mainly on precipitation. In contrast, temperature, sunlight, and evaporation were negatively related to runoff. Among the three parts of the basin, the correlation coefficient of precipitation and runoff in D2 was the least. The reason was that there located the biggest reservoir, the Xinfengjiang Reservoir (Figure 1) in D2. In addition, more than 20 middle and small sized reservoirs and hydraulic projects were constructed in the same period, which made a more significant change of land cover as compared to D1 and D3. Actually, runoff in D2 is mainly controlled by regulation of the reservoir, so precipitation was less correlation with runoff in D2. Table 3 showed that the relevancy between precipitation and runoff was the largest, indicating that precipitation was the major driver for runoff change in the whole study basin, which is similar to the results of correlation analysis.

4.2. Land Use Change Analysis

It has been found that the urban land and farmland increased in the study area [54]. Based on the raster graphics for land use of two periods (1980 and 2000), temporal and spatial variations of land use change in the Dongjiang basin were analyzed with the ArcGIS spatial analysis technique. Figure 2 showed the maps of land use in 1980 and 2000, respectively. Table 4 depicted the detailed information of the land use change in the form of a comparison.

It can be seen from Figure 2 that the dominant land uses in the study area, in 1980, were mainly forest (widely distributed in the whole basin) and farmland (mainly distributed in the upper basin), accounting for 64.56% and 22.33%, respectively. while garden, meadow, and water took a little part with the percentage of 5.94%, 4.71%, and 2.28%, respectively. Urban land mainly distributed in the downstream basin, and the unutilized land took the least percentage in the Dongjiang basin. The area of two types of land use, urban land and farmland, increased from 1980, along with the decrease of garden (2.65%), meadow (1.23%), water (0.49%), and forest (2.38%) in 2000, while the dominant land use was still forest in 2000. In all, since 1980s, the land use in the Dongjiang basin has been characterized with a prominent increase in urban land, a little increase in farmland, and great decrease in forest area, while little change in water area and unutilized land. This land use change due to fast social and economic development had a significant impact on the hydrological response and then contributed to the variability of water resources in the region.

4.3. Runoff Simulation and Runoff Change Analysis
4.3.1. The Value of CN

CN (curve number) can reflect the capacity of runoff yield for the land cover with a continuous spatial distribution. The CN isocline represents the runoff yield capacities in the study basin.

Based on land use maps of two periods (1980 and 2000) and soil type data of the Dongjiang basin, CN distribution maps (see Figure 7) of 3 AMC (antecedent moisture condition) scenarios in the same periods were obtained with spatial interpolation. Figure 8 showed the CN distribution maps for AMC II of the four sub-basins in the two periods.

4.3.2. Model Calibration and Validation

The SCS monthly model was calibrated in the period of 1970–1976 with climatic data and then validated in the period of 1977-1978 by manual calibration method with the acceptable set of parameters after CN value analysis. To measure the performance of the model, we chose Nash-Sutcliffe coefficient of efficiency (NSE) criterion [55] and relative error (RE) as objective function.

The results of simulation were listed in Table 5. The comparisons between simulated and observed runoffs of the calibration and validation periods in Shuntian sub-basin were illustrated in Figures 9 and 10, respectively. The simulated results by the SCS monthly model were comparable with the observations in the four sub-basins. Nash-Sutcliffe coefficients were above 0.8, and RE was within the pale of 0.02 for those all four sub-basins in the calibration period. It can be inferred from Figures 9 and 10 that the simulated hydrograph matched well with the observed one, and the peak flows were coincided in both the calibration and validation periods. Overall, model performance was acceptable within the study domain, and reliable to be extended to reconstruct the natural runoff.

4.3.3. Runoff Simulation in Human Activity Period

Based on the CN value of land use in the first phase (before 1980) and climatic data in the second phase (precipitation and evaporation, from 1980 to 2000) in the four sub-basins, the natural monthly runoff I and its process in the human activity period were simulated by using the SCS monthly model. And the monthly runoff II was obtained under the land use change condition (with the CN value of land use in 2000 as its input) by the same method. The difference between the simulated runoffs in the two phases was the runoff change due to land use change.

The simulated monthly runoff change due to land use change from 1980 to 2000 for each sub-basin was shown in Figure 11. It can be found that the runoff change varied from sub-basin to sub-basin, and the changes in Yuecheng sub-basin and Jiuzhou sub-basin were less than those of in Shuntian and Lantang, which was relevant to the variation scope of the CN value in the basins.

4.3.4. The Quantitative Effect due to Climate and Land Use Changes

The analysis of quantitative impacts on runoff due to human activities particularly land use and climate changes was carried out in the four sub-basins. The runoff change due to land use change was obtained from the simulated runoff by the SCS model in the two periods (human activity and natural periods). The quantitative impacts of climate change and human activities were identified by the method presented in (12), and the results were shown in Table 6.

It can be seen from Table 6 that, take Shuntian sub-basin, for instance, as compared to natural period, the runoff in human activity period increased by 4 × 106 m3. The 24.11% increase of runoff was contributed by the land use change. On the contrary, climate change contributed the 44.6% decrease to runoff, which was superior to that of land use change. A meaningful clue has been found for explaining such a result: as two most important factors of the runoff change, precipitation showed an insignificant increasing trend, which was different from that of temperature (significant increasing trend), and the compound effects led to the reduction of the total runoff. Other factors including natural and human aspects accounted for 31.29% to the runoff change. Overall, all the driving factors can make the runoff change to a different level among which the role played by climate change took nearly the half for each sub-basin. This was in accordance with the result of identifying the quantitative effect of land use and climate change on runoff in the high flow period in the Dongjiang basin [30].

The impact of land use change on the total runoff change in Lantang sub-basin was the highest (29.94%) among the four sub-basins, while the impact ratios of climate change and other human activities were 40.51% and 29.55%, respectively. In summary, human activities contributed 59.49% to runoff change, which revealed that human activities were the major driver for runoff change.

It was obvious that climate change with increased precipitation and decreased evaporation caused the increase of runoff in the study area. Furthermore, changes of inner-annual distribution precipitation also affected the runoff change. Seven pairs of simulated runoff under nearly equal amount of precipitation were shown in Table 7. It can be found from Table 7 that the greater the percentage of precipitation in flood season, the greater the simulated runoff in Shuntian sub-basin, except for 1985 and 1987. It was also implied that precipitation was the major driver for runoff change in Shuntian sub-basin.

To analyze the impact of land use change, we divided annual precipitation into four classes as under 1500 mm, 1500–1800 mm, 1800–2000 mm, and above 2000 mm in the four sub-basins. The total runoff change due to land use change in each class was shown in Figure 12. It can be found that under larger precipitation, the runoff (to be flattened) of flood season was less as compared with that of drought season and the whole year in the four sub-basins. This was attributed to land cover conversion due to construction of reservoirs, which stored much water in flood season for flood prevention and water use. In addition, runoff yield changed to be much quick and more sensitive to precipitation, along with the land cover conversion from forest to urban areas or other vegetation types. The impact of land use change on runoff was more distinctive when the precipitation was larger.

Finally, it can be concluded that the runoff change was affected by many factors, and the contribution ratio of each factor was different. Climate change and human activities (especially the land use change) were the two dominant drivers which contributed about 40%~50% each to the runoff change. Particularly the land use change, whose impact on the runoff change in flood season under large precipitation, should not be neglected. The role played by climate change and human activities including land use change should be considered, respectively, in the analysis of variability and availability of water resources, and then reasonable measures and policies should be taken.

5. Conclusions

Based on land use maps of two time periods in the Dongjiang basin, this study identified the quantitative effects of land use and climate change on the runoff, which revealed some interesting results with the application of the SCS monthly model.

The climate experienced significant changes in the basin over the past 50 years, and the changes were detected using the Mann-Kendall test method with the results of nonsignificant increase in precipitation, significant increase in temperature while decrease in evaporation and sunlight.

Land use has experienced a significant change since the 1980s in the Dongjiang basin with the characteristic of spatial distribution for the CN value under three AMC scenarios in two periods. The average value of CN increased and the CN value varied in each sub-basin due to different land use patterns. In particular, expansion of urban area in the south part of the study area (downstream) and deforestation in the whole area provided contributory factors that affected hydrological processes and subsequently increased runoff.

The SCS monthly model performed well in the four sub-basins of the Dongjiang basin, and the results showed that the runoff change in each sub-basin during two time periods was different. The runoff change in Yuecheng sub-basin was less than that of Shuntian sub-basin and Lantang sub-basin, which was in good correlation to the variation scope of the CN value. The more land use changed, the greater the CN value changed, which resulted in the greater variation of the runoff under the same climatic condition.

The separated quantitative effect of land use and climate change on the runoff showed that climate change and human activities (including land use) contributed half and half, respectively, to runoff change. While land use change independently contributed 20%~30% to the total runoff change. Moreover, the effect of land use on runoff change was in different level under different amounts of precipitation. The effect of climate change on runoff change had a different inner-annual distribution even under a same amount of annual precipitation.

Further research is required to acquire the regional future climate scenarios coupled with the hydrological model of a basin scale under GCMs (general circulation models) with the downscaling technique, so as to further quantify the relations between runoff and climatic variables. In addition, the space-time distribution of floods and droughts resulted from the runoff change should also be studied to provide scientific framework for basin-scale water resources management.

Acknowledgments

The research is financially supported by the National Natural Science Foundation of China (Grant no. 51210013), the Public Welfare Project of Ministry of Water Resources (Grant no. 201301071-02, 201301002-02), the project from Guangdong Science and Technology Department (Grant no. 2010B050300010), the project from Guangdong Water Resources Department (Grant no. 2009-39), Key Project of Chinese National Programs for Fundamental Research and Development (973 program) (Grant no. 2013CB036406, 2010CB951102), and the Pearl-River-New-Star of Science and Technology supported by Guangzhou City (Grant no. 2011J2200051).