Linkage Analysis of Land Use / Cover Patterns and Hydro-Chemical Characteristics in Different Seasons in Ebinur Lake Watershed , China

Ebinur Lake Watershed, with its oases and deserts, is a typically arid and mountainous region on the northern slope of Tianshan Mountains. Along with ever increasing human activities, agricultural and domestic wastes have been directly discharged into river systems around the Watershed, which consequently poses a grave threat to the sustainable development of Xinjiang. Through statistical and spatial analysis, we have determined the relationships between land use/cover (LULC) and hydro-chemical characteristics during rainy and dry seasons in 2014. Spatial patterns of hydro-chemical characteristics as demonstrated by mineralization, pH, electrical conductivity (EC), and SO4, HCO3, Cl−, Mg2+, Na+, Ca2+ and K+ concentrations were examined in 113 sites. Since hydro-chemical characteristics were affected by LULC patterns, this study delineated six zones to analyze the distribution characteristics of hydro-chemical parameters and its relationship with LULC patterns. The results showed that, except for the pH, all variables demonstrated significant spatial differences during dry and rainy seasons. In addition, the relationships between water quality and LULC patterns indicated that the farmland, forest–grassland, water body and salinized land all affected hydro-chemical characteristics during both rainy and dry seasons. Especially, decreased rainfall, irrigation, surface runoff, the area of lake, etc. largely led to the increase in ion content, which had great influence on hydro-chemical characteristics parameters in dry season. Furthermore, we established several stepwise linear multiple regressions models. The results showed that pH, mineralization and Ca2+ were defined by forest-grassland, while the Cl− and Mg2+ were defined by salinized land during the rainy season. The pH and Na+ were estimated based on farmland, whereas Cl− and Ca2+ were estimated based on forest-grassland during the dry season. In conclusion, this research on the relationships between the spatial distribution of hydro-chemical characteristics in Ebinur Lake Watershed and LULC patterns will be significant for the sustainable management of the arid regions in northwest China.


Introduction
Water quality is of great importance to the study of water resources in arid regions.Hydro-chemical characteristics may help understand the environmental and geological conditions, i.e., the controlling factors in which waters are formed [1].Accurate information on the spatial distribution of surface water characteristics is imperative for assessment of environmental monitoring, land surface water management as well as watershed changes [2,3].In particular, there is little or no information on river hydro-chemical characteristics in arid zones in China.Formal references on its hydro-chemical characteristics are also scarce in China [4,5].Generally, stream water chemistry is controlled by numerous natural, e.g., weathering, precipitation inputs, and anthropogenic factors.Hydro-chemical characteristics have deteriorated noticeably in many countries in the past decades due to poor land use practices [6][7][8], which resulted in the material sources and evolution processes of chemical composition and the content of different ions are being very complicate [1].Land use/cover (LULC) patterns have important effects on the quality of river aquatic systems within a watershed [9][10][11]; for example, the development of a watershed from the dominant LULC to an artificial land system frequently reduces the base flow, which is a result of changing the groundwater flow pathways to surface water bodies [12].Numerous problems related to hydro-chemical characteristic are caused by human-driven changes in LULC [13][14][15][16], because these anthropogenic activities are directly reflected in LULC patterns in the river basins [17].Recently, many studies around the world have shown that LULC has a strong impact on the quantity and quality of the stream water in watersheds [18][19][20], and significant correlations exist between LULC patterns and hydro-chemical characteristic parameters [21][22][23].For example, strong relationships have been found between increasing extent of agricultural development in catchments and declining water quality of streams [24].Donohue et al. [25] associated land cover characteristics with the ecological status and identified that urban, pasture lands, and arable were dominant factors impacting the ecological quality of river streams.Lee et al. [9] indicated that water quality of reservoirs in South Korea was closely associated with urban, agricultural, and forest areas.
Meanwhile, numerous researchers in China have examined the relationship between LULC and hydro-chemical characteristic, but few of these studies have included the typically arid areas in Ebinur Lake Watershed.Zhang et al. [26] reported that LULC significantly impacted the hydro-chemical ecosystem in this area.Given seasonal differences, the unbalanced distribution of precipitation amount results in an apparent variation in surface runoff and further imbalance in the spatial distribution of hydro-chemical characteristics in the research zone [27][28][29].During the wet season (May) in Ebinur Lake Watershed, melted water from mountain ice and snow is collected, which promotes flow in Bortala River and Jing River, thereby resulting in a significant increase in surface runoff.During the dry season, however, elements in water can be partially absorbed and purified to a certain extent because of the rise in temperature and the growth of aquatic plants in rivers and lakes.Therefore, a significant change in surface runoff and seasonal change in the research zone are important factors that result in noticeable differences in the spatial distribution characteristics of hydro-chemical characteristics during the wet and dry seasons.In this study, we performed a comparative examination of LULC in the target area using GIS procedures and RS satellite images in May and September of 2014 to assessment of the impacts of LULC on the hydro-chemical characteristics in Ebinur Lake main tributaries.The main objectives of this study are to: (1) Explore the spatial patterns of hydro-chemical characteristics during rainy and dry seasons; (2) Analyze LULC changes during rainy and dry seasons; and (3) Study LULC patterns' effects on hydro-chemical characteristics during dry and rainy seasons.

