Evaluating the Scale Effect of Soil Erosion Using Landscape Pattern Metrics and Information Entropy : A Case Study in the Danjiangkou Reservoir Area , China

The regular patterns of soil erosion tend to change at different scales of observation, affecting the mechanism of soil erosion and its evolution characteristics. This phenomenon has essential scientific significance for the rational allocation of land resources and for studies on sustainable ecosystems. As an important agricultural area in China, Danjiangkou reservoir is threatened by severe soil erosion. In this study, we selected four kinds of landscape pattern metrics, including patch density, fractal dimension, Shannon diversity index, and connectivity, to analyze soil erosion intensity in the Danjiangkou reservoir area at different scales based on landscape ecological principles. In addition, we determine the optimum research scale of the experimental area by calculating the information entropy value of soil patches at different scales. The findings suggest that: (1) the landscape pattern of soil erosion in the experimental area is obviously scale-dependent, and the responses to scale differ from index to index; (2) as the scale of observation increases, the fragmentation of soil patches is weakened, the stability of different landscape components is enhanced, and the soil becomes less vulnerable to erosion; and (3) based on information entropy theory, 60 m is confirmed to be the optimum scale of this study.


Introduction
As a core subject of soil and water conservation research in sustainable ecosystems, soil erosion is affected by various natural and human factors [1,2].Climate, lithology, soil types, and especially intense land use/cover changes make soil vulnerable to runoff generation and water erosion, leading to the loss of nutrients from the soil, reservoir siltation, and devastating floods [3].Climatic parameters such as soil moisture, wind, and precipitation are all prone to influence soil erosion.Soil moisture is especially one of the most important factors controlling erosive processes [4], affecting the partitioning of rainfall into infiltration and runoff and consequently influencing soil erosion [5].Another parameter influencing soil erosion is soil type, as the impact and sensitivity of different soil types to the erosion characteristics arising from land use change are quite different.Purple soil, for instance, was determined to be more sensitive to rainfall intensity in the land-use change from farmland to forest than red soil, black soil, and loess [6].Moreover, some human impacts have altered land use structures and triggered soil erosion [7].These land use changes, along with the climate, steep slope, and the fragility of the soil, make the study region very prone to erosion.How to prevent the soil erosion and recover the ecosystems through reasonable land use has become a challenging issue.
Furthermore, the results of studies on soil erosion could vary due to the changes in the scales of observation [8].Scale is a kind of space or time unit used in studies of certain phenomena, which can be generally divided into a time-scale and space-scale [9,10].Space-scale refers to the resolution of the smallest recognizable unit; that is, the granularity or the cell size, taking a remote sensing (RS) image as an example.The scale of RS images is an essential factor that must be taken into account when analyzing geospatial issues produced from the RS technique [11].Methods that deal with spatial autocorrelations are common in non-geostatistical studies.For example, Moran's 'I' and Geary's 'C' statistics are used in spatial correlations for the analysis of raster images [12,13].Meanwhile, methods based on geo-statistical measures of spatial variability are also widely developed.The most effective method is based on computing a set of textural measures of spatial variability (TMSV) using a moving window [14].However, the choice of an appropriate moving window size is an important and challenging issue when researching scale-related problems in remote sensing, especially in the scale conversion of raster images.The scale analyses in remote sensing have been achieved by many scholars.For instance, Zhang et al. studied the scale effect of remote sensing images using histo-variograms and conducted scale conversion with an aggregation moving window [15].Lausch et al. tested and evaluated the suitability of hyper-spectral remote-sensing techniques in detecting the predisposition of spruce vegetation to beetle damage on different spatial scales [16].The spatial scale effect of mapping plastic-mulched farmland (PMF) was elucidated by Hasituya et al. using high-resolution imagery [17].
Process-based models of ecological function are usually parameterized using measurements at a particular plot, which simplifies only a small fraction of the landscape [18].This limited sampling means that it is difficult to parameterize ecosystem models across space, which hinders our ability to accurately estimate multi-scale landscapes.Landscape pattern is, to some extent, dependent on scale [19], and changes in scale also have certain impacts on the analyses of landscape patterns [20].Therefore, landscape pattern research emphasizes the spatial heterogeneity of landscape features [21] and the determination of scales to quantify such variability [22].Morris et al. quantified the particulate organic matter (POM) cover of a certain area and studied the landscape-scale spatial variation in POM supply to the abyssal seafloor [23].Furthermore, landscape metrics can be used as indicators of the heterogeneity of different landscape components, which have been studied by many researchers.For example, Zhang et al. assessed ecological security characteristics in the Yuan River watershed based on the relationship between landscape patterns and soil erosion [24].Morelli et al. explored several models using the metrics of landscape heterogeneity and altitude as predictors of bird diversity in different spatial scales [25].In another study, the spatial heterogeneity among Greek landscapes was investigated by Plexida and others, where they identified appropriate landscape indicators using different extents [26].However, among all the above-mentioned studies, few have investigated the importance of the spatial-scale heterogeneity and scale effects on soil erosion in a reservoir region, despite some efforts to research forest and hydrological scale effects.Moreover, appropriate approaches to determine the optimum scale for soil erosion studies are also new area to explore.Here, from the perspective of landscape ecology, we propose a study to investigate 'landscape pattern-soil erosion' processes by using landscape metrics and explore the scale effect of soil erosion in an important Chinese agricultural reservoir.
Similarly, few if any examples of scaling in landscape ecology consider the inherent loss or change in information among aggregation process across spatial scales, despite the widespread use of information-theoretic approaches in modeling [27,28].Changes in information depend on how information is quantified, while techniques from information theory can be applied to quantify how information is gained or lost upon increasing or decreasing the spatial scale of observations [29].In addition to the selection of entropy for quantifying the amount of information, it can be also chosen for studies on time-scale effects combined with wavelet analysis [30].Here, we investigate approaches from information theory to determine the optimum research spatial-scale for process-based soil erosion studies.Furthermore, under this optimum cell size, we intend to explore the time-scale effect of soil erosion data combined with wavelet analysis in our future study.
Our research design is based on the Danjiangkou reservoir, China.As an important area in terms of agricultural and economic production in Hubei Province of China, the Danjiangkou reservoir has been threatened by soil erosion.Therefore, this reservoir is taken as a typical area to analyze the relationship between soil erosion and landscape patterns under different spatial scales, as well as to explore the inherent law of changes in soil erosion with spatial scales.In this paper, we aim to detect the scale effect of different landscape pattern metrics for 'landscape pattern-soil erosion' studies by using geographic information system (GIS) and remote sensing (RS) technologies.In addition, we attempt to propose a new idea to determine the optimum spatial scale of soil erosion in the reservoir area based on approaches from information entropy theory.The remainder of this paper is organized as follows.After the introduction, Section 2 describes the methodology and data used in this study.Section 3 presents and discusses the main results of our analysis.Section 4 provides a conclusion of the paper and points out limitations for future study.process-based soil erosion studies.Furthermore, under this optimum cell size, we intend to explore the time-scale effect of soil erosion data combined with wavelet analysis in our future study.Our research design is based on the Danjiangkou reservoir, China.As an important area in terms of agricultural and economic production in Hubei Province of China, the Danjiangkou reservoir has been threatened by soil erosion.Therefore, this reservoir is taken as a typical area to analyze the relationship between soil erosion and landscape patterns under different spatial scales, as well as to explore the inherent law of changes in soil erosion with spatial scales.In this paper, we aim to detect the scale effect of different landscape pattern metrics for 'landscape pattern-soil erosion' studies by using geographic information system (GIS) and remote sensing (RS) technologies.In addition, we attempt to propose a new idea to determine the optimum spatial scale of soil erosion in the reservoir area based on approaches from information entropy theory.The remainder of this paper is organized as follows.After the introduction, Section 2 describes the methodology and data used in this study.Section 3 presents and discusses the main results of our analysis.Section 4 provides a conclusion of the paper and points out limitations for future study.

