Measurement of Ecological Land Use / Cover Change and Its Varying Spatiotemporal Driving Forces by Statistical and Survival Analysis : A Case Study of Yingkou City , China

Detecting the spatial-temporal dynamics mechanisms and driving forces of ecological land change will offer a scientific basis for the sustainable utilization and ecological protection of regions undergoing rapid urbanization. This study examined the variations in ecological land from 2000 to 2015 in Yingkou by spatial statistical analysis with four land use/cover maps (2000, 2005, 2010, and 2015) interpreted by Landsat Thematic Mapper images and Google Earth maps. This study also measured the varying spatiotemporal drivers of typical ecological land use/cover types by survival analysis. The results indicated that ecological land was the main land use/cover type, and that the area decreased over time with a large transformation into agriculture and urban lands with significant temporal dynamics. The driving mechanisms became more complicated over time, and different time-dependent covariates significantly impacted the ecological land, forest land, and wetland losses. However, the distance from the city and different ways were the most important influencing spatiotemporal variables of the loss of ecological land, forest land, and wetlands over time. This study demonstrated the combined influence of a series of factors on ecological land loss. The spatial determinants and their impacts varied over time, especially the accessibility factors.


Introduction
Land use has undergone fundamental changes and it has substantially influenced biogeochemical and hydrological cycles, regional environment, biodiversity, and ecosystem services [1][2][3][4].The world, especially the developing country of China, has experienced dramatic land use/cover changes and growing conflicts between economic development and conservation in the recent decades [5][6][7].Studies on land use/cover change have focused on changes to types of land use and land production, rather than the variations in the land's ecological functions [5,8].Reasonable methodologies of land use classification and monitoring programs have emphasized the notion that ecological functions should be adopted to manage natural resources, monitor the impact of urbanization and human activities on the ecological environment, and coordinate the relationship between production, livelihood, and ecological functions [9,10].Some studies have categorized land use/cover types into agricultural, urban, and ecological lands [5,6].In this frame, ecological lands, including forest land, wetlands, grassland, orchards, and other natural land use/cover types can maintain the ecological process, balance the regional ecosystem, meet reasonable human needs, and improve the quality of human lives [10][11][12].However, the causes, patterns, processes, and consequences of ecological land change still require further study.Identifying the driving forces of ecological land change is still essential to understanding the characteristics of the changes in order to mitigate the negative effects, promote the desired outcome, and protect natural resources in land-use planning.
The driving mechanisms of change to different ecological land types have been explored in some areas, such as forest land [13][14][15], wetlands [16,17], and ecological land [8,11,18].However, the driving mechanisms and factors of change differ for dissimilar ecological land types from a regional perspective, due to land use type features, the natural environment, and human activities [8].Thus, changes in different ecological land types have been rarely discussed.Ecological land, which is an important natural resource, is strongly influenced by the combination of the natural/spatial underlying factors, socio-economic factors, cultural factors, and land-use planning, rather than a single key variable [8].Physical factors, such as elevation and slope have been considered as generally restrictive factors to ecological land, and ecological land with flat slopes and lower elevation, which are more likely to be exploited into becoming agriculture land and urban land [8,11,19,20].Socio-economic factors, such as population and gross domestic product (GDP) have frequently affected ecological land change in two ways [5,8].On one hand, with rapid urbanization and economic development, populations move from the countryside to urban areas, and the conflict between regional construction land and ecological protection increases.On the other hand, people begin to enjoy holidays in rural areas with economic developments, which also affects regional ecological land patterns.Proximity factors (i.e., distance to city, distance to towns, and distance to residential lands) impact on the ecological land through locational effects on a small scale in that ecological land near areas with intensive human activities are most greatly affected by the living and production benefits of human beings [8,11,18].Frequently, accessibility factors (i.e., distance to national ways, province ways, and county ways) guide the expansion of urban and agriculture land and promote the land-use demand of different populations through influence on land price, and then they influence the spatial distribution of regional ecological land [8,21].Land-use policy and urban planning (i.e., conservation area and land use policy) also influence ecological land change [3,13,22,23].Climatic and hydrological conditions (i.e., annual rainfall, evaporation, quantity of inflow water) [16], soil organic matter [8], and other factors have also been linked to ecological changes.These factors might be the main determinants of the ecological land use/cover changes at the regional scale [24], especially in the region of China [8,11,18].Therefore, these driving factors for measuring the varying spatiotemporal forces of ecological land in the study area are the focus of this study.
Many mathematical approaches, such as literature reviews [25,26], reasoning [27,28], regressions [29,30], and spatial statistical models [31] have been developed to examine the relationship between land use/cover and the influencing factors.Mostly regression models are used, such as the logistic model [11], binary logistic regression [18], stepwise logistic regression [20], partial least square regression [5], and the spatial analysis model [8,32].Driving forces are identified at different times and regions by using these models.These approaches identify the main explanatory variables over a specified period without considering the long-term temporal changes in different periods.Whilst the driving forces of ecological land identified in different studies seem to be different, it is natural to develop the question: do the driving forces of ecological land change vary spatially and temporally?Few studies have shown that the driving force of ecological land varies based on the analyses of data in various models used over different time periods [11].However, the measurement of varying spatiotemporal forces of ecological land use/cover change by one model is rarely studied.Therefore, efficient methods must be introduced to determine and to evaluate the dynamic driving forces of ecological land change quantitatively.
Survival analyses have comprehensively considered the occurrence and timing of a research event, which have been beneficial for longitudinal research that handles censoring and time-varying explanatory variables, and that are superior to conventional statistical methods [33][34][35].The Cox hazard regression (a semi-parametric model for survival analysis) can easily incorporate time-dependent covariates that may change over time [35,36].Survival analysis is extremely useful for studying different kinds of events in natural and social sciences, such as the recovery of patients, animal death [37], and deforestation [38,39].An efficient approach is to analyze the varying spatiotemporal forces of land use/cover change [33,34,40] in order to identify the varying drivers and their effects.
To mitigate the negative effects of ecological land change, the varying driving forces of ecological land loss and their spatiotemporal features in Yingkou City were identified.Specifically, this study aimed to (1) investigate the ecological land change in Yingkou from 2000 to 2010, (2) quantify the varying spatiotemporal drivers of ecological land loss by survival analysis using the Cox hazard model with time-dependent variables, and (3) offer a scientific basis for natural resource protection and ecosystem strengthening.