Study Area
Ebinur Lake Watershed is located at 43  1).The total area of the Watershed is 50,621 km 2 , which consists of mountains, plains and water bodies [30].Ebinur Lake Watershed was once fed by 12 branch rivers from the Bortala, Jing and Kuytun River systems.The Watershed is now mainly recharged by alpine glacier melt water and mountain precipitation.Total runoff is 37.46 × 10 8 m 3 /year [31].Due to environmental changes and human activities (i.e., agricultural development in oases), some rivers have gradually lost their hydraulic connections, with only Bortala River and Jing River now feeding Ebinur Lake.Ebinur Lake is a closed saline lake with the water salinity from 80 g/L to 120 g/L; sodium sulfates and chlorides predominate in the composition of dissolved salts [32].At the beginning of the twentieth century, the area of the lake was about 1300 km 2 .At present time (late September 2014), the area of the lake is only about 420 km 2 as seen in the satellite image.The surface from the water evaporation represents a vast bare solonchak with a high amount of salts (more than 40% in the upper 2 cm and 10-30% in the lower layer of 2-5 cm) made up of predominantly sodium sulfates and chlorides [33].
Water 2017, 9, 888 3 of 20 composition of dissolved salts [32].At the beginning of the twentieth century, the area of the lake was about 1300 km 2 .At present time (late September 2014), the area of the lake is only about 420 km 2 as seen in the satellite image.The surface from the water evaporation represents a vast bare solonchak with a high amount of salts (more than 40% in the upper 2 cm and 10-30% in the lower layer of 2-5 cm) made up of predominantly sodium sulfates and chlorides [33].

Satellite Image Acquisition and Processing
Twelve Landsat OLI images were acquired in May and September of 2014 over Ebinur Lake Watershed (Table 1).Data were downloaded from USGS website in level 1G product [34] with a spatial resolution of 30 m. Furthermore, all images were selected with less than 20% cloud cover and geo-referenced images were obtained to generate the entire study area images based on a mosaic tool [35].In this study, the images were acquired in different seasons.Thus, it was crucial to make an atmospheric correction and radiometric normalization to all images to eliminate effects of haze and dust as well as solar angle variations.Geometric rectification using 1:50,000 topographic maps was 4erformed on all remotely sensed images.A second order polynomial fit was applied to correct remotely sensed images using 50 ground control points (GCPs) selected from topographic maps, with a root mean square error (RMSE) less than 0.5 pixels.Radiometric calibration was carried out using the Environment for Visualizing Images software (ENVI Version 5.0 (Exelis Visual Information Solutions, Boulder, CO, USA)) [36].The atmospheric correction for the OLI images was undertaken by using the Fast Line-of-Sight Atmospheric Analysis of Spectral Hypercube (FLAASH) method of ENVI 5.0 [36].

Satellite Image Acquisition and Processing
Twelve Landsat OLI images were acquired in May and September of 2014 over Ebinur Lake Watershed (Table 1).Data were downloaded from USGS website in level 1G product [34] with a spatial resolution of 30 m. Furthermore, all images were selected with less than 20% cloud cover and geo-referenced images were obtained to generate the entire study area images based on a mosaic tool [35].In this study, the images were acquired in different seasons.Thus, it was crucial to make an atmospheric correction and radiometric normalization to all images to eliminate effects of haze and dust as well as solar angle variations.Geometric rectification using 1:50,000 topographic maps was performed on all remotely sensed images.A second order polynomial fit was applied to correct remotely sensed images using 50 ground control points (GCPs) selected from topographic maps, with a root mean square error (RMSE) less than 0.5 pixels.Radiometric calibration was carried out using the Environment for Visualizing Images software (ENVI Version 5.0 (Exelis Visual Information Solutions, Boulder, CO, USA)) [36].The atmospheric correction for the OLI images was undertaken by using the Fast Line-of-Sight Atmospheric Analysis of Spectral Hypercube (FLAASH) method of ENVI 5.0 [36].