Study Area and Data
Danjiangkou Reservoir is located in the Yunxian County of Hubei Province, China.The reservoir is bordered by the Danjiangkou city in the east, Xiling of Yunxian County in the west, Wudangtown in the south, and XichuanCounty of Henan Province in the north.The region is situated between 32°36′-33°48′ N and 110°59′-111°49′ E. Figure 1a,b show the relative location and the enlarged image of the area, respectively.The reservoir's storage capacity reaches 19 billion m 3 , which makes it the second-largest artificial reservoir in China.It lies in the transitional zone between the northern subtropics and the warm temperate zone and is affected by the typical monsoon continental semi-humid climate.The region is a relatively independent natural unit, the terrain of which is high in the northwest and low in the southeast.With complex terrain and diverse natural conditions, the reservoir is suitable for the production of vegetation.However, due to poor natural conditions and the impact of human reclamation, the area has experienced different degrees of soil erosion.
The original data used in this study include; Landsat-5 TM images with a resolution of 30 m for the years 2002 and 2009, Digital Elevation Model(DEM) data (with a resolution of 30 m), and other fundamental materials related to soil erosion such as 'Soil records of Danjiangkou City'.

Data Processing
In our study, data are processed in GIS and the remote sensing platforms ArcGIS and ENVI.The first data processing step is to transform projection, wherein the administrative data and remote sensing images are projected under a unified projection.Then we use ENVI 5.0 to complete the The reservoir's storage capacity reaches 19 billion m 3 , which makes it the second-largest artificial reservoir in China.It lies in the transitional zone between the northern subtropics and the warm temperate zone and is affected by the typical monsoon continental semi-humid climate.The region is a relatively independent natural unit, the terrain of which is high in the northwest and low in the southeast.With complex terrain and diverse natural conditions, the reservoir is suitable for the production of vegetation.However, due to poor natural conditions and the impact of human reclamation, the area has experienced different degrees of soil erosion.
The original data used in this study include; Landsat-5 TM images with a resolution of 30 m for the years 2002 and 2009, Digital Elevation Model(DEM) data (with a resolution of 30 m), and other fundamental materials related to soil erosion such as 'Soil records of Danjiangkou City'.