Study Area and Data Processing
The Yingkou region (39 • 55 N-40 • 56 N and 121 • 56 E-123 • 02 E) is in the Northwestern part of Liaodong Peninsula, China (Figure 1).The region covers 527,500 ha, and mountainous areas account for approximately 50% of the area.Benefiting from these topographic and geography features, the area has abundant forest and wetland resources, making up about 46% and 6% of the total area, respectively [41].Over the past 15 years, the region has undergone remarkably rapid economic growth and urbanization, and a considerable decrease in forest land and wetlands [41].Therefore, the natural resources should be analyzed and monitored in Yingkou, to protect natural environment and ecosystem processes.Survival analyses have comprehensively considered the occurrence and timing of a research event, which have been beneficial for longitudinal research that handles censoring and time-varying explanatory variables, and that are superior to conventional statistical methods [33][34][35].The Cox hazard regression (a semi-parametric model for survival analysis) can easily incorporate time-dependent covariates that may change over time [35,36].Survival analysis is extremely useful for studying different kinds of events in natural and social sciences, such as the recovery of patients, animal death [37], and deforestation [38,39].An efficient approach is to analyze the varying spatiotemporal forces of land use/cover change [33,34,40] in order to identify the varying drivers and their effects.
To mitigate the negative effects of ecological land change, the varying driving forces of ecological land loss and their spatiotemporal features in Yingkou City were identified.Specifically, this study aimed to (1) investigate the ecological land change in Yingkou from 2000 to 2010, (2) quantify the varying spatiotemporal drivers of ecological land loss by survival analysis using the Cox hazard model with timedependent variables, and (3) offer a scientific basis for natural resource protection and ecosystem strengthening.

