Identifying the Relationships between Water Quality and Land Cover Changes in the Tseng-Wen Reservoir Watershed of Taiwan

The effects on water quality of land use and land cover changes, which are associated with human activities and natural factors, are poorly identified. Fine resolution satellite imagery provides opportunities for land cover monitoring and assessment. The multiple satellite images after typhoon events collected from 2001 to 2010 covering land areas and land cover conditions are evaluated by the Normalized Difference Vegetation Index (NDVI). The relationship between land cover and observed water quality, such as suspended solids (SS) and nitrate-nitrogens (NO3-N), are explored in the study area. Results show that the long-term variations in water quality are explained by NDVI data in the reservoir buffer zones. Suspended solid and nitrate concentrations are related to average NDVI values on multiple spatial scales. Annual NO3-N concentrations are positively correlated with an average NDVI with a 1 km reservoir buffer area, and the SS after typhoon events associated with landslides are negatively correlated with the average NDVI in the entire watershed. This study provides an approach for assessing the influences of land cover on variations in water quality.


Introduction
Land use and land cover changes, associated with human activities and natural factors, compromise many ecosystem services in a watershed [1,2]. For example, forestland converted to agricultural or urban land may have increased erosion, runoff, and flooding [3]. Changes in land use and land cover interact with anthropogenic and natural drivers to affect the water quality of watersheds. Studies have used environmental and landscape data to examine the relationships between land use and land cover changes and suspended sediments [4][5][6][7] and nutrients [1,[7][8][9][10]. Ahearn et al. showed that land use and land cover exert the greatest control over water quality in the Cosumnes Watershed, California [7]. The percentage of agricultural coverage had a significant influence on nutrient loading. Sliva and Williams used multivariate analysis to determine whether there was a correlation between water quality and landscape characteristics within the local Southern Ontario watersheds in Canada. They compared the influences of buffer zones and whole catchment landscape characteristics on water quality [9]. Li et al. showed the impact of land use and land cover on the water quality in the Upper Han River basin, China [10]. The correlation and regression analysis indicated that water quality was significantly related to vegetated coverage.
Water quality is controlled by numerous anthropogenic and natural factors [7]. The quality of receiving waters is affected by human activities in a watershed by point sources, such as wastewater treatment facilities, and non-point sources, such as runoff from urban areas and farmland [11]. Understanding non-point source pollution requires an understanding of how particular land covers influence water quality within a watershed. The extent that land covers hierarchically affect water quality at space-time scales is a key question. The most widely used land cover index in this context is the normalized difference vegetation index (NDVI), which is a function of red and near-infrared spectral bands [12]. On a regional scale, multi-temporal NDVI images are practical for monitoring vegetation dynamics. The multi-temporal NDVI is useful for classifying land cover and detecting the dynamics of vegetation [13,14]. However, major changes in the NDVI are noted near landslides that were induced by disturbances in Taiwan [15,16]. For example, a typhoon is one of major natural disturbances to land cover. Sequent typhoons and rainstorms cause abnormal destruction to the vegetation; this destruction is influenced by rainfall distributions and typhoon paths [15].
The NDVI data were derived from SPOT satellite images in the Tseng-Wen Reservoir Watershed, Taiwan, before and after Typhoon Morakot and several other large typhoons (e.g., Typhoon Mindulle in 2004, Haitang in 2005, Sepat in 2007, Kalmaegi in 2008, and Fanapi in 2010) [17] to identify the changes to land cover. To represent land use and land cover change, an evaluation of multiple NDVI spatial scales was conducted. The study identified and delineated the relationships between temporal variations of the NDVI and water quality in the study area.

Study Area
The Tseng-Wen Reservoir is a multipurpose reservoir designed for flood control, hydroelectric power generation, irrigation, water supply, recreation, and flow augmentation. The storage capacity of the Tseng-Wen reservoir is 608 × 10 6 m 3 , but its effective water storage is 490 × 10 6 m 3 , and the hydroelectric plant capacity is 50 MW. The Tseng-Wen Reservoir Basin is located in the upstream area of the Tseng-Wen River system in Chiayi County ( Figure 1). The entire watershed area of this river basin is 1,176 km 2 , in which the Tseng-Wen Reservoir watershed covers 481 km 2 . The average slope of this river basin is approximately 1/57. Average rainfall in this watershed area is approximately 3,000 mm per year and annual average temperature is about 23.4° Celsius. Rich soil in the watershed is suitable for fruit and tea farms. Agriculture has restricted near the reservoir but tourism has developed in recent years.
Numerous  [17]. Especially, Typhoon Morakot struck Taiwan from 7 to 9 August, 2009, and produced record-breaking rainfall and catastrophic damage in Southern Taiwan. The typhoon produced copious amounts of rainfall peaking at 2,777 mm.

Satellite Images
Multi-temporal Système Pour l'Observation de la Terre (SPOT) satellite images obtained after typhoons in 2001, 2003-2005, and 2007-2010 were used to quantify land cover changes in our study. Details for the dates are listed in Table 1. For atmospheric correction, Fast Line of Sight Atmospheric Analysis of Spectral Hypercube (FLAASH) is applied to correct the visible and near-infrared wavelengths in the satellite images [18]. Then, the NDVI maps were derived from the SPOT images taken in 2001 with a 20 m resolution, and with a 10 m resolution in 2003-2010. Moreover, the SPOT images were classified using supervised classification by the software package ERDAS IMAGINE. Land-use types were classified into the following six categories: forested land, built-up land, landslide, grassland, water, and bare land [19]. The reference maps were the aerial photographs by the Aerial Survey Office, Forestry Bureau, Taiwan.

Water Quality Data
Seasonal time series of water-quality data monitored in the reservoir were obtained from Taiwanese EPA Web sites [20]. The water quality data observed at three stations were obtained from 2001 to 2010 and data sampling frequency was three months. The sampling sites are shown in Figure 1. The water-quality variables, such as nitrate-nitrogen (NO 3 -N), suspended sediments (SS), chemical oxygen demand (COD), dissolved oxygen (DO), total phosphorus (TP), and turbidity, were derived for analysis. The variables are used as general indicators of water quality. For example, the COD is commonly used to measure the amount of organic compounds in water. As the DO in water drops below a threshold, aquatic life is under stress. The presence of high nitrates and TP concentrations in water indicates possible pollution of the water. Turbidity is the haziness of a fluid caused by the SS, which are solid particles usually transported by flowing water.

Regression Model
Reservoir water chemistry was sampled at the outlets and downstream of the reservoir every three months. We refer to the water samples after typhoons in the fourth season and acquire annual values each year. We ran a series of models to examine the correlations between land-cover and the water quality variables.
Our basic statistical tool was stepwise multiple linear regression, with backward selection of variables and p = 0.1 to enter or remove variables. Cases with missing data were excluded. Statistical analyses were done using SPSS 10.0. The dependent (response) variables are NO 3 -N and SS concentrations that are selected from the high-NDVI correlated water quality factors. The variable details are shown in Section 3.2. In addition, the independent variables are average NDVI at various spatial scales such as average NDVI in whole watershed, average NDVI in reservoir 1 km, 2 km and 3 km buffer zones, can be represented as: (1) where represent NDVI value varies with space s and time t; is average NDVI in whole watershed or reservoir buffer zone during time t; S is the domain set that defined as whole watershed or reservoir buffer zone; num is the number of the pixels in the set S. .81%, grassland has increased 2.35%, landslide has increased 2.12%, bare land has increased 0.55%, and built-up land has decreased 0.21% (Figure 3). The results matched previous studies [21,22] that many landslides in the Tseng-Wen reservoir watershed were caused by typhoons. Table 1 (Figure 4) are sufficient to present land cover changes induced by disturbances, particularly by spatial structure, variability, and spatial correlation. The disturbances impacted the fragmentation and interspersion of the low NDVI patches and created heterogeneous patterns across the landscape within the affected area [16]. However, land cover change may be different in the spatial scales. The box plot shows that the range between the lower and upper quartiles in the NDVI decreases when the buffer zone increases ( Figure 5). Table 2 shows the correlation coefficients between the average NDVI in the watershed and average water quality factors during the whole year and after typhoon. The average value of the NDVI in the watershed is strongly correlated with NO 3 -N during the entire year and is strongly correlated with SS after typhoon from 2001 to 2010. These water quality factors such as SS and NO 3 -N concentrations are the indices for water quality assessment when considering land cover changes. The average SS concentrations after typhoons and the average annual NO 3 -N in the sites are used by later descriptive statistics and regression analysis.       concentrations are cyclical. Typhoons and heavy rainfalls trigger large sediment discharge into the rivers of Taiwan and cause high-suspended sediment concentrations during these events. Most nitrate concentrations in the water drains from agricultural land. However, Figure 6(b) shows NO 3 -N concentrations vary with a decreasing trend. Table 4 shows the regression models for water quality and land cover changes at various spatial scales. The factors comprise the average NDVI in the watershed, and the average NDVI in 1 km, 2 km, and 3 km buffer zones. Results show that annual nitrate-nitrogen concentrations are positively correlated with the NDVI with a 1 km buffer area. However, SS is negatively impacted by the average NDVI in the watershed, suggesting that typhoons impact land cover change in the watershed. For example, typhoons cause landslides, and these are a major source of soil erosion and sediment yield in the watershed. The average NDVI in the watershed adversely impacts the water quality, and therefore, increases sediments associated with water quality. The average NDVI in the watershed becomes a key factor influencing the SS concentration. Typhoon events are major natural disturbances causing NDVI changes and also cause serious landslides [16,23,24]. Both important factors affecting soil erosion and sediment delivery to river channels are changes in land use and climate. Due to destruction of vegetation and increased soil exposure in the watershed after rainstorms and typhoons, the NDVI values decreased. The rainstorms and typhoons cause divergent destruction of vegetation, and led to an increase in the potential for soil erosion.

Relationship Between NO 3 -N Concentration and Land Cover Change
Results show that the dominant explanatory variables in NO 3 -N cases have an average NDVI with a 1 km buffer zone. In some reservoirs and lakes, the primary indicator of agriculture is dependent on NO 3 -N concentration. The nitrate concentration is correlated with agricultural practices during the high-flow period. Human activity alters the patterns of nitrate concentrations during storm events in the agricultural catchment [25]. Since Taiwan joined the World Trade Organization (WTO) in 2002, imported agricultural products are cheaper than domestic ones, thus negating the need for extensive agriculture areas. This corresponds to the data that tea farms in Chiayi County decreased from 2,292 to 2,189 ha from 2005 to 2011 [26]. Forestry, agriculture, and anthropogenic activities impact the quality of water over short and long periods [27]. The SS and NO 3 -N are typically sensitive in landslide and agriculture land areas. The previous results match that percent agricultural coverage had a significant influence on both SS and nitrate-N loading [7]. NDVI variation results imply that as land cover changes; hence, the multi-scale NDVI, which is one of the indices in the watershed of land cover changes, is associated with water quality and is hard to directly link with agriculture. However, the land use classification of SPOT images is also hard to identify the agriculture land. Further study could consider the high-resolution satellite images in the land use classification.

Conclusions
This study examined the NDVI images from 2001 to 2010 based on SPOT imagery data. The imagery shows that the land cover changes in the study vary with the influences of typhoons and human activities. Satellite image data showed a general decline in the acreage of vegetation cover implying increased landslide and decreased forest pressure on the vegetation resources.
Land cover change had a significant influence on both suspended solid and nitrate-nitrogen loadings. Simple regressions were performed that showed water quality is related to land cover in various spatial scales. Annual NO 3 -N concentration is positively correlated with an average NDVI with a reservoir 1 km buffer area, but SS are negatively correlated with an average NDVI in the watershed after typhoon events. Understanding the relationship between land cover change and water quality is useful for watershed management and pollution prevention plans. Further study should add additional spatial and independent variables to the models, such as percentage of land use type, anthropogenic activities and typhoon precipitation.