Water Sampling Collection
The surface water quantity in March, April and May are greater than that in June, July, August, September and October in Ebinur Lake Watershed.Thus, we regard May as in rainy season and September in dry season [37].Water samples were collected in two field surveys in the Watershed during the rainy (May) and dry (September) seasons in 2014.In the last 30 years, cotton has been the major crop in the Watershed with an annual harvest of 7.5 × 10 7 kg [38].Urban residents around the Watershed live primarily on agriculture and husbandry.No heavy industry is present.Therefore, point-source pollution from industrial wastewater was not considered in this study.Based on on the field survey, sampling points were set in the main entrances and exits of lake, the confluence of rivers, the typical area such as farmland around, etc.Our samples were collected from the farmland in Jinghe County and Tuotuo Village around Ebinur Lake, a national ecological zone in Ebinur Lake Bird Isle, and the Ganjia Lake Haloxylon Forest National Nature Reserve.We obtained altogether 56 and 57 samples for rainy and dry seasons, respectively.The sampling points were separated at approximately 5-km intervals along the horizontal gradient from west to east and samples were taken both from mainstreams and tributaries.Ten representative parameters were chosen from the measurements, namely the pH, EC, SO 4 2− , HCO 3 − , Cl − , Mg 2+ , Na + , Ca 2+ , K + and mineralization, with the last being variable with the content of ions, molecules and compounds.In the field sampling, water was sampled at the depth of 0.5 m below the surface water using previously acid-washed high density polyethylene (HDPE) bottles and the environment of sampling points was recorded.
In situ parameters such as pH, EC, mineralization, etc. were measured with portable electronic meters immediately in the field after sampling.All samples were packed and stored in refrigerated containers, tried to slow the progress of the physical function and chemical action before being sent to the labs for further measurement or data collection.In the laboratory, measurement parameters were SO 4 2− , HCO 3 − , Cl − , Mg 2+ , Na + , Ca 2+ and K + .
The value of pH was measured using pH-40A portable pH acidity meter (Hangzhou Lohand Biological Co., Ltd., Hangzhou, China); EC was determined by the conductivity meter (Sensortechnik Meinsberg, Bad Pyrmont, Germany); Mineralization was determined using gravimetric method; SO 4 2− was determined by EDTA-indirect titrimetry; HCO 3 − was determined with the double indicator neutral method; Cl − was determined by AgNO 3 titration; Mg 2+ and Ca 2+ were determined with the EDTA complexometric method; and Na + and K + were determined by the flame photometry.
The determination of water parameters in the laboratory complied with national standard methods defined in [39,40] .

Classification System for LULC Patterns in Ebinur Lake Watershed
Land use maps interpreted from 2014 Landsat OLI images (30 m resolution) were used to measure the LULC composition in the Watershed using ArcGIS version 10.0.Many land cover classification methods have been proposed in recent years [41,42].Compared with these novel classifiers, the decision tree (DT) classifier has generally been used because of its nonparametric and characteristic tree structure [43,44].DT is made up of a series of binary decisions that are used to determine the correct category for each pixel.The DT method also is sufficiently robust to deal with the loss caused by incomplete remote data.Thus overcome the limitation of the data distribution requirements with the maximum likelihood method [45].Thus, in this paper, DT was applied to Landsat OLI images which had been acquired in May and September 2014.
The LULC was classified into five categories based on the national land resource classification system, on the field survey that had been conducted many times to set up the interpretation signs [46], and on the actual LULC patterns of Ebinur Lake Watershed: (1) Farmland, mostly planted with maize, cotton, grape and medlar; (2) Forest-grassland, including forest land, grassland and pasture for livestock; (3) Water body, including reservoir, river, and pond; (4) Salinized land, including slight saline land, moderate saline land and heavy saline land (statistics indicate that soil salinization in Ebinur Lake area mainly occurs in Bortala River, Jing River, the villages and towns surrounding Ebinur Lake, downstream of Daheyanzi River, and north of Bole); (5) Others, including gravel, bare ground, and bare rock.Google Earth was used to verify LULC patterns.

