Evolution and Prediction of Landscape Pattern and Habitat Quality Based on CA-Markov and InVEST Model in Hubei Section of Three Gorges Reservoir Area (TGRA)

: The spatial pattern of landscape has great inﬂuence on the biodiversity provided by ecosystem. Understanding the impact of landscape pattern dynamics on habitat quality is signiﬁcant in regional biodiversity conservation, ensuring ecological security guarantee, and maintaining the ecological environmental sustainability. Here, combining CA-Markov and InVEST model, we investigated the evolution of landscape pattern and habitat quality


Introduction
As an important part and one of the major causes of global change, landscape pattern has great influences on ecosystems and the goods and services they provide to people ("ecosystem services") [1], the change of which will lead to variations in lots of natural phenomena and ecological processes [2,3].The term "ecosystem services" was used quite broadly to include any ecosystem process or function that contributes to human [4].Biodiversity is intimately linked to the production of ecosystem services.As a proxy for biodiversity, habitat quality refers to the ability of the environment to provide conditions appropriate for individual and population persistence, ranging from low to medium to high, resources available for survival, reproduction, and population persistence, respectively [5].Generally, the spatial landscape pattern can reflect the regional quality of ecological environment, and they affect each other, interact with each other, and are closely related [6].Changes in landscape will cause corresponding changes in the ecosystem services [7][8][9].By studying the spatial and temporal changes of landscape pattern, we can obtain the corresponding changes in biodiversity.For example, landscape transition from woodland to construction land has a great impact on regional biodiversity, which would cut off the spatial connectivity [10,11].The Three Gorges Reservoir Area (TGRA) is considered to be one of the most important biodiversity hotspots in China [12].The construction of the Three Gorges Dam has huge effects on both the environment and society.Evidence from many sources builds an overwhelming picture of pervasive biodiversity decline in this region after the completion of the Three Gorges Project [13][14][15].This serious consequence for biodiversity has prompted a wide-ranging response from both governments and civil society.Therefore, analyzing the spatial and temporal changes of landscape pattern in past, and exploring the impact of landscape pattern dynamics on biodiversity are beneficial to formulate land use management strategy scientifically and effectively, which can guide the sustainable development of regional society and economy, maintain regional biodiversity, and ensure regional ecological security [16].
Changes in landscape pattern can be analyzed and simulated by various models, such as GIS-optimization modeling [17], random prediction model [18], empirical regression model [19], Cellular Automata (CA) [20,21], and so on.CA is one of the most representative geo-models, which is widely used in many fields such as land use, geomorphologic evolution, and urban expansion [22].The CA-Markov model synthesizes the ability of CA model in simulating spatial changes in complex systems and the advantage of Markov model in long-term forecasting, which improves the prediction precision of the land use types transformation, and effectively simulates the spatial variation of land use structure [23].The CA-Markov model, as a scientific and practical geo-model, is broadly applied to modelling and realistic prediction of land-use patterns.Based on cellular automata, Batty et al. [24] proposed a Dynamic Urban Evolution Model (DUEM), and systematically and effectively simulated urban expansion.Nourqolipour et al. [25] simulated spatial patterns of oil palm expansion by spatial and temporal prediction model integrating CA model, multi-criteria evaluation (MCE), and Markov chain (MC) in the Kuala Langat district, Malaysia.The results from above researches and other applications all demonstrated the effectiveness of CA-Markov model, which has a common applicability for modelling and realistic prediction of landscape patterns.
Changes in landscape pattern will lead to corresponding changes in the composition of ecosystem, and biodiversity.The change of habitat quality can directly exhibit the regional change of biodiversity and landscape pattern.The changes ecological processes are of great significance to regional maintenance of ecosystem service function and sustainable management.There were many models developed to evaluate the function of ecosystems services, such as Multiscale Integrated Models of Ecosystem Services (MIMES) [26], Artificial Intelligence for Ecosystem Services [27], Biodiversity module in Integrated Valuation of Ecosystem Services and Tradeoffs model (InVEST) [28] and so on.InVEST was developed as part of the Natural Capital Project [29], a partnership between Stanford University, the University of Minnesota, The Nature Conservancy, and World Wildlife Fund in 2007, whose aim is to provides a powerful tool for simultaneously quantifying and valuing multiple ecosystem services generated by a landscape.InVEST can evaluate various service functions of ecosystems such as biodiversity, soil conservation and carbon reserves, with advantages of less input data, accurate output data and intuitive spatial visualization.InVEST thus was widely used by its low application cost, high evaluation accuracy, and strong spatial analysis functions, which has been successfully applied to comprehensive assessment of the ecosystems in different areas [30][31][32].Patterns in biodiversity are inherently spatial, therefore, which can be estimated by analyzing landscape maps in conjunction with threats.The above researches all obtained good evaluation results, and the InVEST model is proved scientific and effective.
CA-Markov model and InVEST have achieved effective results in their respective fields, less research has been done about the combination of two above models.By varying landscape and evaluating the output from InVEST, we can provide useful information to managers and policy-makers weighing the tradeoffs in ecosystem services, biodiversity conservation, and other landscape objectives.The Three Gorges Project is the largest hydropower project in the world today, the construction of which brings enormous benefits of flood control, electricity generation, shipping and water supply.Likewise, it has a long-term but far-reaching impact on the ecology and environment of the reservoir area and the whole basin.The ecological and environmental impact of the Three-Gorges Project has been the focus of all worlds' attention.As the special ecological functional zone, ecological security status of TGRA concerns social and economic sustainable development of the whole Yangtze River basin even the whole China.Owing to abundant in energy, resources and biodiversity, Hubei Section of TGRA has become the essential area of natural environment construction.Many scholars have studied the status of natural environment and the ecosystem health in this region, including ecological environment sensitivity, ecosystem health condition, soil erosion sensitivity and so on [33,34].
To better understand the impact of spatial and temporal dynamics in landscape pattern on biodiversity in Hubei Section of TGRA, it is important to extract the landscape pattern and calculate the habitat quality.This study therefore aimed to investigate and predict landscape pattern and habitat quality by combining CA-Markov model and InVEST model, and explore biodiversity responses under the influence of landscape pattern dynamics.Within the context of biodiversity decline, this study is of practical implications for regional biodiversity conservation, natural environment protection and economic sustainable development.

Study Area
Hubei section of Three Gorges Reservoir Area (TGRA) is located in the middle reaches of the Yangtze River, which lies at the head of the TGRA.The Hubei section of TGRA covers five parts (Yiling districts, main urban areas of Yichang city, Zigui county, Xingshan county and Badong county), extending from 110.06 • to 111.65 • E longitude and 31.06 • to 31.56 • N latitude (Figure 1).The total area of this land is approximately 11,895 km 2 , with a population of 1.48 million people.Mountains and hills are the dominant terrain, with a summit that reaches 1000 to 2000 m.In addition to Yangtze River, some headwater rivers and streams are located in northern parts of Badong county, southern parts of Xingshan county and southeastern part of Yiling districts.This region has a typical subtropical monsoon climate with an annual average precipitation from 1000 to 1400 mm.Sub-tropical evergreen broad-leaved forest and warm coniferous forest are the main vegetation types.The soil type mainly includes yellow brown soil, rock soil, yellow soil, and purple soil.