Data Processing
In our study, data are processed in GIS and the remote sensing platforms ArcGIS and ENVI.The first data processing step is to transform projection, wherein the administrative data and remote sensing images are projected under a unified projection.Then we use ENVI 5.0 to complete the preprocessing of the remote sensing images of two-year periods, including splicing, cropping, and merging.Through the method of supervised classification, we categorize land use types in the study region into five categories; water, forestland, arable land, construction land, and unused land.Then the land use data is converted into a vector format.The classification of the land-use system is presented in Table 1.
Table 1.The land use classification system.

First Level Classes Second Level Classes Descriptions
Arable land -Cultivated lands for crops, including mature cultivated land, newly-cultivated land, reclamation land, shifting cultivated land; crop land in which a crop is a dominant species and land that is used for intercropping such as crop-fruiter; other lands that have been shifted to cultivation temporarily.

Paddy field
Arable land that has enough water supply and irrigation facilities for the cultivation of paddy rice and other aquatic crops.

Dry land
Arable land with no water supply and irrigation facilities; arable land that is used for planting dry farming crops with natural precipitation.
Forest land -Lands growing trees, including arbor, shrub, bamboo, and for forestry use.

Forest
Natural or planted forests with tree canopy cover >30%.

Shrub
Lands covered by trees under 2 m high with canopy cover ≥40%.

Others
Lands, including sparse woodland and nurseries; lands with the growth of herbs-based land.
Construction land -Lands used for settlements, business use, factories, transportation infrastructures, and so on.

Residential land
Lands used for urban and rural settlements and their ancillary facilities.

Business land
Lands used for business use, including wholesale and retail land, accommodation land, financial land, and so on.

Others
Lands used for industrial and mining use, public administration, and public service and lands used for transportation facilities such as airports.
Water -Lands covered by natural water bodies or lands with facilities for water reservation.

Reservoir
Man-made facilities with a total storage capacity more than 10 million cubic meters.

Lakes and rivers
Lands covered by lakes or rivers, including canals.
Unused land -Lands are not utilized or hard to make use of.

Sandy land
Lands with a surface of sand cover, with no more than 5% vegetation cover.

Bare land
Bare land with exposed soil or rocks, basically with no vegetation cover.

Saline land
Lands with saline accumulation and sparse vegetation.