Accuracy Assessment
In this study, the Kappa coefficient was used for a quantitative assessment, which is one of the most widely used methods in the assessment of classification accuracy [47].Using stratified sampling methods, fifty samples were selected in Ebinur Lake Watershed as the "ground-truthing" pixels by GPS positioning.The LULC classification accuracy was evaluated by the overall accuracy and the Kappa coefficient based on the error matrix [48].The Kappa coefficient was calculated as: where P is the total number of pixels from the reference data; P ii is the total number of correct pixels from the category; P i+ is the total number of pixels for the category derived from the classified data; P +i is the total number of pixels for the category derived from the reference data; and n is the total number of categories.In this study, the total accuracy and Kappa coefficient for classification were equal to 91.2% and 0.8848 in the rainy season, and 94.5% and 0.9335 in the dry season, respectively.

Spatial Analysis
Firstly, we use stratified sampling methods to select the points in each LULC patterns by ArcGIS and collected the ground reflection value of every point to analyze the relationship between the ground reflection value of LULC patterns and the hydro-chemical characteristic parameters.Due to the differences of the spatial distribution of LULC patterns, 30 points are selected in each LULC patterns surrounding the actual water sampling points.Then, the ground reflection value of each LULC pattern was extracted, and the correlation between the ground reflection values of different LULC patterns and hydro-chemical characteristic parameters was analyzed.In addition, this study was based on the DEM and stream networks of the Watershed, and delineated six zones to analyze the distribution characteristics of physicochemical parameters of hydro-chemical characteristics and its relationship with LULC patterns in different zones during dry and rainy seasons.Ebinur Lake Watershed was divided into six zones: (1) Zone 1, including Wenquan County and Bole city, mainly from the upstream to the middle part of Bortala River; (2) Zone 2, including Jinghe County and Jinghe Oasis, mainly farmland; (3) Zone 3, mainly including Ebinur Lake, a national ecological zone in Ebinur Lake Bird Isle and the Ganjia Lake Haloxylon Forest National Nature Reserve; (4) Zone 4, including reservoirs and farmland in Wusu County; (5) Zone 5, mainly including farmland in Wusu County; and (6) Zone 6, including Kuitun River and nearby mountains.Meanwhile, because of the uneven distribution of water sampling points in six zones, using the zone to set up the relationship between various LULC patterns and hydro-chemical characteristics parameters harbors some uncertainties.Thus, this study set up 100 m, 200 m, 300 m, 400 m, 500 m, and 600 m buffers around water sampling points using ArcGIS spatial analysis tools.Preliminary analysis found that the relationship between various LULC patterns and hydro-chemical characteristics parameters were better under the 500 m buffer.Thus, land use/cover patterns and their influences on hydro-chemical characteristic were only analyzed under the 500 m buffer during dry and rainy seasons.

Statistical Analysis
For hydro-chemical characteristic properties, single factor analysis of variance (ANOVA test) with the post-hoc test (least-significance difference, LSD).Li et al. [49] used to analyze the variance in hydro-chemical characteristic parameters in different sites at a significance level of p < 0.05.Factor analysis (FA) was used to identify the change of hydro-chemical characteristic factors that affected hydro-chemical characteristics during rainy and dry seasons.Before FA was performed, the hydro-chemical characteristic data were initially standardized by z-scale transformation to avoid misclassification because of the wide differences in data units and dimensionality [39].Kaiser-Meyer-Olkin (KMO) and Bartlett's sphericity tests were conducted to measure the adequacy of the sampling data for FA.For the FA, a given variable was considered to be moderately important when PC loadings are about 0.50 and 0.75 and to be an important contributor to a PC if its absolute loading exceeds about 0.75 [50].Correlations between LULC patterns and hydro-chemical characteristics parameters were tested using Pearson's correlation coefficients with statistical significance at p < 0.01 and p < 0.05 levels (2-tailed), respectively.We conducted Stepwise linear multiple regressions, an effective approach to identify significant land use patterns to explain hydro-chemical characteristics variation in a watershed.For the correlation and regression analysis, the one-sample Kolmogorov-Smirnov test was used in testing the normal distribution for all variables [51].The results showed asymptotic significance (Asymp.Sig.) (2-tailed) values (p values) ranging from 0.148 (Mg 2+ ) to 0.939 (pH) and greater than 0.05, suggesting normal distribution.All statistics were calculated by using SPSS 16.0.The conceptual flow chart describing the methodology is shown in Figure 2.
hydro-chemical characteristics variation in a watershed.For the correlation and regression analysis, the one-sample Kolmogorov-Smirnov test was used in testing the normal distribution for all variables [51].The results showed asymptotic significance (Asymp.Sig.) (2-tailed) values (p values) ranging from 0.148 (Mg 2+ ) to 0.939 (pH) and greater than 0.05, suggesting normal distribution.All statistics were calculated by using SPSS 16.0.The conceptual flow chart describing the methodology is shown in Figure 2.