Study Area and Data Processing
The Yingkou region (39°55′ N-40°56′ N and 121°56′ E-123°02′ E) is in the Northwestern part of Liaodong Peninsula, China (Figure 1).The region covers 527,500 ha, and mountainous areas account for approximately 50% of the area.Benefiting from these topographic and geography features, the area has abundant forest and wetland resources, making up about 46% and 6% of the total area, respectively [41].Over the past 15 years, the region has undergone remarkably rapid economic growth and urbanization, and a considerable decrease in forest land and wetlands [41].Therefore, the natural resources should be analyzed and monitored in Yingkou, to protect natural environment and ecosystem processes.Mapper images and Google Earth maps using the ENVI 5.1 software.These land use/cover maps represent the same data sources, processing methods, and classification systems, and they are the most appropriate data to identify changes in land use/cover.Overall accuracies were approximately 90% and the Kappa coefficient was 0.8 [41,42].The land use/cover classification system included cropland, forest land, urban land, wetlands, and other land.In this classification, other land comprised grassland, bare land, and so on, which are less disturbed by human activities.In this study, the land use/cover classes were divided into three sub-types (i.e., ecological, agricultural, and urban lands) according to their land functions, and natural or human disturbance characteristics using the ESRI ArcGIS 10.2 software.Ecological lands, including forest land, wetlands, and other natural land provide ecological goods and services, and have natural characteristics [5].The spatial distribution of land use/cover types is shown in Figure 2.
Sustainability 2018, 10, 4567 4 of 18 same data sources, processing methods, and classification systems, and they are the most appropriate data to identify changes in land use/cover.Overall accuracies were approximately 90% and the Kappa coefficient was 0.8 [41,42].The land use/cover classification system included cropland, forest land, urban land, wetlands, and other land.In this classification, other land comprised grassland, bare land, and so on, which are less disturbed by human activities.In this study, the land use/cover classes were divided into three sub-types (i.e., ecological, agricultural, and urban lands) according to their land functions, and natural or human disturbance characteristics using the ESRI ArcGIS 10.2 software.Ecological lands, including forest land, wetlands, and other natural land provide ecological goods and services, and have natural characteristics [5].The spatial distribution of land use/cover types is shown in Figure 2.

Ecological Land Change Assessment
Ecological land use/cover changes were analyzed using four land-use maps with the same land use/cover classification system for the reference years 2000, 2005, 2010, and 2015, by spatial statistical analysis.The spatial statistic (i.e., cross-tabulation analysis) was used to construct the transfer matrix of land use/cover by overlaying and intersecting land use/cover maps [43,44]  Ecosystem service values (ESV) were also calculated to reflect the ecosystem functions change of ecological lands, which can be expressed as follows [45]:

Ecological Land Change Assessment
Ecological land use/cover changes were analyzed using four land-use maps with the same land use/cover classification system for the reference years 2000, 2005, 2010, and 2015, by spatial statistical analysis.The spatial statistic (i.e., cross-tabulation analysis) was used to construct the transfer matrix of land use/cover by overlaying and intersecting land use/cover maps [43,44]  Ecosystem service values (ESV) were also calculated to reflect the ecosystem functions change of ecological lands, which can be expressed as follows [45]: where ESV indicates the total value of the ecosystem services in the ecological land; A i donates the area of the ecological land i (i.e., forest land, wetlands, and other natural land); VC signifies the valuation coefficients of ecosystem services j in ecological land i.The value coefficients for ecosystem services of different ecological land use/cover were calculated by the evaluation method for ecosystem service values based on a per unit area [45].