Undeveloped land
New lands in towns or villages that have not been put into use yet.
Vegetation coverage data is extracted from remote sensing images, and slope data is acquired from the DEM.The first step is to calculate the NDVI (normalized differential vegetation index), which is used for constructing vegetation coverage grid data in ArcGIS from the images; the second step is to generate slope grid data, which is derived from the DEM using the 'slope' function in the software.Finally, these grid data are converted into vectors to be consistent with the land use data.In this paper, we adopt a method to study soil erosion using a remote sensing monitoring method based on the indicator rules proposed by Zhou in 2005 [31], then overlay the vectors of land use, vegetation coverage, and slope data in order to obtain the soil erosion map.
According to the industrial standard of the People's Republic of China for the classification and gradation of soil erosion (SL 190-2007), as well as the different combination of topography, geomorphology, land use, and vegetation, the Beijing Municipal Water Resources Bureau formulates a grading system of soil erosion (Table 2).This system divides the erosion intensity into six degrees: (1) micro (no erosion or erosion is not obvious, including forest land with vegetation coverage >75% and lands with slope <5 • but not including arable land with a slope of 3~5 • ); (2) mild (soil erosion is more obvious, including cultivated lands with a slope of 3~8 • , land with vegetation coverage of 60~70% and a slope of 5~25 • , or land with vegetation coverage of 45~60% and a slope of 5~15 • , etc.); (3) moderate (soil erosion is pretty obvious, including cultivated lands with a slope of 8~15 • , land with vegetation coverage <30% and a slope of 5~15 • , etc.); ( 4) deep (strong erosion, mainly distributed in mountains with little vegetation, including land with vegetation coverage of 45~60% and a slope >35 • , etc.); ( 5) intensive (intensive erosion, mainly distributed in high hilly slope); (6)and severe (severe erosion, including bare land with sparse vegetation or arable land with a slope higher than 35 • ).(3) To test various scale effects of soil erosion, spatial aggregation and scale conversion are processed in ArcGIS to derive a series of data with different spatial resolutions ranging from 30 m to 240 m.We resample the original map (with a resolution of 30 m) and aggregate cells into several larger scale datasets with 30 m increments using an arithmetic average variability weighted (AAVW) algorithm.Then, based on the fractal theory and landscape pattern metrics, we treat the soil erosion unit as a patch from the view of landscape ecology, calculating and comparing the landscape pattern metrics at different spatial scales in order to explore and reveal the inherent law of soil erosion change with scale.

Calculation of Landscape Pattern Metrics
Landscape pattern metrics can be used to reflect and evaluate the spatial characteristics and the distribution of the landscape [32] such as its fragmentation, complexity, or diversity.In this study, we selected four representative pattern metrics, including Patch Density (PD), Perimeter-Area Fractal Dimension (PAFRAC), Shannon's Diversity Index (SHDI), and Connectivity (CONNECT).
PD is the number of patches per unit area.The higher the PD, the higher the degree of landscape fragmentation.It reflects the uniform distribution of landscape patches and is conducive to the comparison between landscapes of different sizes.
It can be calculated as follows: where A is the area of a certain landscape type (m 2 ), n is the number of patches, and '10000 × 100' represents the number of patches per 100 hectares.
PAFRAC is an index to reveal the relationship between the patch shape and its area, reflecting the complexity of patches and landscape patterns at a given scale of observation [33,34].Its values are between one and two, with stronger self-similarity in patches when converging towards one, and the shape of patches is also simpler.The expression is shown as follows: where n ij is the number of patches typed as i, a ij is the jth patch area of patch type i (m 2 ), and p ij is the jth patch perimeter of patch type i (m).SHDI reflects the different landscape patterns or the diversity index of a certain kind of landscape pattern at different times.In some landscapes, the richer the land-use, the higher the degree of fragmentation and the higher the value of SHDI.Its calculation formula is: where m is the number of landscape types and P i is the proportion that type i accounts for in all types.
CONNECT is the connection degree of patch-to-patch distance [35,36].In soil erosion studies, the spatial connectivity of landscape units is very important.For example, the increased connectivity of vegetation patches is conducive to the inhibition of soil erosion [37].The connect index can be expressed as follows: where c ijk represents the connection situation between the same type of patches j and k within a given connection threshold (1 for connection and 0 for no connection), n i is the number of patches for type i, and m is the number of patch types.