Landscape Map
To obtain the landscape pattern information of the study area, three sets of remote sensing images (respectively obtained in1990, 2000 and 2010) need to be processed.Landsat-5 TM data with low cloud cover in summer time were selected as basic remote sensing sources.The Landsat-5 TM products with a spatial resolution of 30m are available from the United States Geological Survey (USGS) [35] and have been calibrated and atmospherically corrected for gas, aerosol, and cirrus-cloud effects.Image mosaic and clip were done according to the administrative boundaries of study area, and remote sensing images of the study area in1990, 2000 and 2010 were obtained.
The landscape was classified on the reference of Land Use Category System of the Chinese Academy of Sciences.Table 1 shows the landscape category system for use in our study, which is divided into 6 classes of the first class and 25 classes of the second class.Based on three sets of TM images, the landscape pattern information was extracted by a combination method of supervised classification and visual interpretation.To validate the accuracy of the landscape classification, 145 field points were sampled in June 2010.The field accuracy of first-order classification reaches 90%, and the second-order classification reaches 88%, which meet the research need.

Landscape Map
To obtain the landscape pattern information of the study area, three sets of remote sensing images (respectively obtained in1990, 2000 and 2010) need to be processed.Landsat-5 TM data with low cloud cover in summer time were selected as basic remote sensing sources.The Landsat-5 TM products with a spatial resolution of 30m are available from the United States Geological Survey (USGS) [35] and have been calibrated and atmospherically corrected for gas, aerosol, and cirrus-cloud effects.Image mosaic and clip were done according to the administrative boundaries of study area, and remote sensing images of the study area in1990, 2000 and 2010 were obtained.
The landscape was classified on the reference of Land Use Category System of the Chinese Academy of Sciences.Table 1 shows the landscape category system for use in our study, which is divided into 6 classes of the first class and 25 classes of the second class.Based on three sets of TM images, the landscape pattern information was extracted by a combination method of supervised classification and visual interpretation.To validate the accuracy of the landscape classification, 145 field points were sampled in June 2010.The field accuracy of first-order classification reaches 90%, and the second-order classification reaches 88%, which meet the research need.According to the status of the study area, both natural and human factors were selected to analyze driving forces for evolution of landscape pattern.Evaluation factors from two aspects of natural environment conditions and location were chose as internal driving forces, including DEM, slope, aspect, precipitation, temperature, distance to waters, distance to railway, distance to highway, distance to national road, distance to provincial road, distance to county road, distance to urban land and distance to residential land.As for external driving forces, population density and regional Gross Domestic Product (GDP) were selected from aspects of social and economic development.
DEM data was obtained from the ASTER Global Digital Elevation Model (ASTER GDEM) product, which is developed by the U.S. National Aeronautics and Space Administration (NASA) and Japan's Ministry of Economy, Trade, and Industry (METI).The ASTER GDEM covers land surfaces between 83 • N and 83 • S, and referenced to the 1984 World Geodetic System (WGS84)/1996 Earth Gravitational Model (EGM96) geoid.The product of ASTER GDEM Version 002 with spatial resolution of 30 m are available from Land Processes Distributed Active Archive Center (LP DAAC) of United States Geological Survey (USGS) [36].The DEM data in 2010 of study region was processed and projected to an Albers equal area conic projection.Based on DEM data, slope and aspect map were created and extracted by using Slope and Aspect function of raster surface from 3D Analyst tools in ArcGIS 10.2.
The precipitation and temperature data were downloaded from the Meteorological Data Sharing Service System of China [37].The annual average precipitation and air temperature data from 72 meteorological stations in 2010 that fully covered the entire region were selected.The annual precipitation and temperature were interpolated from a point to a regional scale using a kriging method with a spatial resolution of 30 m.Then, cross-validation was performed, and the interpolating error were approximately 76.67 mm and 0.68 degree which meet the accuracy requirements.
The spatial distribution of waters and construction land are all derived from the landscape map of the study area.The buffer analysis and assignment were done by distance function in Idrisi 17, and the buffer map of waters, traffic road at different level and town was drawn by setting 1.0 km buffer radius.
Night-time light data is a powerful tool to estimate socio-economic status as well as map gross domestic product (GDP) distribution and urban expansion [38,39].The DMSP-OLS NTL composite has evolved into Version 4, containing global NTL time series 1992-2013.The composite for the year of 1992 and 2010 are available from the NOAA's NGDC [40].Due to lack of data before 1992, the year composite data of 1992 was considered as night-time light data for the year of 1990.The stable light image, which has discarded ephemeral lights and background noise, from the composite was selected to reflect GDP status.The stable light images in 1990 and 2010 of study region was processed and projected to an Albers equal area conic projection.
Furthermore, population data at the level of community in 1990 and 2010 were collected from the fourth and sixth census of China respectively.Based on point data, Tyson polygons were constructed and the population density of each polygons was calculated.Then, the regional population density in 1990 and 2010 was obtained by interpolating points to a regional scale using an ordinary Kriging method with a spatial resolution of 30 m.
In the comprehensive evaluation of multiple factors, different factors often have different dimensions and dimensional units.To eliminate the effect of dimension, dimensionless treatment should be adopted before analysis.The data after dimensionless processing can accurately reflect the information contained in the original data.According to the research needs and data characteristics, equalization method in dimensionless treatment was used to reduce the loss of information in data and preserve differences in variation degree among factors.The equalization method follows that: where, x i is the original value of each factor, and x i is the average value of all original value.

Methodology
There were four main steps in conducting this research.Figure 2 shows the details for each step.Firstly, the temporal and spatial variation characteristics of landscape pattern were analyzed by calculating the transition matrix of Markov model.Secondly, the habit quality was evaluated by using biodiversity module of InVEST model.Thirdly, the driving force for alterations in landscape pattern was analyzed.Finally, Ca-Markov model together with InVEST model were used to simulate future landscape pattern and habit quality in 2020 under different scenarios.

Markov Chain
Markov chain is the simplest of stochastic models which is a transition matrix [41], and has been widely used for land cover change studies at various spatial scales [42].In Markov model, the initial probability and the transition probabilities between different states are calculated, and the trend of landscape pattern changes over time is determined.Finally, the landscape pattern is predicted.Please note that the transition between landscape types is often non-directional, which means that the transition direction for the landscape types of current pixel allows for many possibilities in the next moment [43].The transition probability of each landscape type is different, and the size of each possibilities are affected by many factors.The transition of landscape types is bi-directional, which not only transit from the current types into other types, but also from other types to the current types [44].The types and quantity of landscape change constantly in the process of stochastic transition.The transition process of landscape types accords with the conditions of Markov's research, which is difficult to describe accurately by mathematical methods [45,46].The output of Markov model is the probability of transition, P ij is between state i and j.In a landscape with multiple land covers or land uses, the transition probability P ij would be the probability that a landscape type (pixels) i in time t0 changes to landscape type j in time t1.As the transitions are probabilities, it follows that: where, n is number of landscape types, and Pij should meet two necessary conditions: 0 According to the principle of map algebra, the Markov transition matrix is established to quantitatively analyze the direction and intensity of the transition between different landscape types of second class.This transition matrix also can be used to predict the future landscape at time t2.