Driving Forces
Landscape change is affected by distinct combinations of cultural, political, and natural drivers [46].Many studies have recognized the analytical framework of the influencing factors of ecological land change, the variables of which were considered in the survival analysis in the current study with the available data (Table 1 and Figure 3).( 1) Elevation and slope are generally restrictive factors of ecological land, which may change in the area with flat slope and low elevation [8].(2) Two socio-economic factors, population and GDP [5,8], affect land change because the conflict between regional construction land and ecological protection increases with rapid urbanization and economic development.(3) Three proximity factors; namely, distance to city, distance to towns, and distance to residential lands [8,11,18], affect ecological lands through locational effects on a small scale; that is, ecological lands near intensive human activity areas are expected to be most greatly affected by the living and production benefits of human beings.(4) Three accessibility factors, namely, the distance to national, province, and county ways [21] impact land use/cover in promoting the land use/cover demand of different populations by influencing the land price.All variables were normalized into [0, 1] with a z-score standardization method and used as the data for determining the drivers of ecological land change.

Survival Analysis
A survival analysis can be used to investigate the occurrence and timing of events such as land use/cover change for a specific period [47].For ecological land loss, the conversion of land units from ecological land to other types is equivalent to an individual death event, which indicates that the survival time of an ecological land cell is the time duration until its state conversion [40].
In the survival analysis, the unchanged probability of ecological land can be calculated using the following formula [34]: where T represents the time for ecological land unchanged, represents the cumulative unchanged rate of ecological land, which implies that the ecological land remains unchanged after the observation time t, r P represents the probability of ecological land that is unchanged.The derivative of ) (t S is equivalent to the hazard function:

Survival Analysis
A survival analysis can be used to investigate the occurrence and timing of events such as land use/cover change for a specific period [47].For ecological land loss, the conversion of land units from ecological land to other types is equivalent to an individual death event, which indicates that the survival time of an ecological land cell is the time duration until its state conversion [40].
In the survival analysis, the unchanged probability of ecological land can be calculated using the following formula [34]: where T represents the time for ecological land unchanged, S(t) represents the cumulative unchanged rate of ecological land, which implies that the ecological land remains unchanged after the observation time t, P r represents the probability of ecological land that is unchanged.The derivative of S(t) is equivalent to the hazard function: where h(t) represents the risk rate of ecological land loss between t and t + ∆t( ∆t → 0 ), given that the ecological land is unchanged over time t.
In measuring the time-varying characteristics of the spatial driving factors, the Cox model with time-independent covariates was used to obtain h(t) [48].The proportional hazard model with time-dependent covariates, as an extension of the Cox model, is usually written as [35]: The formula can also be expressed as follows: where h i (t) represents the probability of ecological land loss over time t, λ 0 (t) denotes the baseline hazard function, which is related to t but is not the variable, x ik (t) represents the time-dependent variables, and it is sometimes expressed as The changing state of ecological land loss was coded as a chain: 0000, 1000, 1100, 1110, and 1111.The model requires the observations to be censored, such as left, interval, and right censored data [34].The forms of censored observations depend on whether and when the ecological land loss occurs (Figure 4).A combination of systematic and random sampling schemes was used to minimize autocorrelation, and to maintain the representativeness of the samples [33,49].), given that the ecological land is unchanged over time t.
In measuring the time-varying characteristics of the spatial driving factors, the Cox model with timeindependent covariates was used to obtain ) (t h [48].The proportional hazard model with time-dependent covariates, as an extension of the Cox model, is usually written as [35]: ( ) The formula can also be expressed as follows: where represents the probability of ecological land loss over time t, ) ( 0 t λ denotes the baseline hazard function, which is related to t but is not the variable, according to land use/cover, due to the difficulty in obtaining the time of land use/cover change.The changing state of ecological land loss was coded as a chain: 0000, 1000, 1100, 1110, and 1111.The model requires the observations to be censored, such as left, interval, and right censored data [34].The forms of censored observations depend on whether and when the ecological land loss occurs (Figure 4).A combination of systematic and random sampling schemes was used to minimize autocorrelation, and to maintain the representativeness of the samples [33,49].The Cox regression model was implemented on the basis of statistical software SAS9.4.Three formulas (i.e., BRESLOW, EFRON, and EXACT methods) can be alternated to handle the tied data in which two events occur simultaneously.Schoenfeld residuals were used to detect possible departures from the proportional hazard assumption [50].If the assumption test is significant for whether the Schoenfeld residuals are correlated with time or some function of time, then the Schoenfeld residuals should be independent of time.Therefore, the correlation analysis between the Schonefeld residuals of spatial variable and time variable t was used to identify the time-dependent variable xik (t) for different types of ecological land.Three alternative variables, namely, a likelihood-ratio test, a score test and a Wald test, were used to display the null hypothesis that all coefficients were 0 [35].The −2 log likelihood, Alaile's The Cox regression model was implemented on the basis of statistical software SAS9.4.Three formulas (i.e., BRESLOW, EFRON, and EXACT methods) can be alternated to handle the tied data in which two events occur simultaneously.Schoenfeld residuals were used to detect possible departures from the proportional hazard assumption [50].If the assumption test is significant for whether the Schoenfeld residuals are correlated with time or some function of time, then the Schoenfeld residuals should be independent of time.Therefore, the correlation analysis between the Schonefeld residuals of spatial variable and time variable t was used to identify the time-dependent variable x ik (t) for different types of ecological land.Three alternative variables, namely, a likelihood-ratio test, a score test and a Wald test, were used to display the null hypothesis that all coefficients were 0 [35].The −2 log likelihood, Alaile's information criterion (AIC), and Schwartz's Bayesian criterion statistics were used to measure the fit of the models.A comparison of the fitted statistics can be used to choose an improved performance model with a low value in order to explore the driving forces of ecological land loss effectively.

Ecological Land Change
The amount of land use/cover changed between ecological land and non-ecological land during 2000-2005, 2005-2010, and 2010-2015 are shown in Tables 2 and 3. Ecological land was the main land use/cover type and it accounted for over 50% of the total area in four years.Forest land was the main component of ecological land, followed by wetlands.It accounted for over 80% of the ecological land.The proportional area of ecological land decreased from 55.02% to 51.76%, with an area of 17,110.6 ha.Each type of ecological land showed a decreasing trend except for water in 2005 and other lands in 2010.3).Most of the ecological lands in the Western area were transformed into urban lands from 2000 to 2010, whereas most of the scattered ecological lands were transformed into agricultural lands in 2010-2015 (Table 3 and Figure 5).The lost wetlands, which were concentrated in the Northwest coastal regions, were transformed into new urban lands through the reclamation of sea and water regions from 2005 to  4).The trend was mainly dominated by forest land, due to its occupation of a higher proportion of area and the coefficient ecosystem services.This trend indicated that the ecosystem service values of Yingkou had a fluctuating trend, which was not only impacted by the decreased area of ecological lands, but also influenced by the conditions of the environment features.