Grading Evaluation of Soil Erosion
Different land-use types, vegetation coverage, topographic gradients, and their combinations correspond to different soil erosion intensity grades.By using the spatial analysis function of ArcGIS, we first overlay the land-use classification map, the vegetation coverage gradation map, and the slope gradation map; then we count the number of cells of soil erosion under each intensity and aggregate the cells with the same erosion degree to obtain the classification tables of soil erosion intensity (Tables 3  and 4) and soil erosion intensity grading diagrams (shown in Figures 2 and 3) in Danjiangkou reservoir in the years 2002 and 2009, respectively.From the charts, we can conclude that the main degrees of soil erosion in the Danjiangkou reservoir area between 2002 and 2009 are micro-erosion and mild-erosion.The total area of erosion decreases to a certain extent in 2009 compared with that in 2002.In the year 2002, micro and mild degrees of soil erosion account for 87.01% of the total erosion area, while moderate erosion and above accounts for 12.99%.However, in 2009, the proportion of micro-and mild-erosion is 90.01%and the percentage of moderate erosion and above is only 9.99%.The results show that, during this period, the intensity of soil erosion in the Danjiangkou reservoir is alleviating, while those areas of more intensive erosion gradually transformed into slighter ones.The total area of micro-erosion increased significantly and the overall level of soil erosion experience a declining trend.
The change in soil erosion can be attributed to natural and human factors.In addition to changes in land use structure, other aspects such as the 'grain for green project' and local residents' promoted awareness of environmental protection, have all contributed to the alleviation of soil erosion.Moreover, the change in agricultural tillage method, wherein the high-load and intensive farming style has been transformed to new eco-friendly one, is also beneficial to lighten the erosion.

Scale Effects of Soil Erosion Patches
Scale effects of soil erosion patches can be explored through analyzing the differences and From the charts, we can conclude that the main degrees of soil erosion in the Danjiangkou reservoir area between 2002 and 2009 are micro-erosion and mild-erosion.The total area of erosion decreases to a certain extent in 2009 compared with that in 2002.In the year 2002, micro and mild degrees of soil erosion account for 87.01% of the total erosion area, while moderate erosion and above accounts for 12.99%.However, in 2009, the proportion of micro-and mild-erosion is 90.01%and the percentage of moderate erosion and above is only 9.99%.The results show that, during this period, the intensity of soil erosion in the Danjiangkou reservoir is alleviating, while those areas of more intensive erosion gradually transformed into slighter ones.The total area of micro-erosion increased significantly and the overall level of soil erosion experience a declining trend.
The change in soil erosion can be attributed to natural and human factors.In addition to changes in land use structure, other aspects such as the 'grain for green project' and local residents' promoted awareness of environmental protection, have all contributed to the alleviation of soil erosion.Moreover, the change in agricultural tillage method, wherein the high-load and intensive farming style has been transformed to new eco-friendly one, is also beneficial to lighten the erosion.

Scale Effects of Soil Erosion Patches
Scale effects of soil erosion patches can be explored through analyzing the differences and changing regulations of landscape pattern metrics under different spatial scales [38,39].Based on the original eroded raster images (with a resolution of 30 m), we generated the soil erosion images at different scales using GIS technology according to an incremental size of 30 m.The maximum raster scale is set as 240 m.Then we calculate and compare the scale response characteristics of each landscape index in two different years.The results are listed in Tables 5 and 6.It can be noted from the tables that the index values of soil erosion patches have certain responses to the change of scale.The metrics of PD, PAFRAC, and CONNECT are more sensitive to the changes of scale, while the SHDI index has no fluctuation at all. Figure 4 shows the responses of four different metrics to the varying scales.
In terms of these metrics, PD has the largest variation with scale (Figure 4a).With the increase of the scale of observation, the patch density (PD) of the soil erosion unit shows a decreasing trend.The number of patches per unit area of soil is a characterized parameter of landscape fragmentation [40], and the degree of landscape fragmentation is related to natural geographical factors as well as the interference of human economic activities [41].The higher the soil ecological vulnerability grade is, the less vulnerable soil becomes to erosion and therefore the less intensive the fragmentation of patches is likely to be.
Figure 4b,c show the responses of the PAFRAC and CONNECT metrics due to the increase of scale, respectively, indicating that both of these metrics respond significantly to spatial scales.It is widely acknowledged that the area-perimeter fractal dimension (PAFRAC) is used to describe the complexity of the shape of the patch boundary [31,42].Therefore, the shape of patches is closely related to their boundary lines, and the shape is supposed to become more complex if the boundary line does.In this study, the growing trend of the PAFRAC index shows that the structure of patches has a moderate scale effect and is projected to suffer higher complexity with the increase of scale.
As for the CONNECT index, the degree of connectivity in landscapes can reflect the diversity, stability, and integrity of different landscape patterns [43].Therefore, stronger connectivity is also beneficial to the suppression of soil erosion and the maintenance of an ecosystem's integrity [44].Figure 4c indicates that the CONNECT index represents scale-dependency to some extent and peaks at the scale of 150 m.This Sustainability 2017, 9, 1243 9 of 15 is due to the fact that the stability of different landscape components is enhanced when the scale of the observation is smaller than 150 m, after which the connection of erosion patches is weakened.
Last but not least, the diversity index (SHDI) does not show such significant fluctuations (shown in Figure 4d) as the previous three indices do, which illustrates that the scale effect of the SHDI index is relatively weaker compared to the other metrics, as it responds faintly according to the increase of scale.Therefore, within the selected range of scales, no patch types disappear.
As indicated above, based on the 'scale-pattern-process' principle of landscape ecology, we can conclude that, with the increase of scale, the scale effect of PD shows that the fragmentation of soil patches is weakened, and the less vulnerable soil succumbs to erosion.The scale response of PAFRAC indicates that the structure of patches tends to suffer higher complexity; the change of CONNECT suggests that the stability of different landscape components is enhanced when the scale is below 150 m, after which the connection of erosion patches is weakened; and the weak scale effect of the SHDI index illustrates that, within the selected range of scales, nearly no disappearance of patch types occurs.