Habitat Quality
The biodiversity module in InVEST model maps the quality of habitat for a target conservation objective [47].Landscape maps are converted into habitat maps by defining what landscape counts as habitat for various species [48].Habitat quality is a function of the landscape type in a grid cell, the landscape in surrounding grid cells, and the sensitivity of the habitat in the grid cell to the threats posed by the surrounding landscape [49].Landscape in surrounding grid cells can modify habitat quality in a grid cell.Each landscape type is given a habitat suitability or quality score of 0 to 1 for each particular measure of biodiversity with perfectly suitable habitat scored as 1 and non-habitat scored as 0 [50].
Sources of degradation were considered as those human modified landscape types (e.g., urban, agriculture, and roads) that cause edge effects [51,52].Edge effects refer to changes in the biological and physical conditions that occur at a patch boundary and within adjacent patches (e.g., facilitating entry of predators, competitors, invasive species, toxic chemicals and other pollutants).The sensitivity of each habitat type to degradation is general principles of landscape ecology and conservation biology [53,54], and specific to each measure of biodiversity.The habitat quality was calculated combining landscape types information and threats of biodiversity, and the equation is shown in formula (3): where the quality of habitat in grid cell x that is in landscape type j be given by xj +k 2 = 0.5.D xj is the total threat level in grid cell x with landscape or habitat type j, and the equation is shown in formula (4): where y indexes all grid cells on r's raster map and Y r indicates the set of grid cells on r's raster map.S jr is the sensitivity of habitat type j to threat r where values closer to 1 indicate greater sensitivity.Please note that each threat map can have a unique number of grid cells due to variation in raster resolution If S jr = 0 then D xj is not a function of threat r.Also note that threat weights are normalized so that the sum across all threats weights equals 1.
The impact of threat r that originates in grid cell y, r y , on habitat in grid cell x is given by i rxy and is represented by the following equations, where d xy is the linear distance between grid cells x and y and d r max is the maximum effective distance of threat r's reach across space.For example, if an exponential decline is selected and the maximum impact distance of a threat is set at 1 km, the impact of the threat on a grid cell's habitat will decline by ~50% when the grid cell is 200 m from r's source.If i rxy > 0 then grid cell x is in degradation source r y 's disturbance zone.ω r is the threat weight.β x indicates the level of accessibility in grid cell x where 1 indicates complete accessibility.Some input parameters in model were determined from the literatures [32,49] and 12 experts, who are all in the field of regional ecological assessment and environmental simulation from the key research and development project of typical vulnerable ecological restoration and protection in China.The following parameters were input into the model, including landscape maps (three periods), threats factor layers (including paddy field, dry land, urban land, rural resident land, industrial and traffic land and bare land), index of threats and sensitivity (Table 2), weight of threat factor (Table 3), sensitivity of landscape types to each threat, and accessibility to threats, and half-saturation constant.The habitat quality maps were obtained after running the biodiversity module, and the grid resolution was set to 30 m according to remote sensing data requirements.Habitat quality scores should be interpreted as relative scores with higher scores indicating landscapes more favorable for the given conservation objective.The score of landscape habitat quality cannot be interpreted as a prediction of species persistence on the landscape or other direct measure of species conservation.The InVEST habitat model does not convert habitat quality measures into monetary values.[55,56].The principle of LR rests on the analysis of a problem, in which a result measured with dichotomous variables such as 0 and 1 or true and false, is determined from one or more independent factors [57].The goal of LR is to find the best fitting (yet reasonable) model to describe the relationship between a dependent variable (the presence or absence of event) and a set of independent variables [58,59].Compared with linear regression analysis models, the advantage of LR is that the variables may be either continuous or discrete, the value of variable can be any value in the real number range, and they do not necessarily follow normal distributions.Maximum likelihood estimation is used in the algorithm of LR model after transforming the dependent variable into a logit variable [55,60].Logistic regression has been widely applied for driving force analysis by many researchers [58][59][60][61][62][63], which can generate the regression relation with dichotomous variables such as 0 and 1 or true and false [64].The equation of LR model follows (7): where, P i is the probability that a certain landscape type i may appear or disappear in each grid.Function Y is represented as logit(P i ), i.e., the log (to base e) of the odds or likelihood ratio that the dependent variable.X n denotes each influencing factor.∝ represents the constant term.β 1 , β 2 , • • • , β n are the partial regression coefficient of LR model, representing the impact of the independent variable X n on Y.When the coefficient is positive and statistically significant, which means that the event of Y event is more likely to occur as the value of corresponding independent variables increases under the control of other independent variables.When the coefficient is negative and statistically significant, which means that the odds of the event Y occurring decreases as the value of corresponding independent variables increases under the control of other independent variables.A coefficient of 0 does not change the odds one way or another.
In this study, LR establishes a functional relationship between the binary coded landscape locations (presence or absence of a landscape) and different factors that are recognized as playing a role in landscape development.For each landscape type, when the landscape type is changed, its value is set as Y = 1; otherwise, Y = 0. LR generates the model statistics and coefficients of a formula useful to describe the relationship between the binary coded landscape locations (presence or absence of a landscape) and different factors.Relative Operating Characteristics (ROC) is usually used to test the regression effect [65].
LR model estimates the probability of a certain event occurring, therefore, selecting influencing factors that have a significant impact on the landscape pattern, and determining the quantitative relationship between them.According to the status of the study area, both natural and human factors were selected to construct the index system of driving forces for landscape evolution, and each factor was spatialized.Before the analyzing work of LR model, it is necessary to normalize the data of different measuring scales.The results should be interpreted with caution if the data fail to normalized in a manner LR model needs.In the application for exploring influencing factors on landscape variation, the common solution is to create layers of binary values for each class of an independent parameter [59,66,67].Here, we used binary landscape data (when the landscape type is changed, its value is 1; otherwise, 0.) from all over the area.Then, the driving mechanism of landscape evolution from 1990 to 2010 can be explored by using LR model.