Characteristics of Hydro-Chemical Characteristics
Given seasonal differences, the unbalanced distribution of precipitation results in an apparent variation in surface runoff and further imbalance in the spatial distribution of hydro-chemical characteristics in the research zone [27][28][29].Overall, hydro-chemical characteristic parameters changed significantly during rainy and dry seasons (Table 2).In the rainy season, the coefficient of variation (CV) of EC, mineralization, HCO3 − , SO4 2− , Ca 2+ , Mg 2+ , Na + and K + exceeded 50%.In the dry season, the CV of Cl − , SO4 2− , Ca 2+ and Na + exceeded 50%.The results showed that the ion

Characteristics of Hydro-Chemical Characteristics
Given seasonal differences, the unbalanced distribution of precipitation results in an apparent variation in surface runoff and further imbalance in the spatial distribution of hydro-chemical characteristics in the research zone [27][28][29].Overall, hydro-chemical characteristic parameters changed significantly during rainy and dry seasons (Table 2).In the rainy season, the coefficient of variation (CV) of EC, mineralization, HCO 3 − , SO 4 2− , Ca 2+ , Mg 2+ , Na + and K + exceeded 50%.In the dry season, the CV of Cl − , SO 4 2− , Ca 2+ and Na + exceeded 50%.The results showed that the ion concentration distributes unevenly in the Watershed.During dry and rainy seasons, the results of ANOVA (Table 3) showed that the significant (Sig.)value of pH was greater than the significance level of 0.05, so the distribution of pH was uniform in each zone.pH range was mainly from 7.8 to 8.1.Overall, most of the variables showed significant spatial differences (p < 0.05) among different zones in rainy season and less significant differences in dry season.Samples and their physicochemical parameters during dry and rainy seasons in each zone of the Watershed are presented in Tables 4 and 5 and Figure 3.The values of EC, HCO 3 − , Cl − , Na + , Ca 2+ and K + were higher in Zone 3, which mainly included Ebinur Lake, a national ecological zone in Ebinur Lake Bird Isle and the Ganjia Lake Haloxylon Forest National Nature Reserve.Zone 3 was a seriously salinized area in the Watershed.The surrounding area of Ebinur Lake was particularly affected by human activities.As a result, the hydro-chemical characteristics in this region were more seriously affected.In addition, the highest value of mineralization, SO 4 2− and Mg 2+ were located in Zone 2, which mainly included the farmland of Jinghe oasis.The low concentration of EC, SO 4 2− , Cl − , Mg 2+ , Na + , Ca 2+ and K + were recorded in Zones 1 and 5.The value of mineralization in Zone 5 was higher than that in Zone 4. The low concentration of HCO 3 − was found in Zone 5, which mainly included the farmland in Wusu County.
Water 2017, 9, 888 activities.As a result, the hydro-chemical characteristics in this region were more seriously affected.
In addition, the highest value of mineralization, SO4 2− and Mg 2+ were located in Zone 2, which mainly included the farmland of Jinghe oasis.The low concentration of EC, SO4 2− , Cl − , Mg 2+ , Na + , Ca 2+ and K + were recorded in Zones 1 and 5.The value of mineralization in Zone 5 was higher than that in Zone 4. The low concentration of HCO3 − was found in Zone 5, which mainly included the farmland in Wusu County.Through the factor analysis (FA), we have identified the changes of hydro-chemical characteristic factors that affect hydro-chemical characteristics during dry and rainy seasons.During the rainy season, for instance, the KMO and Bartlett's sphericity test results were 0.803 and 434.174 (df = 45, Sig.< 0.001), respectively.During the dry season, they became 0.818 and 973.172 (df = 45, Sig.< 0.001), respectively.The results show that the water sample data during rainy and dry seasons are suitable for factor analysis.
During the rainy season, the first two rotated factors with eigenvalue of 1 or greater were extracted using Varimax with Kaiser Normalization, explaining 78.9% of the total variance in the hydro-chemical characteristics data set (Table 6).For Factor 1 in the FA study, the important variables were EC, mineralization, SO 4 2− , Ca 2+ , Mg 2+ , K + and Na + .The moderately important variable was Cl − .Similarly, for Factor 2 in the FA study, the results showed that the Watershed was to be of alkaline hydro-chemical characteristics with high loadings on pH, Cl − and Mg 2+ .During the dry season, the same two rotated factors with eigenvalue of 1 or greater were extracted using Varimax with Kaiser Normalization, explaining 86.5% of the total variance in the hydro-chemical characteristics data set.For Factor 1 in the FA study, the results showed that, except the pH, the other water parameters were important variables.In Factor 2, the results showed that the Watershed was with high loadings on pH.Overall, the majority of water parameters were moderately important variables in FA during rainy and dry seasons.In recent years, with increasing cultivation in the upstream of Ebinur Lake Watershed, irrigation water usage increased, causing the lake water levels to fall and the lake area to shrink.At the same time, the salt flowing with the river water into Ebinur Lake exacerbated the latter's salinization, resulting in the increase of ion content.Especially, Factor 1 in the FA study during rainy and dry seasons showed that the Ca 2+ and Mg 2+ were higher than other cationic contents.Ca 2+ is easy to form into CaCO 3 and then precipitate.With increasing salinization of water, the content of Mg 2+ greatly exceeds the Ca 2+ .With the increase of salt content, the ratio of Ca 2+ and Mg 2+ becomes smaller and smaller, eventually causing the imbalance of Ca 2+ and Mg 2+ , which is harmful to fishery in the Watershed [52].