Optimal Scale Selection of Soil Erosion Based on Information Entropy
At present, many researchers have explored how to determine the best scale to conduct soil erosion research.For example, Lv et al. [45] adopt a method of drawing diagrams of landscape metrics to assign the scale threshold and identify the turning points of landscape metrics to determine the optimal scale.However, this method is significantly affected by subjective factors, which undermine the precision of the results.Other scholars have studied the spatial transfer of soil

Optimal Scale Selection of Soil Erosion Based on Information Entropy
At present, many researchers have explored how to determine the best scale to conduct soil erosion research.For example, Lv et al. [45] adopt a method of drawing diagrams of landscape metrics to assign the scale threshold and identify the turning points of landscape metrics to determine the optimal scale.However, this method is significantly affected by subjective factors, which undermine the precision of the results.Other scholars have studied the spatial transfer of soil erosion parameters in agricultural loess areas and found that erosion parameters optimized at smaller scales could be transferred to higher scales without the need for calibration [46].In this paper, based on the theory of information entropy, we propose a new method to identify the optimal scale of soil erosion research by calculating the information entropy values of erosion patches at different scales, which provides a new way to study scale effects.The basic theoretical framework of this method is as follows.
Land use structure can be analyzed at different scales using entropy [47].The process of a soil erosion system embodies the proceeding of time and space [48].In addition to being subject to more or less external constraints, any physical system always has certain internal degrees of freedom, therefore causing the elements within the system to be in different states.If the entropy value is larger on the scale, then the sequence is less self-similar on the scale and the structure is less simple; on the contrary, it is more complex.The maximum entropy occurs when the abundance of the state of things automatically reaches the maximum.In other words, things always strive for (or show) the greatest freedom under constraints, which is the fundamental principle of nature.In most cases, the change of scale will affect the accuracy of the classification of remote sensing images, thereby affecting the amount of information contained in soil erosion patches and the complexity of the erosion system itself.Thus, the assumption is that there must be an optimal scale under which erosion patches contain the most abundant information.In this regard, such scale can be considered the most appropriate.
If the proportion of each erosion category is equal, the information entropy value can reach its maximum value of 1.At this time, the soil erosion system is supposed to be the most complicated, and the information (i.e., diversity) is also the richest in each grid of the patch.However, some types of erosion tend to change to other types over time (for example, areas of moderate erosion and more intensive erosion have transformed into areas of micro and mild erosion over the years in our study area).Then, the order degree of the system will increase and the information entropy value will decrease.The information entropy value will drop to 0 when there is only one type left.Therefore, there must be an extreme information entropy value (closest to 1) when the value of soil erosion changes with scale.The optimal scale selection based on the theory indicated above is calculated as follows.
In information theory, entropy is a measure of the uncertainty of a system.The greater the entropy, the greater the degree of disorder a system contains.Its definition can be expressed as: where SH is the information entropy value, n is the number of all events, i represents a random event that occurs in the system, and P i is the occurrence probability of i. Obviously, ∑P i =1 satisfies the normalization condition.For the experimental area, define: , where I min = 30m, I max = 240m.A i is the eroded area (m 2 ) of different types of soil erosion (classified in Table 1), i ∈ 1, 2, • • • , 6, and the total area of erosion is defined as A.
P i represents the percentage that each erosion type accounts for in the total area, which can be calculated as follows: When A 1 = A 2 = . . .= A 6 , SH max = 1.
When there is only one A i and A i = 0, other erosion types have no area at all; then SH min = 0. Thus, there must be I m ; when I k ≤ I m , SH increases, and when I k ≥ I m , SH decreases.The scale corresponding to the extreme point of the SH value (closest to 1) is supposed to be the optimal scale of soil erosion in this study area.
In this paper, we take the data from the year 2009 as an example.The soil erosion images of Danjiangkou reservoir with different resolutions are graded at each scale (using the method described in Section 3.1.1).Then, the values of A i and P i and the information entropy value are calculated at each I k .Table 7 shows the results (where P 1 ~P6 is related to the proportion of the areas of micro, mild, moderate, deep, intensive, and severe erosion to the total erosion area).Based on the results in Table 7, we can draw the response of the information entropy value of soil erosion to the change of scale in our study area (Figure 5).The results show that the local maximum at a scale of 60 m is highlighted in the charts, and there is an apparent decline in SH beyond the scale of 120 m, which indicates an information loss.Thus, we can reasonably assume that the optimal scale threshold in this study is 60-120 m and that the best scale is located in this threshold interval.In the threshold, the value of SH decreases progressively from the scale of 60 m, which would, therefore, highlight continuous information loss with the increase of scale.Thus, we can conclude that 60 m is the optimal scale in our study area.
As a reflection of the state of open systems [49] such as soil erosion systems, the change of entropy can be used not only to explore the optimal scale, but also to indicate the progress of the erosion system.The evolution of soil erosion follows thermodynamic laws of open systems, wherein a local decrease of entropy is possible and free energy is used to build complex structures while global entropy increases.The fluctuations in open systems decrease the entropy, whereas the irreversible processes restore the system to its initial state.Therefore, the entropy in an isolated system will not reduce continuously (according to the Second Law of Thermodynamics) [50].A part of the entropy decrease is always accompanied by another part of the entropy increase so as to maintain the stability of open systems.As an open system, soil erosion system follows these dynamic laws with no exception.In the scale response analysis of the information entropy value of soil erosion in the Danjiangkou reservoir area, the entropy increases between the scale of 30 and 60 The results show that the local maximum at a scale of 60 m is highlighted in the charts, and there is an apparent decline in SH beyond the scale of 120 m, which indicates an information loss.Thus, we can reasonably assume that the optimal scale threshold in this study is 60-120 m and that the best scale is located in this threshold interval.In the threshold, the value of SH decreases progressively from the scale of 60 m, which would, therefore, highlight continuous information loss with the increase of scale.Thus, we can conclude that 60 m is the optimal scale in our study area.
As a reflection of the state of open systems [49] such as soil erosion systems, the change of entropy can be used not only to explore the optimal scale, but also to indicate the progress of the erosion system.The evolution of soil erosion follows thermodynamic laws of open systems, wherein a local decrease of entropy is possible and free energy is used to build complex structures while global entropy increases.
The fluctuations in open systems decrease the entropy, whereas the irreversible processes restore the system to its initial state.Therefore, the entropy in an isolated system will not reduce continuously (according to the Second Law of Thermodynamics) [50].A part of the entropy decrease is always accompanied by another part of the entropy increase so as to maintain the stability of open systems.As an open system, soil erosion system follows these dynamic laws with no exception.In the scale response analysis of the information entropy value of soil erosion in the Danjiangkou reservoir area, the entropy increases between the scale of 30 and 60 m and reaches its peak at 60 m, before decreasing dramatically with the further increase of scale.