CA-Markov Model
CA-Markov model consists of Markov chain and Cellular Automata (CA) model.As a method of studying nonlinear science, CA model is a spatially dynamic model, whose time, space and state are discrete.CA model consists of four parts, including cellular and its state, cellular space, cellular neighborhood and transition rules.In cellular space, every cell has its limited and particular state, and will be update according to the defined local rules.The interaction of these local rules forms a dynamic evolution system.In a CA model, the transition of a cell from one landscape to another depends on the state of the neighborhood cells [68].A cell will have a higher probability to transit to land-use type 'A' than to a land-use type 'B' if the cell is in closer proximity to land-use type 'A'.Thus, the CA model not only uses the information of the previous state of a landscape type as done by a Markov model but also uses the state of neighborhood cells for its transition rules.
CA model has been widely used for land use and land cover (LULC) change analysis, which can be expressed by the following formula: where, S is the set of finite and discrete cellular states, t and t + 1 represent different moment in time, N is cellular neighborhood, and f is the cellular transition rules in local space.CA model adds spatial character to a Markov model, thus combining CA-Markov model.In CA-Markov model, each landscape pixel is a cellular, and the landscape type of each cellular refers to a cellular state.With the support of GIS software, the transition of cellular state is determined by calculating the transition area matrix and conditional transition probability maps, and the change of landscape pattern is simulated [23].CA-Markov module of IDRISI 17 was used to predict future landscape.The specific implementation process is as follows: In the first step, a Markov model is run using the landscape maps of different period to predict for future.In doing so, the CA Markov model makes the assumption that the driving factors affecting change in all time periods are same.The outputs of this step are a transition area matrix, a transition probability matrix and a series of transition probability maps (probability of each pixel to each landscape types).
In the second step, CA analysis is run using the landscape maps and outputs from the Markov module, i.e., transition area matrix, transition probability maps and contiguity.The combined cellular automata and Markov chain method of the CA Markov module adds an element of spatial contiguity as well as knowledge of the likely spatial distribution of transitions to Markov chain analysis.The inherent suitability of each pixel and the restriction factors for each landscape types are established.The inherent suitability of each pixel is obtained by using the transition probability maps in an iterative process of CA Markov module.It is considered that within a certain threshold value, the suitability is the highest, which is set for 1.When exceeding the threshold, the Sigmoid function decreases with the distance until the extremely unsuitable distance is reached, and the adaptability is set for 0. After exceeding this distance, the suitability will not change.Next, an Analytic Hierarchy Process (AHP) method is used to weight sum the suitability layers the results of LR analysis.Then, OWA aggregation is done by selecting Location 4. In addition, the restriction factors obtained by LR analysis also were selected to put in model.
In the third step, the start time and the number of CA cycles were determined.One year represents one iteration, the iteration is taken as 10.For our present study, the time period of 1990 and 2000 should be used as reference map for validating the corresponding dates predicted maps of 2010, because the data were not used to build the respective models.If the predicted 2010 map is similar to reference 2010 map, we considered that the model performed well.
In the final step, validation was done simply by cross tabulating the predicted and observed landscape map of respective years.In present study, the predicted landscape map of 2010 was used to compare with reference 2010 map.The Kappa index was used to test the prediction accuracy of model, which can be expressed by the following formula: Kappa = (P 0 − P c )/(P 0 + P c ) (10) where P 0 is the correct ratio of simulation, P c is the correct ratio of random simulation, and P p is the correct ratio of ideal classification.

Spatial Landscape Pattern Analysis
Woodland, farmland and grassland are the main landscape types in the study area, accounting for 95% of the total area, of which 75% is woodland (Figure 3).From 1990 to 2010, the landscape pattern of the study area has changed slightly (Figure 4).Among the landscape types, woodland, waters and construction land had the greatest change in proportion.The area of woodland, farmland and grassland has been decreasing gradually, and the area of waters and construction land increased in the past two decades.The area of unused land slightly declined.The area of woodland, farmland and grassland in 1990 were 9631.88 km 2 , 1598.50 km 2 and 815.38 km 2 .By 2010, the area of above landscape types was 9539.11km 2 , 1579.64 km 2 and 813.85 km 2 respectively.The area of waters and construction land in 1990 were 131.02 km 2 and 67.71 km 2 , and the area of above landscape types was 190.50 km 2 and 121.51 km 2 in 2010.These landscape variations fit the facts of the Three Gorges Dam construction, impoundment of the Three Gorges Reservoir, immigration and expansion of construction land in study area.
Figure 5 exhibits the landscape transition map from 1990 to 2010 by overlaying landscape maps of before and after periods.The landscape transition matrix was statistically sorted.The results showed that 98.62% of the landscape types did not change within 20 years, and only about 1.38% of the landscape types changed (Table 4).The most varied landscape types were woodland, farmland, waters and construction land.The land occupancy of Three Gorge Dam construction and the impoundment of TGR could be the reason for the decreasing area of woodland and farmland located with the average elevation of 140 m.The increasing area of waters and construction land could be explained by the impoundment of TGR, migration and resettlement, and urbanization.Figure 6 illustrates the transfer in and out rate of landscape types.From 1990 to 2010, shoaly land and unused land had relatively high transfer-out rate, with a decreasing ratio of 41.86% and 27.05% respectively.Most of shoaly land mainly transferred into canal and reservoir or pond.Unused land mainly transferred into urban land, industrial and traffic land.Waters (canal, lake, reservoir or pond) and industrial and traffic land had relatively high transfer-in rate.Industrial and traffic land has highest transfer-in rate.The area of industrial and traffic land in 2010 was six times bigger than in 1990, and they were mainly transferred from woodland and farmland.Lake has the second-highest transfer-in rate with an increasing ratio of 136.36%, mainly transferring from woodland.
In the whole study area, waters and construction land have been the most dramatic changed landscape types.These variations mainly concentrated in main urban areas of Yichang city, extending downstream along the Yangtze River.Some of the variation located in Badong County.The expansion of waters was distributed around Three Gorges dam site in Zigui County and along the main and tributary banks of the Yangtze River.The strong increasing trend of construction land could be explained by the continuous improvement of settlement scale, the demand for supporting and improving facilities after three periods of immigration of the pilot project.

Analysis of Habitat Degradation
The value of habitat degradation index ranging from 0 to 1 represents the relative habitat degradation level of the current landscape, where 1 denotes high degradation and 0 denotes low degradation.Figure 7 shows the spatial distribution map of habitat degradation.The maximum value of habitat degradation index were 0.1188, 0.1194 and 0.1195 respectively, showing a gradually increasing trend from 1990 to 2010.In terms of spatial pattern, the areas with higher habitat degradation index were mainly distributed in main urban areas of Yichang city, southeastern parts of Yiling districts and south-central parts of Xingshan county, and the corresponding landscape types were urban land, rural residential land and canal.The areas with relatively low habitat degradation index were mainly distributed in the central and southwest of study area, such as Zigui county and Badong county, and the corresponding landscape type was woodland.Overall, the TGR area (Hubei section) developed rapidly from 1990 to 2010, causing a causing a continuously increasing trend of habitat degradation.