Influencing Factors of Ecological Land Loss by Survival Analysis
The results of the Cox proportional hazard regression using the three tied handling methods (EXACT, BRESLOW, and EFRON) with 10 variables, and the results of the correlation analysis between the Schonefeld residuals of the spatial variable and the time variable t are shown in Tables 5, A1 and A2.The different results of the correlations in the three tied handling methods means that the time-dependent variables xik (t) in these methods were different.Therefore, the EXACT model could explore the driving forces of ecological land loss effectively in Yingkou, as in the study by Allison (2010).For the EXACT model, all of the time-dependent variables were implemented into the model for ecological land.xgdp (t) and xprovince (t) were not implemented for forest land and xslope (t) was excluded for wetlands.
Table 5. Correlation between the Schoenfeld residuals and the time variable t using the EXACT method.4).The trend was mainly dominated by forest land, due to its occupation of a higher proportion of area and the coefficient ecosystem services.This trend indicated that the ecosystem service values of Yingkou had a fluctuating trend, which was not only impacted by the decreased area of ecological lands, but also influenced by the conditions of the environment features.

Influencing Factors of Ecological Land Loss by Survival Analysis
The results of the Cox proportional hazard regression using the three tied handling methods (EXACT, BRESLOW, and EFRON) with 10 variables, and the results of the correlation analysis between the Schonefeld residuals of the spatial variable and the time variable t are shown in Tables 5, A1 and A2.The different results of the correlations in the three tied handling methods means that the time-dependent variables x ik (t) in these methods were different.Therefore, the EXACT model could explore the driving forces of ecological land loss effectively in Yingkou, as in the study by Allison (2010).For the EXACT model, all of the time-dependent variables were implemented into the model for ecological land.x gdp (t) and x province (t) were not implemented for forest land and x slope (t) was excluded for wetlands.The results of the Cox hazard model with time-dependent variables using tied handling methods for different ecological land types from 2000 to 2015 are shown in Tables 6 and A3.The test of all models in the "Test Global Null Hypothesis: BETA = 0" section suggested that the null hypothesis (all coefficients of variables were 0) should be rejected, given a very strong statistical significance for each statistic (likelihood ratio, score, and Wald) (p < 0.001).Comparison of the fitted statistics of the three tied handling models showed that the Cox model using EXACT methods could measure the varying spatiotemporal forces of ecological land change effectively.The fitted statistics of the Cox model were higher than those generated from the Cox model with time-dependent covariates.For example, the fitted statistics AIC of the time-independent model for ecological land was 55,164.343,whereas that estimated for the time-dependent hazard models was 51,539.740.Evidently, use of time-dependent variables is necessary.x elevation (t) was not used to calculate h i (t) for ecological land because the coefficients were insignificant.The ecological land loss was caused by the combination of physical, socio-economic, proximity, and accessibility variables and the impact varied over time.The suggestion that the driving mechanisms of ecological land loss became more complex and diverse over time, and that multiple drivers exist rather than a single key driver, was consistent with the results of Wang et al. (2018), which analyzed the socio-economical drivers of ecological land from 1984-2012 in China.
The coefficients of the physical variable x slope and the interaction x slope (t) were −0.875 and 0.747, respectively, which could be combined (−0.875 + 0.747 × time) × x slope .Therefore, ecological lands in flat areas were more likely to be transformed but this trend was reversed in 2010.A possible reason for this was that the expansion of the urban and agricultural lands in the flat areas would save costs and promote ecological protection [11], and that the newly built areas and agriculture lands might therefore be located in areas with slightly steeper topographies due to the shortage of land supply in flat areas over time.Meanwhile, elevation had a stable positive effect on ecological land loss because the ecological land distribution in higher mountains would be affected by rainfall and soil erosion, as in the study by Xie et al. (2017).
Socioeconomic variables, economy, and population had different significant impacts on ecological land loss.The relative risk of ecological land loss greatly decreased during economic development, and the decreasing trend was mitigated over time.The increase in population posed a considerable threat to ecological land loss, as seen in the study result by Xie et al. ( 2014) and the impact reserve in 2015.This condition might be due to regions with relatively slow economic development improving the level of regional economic development by occupying ecological lands, which would enhance the contradiction between development and protection.In addition, the available resources were insufficient to meet the demand in densely populated regions.As a result, natural resources were at great risk, and they gradually decreased due to the existing demographic pressures.
The probability of ecological land loss increased near cities, towns, and residential lands in the previous periods [11].This finding indicated that high intensity human activities would have significantly negative impact on ecological protection.However, this trend was gradually reversed toward the end of the study period, which indicated that ecological lands near cities and towns must meet the supply-demand balance of ecological lands for human beings.
The coefficients of distance of ecological lands to roads suggested that the ecological lands farther from national, provincial, and county ways had higher chances of disappearing.This trend weakened over time, in that the ecological land had a higher loss risk at the national, provincial, and county levels during the period 2010-2015.This negative affection was consistent with the study by Xie et al. (2017).These varying relationships did not follow a monotonic pattern, due to the main function of national and provincial ways being external transportation, and because they are always far from urban development areas, as shown in Figure 3.The value of these accessibility variables may vary over time as the new road infrastructure is constructed.
Compared to ecological land loss, the variables x population , x town , x nation , x elevation (t), x gdp (t), x population (t), x town (t), and x province (t) were not the important drivers of forest land loss.In other words, the impact of physical, socio-economic, and proximity factors on forest land loss was stable over time, except for slope, distance to city, and residents.Moreover, the distance to different roads had an intensive positive impact on forest land, which was not consistent with the trend for ecological land.This finding might be due to the growth characteristics of forest land, and that forest land near the roads might be protected in order to maintain regional environments and ecosystems.Slope, distance to city, and residents were the most important spatial determinants of forest land loss, where the impact varied significantly over time.
The influencing variables of wetland loss evidently differed from those of ecological and forest land losses (Table 4) especially for topographic elements.The time-dependent variables x elevation (t) and x slope (t) which were removed from the model, showed that the relative risk of wetland loss increased in the region within steep-and low-elevation areas, and that the impact was stable.This spatial relationship between wetland loss and terrain features was related to the distribution characteristics and evolution of wetlands.The loss risk of wetlands was significant when the distance to different ways increased.This effect increased over time.However, socio-economic factors and distance to city and ways, were the important influencing variables of wetlands loss over time.
In summary, the drivers of different ecological land types and their effects varied spatially and temporally.Distance to city and different ways were the most influential spatiotemporal variables of the loss of ecological land, forest land, and wetlands over time.However, there were also different effects on the different ecological land types.For example, socio-economic variables (i.e., economy and population) had varying spatiotemporal impacts on the ecological land and wetlands, but they had a stable effect on forest land.The slope had a negative impact on ecological land and forest land during the early period, which converted to a positive effect during the final periods, whereas it always had a positive effect on wetlands.The differences between the different types were always due to spatial distribution, function, and human demands.There are however, other factors linked to these base features.Therefore, future research based on big data could be developed to analyze the dynamic driving mechanisms of ecological land change, as well as land use/cover change.