Conclusions
In this paper, we use RS and GIS technologies to explore the scale effect of soil erosion in the Danjiangkou reservoir area between 2002 and 2009.We also use landscape pattern metrics to reveal the scale effect and the evolution law of soil erosion at different scales.In addition, we analyze and identify the optimal scale for the soil erosion of this agricultural area based on information entropy theory.The conclusions are as follows: (1) The main degrees of soil erosion in the Danjiangkou reservoir area are micro-and mild-erosion.
Compared to the level in 2002, the overall degree of erosion has decreased during this period.
The micro-and mild-erosion are concentrated in the arable land with a shallow slope and construction land, while the more severe erosion is clustered in the forestland with little vegetation coverage and a steep slope.(2) At the holistic level of landscape patterns, the density of soil erosion patches (PD) presents a declining trend with the increase of scale, which demonstrates that the fragmentation of soil patches is weakened and that the less vulnerable soil succumbs to erosion; the area-perimeter fractal dimension (PAFRAC) and the connectivity (CONNECT) fluctuate considerably with scale, showing that the structure of patches tends to suffer higher complexity and that the stability of different landscape components is enhanced when the scale is below 150 m.The response of the diversity index (SHDI) to scale is relatively weak, indicating that nearly no disappearance of patch types occurs within the selected range of scale.(3) Information entropy theory is used to reveal the information richness of soil erosion patches.
Using this correlation, we identify that the optimal scale threshold in this study is 60-120 m, the value of SH decreases progressively from the scale of 60 m, and there is an obvious decline in SH beyond the scale of 120 m, highlighting an information loss.Thus, the best scale for this reservoir area is supposed to be 60 m.
There are two areas we can improve in the study presented in this paper.First, the landscape pattern metrics we chose in this paper have some limitations if more land use types are taken into consideration.Therefore, it is necessary to include more representative metrics into future research.Secondly, in this paper we investigated approaches from information theory to determine the optimum research spatial-scale for process-based soil erosion studies.Nevertheless, we are also fully aware that time-scale effects of soil erosion can be considered if combined with other mathematical models.Therefore, on top of the optimum spatial cell size method we propose in this paper, we will explore the time-scale effect of soil erosion data using entropy and wavelet analysis in our future work.