Analysis of Habitat Quality
Higher values indicate better habitat quality vis-a-vis the distribution of habitat quality across the rest of the landscape.To facilitate the observation of changes in habitat quality, the division method of natural breaks was adopted to assign grades according to the output values of habitat quality.The habitat quality was divided into five grades, they are low, relatively low, medium, relatively high, and high, and these grades are on behalf of poor habitat quality, relatively poor habitat quality, medium habitat quality, relatively high habitat quality, and high habitat quality respectively (Table 5).Based on divisions, the spatial distribution map of habitat quality grade in three phases was obtained (Figure 8), and the variations of habitat quality were statistically analyzed (Table 6).Overall, woodland, farmland and grassland occupied for 95% of the total area, thus the habitat quality of whole area was at a relatively high grade.The areas with high habitat quality grade was mainly distributed in Badang, Zigui and Xingshan county where the coverage of woodland.The areas with low habitat quality grade was mainly distributed in main urban areas of Yichang city, southeastern parts of Yiling districts, southern parts of Xingshan county and shoaly areas along the banks of the Yangtze River, where the corresponding landscape types are construction land and canals with higher elevation.
The overall habitat quality score declined slightly from 1990 to 2010.Habitat quality achieved a score of 0.8653 in 1990, 0.8259 in 2000 and 0.8155 in 2010.Most of areas was in a relatively high grade of habitat quality, the overall score of habitat quality decreased over the past twenty years.The magnitude of the declines for overall habitat quality was likely to be slight, but obvious in some local areas.The area of landscape with high habitat quality grade decreased obviously.The area of landscape with low and relatively low habitat quality grade increased.In 1990, the areas with high grade of habitat quality occupied for 84.78% of the total area, and accounted for 84.02% in 2010 with a decline area of 0.76%.The area of low and relatively low habitat quality grade increased by 0.3% and 0.44% respectively.
The overall habitat quality continues to decrease and the urbanization process continues to develop combining with the analysis of the habitat degradation.From 1990 to 2010, the overall fluctuation of habitat quality was relatively small, and the variation areas mainly were distributed in urban areas, canal areas with higher elevation and areas along the banks of the Yangtze River.In general, the fluctuation of habitat quality tended to be stable as a whole.
Based on habitat quality maps of 1990 and 2010, spatial changes of habitat quality were obtained in the past two decades (Figure 9).The landscape transition showed a decreasing trend of habitat quality grade from high to low and relatively low grade in areas of study region.The changes of habitat quality grade from high to low grade were that large areas of woodland and grassland have been developed for residential land, industrial and traffic land.The changes of habitat quality grade from high to relatively low grade were the transitions from woodland and grassland to canals.The construction of the Three Gorges Dam, impoundment in the TGR area and environmental migration contributed the increasing area of construction land and waters, decreasing area of woodland and grassland, and a decline of habitat quality as a whole.

Driving Force for Spatial Variation of Landscape Pattern
With the construction of Dam, impoundment and reservoir resettlement in the TGR area, the landscape pattern has been changed, leading to the decline of habitat quality.According to statistics of landscape pattern, the most varied landscape types were woodland, farmland, waters and construction land.Hence, both natural factors and human factors were selected to explore the driving mechanism of spatial variation of landscape pattern.To investigate the most dominant influencing factors, the following 15 factors are listed, including Night-time light (X1), population density (X2), DEM (X3), slope (X4), aspect (X5), precipitation (X6), temperature (X7), distance to waters (X8), distance to railway (X9), distance to highway (X10), distance to national road (X11), distance to provincial road (X12), distance to county road (X13), distance to urban land (X14), and distance to rural residential land (X15).Table 7 shows the results of LR analysis.

Driving Force for Spatial Variation of Landscape Pattern
With the construction of Dam, impoundment and reservoir resettlement in the TGR area, the landscape pattern has been changed, leading to the decline of habitat quality.According to statistics of landscape pattern, the most varied landscape types were woodland, farmland, waters and construction land.Hence, both natural factors and human factors were selected to explore the driving mechanism of spatial variation of landscape pattern.To investigate the most dominant influencing factors, the following 15 factors are listed, including Night-time light (X 1 ), population density (X 2 ), DEM (X 3 ), slope (X 4 ), aspect (X 5 ), precipitation (X 6 ), temperature (X 7 ), distance to waters (X 8 ), distance to railway (X 9 ), distance to highway (X 10 ), distance to national road (X 11 ), distance to provincial road (X 12 ), distance to county road (X 13 ), distance to urban land (X 14 ), and distance to rural residential land (X 15 ).Table 7 shows the results of LR analysis.The results of correlation analysis denoted that, compared with human factors, natural factors played a more important role in the spatial variation of each landscape types.From 1990 to 2010, the most influencing factors were DEM, slope, precipitation and temperature.The mean value of ROC reached 0.91, and the ROC value of each type was more than 0.8, indicating the good performance of LR model in analyzing the driving factors of landscape pattern variation.
For spatial variation of farmland, the average ROC value was 0.88.DEM, slope and temperature had the relatively high correlation value, which means that these three factors were the crucial determinants for spatial evolution of farmland.For spatial variation of dryland, precipitation and temperature were considered special factors.For spatial variation of paddy field, DEM and slope were the most important explanatory variables.The occurrence of paddy fields is often closely correlated with terrain.The exits of lower elevation and gentle slope greatly impact the occurrence of paddy field.
For spatial variation of woodland, the average ROC value was 0.92.DEM, slope, temperature and precipitation had the relatively high correlation value, illustrating that climate and terrain restricted the occurrence of woodland.
For spatial variation of grassland, the average ROC value was 0.86.Nighttime light, temperature and precipitation had the relatively high correlation value, indicating that human activities and climate affected the presence or absence of grassland.Among them, nighttime light had the highest correlation value, illustrating that the intensification of human economic activities leads to the transition of grassland to artificial surface.
For spatial variation of waters, the average ROC value was 0.94.DEM, slope, temperature and precipitation had the relatively high correlation value.Specially, the coefficient that belongs to the variable nighttime light strongly departs from 0 and led to the inference that human economic activities has a higher effect on the development of canal than any other variable.Along with the increasing activities of human beings, the demand for canal increased as well as the awareness of canal conservation.Terrain and climate were the dominant explanatory variables for the occurrence of lake, reservoir or pond, and shoal land.There were significant negative correlations between the development of those waters and elevation, which concluded that the higher elevation, the less probability of occurrence of lake, reservoir or pond, and shoal land would be.Those waters mainly occurred in relatively gentle terrain, and has significant correlation with precipitation and temperature.
For occurrence of construction land, the average ROC value was 0.92.DEM, temperature and precipitation, and nighttime light showed the relatively high correlation value.Terrain, climate and human economic activities were the dominant explanatory variables for the occurrence of construction land.For development of urban land, temperature and precipitation had relatively high correlation value, which indicated that climate has a higher effect on the development of urban land.For development of rural residential land, nighttime light had the highest correlation value.The impact of nighttime light variation on the occurrence of rural residential land was quite remarkable.Frequent and intensive human activities have caused increasing rural residential areas.For occurrence of industrial and traffic land, DEM had highest correlation value, which concluded that the lower elevation, the more likely the occurrence of industrial and traffic land would cause.
For occurrence of bare land, the ROC value was 0.87.Nighttime light, DEM, Slope, temperature, precipitation and distance to urban land had relatively high correlation value, illustrating that human activities, terrain, and climate all affected the occurrence of bare land.As human activities intensified, there will be a certain damage to vegetation, thus leading to the increasing probability of occurrence of bare land.
Although natural factors were highly correlated with presence or absence of landscape, they were not the crucial factors to drive the evolution of landscape distribution.Natural factors (e.g., Terrain and climate) evolve relatively slow, and usually will not change much in a shorter period (<20 years).To a certain extent, the driving force of socio-economic factors plays a role in evolution of landscape pattern, which cannot be neglected.The rapid development of central China brings out the economic rise in recent years, which can be reflected through the changes of nighttime light data.In our study, changes in nighttime light data, as a restriction input factor, was selected to add into the forecast model of landscape pattern, thus the potential evolution direction could be predicted.