LULC patterns in Different Zones
According to the classification results (Figure 4), we can clearly see that the farmland area increased from May to September as the crops grow in rainy season and become ripe for harvest in dry season.Meanwhile, with a shrinking water area, the salinized land increased in September.Since Ebinur Lake Watershed was characterized by little precipitation, strong evaporation, and much wind [26], higher temperatures cause more evaporation and worsen the arid conditions in the Watershed.Thus, in dry season, combined with a large evaporation, the salt was brought to the surface, causing salinization in the surroundings of Ebinur Lake.
Based on classification, the area of various LULC patterns showed obvious differences in six zones.In Figure 5 and Table 7, compared with the rainy season, the area of farmland obviously increased in each zone in the dry season; the area of forest-grass land and water body in each zone decreased.While the salinized land showed an increasing trend in each zone, the area of bare ground and bare rock in rainy season was less than that in the dry season in Zones 1, 2 and 6, and the area of others was greater than that in the dry season in Zones 3-5.Overall, the area of bare ground and bare rock changed little.Based on classification, the area of various LULC patterns showed obvious differences in six zones.In Figure 5 and Table 7, compared with the rainy season, the area of farmland obviously increased in each zone in the dry season; the area of forest-grass land and water body in each zone decreased.While the salinized land showed an increasing trend in each zone, the area of bare (i) (j)   Overall, LULC patterns are conspicuously different in the six zones (Figure 6).In the rainy season, large forest-grassland areas, over 50% of total land area, are distributed in Zone 4. Water areas have extensive distributions from 3.9% (Zone 5) to 8.1% (Zone 3).Around Ebinur Lake and in the Bohe and Jinghe tributaries, salinized land area ranges from 8.1% (Zone 2) to 15.7% (Zone 3).Gravel, bare ground and bare rock also represent a big proportion in Zones 1, 4 and 6, ranging from 53.1% (Zone 6) to 66.5% (Zone 4) of their respective drainage areas.
In dry season, farmland area ranges from 8.1% (Zone 2) to 22.9% (Zone 4).Most forest-grassland areas, over 20% of the total land area, are distributed in Zones 1, 2, 4 and 5. Water body is mainly found in Zones 2 and 3. Due to human activities such as farming irrigation, salinized land was concentrated in large farmland areas.Meanwhile, reduced precipitation in dry season leads to the contraction of Ebinur Lake and expansion of the salinized land around it.Large areas of salinized land are mainly found in Zone 2 and 3.In dry season, the phenomenon of salinization was more serious than in rainy season.Salinized land mainly contained the sodium, calcium, magnesium carbonate, sulfate and chloride.The proportion of salinized land is the primary predictor for mineralization values during the dry season.Bare ground and bare rock are at a high percentage and over 50% are distributed in Zones 1 and 6.In dry season, farmland area ranges from 8.1% (Zone 2) to 22.9% (Zone 4).Most forest-grassland areas, over 20% of the total land area, are distributed in Zones 1, 2, 4 and 5. Water body is mainly found in Zones 2 and 3. Due to human activities such as farming irrigation, salinized land was concentrated in large farmland areas.Meanwhile, reduced precipitation in dry season leads to the contraction of Ebinur Lake and expansion of the salinized land around it.Large areas of salinized land are mainly found in Zone 2 and 3.In dry season, the phenomenon of salinization was more serious than in rainy season.Salinized land mainly contained the sodium, calcium, magnesium carbonate, sulfate and chloride.The proportion of salinized land is the primary predictor for mineralization values during the dry season.Bare ground and bare rock are at a high percentage and over 50% are distributed in Zones 1 and 6.

