Responses of Sediment Yield to Vegetation Cover Changes in the Poyang Lake Drainage Area, China

Ascertaining the relationships between sediment transport processes and vegetation cover is essential for watershed soil and water conservation. However, it is not easy to realize this target on a large scale. In this study, the location-weighted landscape contrast index (LCI) based on the " source-sink " theory of ecological processes was introduced to unravel the response of sediment yield to vegetation cover changes in the Poyang Lake drainage area (the largest freshwater lake in China). A modified location-weighted landscape contrast index (mLCI) was proposed to improve the accuracy and efficiency of the LCI. The average vegetation fraction of the study area significantly increased from 46. the area-specific sediment yield (ASY) sharply decreased by 55%. In Period II, the " source " vegetation patches were found further away from the water bodies, and were more likely to be located on gentler slopes than in Period I. From Period II to Period III, the ASY sharply increased by 83%. In Period III, " source " patches were found closer to the water bodies than in Period II, whereas the " sink " patches were found further away from the water bodies. The high statistical correlation between LCIs/mLCIs and ASY indicated a sensitive response of the sediment yield to vegetation cover changes that significantly altered the sediment transport processes in the study area. The ASY was better correlated with the mLCI than with the LCI. Three key harnessing sub-watersheds of the Ganjiang watershed were identified by calculating the mLCIs: Qingfengshan, Yuanhe, and Jinjiang. This study supported watershed ecological management in the Poyang Lake drainage area and provided a methodology reference for future sediment transport process studies.