Biodiversity Variation under the Influence of Landscape Pattern
In this section, CA-Markov model was used to adjust the future landscape pattern under different scenarios.Based on simulated landscape patterns, future habitat quality maps were evaluated through InVEST model.Therefore, the impact of landscape pattern dynamics on biodiversity was evaluated.

Forecast of Landscape Pattern
To verify the accuracy of the prediction, the observed and predicted landscape map of 2010 were input into CROSSTAB module in IDRISI.The Kappa coefficient reached 0.89, which showed the well performance of prediction model.Using the validated CA-Markov models will therefore allow us to predict the landscape map of 2020 under the business as usual and the ecological protection scenario.
Under the business as usual scenario, the landscape pattern will not be affected by policy adjustments, and it would change according to the transition probability matrix of each landscape types from 2000 to 2010.In addition, nighttime light together with other restriction factors were used to put into CA-Markov model.Figure 10a exhibits the landscape distribution map of 2020 under the business as usual scenario.By 2020, the variation trend of landscape pattern was basically consistent with the previous decades (Figure 11).The area of farmland, woodland and grassland reduced continuously, and the area of construction land increased continuously.The urban land has the largest changing rate among all the landscape types, with a changing rate of 19.05% from 2010 to 2020.The area of urban land increases rapidly, reaching an area of 60.97 km 2 by 2020.The area of canal has decreased to a certain extent, with a changing rate of 3.25%.Following this business as usual growth pattern, the area of construction land will continue to grow rapidly, and the regional ecological security will be endangered, which is no longer sustainable.
Sustainability 2018, 10, x FOR PEER REVIEW 20 of 28 In our study, changes in nighttime light data, as a restriction input factor, was selected to add into the forecast model of landscape pattern, thus the potential evolution direction could be predicted.

Biodiversity Variation under the Influence of Landscape Pattern
In this section, CA-Markov model was used to adjust the future landscape pattern under different scenarios.Based on simulated landscape patterns, future habitat quality maps were evaluated through InVEST model.Therefore, the impact of landscape pattern dynamics on biodiversity was evaluated.

Forecast of Landscape Pattern
To verify the accuracy of the prediction, the observed and predicted landscape map of 2010 were input into CROSSTAB module in IDRISI.The Kappa coefficient reached 0.89, which showed the well performance of prediction model.Using the validated CA-Markov models will therefore allow us to predict the landscape map of 2020 under the business as usual and the ecological protection scenario.
Under the business as usual scenario, the landscape pattern will not be affected by policy adjustments, and it would change according to the transition probability matrix of each landscape types from 2000 to 2010.In addition, nighttime light together with other restriction factors were used to put into CA-Markov model.Figure 10a exhibits the landscape distribution map of 2020 under the business as usual scenario.By 2020, the variation trend of landscape pattern was basically consistent with the previous decades (Figure 11).The area of farmland, woodland and grassland reduced continuously, and the area of construction land increased continuously.The urban land has the largest changing rate among all the landscape types, with a changing rate of 19.05% from 2010 to 2020.The area of urban land increases rapidly, reaching an area of 60.97 km 2 by 2020.The area of canal has decreased to a certain extent, with a changing rate of 3.25%.Following this business as usual growth pattern, the area of construction land will continue to grow rapidly, and the regional ecological security will be endangered, which is no longer sustainable.Under the ecological protection scenario, the area of woodland and waters were strictly protected, and the transition of farmland and construction land into above two landscape types was rate of 17.72% from 2010 to 2020.The area of grassland reached 957.87 km 2 .The area of urban land increases slightly, with a changing rate of 3.42% and reaching an area of 52.97 km 2 by 2020.The area of shoaly land has decreased sharply, with a changing rate of 16.12%.Following this ecological protection pattern, the quality habitat will increase to the 2000's level and the natural environment will improve, leading to harmonious development between human and nature.

Prediction of Habitat Quality
Based on predicted landscape maps, the habitat quality maps under business as usual scenario and ecological protection scenario were obtained using InVEST model, and the variations of habitat quality were statistically analyzed.
Figure 12 exhibits the future habitat quality map under the business as usual scenario.The variation trend of habitat quality in this scenario was consistent with the previous decades.The overall habitat quality score declined slightly, with an increasing area of low habitat quality in some local regions (Table 8).Habitat quality achieved a score of 0.8051 in 2020, and declined by 0.0104 when compared with 2010.The area of landscape with low habitat quality grade continuously increased and the area of landscape with medium and above habitat quality grade continuously decreased.The areas with relatively low and low grade of habitat quality occupied for 15.30% of the total area, with an increasing rate of 0.14% when compared with 2010.From 2010 to 2020, the transition of landscape pattern from high to low habitat quality grade occurred in some parts of the study area, and most of areas with low-grade habitat quality were construction land.The main change process of high-grade habitat quality was that large area of woodland and grassland was developed for construction land, which concentrated in main urban areas of Yichang city, southeastern parts of Yiling districts and northern parts of Xingshan county.The increasing urbanization and the construction of industrial and traffic land, decreasing of woodland and grassland would contribute the continuously decline of habitat quality as a whole.Under the ecological protection scenario, the area of woodland and waters were strictly protected, and the transition of farmland and construction land into above two landscape types was required to strengthen.The transition setting was as follows: the transition probability from farmland and woodland to construction land reduced by 3% and 1% respectively, from waters to construction land decreased by 2%, and from farmland to woodland increased by 1%. Figure 10b illustrates the simulated landscape map of 2020 under the ecological protection scenario.By 2020, the variation trend of landscape pattern was quite different from the previous decades, showing a decreasing area of farmland, increasing area of woodland and grassland, and relatively constant area of waters (Figure 11).Grassland has the largest changing rate among all the landscape types, with a changing rate of 17.72% from 2010 to 2020.The area of grassland reached 957.87 km 2 .The area of urban land increases slightly, with a changing rate of 3.42% and reaching an area of 52.97 km 2 by 2020.The area of shoaly land has decreased sharply, with a changing rate of 16.12%.Following this ecological protection pattern, the quality habitat will increase to the 2000's level and the natural environment will improve, leading to harmonious development between human and nature.