Policy Implications of Ecological Land Protection
Land use policy is the administrative approach that is adopted by government departments to ensure the rational and effective use of land resources.It is the guiding factor of land use/cover, and it guides and regulates the quantity and spatial distribution of urban land use types [22].Therefore, on the basis of the results of this study and literature reviews, the following policy proposals are made: (1) Relevant policies on ecological land in China should focus on ecological land protection.
For example, the natural forest protection project, the conversion of sloping farmland to forest, and other policies that influence the spatial distribution of ecological lands [23].Despite the numerous laws and policies on the protection of ecological land and urban restrictive development, implementation and supervision at this stage is insufficient.The decreasing trend of ecological land remains evident (Table 2) and much attention should be paid to the functions of law enforcement and supervising departments.(2) The driving forces of ecological land loss are more complex than they previously were.Thus, a large number of policy groupings, such as a combination of policies that consider physical, socio-economic, proximity, and accessibility factors, should be implemented to protect natural resources.For example, population and distance to ways (i.e., the main driving forces of all ecological land changes), scientific and rational population distribution, and traffic construction are required to protect important ecological lands.However, these drivers vary over time.Thus, policies considering varying drivers also require flexibility or timeliness.(3) Spatial factors have variable effects on dissimilar ecological land types.These effects vary over time.For example, distance to roads has a varying effect on the different ecological land types.Therefore, some policies that optimize road networks, guide human activities rationally, and that coordinate the conservation mechanism of different ecological lands may be effective in constructing regional ecological spatial patterns.(4) Many studies have shown that the implementation of land use policies on different natural ecological resources in China has not achieved satisfactory results in recent decades [51].
The reason is that different lands in China are managed and controlled by various and multiple functional departments, and that this land management system results in conflicts or gaps in land use planning [51].The construction of the Ministry of Natural Resources of the People's Republic of China provides a large opportunity to modify this situation.Therefore, some policies that comprehensively consider different resource utilization schemes and that coordinate different functions play an important role in balancing the development of important ecological land protection in the region.