Linkage between LULC Patterns and Hydro-Chemical Characteristics
First, we extracted the ground reflectance value of various LULC patterns and analyzed the correlation between the ground reflectance values and hydro-chemical characteristics (Table 8).The result showed that the ground reflectance values of various LULC patterns had good correlation with hydro-chemical characteristic parameters.Especially in rainy season, the result showed that the ground reflectance value of farmland had was negatively correlated with Cl − .Usually, Cl − was affected by fertilizer and organic fertilizer around the farmland.The soil salinity is higher in Xinjiang, chlorine salt in soil salt is greater.In recent years, with the development of the soil experiments on fertilizer formula, the chloride ion content in the fertilizer is generally high [53].In rainy season, the crops in the study area had not yet started to grow, so the fertilizer had less influence on Cl − .At the same time, the ground reflectance value of water body is related to pH under the significance level of 0.05.Meanwhile, the relation between the ground reflectance value of salinized land and hydro-chemical characteristics parameters were not significant, because the salinized land was not the main factor affecting hydro-chemical characteristics parameters in rainy season.In dry season, the ground reflectance values of salinized land had great impact on hydro-chemical characteristic parameters.With the decrease of the water area, the salinized land increased in dry season.Because it had higher reflectivity, the ground reflectance values of salinized land had great influence on hydro-chemical characteristic parameters.
Then, we discussed and analyzed the correlation between hydro-chemical characteristics and LULC patterns at different zones during rainy and dry seasons.The results are presented in Table 9.
During the dry season, farmland predicts more hydro-chemical characteristic variables than during the rainy season because of intensive fertilization and agricultural runoff [10,13,54].The results show that the proportion of farmland is positively correlated with pH (p < 0.01) and negatively related with EC, mineralization, Cl − , Ca 2+ and Na + (p < 0.01).Forest-grassland displays significant positive correlations with pH, EC, mineralization, Cl − and Na + (p < 0.01), and negatively correlations with Ca 2+ (p < 0.01).Water body in different zones is positively correlated with EC, mineralization, SO 4 2− , Ca 2+ , Mg 2+ and K + (p < 0.01), and negatively related with pH.Except the HCO 3 − , SO 4 2− , Ca 2+ , Mg 2+ and K + , the salinized land shows positive correlations with other hydro-chemical characteristics parameters.Stepwise linear multiple regressions models create a "goodness of fit" (R 2 values > 0.50) [55].EC, HCO 3 − , and SO 4 2− during the dry and rainy seasons are predicted based on the proportions of water body (Table 11).During the rainy season, pH, mineralization and Ca 2+ are defined by forest-grassland, while the Cl − and Mg 2+ are defined by salinized land.During the dry season, pH and Na + are estimated based on farmland, whereas Cl − and Ca 2+ are estimated based on forest-grassland.Specifically, mineralization is explained by salinity during the dry season.Yan et al. [56] discussed the formation of Ebinur Lake and its evolution process from the point of farmland, livestock, population and industry.His study gave preliminary analysis on the shrinking of the area of Ebinur Lake and the changing of watershed ecological environment over the recent 40 years.The results showed that the water consumption in Ebinur Lake area was negatively related to human activities, increasing vegetation damage caused by human activities, grassland degradation, desertification and salinization.Compared with those a decade ago, the mineralization and the ion content of Ebinur Lake both increased greatly [52].The phenomenon of water salinization should not be ignored.
Hydro-chemical characteristics is generally linked to LULC in the watershed [57].Human activities on LULC influence the types and degree of pollution.Therefore, measuring the proportions of certain LULC patterns in a watershed might enable us to conveniently predict hydro-chemical characteristics.