Prediction of Habitat Quality
Based on predicted landscape maps, the habitat quality maps under business as usual scenario and ecological protection scenario were obtained using InVEST model, and the variations of habitat quality were statistically analyzed.
Figure 12 exhibits the future habitat quality map under the business as usual scenario.The variation trend of habitat quality in this scenario was consistent with the previous decades.The overall habitat quality score declined slightly, with an increasing area of low habitat quality in some local regions (Table 8).Habitat quality achieved a score of 0.8051 in 2020, and declined by 0.0104 when compared with 2010.The area of landscape with low habitat quality grade continuously increased and the area of landscape with medium and above habitat quality grade continuously decreased.The areas with relatively low and low grade of habitat quality occupied for 15.30% of the total area, with an increasing rate of 0.14% when compared with 2010.From 2010 to 2020, the transition of landscape pattern from high to low habitat quality grade occurred in some parts of the study area, and most of areas with low-grade habitat quality were construction land.The main change process of high-grade habitat quality was that large area of woodland and grassland was developed for construction land, which concentrated in main urban areas of Yichang city, southeastern parts of Yiling districts and northern parts of Xingshan county.The increasing urbanization and the construction of industrial and traffic land, decreasing of woodland and grassland would contribute the continuously decline of habitat quality as a whole.Figure 13 exhibits the future map of habitat quality under the ecological protection scenario.The variation of habitat quality in this scenario represented reverse trend to the previous decades.The overall score of habitat quality increased slightly with an increasing area of high and relatively high habitat quality grade in local regions (Table 9).Habitat quality achieved a score of 0.8155 in 2020, and increased by 0.001 when compared with 2010.The area of landscape with high and relatively high habitat quality grade continuously increased and the area of landscape with low habitat quality grade continuously decreased.By 2020, the areas with relatively high and relatively high grade of habitat quality occupied for 84.75% of the total area, with an increasing rate of 0.2% when compared with 2010.From 2010 to 2020, the transition of landscape pattern from low to high habitat quality grade occurred in some parts of the study area, with the corresponding transitions from bare land to woodland, grassland and waters.The transition of landscape pattern from low to relatively high habitat quality grade occurred with the corresponding transitions from farmland and shoaly land to woodland and grassland.The changes were mainly concentrated in the northwestern areas of Zigui county, the northern areas of Badong county, and the shoaly areas along the Yangtze River in main urban areas of Yichang city.From 2010 to 2020, the policy implementation of returning farmland to forest may lead to the transition of paddy field and dry land to woodland.The enhancement of ecological protection awareness and series policies on vegetation protection would bring about the transition of shoaly land and bare land to woodland.The policy implementation of  Figure 13 exhibits the future map of habitat quality under the ecological protection scenario.The variation of habitat quality in this scenario represented reverse trend to the previous decades.The overall score of habitat quality increased slightly with an increasing area of high and relatively high habitat quality grade in some local regions (Table 9).Habitat quality achieved a score of 0.8155 in 2020, and increased by 0.001 when compared with 2010.The area of landscape with high and relatively high habitat quality grade continuously increased and the area of landscape with low habitat quality grade continuously decreased.By 2020, the areas with relatively high and relatively high grade of habitat quality occupied for 84.75% of the total area, with an increasing rate of 0.2% when compared with 2010.From 2010 to 2020, the transition of landscape pattern from low to high habitat quality grade occurred in some parts of the study area, with the corresponding transitions from bare land to woodland, grassland and waters.The transition of landscape pattern from low to relatively high habitat quality grade occurred with the corresponding transitions from farmland and shoaly land to woodland and grassland.The changes were mainly concentrated in the northwestern areas of Zigui county, the northern areas of Badong county, and the shoaly areas along the Yangtze River in main urban areas of Yichang city.From 2010 to 2020, the policy implementation of returning farmland to forest may lead to the transition of paddy field and dry land to woodland.The enhancement of ecological protection awareness and series policies on vegetation protection would bring about the transition of shoaly land and bare land to woodland.The policy implementation of returning farmland to forest and the enhancement of ecological protection awareness would contribute the continuously increase of habitat quality as a whole.

Discussion
This study provides an effective way to monitor the evolution of landscape pattern and habitat quality, and understand the influence of landscape pattern dynamics on biodiversity in Hubei section of TGRA.However, it should be noted that the results must be interpreted with caution.
The world-famous Three Gorges Project (TGP) is the largest hydropower station in the world and also the largest water resources and hydropower project constructed in China.The construction of TGP will generate enormous benefits in flood control, power generation, shipping and water supply, however, also may have a long-term and far-reaching impact on the ecology and environment of TGRA.The construction of the 180-m-tall dam was officially started in 1994 and the water level rise to 175 m after the three Gorges Reservoir impoundment in 2009.The study period from 1990 to 2010 we selected can fully cover the time span of Three Gorges Dam project.From 1900 to 2010, the variation of the landscape and the value of habitat degradation index seem minimal, while the value of habitat quality exhibits relatively high scores.Given that a small variation of landscape pattern occurred, however, the area of construction land has increased more than 79.46%,

Discussion
This study provides an effective way to monitor the evolution of landscape pattern and habitat quality, and understand the influence of landscape pattern dynamics on biodiversity in Hubei section of TGRA.However, it should be noted that the results must be interpreted with caution.
The world-famous Three Gorges Project (TGP) is the largest hydropower station in the world and also the largest water resources and hydropower project constructed in China.The construction of TGP will generate enormous benefits in flood control, power generation, shipping and water supply, however, also may have a long-term and far-reaching impact on the ecology and environment of TGRA.The construction of the 180-m-tall dam was officially started in 1994 and the water level rise to 175 m after the three Gorges Reservoir impoundment in 2009.The study period from 1990 to 2010 we selected can fully cover the time span of Three Gorges Dam project.From 1900 to 2010, the variation of the landscape and the value of habitat degradation index seem minimal, while the value of habitat quality exhibits relatively high scores.Given that a small variation of landscape pattern occurred, however, the area of construction land has increased more than 79.46%, showing the increasing effect of human activity.The rising value of habitat degradation index and decreasing value of habitat quality can explain the degrading trend of habitat quality in study area.
Generally, the impacts of Three Gorges Dam project on biodiversity and habitat quality are far-reaching and lasting, which cannot be ignored.Compared to other similar study results [69][70][71], dam construction has serious consequences for aquatic ecosystems, especially in fish species and habitat area, leading to an inevitable extirpation of small, dammed-off populations, habitat and fragmentation of landscape.As for biodiversity, many scholars showed evidences of pervasive biodiversity decline in this region [12][13][14].The declined value of habitat quality index over the past 20-year could exhibited the impacts of Three Gorges Dam project on habitat quality, which agreed with these above conclusions.Biodiversity measurement over large areas in traditional way by calculating the number of a species in a given area [72] is far more costly and time consuming.Thus, this study could provide possibility for measuring regional alterations in biodiversity by analyzing the variation of habitat quality through remote sensing technology.
When monitoring landscape pattern dynamics, precise spatial and temporal information is required.Landsat data has been used for land-use or landscape pattern mapping in many previous studies [22,45,73,74].The remote sensing data source in this study are mainly from Landsat-5 TM in summer time.However, a spatial resolution of 30 m and inconsistent acquisition time would affect the accuracy of the extraction results for landscape pattern.The acquisition time of remote sensing images would influence the extraction accuracy of waters.The water level of reservoirs and canals variated seasonally in a year.The extraction areas of waters from remote sensing data in flood period is generally bigger than in dry period.The results could change depending on the spatial-temporal resolution and the acquisition time of remote sensing images.
Given the lack of consideration of policy intervention impact on the landscape pattern, the future pattern prediction simulated by CA-Markov model obtained non-ideal performances under the ecological protection scenario.It is hard to quantify the effect of the policies [75][76][77], such as returning farmland to forestry, eco-migration policies [78], and check the coverage of all kinds of nature protection zones in space.Thus, the influence of policy intervention on landscape variation were not considered as a restriction factor.Similar studies of the influence of policy intervention on landscape variation were also found using a CA Markov model [79,80], which showed a positive impact of the National forest policy on the forest cover.If the quantified policy intervention added into the simulation model and succeed in reality, then post 2020 we expect more accurate landscape prediction and regional natural environment optimization.
As a proxy for biodiversity, habitat quality in InVEST model refers to the ability of the environment to provide conditions appropriate for individual and population persistence [81][82][83].The model assumes that areas with high habitat quality grade have higher biodiversity.Once similar habitats are destroyed, biodiversity is lost [5].However, areas with good habitat quality do not necessarily have high biodiversity.Moreover, the principle of this module is more inclined to vegetation diversity, which has certain limitations in the study area.In addition, we should recognize that a landscape has an artificial boundary where the habitat threats immediately outside of the study boundary have been clipped and ignored.Consequently, the threat intensity on the edges of a given landscape will always be less.All threats in InVEST model are additive, in some cases, the collective impact of multiple threats may be much greater than the sum of individual threat levels would suggest, which collective impact is not considered in InVEST model.