Danjiangkou
Reservoir is located in the Yunxian County of Hubei Province, China.The reservoir is bordered by the Danjiangkou city in the east, Xiling of Yunxian County in the west, Wudangtown in the south, and XichuanCounty of Henan Province in the north.The region is situated between 32 • 36 -33 • 48 N and 110 • 59 -111 • 49 E. Figure 1a,b show the relative location and the enlarged image of the area, respectively.Sustainability 2017, 9, 1243 3 of 15

Figure 1 .
Figure 1.Location of Danjiangkou Reservoir in China: (a) the relative location; and (b) the enlarged remote sensing image of the reservoir.

Figure 1 .
Figure 1.Location of Danjiangkou Reservoir in China: (a) the relative location; and (b) the enlarged remote sensing image of the reservoir.

Figure 3 .
Figure 3. Intensity grading diagram of soil erosion in the reservoir area in 2009.

Sustainability 2017, 9 , 1243 10 of 15 Figure 4 .
Figure 4. Scale effects of four different metrics at the varying spatial scales.(a) The figure of the patch density index; (b) figure of the perimeter-area fractal dimension index; (c) figure of the connectivity index; and (d) figure of the Shannon's diversity index.

Figure 4 .
Figure 4. Scale effects of four different metrics at the varying spatial scales.(a) The figure of the patch density index; (b) figure of the perimeter-area fractal dimension index; (c) figure of the connectivity index; and (d) figure of the Shannon's diversity index.

Figure 5 .
Figure 5. Changes of soil erosion information entropy values at different scales.

Figure 5 .
Figure 5. Changes of soil erosion information entropy values at different scales.

Table 2 .
Grading index of soil erosion intensity.
Note: Intensity classification is judged by land type and erosion factors.For forestland, this is done by using vegetation coverage factor and slope factor; for arable land, by using slope factor directly to determine the erosion intensity; and for other types (water, construction land, and unused land identified as micro-erosion.

Table 3 .
Classification of soil erosion intensity in 2002.

Table 4 .
Classification of soil erosion intensity in 2009.
Figure 2. Intensity grading diagram of soil erosion in the reservoir area in 2002.

Table 4 .
Classification of soil erosion intensity in 2009.
Grade Number Soil Erosion Intensity Number of Cells (30 m × 30 m) Percentage

Table 5 .
Landscape pattern index values at different scales (in 2002).

Table 6 .
Landscape pattern index values at different scales (in 2009).

Table 7 .
Information entropy values of soil erosion under different scales.

Table 7 .
Information entropy values of soil erosion under different scales.