Conclusions
This study analyzed the spatiotemporal characteristics and varying drivers of ecological land loss in Yingkou, China from 2000 to 2015.On the basis of spatial statistical and survival analyses, the following main conclusions are drawn.(1) Ecological land was the main land use/cover type in the study area, and it decreased over time.The loss of ecological land was mostly observed in Northwestern coastal regions with high urbanization.(2) The influencing spatial factors of the dynamic changes in ecological land varied among different ecological land types.For ecological land, population, GDP, and distance to the city and ways were the most important spatiotemporal variables.For forest land, distance to city and slope were the dominant spatial determinants of the loss over time.For the wetlands, population, and distance to the city and ways were the important influencing variables of the loss, and the impact of accessibility factors varied significantly over time.However, distance to city and different ways were the most important influencing spatiotemporal variables for the loss of ecological land, forest land, and wetlands over time.(3) The regulatory and enforcement efforts in constructing the Ministry of Natural Resources of the People's Republic of China should be improved, in order to balance the expansion of urban and agricultural lands, and the protection of ecological land.

Figure 1 .
Figure 1.Location of the study area in China.

Figure 1 .
Figure 1.Location of the study area in China.The land use/cover data with a pixel size of 30 × 30 m in 2000, 2005, 2010, and 2015 were sources from the study of Liu et al. (2016) and Zhang et al. (2017), interpreted based on Landsat Thematic
to reflect the area change between the ecological land use/cover and non-ecological land during the 2000-2005, 2005-2010, and 2010-2015 periods in Yingkou.
to reflect the area change between the ecological land use/cover and non-ecological land during the 2000-2005, 2005-2010, and 2010-2015 periods in Yingkou.