Conclusions
The results of this study demonstrate the relationships between LULC patterns and water quality of Ebinur Lake Watershed during dry and rainy seasons.The results show that: (1) During dry and rainy seasons, the distribution of pH was uniform in each zone and the range was mainly from 7.8 to 8.1.Overall, most of the variables showed significant spatial differences (p < 0.05) among different zones in rainy season and less significant differences in dry season.(2) Compared with the rainy season, as crops become ripe in the dry season, the area of farmland significantly increases in each zone; the area of forest-grassland and water body in each zone decreases.While the salinized land shows an increasing trend in each zone, the overall area of others changed little.(3) During the rainy season, pH, mineralized degree and Ca 2+ are defined by forest-grassland, while the Cl − and Mg 2+ are defined by salinized land.During the dry season, pH and Na + are defined based on farmland, whereas Cl − and Ca 2+ are defined based on forest-grassland.Especially, mineralized degree is explained by salinized land during the dry season.
From the results, the salinized land had great influence on hydro-chemical characteristics.In recent years, a serious salinization phenomenon existed in study area and the salt as rivers into pond, reservoir, Ebinur Lake, etc., resulting in the increase of ion content.Statistics indicate that salinized land in Ebinur Lake Watershed mainly distributed in Bortala River, Jing River, the surrounding villages and towns of Ebinur Lake, and north of Bole City.Furthermore, salinized land had seriously affected the farming, farmers had to use more fertilizer than before to increase yield.Thus, human activities to some extent affected ion concentration.
Most rivers in Xinjiang are characterized by low water yield, short flow, small water environmental capacity, and poor self-cleaning capability.Thus, scientifically utilizing and protecting the water resources are important issues that could help to achieve the sustainable development in Xinjiang.Thus, scientifically utilizing and protecting the water resources are important issues that could help to achieve the sustainable development in Xinjiang.
This study aims to explore the relationship between LULC and hydro-chemical characteristics.Due to the problem of data quality in the watershed, the data of 2014 were selected to explore the relationship between them.In the follow-up work, we hope to further explore the relationship between LULC and hydro-chemical characteristics in recent years, and compare the effects of land use change on hydro-chemical characteristics.

Figure 1 .
Figure 1.Location and Landscape Photos of the Study Area: (a) The water environment of Ebinur Lake Watershed; (b) The salinized land of Ebinur Lake Watershed; and (c) Ganjia Lake Haloxylon Forest National Nature Reserve.

Figure 1 .
Figure 1.Location and Landscape Photos of the Study Area: (a) The water environment of Ebinur Lake Watershed; (b) The salinized land of Ebinur Lake Watershed; and (c) Ganjia Lake Haloxylon Forest National Nature Reserve.

Figure 2 .
Figure 2. Conceptual model for the methodology.

Figure 2 .
Figure 2. Conceptual model for the methodology.

Figure 4 .
Figure 4.The change of Land use/cover (LULC) in Ebinur Lake Watershed during the rainy and dry seasons in 2014: (a) May 2014 remote sensing image; (b) LULC classification in May 2014; (c) September 2014 remote sensing image; (d) LULC classification in September 2014.

Figure 4 .
Figure 4.The change of Land use/cover (LULC) in Ebinur Lake Watershed during the rainy and dry seasons in 2014: (a) May 2014 remote sensing image; (b) LULC classification in May 2014; (c) September 2014 remote sensing image; (d) LULC classification in September 2014.

Figure 6 .
Figure 6.Composition (%) of land use patterns in 6 different zones in Ebinur Lake Watershed: (a) Composition (%) of land use patterns in May 2014; (b) Composition (%) of land use patterns in September 2014.

Figure 6 .
Figure 6.Composition (%) of land use patterns in 6 different zones in Ebinur Lake Watershed: (a) Composition (%) of land use patterns in May 2014; (b) Composition (%) of land use patterns in September 2014.

Table 1 .
Landsat OLI Coverage of Ebinur Lake Watershed in May and September of 2014.

Table 1 .
Landsat OLI Coverage of Ebinur Lake Watershed in May and September of 2014.

Table 2 .
Summary of hydro-chemical characteristics parameters observations of Ebinur Lake Watershed during rainy and dry seasons.

Table 2 .
Summary of hydro-chemical characteristics parameters observations of Ebinur Lake Watershed during rainy and dry seasons.

Table 3 .
The statistical results of hydro-chemical characteristic parameters from ANOVA.

Table 4 .
The sample quantity in different zones.

Table 5 .
The statistical results of mean ± S.E. in different zones.

Table 6 .
Factor loadings of the 10 variables on rotated component matrix for all hydro-chemical characteristic data during dry and rainy seasons in Ebinur Lake Watershed.

Table 7 .
The area of LULC patterns in six zones of Ebinur Lake Watershed.

Table 8 .
Correlation between the ground reflectance and hydro-chemical characteristics during rainy and dry seasons.

Table 10 .
Correlation between LULC patterns and hydro-chemical characteristics in the 500 m buffers during rainy and dry seasons (56 points in rainy season and 57 points in dry season).

Table 11 .
Stepwise linear multiple regressions models for hydro-chemical characteristics parameters and LULC patterns during dry and rainy seasons in Ebinur Lake Watershed.