Conclusions
The objective of this research was to investigate the evolution of landscape pattern and habitat quality, and explore biodiversity responses under the influence of landscape pattern dynamics in Hubei section of TGRA.Several conclusions were made as follow.The overall landscape pattern has a slight but obvious change from 1990 to 2010.Terrain and climate were the most influencing natural factors compared with human factors.CA-Markov model performed well in future landscape prediction by setting socio-economic factors as restriction factor.The quality habitat has been declining over the past 20 years.Simulated habitat quality maps exhibited two opposite trends to the previous decade, indicating the impact of landscape pattern dynamics on biodiversity.The findings from this article offer the potential to analyze the impact of landscape pattern dynamics on habitat quality, and provide scientific basis for optimizing regional natural environment, as well as effective decision-making support to local government for landscape planning and management, regional biodiversity conservation, and economic sustainable development.

Sustainability 2018 ,
10,  x FOR PEER REVIEW 4 of 28 tropical evergreen broad-leaved forest and warm coniferous forest are the main vegetation types.The soil type mainly includes yellow brown soil, rock soil, yellow soil, and purple soil.

Figure 1 .
Figure 1.Geographic location of study area.

Figure 1 .
Figure 1.Geographic location of study area.The study region is situated in the transition zone from the second to the third terrain ladder in China, with abundant relic plants left over from quaternary glacier.As a rare plant gene bank in China, there are many national forest parks and natural reserves in this region.Abundant water and ore resources have made it one of the rarest eco-landscapes in China.Hence, the natural environment protection is of great importance to the construction of ecological civilization in Hubei province and even in China.China Western Development, water conservancy and hydropower projects, reservoir migration and other factors have exerted a profound impact on the ecological system of TGRA, leading to the increased risk of ecological degradation for a long time.Construction of the 180-m-tall dam was officially started in 1994.After the Three Gorges Reservoir impounding, the water level of dam has reached 135 m in 2003, 156 m in 2006, and 175 m in 2010.The construction of the Three Gorges Dam may have serious consequences for plant and animal species.Meanwhile, many positive implementations of policies and projects on vegetation protection promote the restoration for ecosystem structure and function of TGRA, such as Yangtze River protection forest engineering, natural forest protection project, and returning farmland to forest program.

Figure 2 .
Figure 2. Flow chat for evolution and prediction of landscape pattern and biodiversity process.
in H j and decreases in D xj .The k constant is the half-saturation constant and is set by the user.A grid cell's degradation score is translated into a habitat quality value using a half saturation function where the user determines the half-saturation value in InVEST model.The parameter k is equal to the D value where 1 − and grassland in 1990 were 9631.88 km , 1598.50 km and 815.38 km .By 2010, the area of above landscape types was 9539.11km 2 , 1579.64 km 2 and 813.85 km 2 respectively.The area of waters and construction land in 1990 were 131.02 km 2 and 67.71 km 2 , and the area of above landscape types was 190.50 km 2 and 121.51 km 2 in 2010.These landscape variations fit the facts of the Three Gorges Dam construction, impoundment of the Three Gorges Reservoir, immigration and expansion of construction land in study area.

Figure 4 .
Figure 4. Landscape area of first class in 1990, 2000 and 2010.

Figure 5 .
Figure 5. Spatial variation map of landscape pattern between 1990 and 2010.

Figure 6 .
Figure 6.Transfer rate of landscape types between 1990 and 2010.

Sustainability 2018 ,
10, x FOR PEER REVIEW 14 of 28 degradation index were mainly distributed in main urban areas of Yichang city, southeastern parts of Yiling districts and south-central parts of Xingshan county, and the corresponding landscape types were urban land, rural residential land and canal.The areas with relatively low habitat degradation index were mainly distributed in the central and southwest of study area, such as Zigui county and Badong county, and the corresponding landscape type was woodland.Overall, the TGR area (Hubei section) developed rapidly from 1990 to 2010, causing a causing a continuously increasing trend of habitat degradation.

Figure 9 .
Figure 9. Spatial variation map of habitat quality between 1990 and 2010.

Figure 9 .
Figure 9. Spatial variation map of habitat quality between 1990 and 2010.

Figure 10 .
Figure 10.Forecast maps of landscape in 2020 under (a) Business as usual scenario and (b) Ecological protection scenario.

Figure 10 .
Figure 10.Forecast maps of landscape in 2020 under (a) Business as usual scenario and (b) Ecological protection scenario.

Figure 11 .
Figure 11.The comparison of landscape area at first class in 2010, under business as usual scenario and ecological protection scenario in 2020.

Figure 11 .
Figure 11.The comparison of landscape area at first class in 2010, under business as usual scenario and ecological protection scenario in 2020.

Figure 12 .
Figure 12.Forecast map of habitat quality in 2020 under the business as usual scenario

Figure 12 .
Figure 12.Forecast map of habitat quality in 2020 under the business as usual scenario

Sustainability 2018 ,
10, x FOR PEER REVIEW 23 of 28returning farmland to forest and the enhancement of ecological protection awareness would contribute the continuously increase of habitat quality as a whole.

Figure 13 .
Figure 13.Forecast map of habitat quality in 2020 under the ecological protection scenario

Figure 13 .
Figure 13.Forecast map of habitat quality in 2020 under the ecological protection scenario

Table 1 .
Landscape category system.

Table 2 .
Landscape types and its sensitivity to each threat.

Table 3 .
Attributes of threat data.
2.3.3.Logistic Regression ModelLogistic Regression (LR) model is a nonlinear model, forming a multivariate regression relation between a dependent variable and a set of independent variables

Table 4 .
Statistics information of main landscape changing types from 1990 to 2010 (km 2 , %)

Table 5 .
Classification value of habitat quality in study area

Table 6 .
Area and percentage of habitat quality at all grades (km 2 , %)

Table 6 .
Area and percentage of habitat quality at all grades (km 2 , %)

Table 7 .
Results of Logistic regression analysis.

Table 8 .
Area and percentage of habitat quality at different grades in 2020 under business as usual scenario (km , %)

Table 8 .
Area and percentage of habitat quality at different grades in 2020 under business as usual scenario (km 2 , %)

Table 9 .
Area and percentage of habitat quality at different grades in 2020 ecological protection scenario (km 2 , %)

Table 9 .
Area and percentage of habitat quality at different grades in 2020 under ecological protection scenario (km 2 , %)