Figure 3 .
Figure 3. Driving factors of spatial variables for the implementation of models (2000).

Figure 3 .
Figure 3. Driving factors of spatial variables for the implementation of models (2000).
and β 1 , • • • , β k are the regression coefficients for the variables.β 1 , • • • , β k are the regression coefficients for the time-dependent variables.The time variable t was discretized into several periods (2000-2005, 2005-2010, and 2010-2015) according to land use/cover, due to the difficulty in obtaining the time of land use/cover change.
represents the time-dependent variables, and it is sometimes expressed as t x ik × [40], and 1 coefficients for the time-dependent variables.The time variable t was discretized into several periods(2000-2005, 2005-2010, and 2010-2015)

Figure 4 .
Figure 4. Illustration of the forms of observations.

Figure 4 .
Figure 4. Illustration of the forms of observations.

Figure 5 .
Figure 5. Sampling points for ecological loss between 2000 and 2015 in Yingkou.The ecosystem services in ecological land in Yingkou were 37,645 million in 2000 and they reduced to 22,318 million in 2005, increased to 39,303 million in 2010, and reduced to 27,846 million in 2015 (Table4).The trend was mainly dominated by forest land, due to its occupation of a higher proportion of area and the coefficient ecosystem services.This trend indicated that the ecosystem service values of Yingkou had a fluctuating trend, which was not only impacted by the decreased area of ecological lands, but also influenced by the conditions of the environment features.

Table 1 .
Selected preliminary variables for studying the driving forces of ecological land change.
nation Distance to national ways (m) Road maps from the historical traffic maps of Liaoning and Google Earth, computed using near analysis in ArcGIS x province Distance to province ways (m) x county Distance to county ways (m) Sustainability 2018, 10, 4567 6 of 18

Table 3 .
Transfer matrix between ecological land and non-ecological in Yingkou (ha).

Table 5 .
Correlation between the Schoenfeld residuals and the time variable t using the EXACT method.

Table 6 .
Results of the Cox hazard regression with time-dependent variables using the EXACT method.

Table A1 .
Results of the Cox proportional hazards regression using the tied handling methods.