Introduction
Sediment loss by water is a serious environment problem that is threatening the sustainability of ecological system [1].The Poyang Lake drainage area (covering 122,100 km 2 ) has contributed greatly to the socioeconomic development of China.However, this region has a fragile ecosystem [2].The dominant soil type in this region belongs to red soil with a low infiltration rate, which can prompt surface runoff and sediment fluxes [3].Previous studies indicated that serious soil erosion by water caused severe environmental degradation in the Poyang Lake drainage area [4,5].For decades, this region has exhibited the severest soil erosion in south China [4].The area of soil erosion occupies approximately 20.03% of the total territory, and 99.45% of the soil erosion is caused by water erosion [6].In recent years, the annual amount of sediment yield discharged into Poyang Lake reached 12 million metric tons, with a mean bottom accretion rate of up to 3 mm/a [7].The regulation and storage functions of Poyang Lake have gradually reduced.
Vegetation cover exerts a significant influence on runoff-sediment transportation processes [8][9][10][11][12].Over the past few decades, the Poyang Lake drainage area has experienced a significant change in vegetation coverage [13].Developing an accurate and highly efficient method to elucidate the relationship between vegetation cover changes and sediment transportation processes is of practical significance for watershed management in this region.Manual field monitoring can provide fundamental data to reveal ecological processes associated with landscape changes.However, the use of this approach on a large scale requires a great expenditure of manpower and time.
During the last two decades, ecologists have achieved great progress in landscape spatial data acquisition and improved methods for quantitative analysis using readily accessible remote sensing data and significant advances in GIS technology.From a landscape ecological perspective, the ecological function of landscape patterns is often studied by investigating the statistical relationships between landscape indices and ecological conditions (e.g., exposure to contaminants, diversity and biological balance and physical/chemical attributes) [10,14,15].The patch area, patch shape index, Shannon-Weaver diversity index, dominance index, and fragmentation index are commonly used landscape indices at both the slope and catchment (or regional) scales [8,[16][17][18][19][20].However, with the development of landscape ecology, elucidating ecological processes and their interactions with landscape patterns has become the focus of recent landscape ecology studies [21][22][23].These landscape indices only take into account the composition and spatial configuration of the land cover and do not provide insights into the ecological processes [8,10].Some landscape indices have been established from an ecological process perspective.For example, Borselli et al. [24] established an index of connectivity (IC) that represented a connectivity assessment based on information from the landscape.Mayor et al. [25] derived a spatial index (Flowlength) to quantify the connectivity in runoff source areas considering both the topography and vegetation pattern.These indices effectively address the runoff-sediment process.However, they are both heavily dependent on high-resolution spatial data and thus are not suitable for large-scale studies.
Chen et al. [8,26] established the location-weighted landscape contrast index (LCI) (see the methodology section for a thorough description of this method) to represent the ecological functions of landscape patterns based on the "source-sink" theory which has been commonly used to address the response of ecological process to landscape pattern changes [27,28].Remarkably, the LCI is scale-independent (suitable for different sized watersheds) and regards not only the spatial configuration of the landscape but also the topographical characteristics of each land patch.Its advantages could also benefit from the readily available input data (low cost and wide coverage remote sensing images).Additionally, the series software ArcGIS provides a convenient environment for data processing.Previous studies employing the LCI primarily focused on NPS processes [29].For example, Yue et al. [30] applied the LCI to analyze the influence of "source" and "sink" landscape spatial changes based on the total nitrogen (TN) and total phosphorus (TP) contents to assess water body quality in Shenzhen city, China.Liu et al. [31] studied the effects of land cover patterns on TN and TP in nine typical districts of the upper reaches of the Yangtze River using the LCI.These two studies both found that the LCI effectively quantified the contributions of landscape patterns to NPS processes at the watershed scale.
In a watershed, sediments generated from "source" land patches under the exogenic force (e.g., water flow) will be transported into water body (eventually exported out of the watershed) with part or all of them retarded by "sink" land patches.If the pattern of "source" and "sink" patches changes, the sediment yield observed at the outlet of the watershed should respond.Under this hypothesis, we applied the LCI to investigate the response of sediment yield to vegetation cover in the Poyang Lake drainage area, and addressed the following research questions: (1) What are the change characteristics of vegetation patterns related to topography in the study area?(2) How does the sediment yield respond to vegetation patterns changes?(3) Are there any critical sediment loss sub-watersheds that required special attention in the study area?

Study Area Description
The Poyang Lake drainage area consists of five river watersheds and the lake region (Poyang Lake and its vicinity areas).We selected five river watersheds (Ganjiang, Xiushui, Fuhe, Xinjiang and Raohe) as the study area.Each watershed is named after its chief river.As shown in Figure 1, the study area is distributed between 113 ˝45 1 and 118 ˝35 1 E longitude and 24 ˝32 1 and 29 ˝36 1 N latitude.These five rivers all flow into the largest freshwater lake in China (Poyang Lake), and eventually into the Yangtze River (Changjiang).The hydrological station names and locations can be identified by combining the information in Figure 1 and Table 1.The Ganjing watershed contains fifteen sub-watersheds.Nevertheless, these sub-watersheds lack hydrology measurements and the sediment yield data are inaccessible.The whole study area covers 122,100 km 2 and is characterized by a typical moist, subtropical climate with an average annual temperature of 17.7 ˝C, an average annual rainfall of 1786 mm and a maximum daily temperature of approximately 40 ˝C.The distribution of annual precipitation is uneven, with approximately 50% occurring in the summer when the temperature is relatively high.The north and min-north part of the study area has the lowest elevation, with smoothly undulating hilly lands and flat alluvial plains as the predominate relief accounting for approximately 35% of the total area.With mountains on the other three sides, the whole shape of the territory resembles a "dustpan".The whole study area has symmetrical vegetation coverage.The vegetation fraction in the flat plains with high human disturbance is low.In contrast, the mountain areas are usually covered by thick forests.The dominant soil that is widely distributed in the region belongs to red soil developed from quaternary sediments.This type of soil has a low precipitation infiltration rate, and is degraded and erosion-prone [3].Due to the possibility of severe soil loss, several soil and water conservation measures have been implemented.The most famous is the Grain-for-Green Program that promotes ecological reconstruction via the conversion of steeply sloped farmland into forest and grass [13].At present, the sediment loss in the mountain areas covered by a high vegetation fraction is under control [4].

Study Area Description
The Poyang Lake drainage area consists of five river watersheds and the lake region (Poyang Lake and its vicinity areas).We selected five river watersheds (Ganjiang, Xiushui, Fuhe, Xinjiang and Raohe) as the study area.Each watershed is named after its chief river.As shown in Figure 1, the study area is distributed between 113°45′ and 118°35′ E longitude and 24°32′ and 29°36′ N latitude.These five rivers all flow into the largest freshwater lake in China (Poyang Lake), and eventually into the Yangtze River (Changjiang).The hydrological station names and locations can be identified by combining the information in Figure 1 and Table 1.The Ganjing watershed contains fifteen sub-watersheds.Nevertheless, these sub-watersheds lack hydrology measurements and the sediment yield data are inaccessible.The whole study area covers 122,100 km 2 and is characterized by a typical moist, subtropical climate with an average annual temperature of 17.7 °C, an average annual rainfall of 1786 mm and a maximum daily temperature of approximately 40 °C.The distribution of annual precipitation is uneven, with approximately 50% occurring in the summer when the temperature is relatively high.The north and min-north part of the study area has the lowest elevation, with smoothly undulating hilly lands and flat alluvial plains as the predominate relief accounting for approximately 35% of the total area.With mountains on the other three sides, the whole shape of the territory resembles a "dustpan".The whole study area has symmetrical vegetation coverage.The vegetation fraction in the flat plains with high human disturbance is low.In contrast, the mountain areas are usually covered by thick forests.The dominant soil that is widely distributed in the region belongs to red soil developed from quaternary sediments.This type of soil has a low precipitation infiltration rate, and is degraded and erosion-prone [3].Due to the possibility of severe soil loss, several soil and water conservation measures have been implemented.The most famous is the Grain-for-Green Program that promotes ecological reconstruction via the conversion of steeply sloped farmland into forest and grass [13].At present, the sediment loss in the mountain areas covered by a high vegetation fraction is under control [4].The LCI does not consider the effects of rainfall variation on sediment loss.Therefore, to evaluate the relationship between vegetation cover changes and sediment loss, rainfall variation needed to be eliminated prior to the implementation of our study.We used rainfall erosivity to represent the ability of rainfall to cause soil erosion by considering the rainfall amount and intensity [32].Based on the annual spatial data products of rainfall erosivity (1980-2014) developed by Xiao et al. [33] using the daily precipitation dataset provided by the China Meteorological Data Sharing Service System (http://cdc.nmic.cn/home.do),we investigated the annual rainfall erosivity trends in the five watersheds (Ganjiang, Xiushui, Fuhe, Xinjiang and Raohe) during the 34 year period.As shown in Table 1, three time periods (1992-1994, 2004-2006, and 2011-2013) with similar mean annual rainfall erosivities (MAEs) were selected.The MAE values exhibited in Table 1 had an average of 10,665.6 (MJ¨mm)/(km 2 ¨h¨a) and a standard deviation (SD) of 753.0 (MJ¨mm)/(km 2 ¨h¨a).The relative change between each MAE value and the average value did not exceed ˘10% with the exception of Raohe in Period III (14%).In contrast, the MAE values in the other periods (one period = three years) that were not selected varied greatly (SD = 1404.1 (MJ¨mm)/(km 2 ¨h¨a)) and were not close to the MAE in the selected periods.We believed that the sediment yield changes among these three study periods were primarily caused by changes in the underlying surface rather than rainfall.To increase the representativeness and eliminate the random deviation, the area-specific sediment yield (ASY, the ratio of total sediment yield to the watershed area) and rainfall erosivity for the time Periods I, II and III were included as three-year average values.For example, we used the mean annual ASY over 1992-1994 to represent the situation in Period I.The corresponding sediment yield data from the five watersheds were extracted from the Changjiang sediment bulletins [34-39] and the study of Sun et al. [40].This sediment yield data were observed by the hydrological stations located at the outlet of each watershed.

Remote Sensing Products
The MODIS 16-day vegetation index (VI) products (MOD13Q1 with a spatial resolution of 250 m) of Period II and Period III were obtained from the Goddard Space Flight Center (GSFC) (https://ladsweb.nascom.nasa.gov/index.html).All of the VI products used in our study were collected on the same Julian day (257) when the vegetation was at its utmost luxuriance.Landsat-5 Thematic Mapper (TM) images (with a spatial resolution of 30 m) were obtained from the Institute of Remote Sensing and Digital Earth (IRSD) (http://ids.ceode.ac.cn/query.html)to produce the VI for Period I.The acquisition day for all of the TM images applied in our study was near the Julian day of the MOD13Q1 products.The two sub-pixel model was applied to produce the vegetation fraction (VF) based on the VI products.The details are available in study of Miao et al. [41].We calculated the three-year average VF value for each time period.For convenience, the VF of the study area was classified into five levels: high, moderately high, moderate, moderately low, and low (ě80%, 60%-80%, 40%-60%, 20%-40%, and <20%, respectively).The Shuttle Radar Topography Mission (SRTM) version 2 data were obtained from the Jet Propulsion Laboratory (http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp).All spatial raster data were resampled into 250 m to unify the spatial resolution (Figure 2).

The "Source-Sink Theory" of Sediment Transport
The "source-sink" theory refers to the function of the "source" and the "sink" exerted by the landscape [8].When integrated with a specific ecological process (sediment transport), "source" landscapes are those that can promote sediment production and transport, whereas "sink" landscapes refer to those that can retard or delay sediment transport.Naturally, "sink" landscapes belong to land patches with high vegetation fractions (e.g., forestland) that can effectively reduce sediment loss [42].Notably, "sink" landscapes may not arrest all sediment.In other word, sediments may pass the "sink" landscapes with only part of them arrested.As shown in Figure 3, sediments generated from "source" land patches will be transported into the river (eventually exported out of the watershed) with part or all of them arrested by "sink" land patches.Theoretically, if the spatial distribution of the "source" and "sink" landscapes is balanced and reasonable, less sediment will be transported into the river, and a lower sediment yield will be observed by the hydrological station at the outlet of the watershed.Based on this theory, the LCI index could be suitable to evaluate the rationality of the landscape distribution aimed at sediment conservation by considering the area of the "source" and "sink" landscapes and their topography-related patterns as well as their efficiency of generating or arresting sediments.If the correlation coefficient of the LCI and sediment yield at the outlet of the watershed passes the significance testing, the LCI index possesses practicability.

The "Source-Sink Theory" of Sediment Transport
The "source-sink" theory refers to the function of the "source" and the "sink" exerted by the landscape [8].When integrated with a specific ecological process (sediment transport), "source" landscapes are those that can promote sediment production and transport, whereas "sink" landscapes refer to those that can retard or delay sediment transport.Naturally, "sink" landscapes belong to land patches with high vegetation fractions (e.g., forestland) that can effectively reduce sediment loss [42].Notably, "sink" landscapes may not arrest all sediment.In other word, sediments may pass the "sink" landscapes with only part of them arrested.As shown in Figure 3, sediments generated from "source" land patches will be transported into the river (eventually exported out of the watershed) with part or all of them arrested by "sink" land patches.Theoretically, if the spatial distribution of the "source" and "sink" landscapes is balanced and reasonable, less sediment will be transported into the river, and a lower sediment yield will be observed by the hydrological station at the outlet of the watershed.Based on this theory, the LCI index could be suitable to evaluate the rationality of the landscape distribution aimed at sediment conservation by considering the area of the "source" and "sink" landscapes and their topography-related patterns as well as their efficiency of generating or arresting sediments.If the correlation coefficient of the LCI and sediment yield at the outlet of the watershed passes the significance testing, the LCI index possesses practicability.

Determination of the "Source" and "Sink" Patches
Land patches with high VF can reduce the sediment transport process, whereas those with low VFs are more likely to be eroded.Therefore, vegetation coverage can highly represent the land patch conditions that influence sediment transport process [10][11][12].In our study, the VF produced in Section 2.2.2 was used to distinguish the "source" and "sink" land patch types.In October 2013, we performed a field survey of the study area.An accumulated survey length of 3063 km that included 1336 photos and a total of 93 patches with latitude and longitude information were recorded in the field survey database (Figure 1).Based on this data, we determined that a VF of approximately 60% could be considered the threshold at which the sediment transport process was effectively reduced.This finding coincided with the studies of Yuan et al.
[43] and Wu et al. [44].Hence, land patches with high and moderately high VFs were categorized as "sink" types, and land patches with low, moderately low and moderate VFs were categorized as "source" types.As shown in the Results Section, this threshold proved to be rational.However, it is essential to emphasize that the threshold is related to the specific vegetation cover and the sediment transport processes in the Poyang Lake drainage area, and cautions should be taken when extending this threshold to other regions.

Description of the Topographical-Related Spatial Patterns of "Source" and "Sink"
The topography characteristics of land patches can be important factors that influence sediment transportation on terrestrial surfaces (e.g., elevation, slope, and flow length) [10][11][12]45,46].Theoretically, sediments generated from "source" land patches with a shorter flow length have a greater chance of being transported into water body, thereby contributing to a higher sediment yield at the watershed outlet.The "sink" patches that are distributed near the water body have greater chance to arrest sediment than ones far from water body.Elevation and distance to water have similar functions with regards to flow length.Nevertheless, a steeper slope can enhance sediment transportation [47] and have an opposing function to elevation and distance to water.Simply speaking, a steeper slope, lower elevation and shorter flow length and distance to water result in a greater influence of "source" or "sink" land patches on the promotion or reduction of the sediment yield loss.
The Lorenz curve was applied to depict the topographical distribution characteristics of the "source" and "sink" land patches in a watershed (Figure 4).Using the distance to water as an example, O indicates the water body (e.g., river; distance = 0), C indicates the maximum distance, OC represents the increasing process of distance to water, OA is the corresponding accumulated area percentage of a type of "source" or "sink" land patch, ODBC refers to a Lorenz curve of a type of "source" or "sink" land patch that is primarily distributed near the water body and thus has more influence of the promotion or reduction of the sediment yield in water, OFBC is a case opposite to ODBC, and OEBC represents an even distribution case.Consequently, the area of the irregular polygon (A) encompassed by the Lorenz curve of a type of "source" or "sink" land patch is proportional to its influence on the sediment yield observed at the hydrological station at the outlet of the watershed.

Determination of the "Source" and "Sink" Patches
Land patches with high VF can reduce the sediment transport process, whereas those with low VFs are more likely to be eroded.Therefore, vegetation coverage can highly represent the land patch conditions that influence sediment transport process [10][11][12].In our study, the VF produced in Section 2.2.2 was used to distinguish the "source" and "sink" land patch types.In October 2013, we performed a field survey of the study area.An accumulated survey length of 3063 km that included 1336 photos and a total of 93 patches with latitude and longitude information were recorded in the field survey database (Figure 1).Based on this data, we determined that a VF of approximately 60% could be considered the threshold at which the sediment transport process was effectively reduced.This finding coincided with the studies of Yuan et al.
[43] and Wu et al. [44].Hence, land patches with high and moderately high VFs were categorized as "sink" types, and land patches with low, moderately low and moderate VFs were categorized as "source" types.As shown in the Results Section, this threshold proved to be rational.However, it is essential to emphasize that the threshold is related to the specific vegetation cover and the sediment transport processes in the Poyang Lake drainage area, and cautions should be taken when extending this threshold to other regions.

Description of the Topographical-Related Spatial Patterns of "Source" and "Sink"
The topography characteristics of land patches can be important factors that influence sediment transportation on terrestrial surfaces (e.g., elevation, slope, and flow length) [10][11][12]45,46].Theoretically, sediments generated from "source" land patches with a shorter flow length have a greater chance of being transported into water body, thereby contributing to a higher sediment yield at the watershed outlet.The "sink" patches that are distributed near the water body have greater chance to arrest sediment than ones far from water body.Elevation and distance to water have similar functions with regards to flow length.Nevertheless, a steeper slope can enhance sediment transportation [47] and have an opposing function to elevation and distance to water.Simply speaking, a steeper slope, lower elevation and shorter flow length and distance to water result in a greater influence of "source" or "sink" land patches on the promotion or reduction of the sediment yield loss.
The Lorenz curve was applied to depict the topographical distribution characteristics of the "source" and "sink" land patches in a watershed (Figure 4).Using the distance to water as an example, O indicates the water body (e.g., river; distance = 0), C indicates the maximum distance, OC represents the increasing process of distance to water, OA is the corresponding accumulated area percentage of a type of "source" or "sink" land patch, ODBC refers to a Lorenz curve of a type of "source" or "sink" land patch that is primarily distributed near the water body and thus has more influence of the promotion or reduction of the sediment yield in water, OFBC is a case opposite to ODBC, and OEBC represents an even distribution case.Consequently, the area of the irregular polygon (A) encompassed by the Lorenz curve of a type of "source" or "sink" land patch is proportional to its influence on the sediment yield observed at the hydrological station at the outlet of the watershed.
We can calculate the total A of the "source" patches (three types) in a watershed as follow: and that of the "sink" (two types) as follow: where C is the maximum value of the topographical index.Using the distance to water between the land patches and water body as an example; R xi and S xj indicate the accumulated areal percentages of the Lorenz curves of the "source" and "sink" land patch types, respectively; and M and N represent the numbers of the "source" and "sink" land patch types in the watershed, respectively.A contrast index (CI) can be established to assess the balanced states of the "source" and "sink" patches in a watershed: When CI > 1, the "source" patches are more likely to be distributed near the water than the "sink" patches, and thus contribute more to the sediment loss.When CI < 1, the effect is opposite to the case of CI > 1.When CI = 1, the spatial distribution of the "source" and "sink" landscapes is balanced.We can calculate the total A of the "source" patches (three types) in a watershed as follow: and that of the "sink" (two types) as follow: where C is the maximum value of the topographical index.Using the distance to water between the land patches and water body as an example; Rxi and Sxj indicate the accumulated areal percentages of the Lorenz curves of the "source" and "sink" land patch types, respectively; and M and N represent the numbers of the "source" and "sink" land patch types in the watershed, respectively.A contrast index (CI) can be established to assess the balanced states of the "source" and "sink" patches in a watershed: When CI > 1, the "source" patches are more likely to be distributed near the water than the "sink" patches, and thus contribute more to the sediment loss.When CI < 1, the effect is opposite to the case of CI > 1.When CI = 1, the spatial distribution of the "source" and "sink" landscapes is balanced.

Location-weighted landscape contrast index
In addition to the topography-related patterns, the "source" and "sink" patch areas in a watershed and their efficiency in generating or arresting sediments are important for the sediment transport process.By considering the latter two factors, the final LCI calculation formula aimed at the sediment transport process in a watershed can be written as: where and are the effectiveness parameters of promoting and preventing sediment transport by the "source" and "sink" land patch types, respectively; and and represent the area percentages of the "source" and "sink" land patch types in the watershed, respectively.When LCI = 0, the "source" and "sink" patches reach a balance.When LCI > 0, the contribution of the "source" is greater than the contribution of the "sink", and the vegetation cover in the watershed tends to generate and transport sediments.When LCI < 0, the situation is reversed.
The accurate determination of and can be difficult.We used the vegetation coverage and management factor (C) that were usually implemented in USLE and its revised form RUSLE to

Location-Weighted Landscape Contrast Index
In addition to the topography-related patterns, the "source" and "sink" patch areas in a watershed and their efficiency in generating or arresting sediments are important for the sediment transport process.By considering the latter two factors, the final LCI calculation formula aimed at the sediment transport process in a watershed can be written as: where ϕ i and ϕ j are the effectiveness parameters of promoting and preventing sediment transport by the "source" and "sink" land patch types, respectively; and p i and p j represent the area percentages of the "source" and "sink" land patch types in the watershed, respectively.When LCI = 0, the "source" and "sink" patches reach a balance.When LCI > 0, the contribution of the "source" is greater than the contribution of the "sink", and the vegetation cover in the watershed tends to generate and transport sediments.When LCI < 0, the situation is reversed.
The accurate determination of ϕ i and ϕ j can be difficult.We used the vegetation coverage and management factor (C) that were usually implemented in USLE and its revised form RUSLE to approximately represent the effectiveness of the "source" land patch type parameter.The ϕ i of the moderate, moderately low and low VFs was set at 0.15, 0.2 and 0.4, respectively, according to the relevant achievements of the C values for the different VF levels [48,49].Considering the opposing meanings of ϕ i and ϕ j , we assumed that the functions of a high VF in reducing sediment and a low VF in causing sediment could be approximately counterbalanced; a similar assumption was made concerning a moderately high VF versus a moderately low VF.Thus, the ϕ j values of a high VF and a moderately high VF were set at 0.4 and 0.2, respectively.

A modified Topographic Index
In previous ecological process studies using the LCI, the commonly considered topographical indices (e.g., elevation, slope, flow length, and distance to water) were individually introduced into the LCI calculation.However, some indices, such as elevation, flow length, and distance to water, are highly correlated in most cases (Figure 5), which always leads to similar LCI results [30,31,50].This factor may result in redundancy and reduce the research efficiency.In this study, we proposed a modified topographical index that combine the two topographical indices with relatively low spatial correlativity in the hope of improving the LCI and the research efficiency.The Euclidean distance (ED) to water was applied to replace the elevation and flow length.The slope had an opposing function compared to the sediment transport process to ED and was processed into the modified index.The modified ED (mED) index can be expressed as: where D is the Euclidean distance from land patches to water; D max is the farthest Euclidean distance in a given watershed; S is the slope; and S max is the maximum slope in a given watershed.
We introduce mED into Equation ( 4) to calculate the modified location-weighted landscape contrast index (mLCI), and test its applicability by comparing it with LCI-elevation, LCI-slope, LCI-flow length, and LCI-Euclidean distance.All of the topographical indices utilized in this study were produced in the GIS environment based on the SRTM-DEM data.approximately represent the effectiveness of the "source" land patch type parameter.The of the moderate, moderately low and low VFs was set at 0.15, 0.2 and 0.4, respectively, according to the relevant achievements of the C values for the different VF levels [48,49].Considering the opposing meanings of and , we assumed that the functions of a high VF in reducing sediment and a low VF in causing sediment could be approximately counterbalanced; a similar assumption was made concerning a moderately high VF versus a moderately low VF.Thus, the values of a high VF and a moderately high VF were set at 0.4 and 0.2, respectively.

A modified Topographic Index
In previous ecological process studies using the LCI, the commonly considered topographical indices (e.g., elevation, slope, flow length, and distance to water) were individually introduced into the LCI calculation.However, some indices, such as elevation, flow length, and distance to water, are highly correlated in most cases (Figure 5), which always leads to similar LCI results [30,31,50].This factor may result in redundancy and reduce the research efficiency.In this study, we proposed a modified topographical index that combine the two topographical indices with relatively low spatial correlativity in the hope of improving the LCI and the research efficiency.The Euclidean distance (ED) to water was applied to replace the elevation and flow length.The slope had an opposing function compared to the sediment transport process to ED and was processed into the modified index.The modified ED (mED) index can be expressed as: where D is the Euclidean distance from land patches to water; D is the farthest Euclidean distance in a given watershed; S is the slope; and S is the maximum slope in a given watershed.We introduce mED into Equation ( 4) to calculate the modified location-weighted landscape contrast index (mLCI), and test its applicability by comparing it with LCI-elevation, LCI-slope, LCI-flow length, and LCI-Euclidean distance.All of the topographical indices utilized in this study were produced in the GIS environment based on the SRTM-DEM data.

Temporal and Spatial Variations in Vegetation Covers
Figure 2 exhibits the VF distribution of the study area in the three time periods (1992-1994, 2004-2006, and 2011-2013).The average VFs of the study area in these three periods were 46.2%, 76.5%, and 73.2%, respectively.All of the watersheds (Ganjiang, Xiushui, Fuhe, Xinjiang and Raohe) experienced a sharp increase and then a slight decrease in the average VF (Table 2).We investigated the spatial distribution characteristics of the five VF levels (high, moderately high, moderate, moderately low, and low) in all watersheds in the three periods.The Xinjiang watershed was selected as an example due to its strong representativeness.Taking into account the high spatial correction among the Euclidean distance, elevation, and flow length (Figure 5), slope and Euclidean distance were selected for the spatial distribution investigation of the VFs.As shown in Figure 6, there was no high VF level in Period I. "Source" VFs were more likely to be located on gentler slopes in Period II than in Period I.In Period II, "sink" VFs were also more likely to be located on gentler slopes than in Period I; however, the "sink" VFs were not as evident as the "source" VFs.In Period III, both "source" and "sink" patches (excluding high level) were more likely to be located on steeper slopes.With regards to the Euclidean distance to water, in Period II, the "source" VFs were found further away from the water bodies than in Period I.The curve of the moderately high VFs experienced almost no change from Period I to II.In Period III, the low VFs were found closer to the water bodies, whereas, the high VFs found further away from the water bodies than in Period II.Generally, the change in the spatial distribution of VFs increased the likelihood for sediment entrapment from Period I to Period II.Nevertheless, it decreased the likelihood for sediment entrapment to some extent during Period II to III.

Temporal and Spatial Variations in Vegetation Covers
Figure 2 exhibits the VF distribution of the study area in the three time periods (1992-1994, 2004-2006, and 2011-2013).The average VFs of the study area in these three periods were 46.2%, 76.5%, and 73.2%, respectively.All of the watersheds (Ganjiang, Xiushui, Fuhe, Xinjiang and Raohe) experienced a sharp increase and then a slight decrease in the average VF (Table 2).We investigated the spatial distribution characteristics of the five VF levels (high, moderately high, moderate, moderately low, and low) in all watersheds in the three periods.The Xinjiang watershed was selected as an example due to its strong representativeness.Taking into account the high spatial correction among the Euclidean distance, elevation, and flow length (Figure 5), slope and Euclidean distance were selected for the spatial distribution investigation of the VFs.As shown in Figure 6, there was no high VF level in Period I. "Source" VFs were more likely to be located on gentler slopes in Period II than in Period I.In Period II, "sink" VFs were also more likely to be located on gentler slopes than in Period I; however, the "sink" VFs were not as evident as the "source" VFs.In Period III, both "source" and "sink" patches (excluding high level) were more likely to be located on steeper slopes.With regards to the Euclidean distance to water, in Period II, the "source" VFs were found further away from the water bodies than in Period I.The curve of the moderately high VFs experienced almost no change from Period I to II.In Period III, the low VFs were found closer to the water bodies, whereas, the high VFs found further away from the water bodies than in Period II.Generally, the change in the spatial distribution of VFs increased the likelihood for sediment entrapment from Period I to Period II.Nevertheless, it decreased the likelihood for sediment entrapment to some extent during Period II to III.When the LCIs/mLCIs, average VF and ASY of the five watersheds in the three periods were calculated, we investigated their coupling relationships using a scatter map (Figure 7).Theoretically, there should be a negative correlation between the average VF and ASY.As shown in Figure 7, the coefficient of determination (R 2 ) yielded using linear regression analysis was ´0.24, and did not pass the significance testing (p > 0.05).In contrast, the LCIs/mLCIs (with the exception of LCI-slope) exhibited a highly significant correlation with the ASY.Based on the LCI concept, the ASY not only responded to the change of proportion in the VF levels but also the change in their spatial patterns.Using the Xinjiang watershed as an example, from Period I to Period II, the ASY was sharply reduced from 167 to 53.1 t/(km 2 ¨a) (Figure 7).Accordingly, the watershed experienced a great improvement in vegetation coverage in these two time periods.The areal percentage of high and moderately high levels increased from 3.24% and 0.01% to 55.5% and 21.01%, respectively (Table 2).During the same time period, the changes in vegetation pattern of the Xinjiang watershed increased the likelihood for sediment entrapment, and thus prompted a sharp reduction in the ASY.From Period II to Period III, there was no evident change in the VF (Table 2).However, the increase in the ASY from 53.1 to 75.9 t/(km 2 ¨a) (Figure 7) was relatively remarkable.Therefore, the change in the vegetation patterns (which decreased the likelihood for sediment entrapment as seen in Section 4.1) should be awarded the most responsibility (refer to the Discussion Section for details).As mentioned above, a steeper slope can enhance the sediment generation of the "source" landscape.Theoretically, LCI-slope should have a negative correlation with the ASY.However, the vegetation patches with steeper slope were usually further away from the water bodes (Figure 5).Therefore, the negative correlation could be weakened (p > 0.05).
Generally, the ASY, which was highly related to the LCI/mLCI, showed a sensitive response to vegetation cover changes in the Poyang Lake drainage area.In comparison, we found that the mLCI considering the Euclidean distance modified by slope information had the best correlation with the ASY (R 2 = 0.79, p < 0.01).This result illustrates the advancement of the topographical index of the mED (Figure 8) established in the present study.We can also conclude that, compared to other LCIs (LCI-elevation, LCI-slope, LCI-flow length, and LCI-Euclidean distance), mLCI is the optimum index to indicate the sediment loss risk in the study area.

The mLCIs of the Sub-Watersheds of the Ganjiang Watershed
We applied the mLCIs as the judgment indices of sediment loss risk for the sub-watersheds of Ganjiang, where measurements of sediment yield are sparse.As shown in Table 3, Qingfengshan had the highest mLCI in the three periods (0.45, ´0.07, and 0.04, respectively), followed by Yuanhe (0.23, ´0.28, and ´0.06, respectively) and Jinjiang (0.17, ´0.22, and ´0.07, respectively).The latter two watersheds showed a similar change trend in their mLCIs.These three sub-watersheds of Ganjiang should receive special attention for soil conservation.Sub-watersheds possessing relatively excellent soil conservation statuses were also identified.The mLCIs of Suichuanjiang, Meijiang, and Lian-Mian were lower than other sub-watersheds, and Suichuanjiang was characterized by the lowest mLCIs (0.14, ´0.66, and ´0.45 in Period I, II, and III, respectively).

Discussion
It is conceivable that vegetation cover can strongly influence the runoff process that transports the sediments into water (e.g., river) [10][11][12].Indeed, the vegetation coverage conditions in the Poyang Lake drainage area changed significantly from Period I to Period III.This result may be partially due to the implementation of the "Grain for Green" program and the cultivated land balance policy [13,51,52].Nevertheless, studying the response of sediment yield to vegetation cover change can promote the development of watershed management in the Poyang Lake drainage area.
The LCI index applied in our study is a comprehensive landscape index based on the "source-sink" theory of ecological processes.When the LCI is calculated, the contribution of "source" and "sink" land patches to the sediment loss in watershed is addressed.Therefore, the LCI could be directly related to the sediment yield observed at the outlet of the watershed.Our study confirmed the high accuracy of the LCI in reflecting the sediment yield.For example, there was a remarkable increase in the ASY (from 53.1 to 75.9 t/(km 2 ¨a), representing a 42.9% increase) in the Xinjiang watershed from Period II to Period III.However, the average vegetation fraction decreased by only 3.7%.Thus, the average vegetation fraction change provided a very limited explanation for the ASY change.We observed a remarkable topography-related vegetation pattern change: in Period III, the "source" patches were found closer to the water bodies than in Period II, whereas, the "sink" patches were found further away from the water bodies.This finding may be an important promoting factor for the increase in the ASY.The LCIs took into account vegetation pattern changes and thus presented a remarkable increase (all more than 44%) that responded sensitively to the increase in the ASY.The prominent coupling relationship between the LCIs/mLCIs and ASY confirmed our success in introducing the LCI to study sediment yield variations in response to vegetation cover changes in the Poyang Lake drainage area, China.The mLCIs proposed in this study (which performed better than the LCIs) were used to identify the critical sediment loss sub-watersheds.
Nevertheless, the limitations of the mLCI should be known prior to its application.As mentioned in Section 2.2.1, the mLCI does not consider the effects of rainfall variation on sediment transport processes.Thus, the mLCI may only be applicable for regions with similar precipitation conditions when studying the relationship between vegetation cover and sediment yield.For example, in this study, we selected three periods with similar rainfall erosivities.However, the erosivity of the Raohe watershed in Period III was relatively higher than of the other watersheds in the different periods.This discrepancy is why the LCI/mLCI and ASY did not correlate well in the Raohe watershed in Period III.Another issue is that sediment yield can sometimes be high in periods following great floods even if the rainfall in the given year is not that large.To settle this problem, we let one period consist of three years and used the three-year average value of the area-specific sediment yield and annual rainfall erosivity to represent the situation during the given period.Nevertheless, we believe that more solutions need to be explored in the future.
Many landscape indices have been established to reveal the relationships between sediment loss processes and landscape patterns.The patch area, patch shape index, Shannon-Weaver diversity index, dominance index, and fragmentation index are commonly applied landscape indices.However, their ability to present landscape patterns is limited by certain angles.For example, the patch shape index measures the complexity of the patch shape.Moreover, these landscape indices could not provide insights into the ecological processes [8,10].Comprehensive landscape indices that integrate multiple factors (e.g., topography, vegetation) from an ecological process perspective also exist, such as the index of connectivity (IC) proposed by Borselli et al. [24] and the LCI introduced in this article.Both indices can effectively link landscape patterns to sediment transport processes.However, they each have their own characteristics.As one important advantage, the IC index can estimate the connectivity between the sediment and water body in a continuous dynamic form (each pixel has own contribution to connectivity).However, the IC has a strict precision requirement for the input spatial data.Using the watershed as the response unit, the LCI can efficiently provide recognition of the entire watershed, but the loss of detailed spatial characteristics is inevitable.A comparison study between the LCI/mLCI and other landscape indices should be performed in the future.

Conclusions
The goal of this study was to reveal the responses of sediment yield to vegetation cover changes in the Poyang Lake drainage area by introducing the LCI based on the "source-sink" theory.The mLCI was proposed in this study to improve the accuracy and efficiency.We found that, all of the river watersheds (Ganjiang, Xiushui, Fuhe, Xinjiang and Raohe) experienced a sharp increase in the average VF and the changes in the vegetation pattern increased the likelihood for sediment entrapment from Period I to Period II.From Period II to Period III, the average VF slightly decreased and at the same time the changes in the vegetation patterns decreased the likelihood for sediment entrapment to some extent.The LCIs associated with elevation, flow length, and Euclidean distance to water had significant positive correlations with the ASY.Compared with the LCIs, the mLCIs showed a higher correlation with the ASY.Using the mLCIs as the judgment indices of sediment loss risk, three sub-watersheds (Qingfengshan, Yuanhe, and Jinjiang) of the Ganjiang watershed were identified as key harnessing areas of sediment loss.Rationalizing the topography-related distribution patterns of vegetation covers was strongly suggested in these three sub-watersheds.We extended the original use of the LCI based on the "source-sink" theory of ecological processes.The methodology applied in this study effectively revealed the responses of sediment yield to vegetation cover from the perspective of the sediment transport process.This methodology should be the most practical approach for studies in large-scale watersheds where sediment surveillance is sparse.

Figure 1 .
Figure 1.Geographic location of the study area.The detail shows the digital elevation model of the study area.

Figure 1 .
Figure 1.Geographic location of the study area.The detail shows the digital elevation model of the study area.

Figure 2 .
Figure 2. Vegetation fraction (VF) distribution of the study area in three time periods.

Figure 2 .
Figure 2. Vegetation fraction (VF) distribution of the study area in three time periods.

Figure 3 .
Figure 3.The "source-sink theory" of sediment transport: (a) a watershed unit; and (b) the sediment transport process on the ground.

Figure 3 .
Figure 3.The "source-sink theory" of sediment transport: (a) a watershed unit; and (b) the sediment transport process on the ground.

Figure 4 .
Figure 4. Spatial distribution of the Lorenz curve of "source" and "sink" landscapes.

Figure 4 .
Figure 4. Spatial distribution of the Lorenz curve of "source" and "sink" landscapes.

Figure 5 .
Figure 5. Relationships among Euclidean distance to water, elevation, slope and flow length.The abscissa represents buffer zones of the Euclidean distance to water, named C1, C2 …… C32; the ordinate represents the mean of the normalized elevation, slope and flow length values in each buffer zone.

Figure 5 .
Figure 5. Relationships among Euclidean distance to water, elevation, slope and flow length.The abscissa represents buffer zones of the Euclidean distance to water, named C1, C2 . . .C32; the ordinate represents the mean of the normalized elevation, slope and flow length values in each buffer zone.

Figure 6 .
Figure 6.Association of the spatial distribution of the five VF levels with the slope and Euclidean distance in the Xinjiang watershed.Figure 6. Association of the spatial distribution of the five VF levels with the slope and Euclidean distance in the Xinjiang watershed.

Figure 6 .
Figure 6.Association of the spatial distribution of the five VF levels with the slope and Euclidean distance in the Xinjiang watershed.Figure 6. Association of the spatial distribution of the five VF levels with the slope and Euclidean distance in the Xinjiang watershed.

Figure 7 .
Figure 7.The scatter map of the area-specific sediment yield (ASY) and the average vegetation fraction, the location-weighted landscape contrast index (LCI) and the modified location-weighted landscape contrast index (mLCI).

Figure 8 .
Figure 8. Spatial distribution of the Med.

Figure 7 .
Figure 7.The scatter map of the area-specific sediment yield (ASY) and the average vegetation fraction, the location-weighted landscape contrast index (LCI) and the modified location-weighted landscape contrast index (mLCI).

Water 2016, 8 , 114 11 of 16 Figure 7 .
Figure 7.The scatter map of the area-specific sediment yield (ASY) and the average vegetation fraction, the location-weighted landscape contrast index (LCI) and the modified location-weighted landscape contrast index (mLCI).

Figure 8 .
Figure 8. Spatial distribution of the Med.

Figure 8 .
Figure 8. Spatial distribution of the Med.

Table 1 .
Mean annual area-specific sediment yield (ASY, t/(km 2 ¨a); a indicates the year) and the mean annual rainfall erosivity (MAE, (MJ¨mm)/(km 2 ¨h¨a)) in each watershed.

Table 2 .
Areal percentages of the five vegetation fraction (VF, %) levels in each watershed.