Quantitative Evaluation of Ecosystem Health in a Karst Area of South China

The purpose of this study is to propose a GIS-based mechanism for diagnosing karst rocky desertification (KRD) ecosystem health. Using the Huajiang Demonstration Area in Guizhou Province as a case study, this research offers a multi-factor indicator system for diagnosing KRD ecosystem health. A set of geologic, environmental, and socio-economic health indicators were developed based on remote sensing images from field-investigation, hydrological, and meteorological monitoring data. With the use of grid GIS technology, this study gives an indicator for diagnosing the spatial expression of desertification at a 5 m × 5 m grid scale. Using spatial overlaying technology based on grid data, the temporal and spatial dynamics of ecosystem health in the Huajiang Demonstration Area were tracked over a 10 year time span. The results of the analysis indicate that ecosystem health in the Huajiang Demonstration Area varies regionally, and has overall improved over time. The proportion of healthy area increased from 3.7% in 2000 to 8.2% in 2010. However, unhealthy and middle-health areas still accounted for 78.7% of the total area by 2010. The most obvious improvement of ecosystem health was in an area where comprehensive control measures for curbing KRD were implemented. These results suggest that comprehensive control of KRD can effectively mitigate ecosystem deterioration and improve ecosystem health in karst regions of South China.


Introduction
Karst areas are one of the most ecologically fragile ecosystems in the world [1].Of the world's three largest karst-concentrated areas (European Mediterranean coast, Eastern United States, and Southeast Asia), Southern China's karst area has the largest area of successive bare carbonate rock [2].Guizhou is located in the center of South China's karst area, and has the world's most complex, concentrated, and diverse karst landscape [3].The carbonate rocks are well developed and the outcropped area of karst occupies 61.9% of the province's territory [4].
The karst area in Guizhou Province is characterized by thin and discontinuous surface soil, high connectivity of surface water and groundwater, rapid hydrological cycling, uneven distribution of water and soil resources, high heterogeneity of hydrothermal factors over time and space, and soil rich in calcium (Ca) and magnesium (Mg) and deficient in nitrogen (N), phosphorus (P), and potassium (K).These characteristics result in a very fragile ecosystem, which, combined with the aridity of the climate, makes ecological health in this area particularly challenging to sustain.The fragility of this karst ecosystem is exacerbated by the large population of local inhabitants who have built settlements and established unsustainable agricultural practices in these areas.These human and non-human factors have resulted in karst rocky desertification (KRD), which is marked by reversed ecological succession and the decline of ecological health.KRD in South China has been recognized as a serious issue [5,6].Fortunately, the Chinese government has identified ecological restoration in karst areas as an important strategy for poverty alleviation and sustainable development in China.Policy measures have been adopted to help mitigate KRD, including the "Protection Forest System Construction Program" and "Sloping Land Conversion Program/Grain for Green Project".Between 2008 and 2010, the Chinese government spent 3 billion RMB to select 100 county-level demonstration areas (including 55 counties in Guizhou) for implementing comprehensive control measures to restore ecological health of KRD-affected areas.In 2011 through to 2016, the program expanded from 100 counties to 200 counties (including 78 counties in Guizhou) and a series of major projects have been implemented.As comprehensive control mechanisms of KRD are still in the initial stages, the impact of these mechanisms on the ecological health of karst areas needs to be explored.Therefore, it is necessary to set up a diagnostic framework for analyzing the spatial heterogeneity and temporal responses of ecosystem health in KRD-affected areas.
Assessment of ecosystem health in karst areas is a prerequisite for ecological restoration and sustainable development.However, few researchers have focused on ecosystem health of KRD-affected areas, especially those which have implemented comprehensive control measures.Cao and Su [29] used fuzzy comprehensive evaluation methods to evaluate and compare ecosystem health in three different types of karst area.Zhang et al. [30] applied the ecological footprint model to evaluate the ecosystem health of the Karst Plateau.These previous studies have two short-comings.First, there is no unified standard framework for diagnosing KRD-affected areas.Second, the frameworks used in these studies fail to consider spatial patterns and comprehensive changes of ecosystem health over time.
The purpose of this study is to propose a comprehensive and replicable method for diagnosing the ecosystem health of KRD-affected areas.Using the Huajiang Demonstration Area in Guizhou Province as a case study, this research offers a multi-factor indicator system based on grid GIS for diagnosing the health of KRD-affected ecosystems over time and space.This indicator system can be used to assess and inform the management, restoration, and reconstruction efforts of KRD-affected regions.

Study Area
The Huajiang Demonstration Area (HDA) is a typical karst dry-hot valley environment.It is located in the Huajiang Canyon of the Guizhou Plateau, between Guanling and Zhenfeng Counties.Its geographical coordinates are 105 • 36 30 E-105 • 46 30 E and 25 • 39 13 N-25 • 41 00 N (Figure 1).The total area of HDA is 51.2 km 2 and the total population of the demonstration area was 8612 by the end of 2010.The mean population density was 1.7 persons per km 2 .Carbonate rock occupies > 95% of the total demonstration area.The elevation ranges from 450-1450 m and the average elevation is 1000 m above sea level.HDA is in a subtropical humid monsoon climate zone, with an annual mean temperature of 18.4 • C and average annual precipitation of 1100 mm.Lime soil is distributed across this area, resulting in thin soil layers, discontinuous soil distribution, low forest coverage, poor water retention, and poor drought tolerance.With the exception of steep peaks with a few shrubs, desertified karst areas occupy > 70% of the area [4].   .The images were rectified to the Albers Equal Area Conic projection.Lithologic information was gained from the 1:200,000 scale regional hydrology geological map of the Guanling-Zhenfeng area (map number: G-48-XXII, Figure 2b).Altitude (DEM) data were taken from a 1:10,000 scale regional topographic map.
The map image was converted into a vector image.DEM data were then generated using the 3D Analyst tool in ArcGIS 10 which was released by ESRI (http://www.esri.com/) .Triangular Irregular Network (TIN) is a commonly used algorithm to construct Digital Elevation Mode (DEM).TIN data were created using the "Create TIN From Features" tool, and the "Surface Analysis" tool was then applied to generate a slope map from the TIN surface.Precipitation data in HDA was procured from the Guizhou Province Meteorological Information Center and a portable automatic weather station (model (DAVIS-Vantage Pro26162).Using the supervised classification method to classify land use types, field investigations and sample monitoring were conducted to verify the accuracy of the land use/land cover (LULC).LULC was divided into 13 types, including paddy field, dry land, garden, forest, woodland, grass, river and lake, residential land, industrial and mining land, roadway, and rock [31].We extracted the normalized difference vegetation index (NDVI) from remote sensing images using ENVI5.1 software, and then calculated vegetation coverage (FV) based on this figure [32].The landscape diversity index was extracted using FRAGSTATS 3.3 [33].The research method of soil erosion followed the national standard (SL 461-2009) [34], which divides soil erosion into five levels, including no erosion, mild erosion, moderate erosion, intense erosion, and very intense erosion.Building on the research of Xiong and Zhou [35], rocky desertification data were classified into non-karst areas, no rocky desertification, potential rocky desertification, mild desertification, moderate desertification, and strong rocky desertification.The cultivated land data were derived from the Second Chinese National Land Survey.By referring to data collected by  .The images were rectified to the Albers Equal Area Conic projection.Lithologic information was gained from the 1:200,000 scale regional hydrology geological map of the Guanling-Zhenfeng area (map number: G-48-XXII, Figure 2b).Altitude (DEM) data were taken from a 1:10,000 scale regional topographic map.
The map image was converted into a vector image.DEM data were then generated using the 3D Analyst tool in ArcGIS 10 which was released by ESRI (http://www.esri.com/).Triangular Irregular Network (TIN) is a commonly used algorithm to construct Digital Elevation Mode (DEM).TIN data were created using the "Create TIN From Features" tool, and the "Surface Analysis" tool was then applied to generate a slope map from the TIN surface.Precipitation data in HDA was procured from the Guizhou Province Meteorological Information Center and a portable automatic weather station (model (DAVIS-Vantage Pro26162).Using the supervised classification method to classify land use types, field investigations and sample monitoring were conducted to verify the accuracy of the land use/land cover (LULC).LULC was divided into 13 types, including paddy field, dry land, garden, forest, woodland, grass, river and lake, residential land, industrial and mining land, roadway, and rock [31].We extracted the normalized difference vegetation index (NDVI) from remote sensing images using ENVI5.1 software, and then calculated vegetation coverage (FV) based on this figure [32].The landscape diversity index was extracted using FRAGSTATS 3.3 [33].The research method of soil erosion followed the national standard (SL 461-2009) [34], which divides soil erosion into five levels, including no erosion, mild erosion, moderate erosion, intense erosion, and very intense erosion.Building on the research of Xiong and Zhou [35], rocky desertification data were classified into non-karst areas, no rocky desertification, potential rocky desertification, mild desertification, moderate desertification, and strong rocky desertification.The cultivated land data were derived from the Second Chinese National Land Survey.By referring to data collected by Xiong et al. [4], as well as primary investigation, a comprehensive map of the KRD control mechanisms in HDA was generated (Figure 2c).Social and economic data-including per capita GDP, population density, and population characteristics, were obtained through a combination of local census data and interviews and surveys with local residents.Hemeroby index (HI) data were derived from the geomorphology and land use types of HDA.HI was divided into five levels, including extremely high, high, moderate, low, and extremely low levels of disturbance.
Sustainability 2016, 8, 975 4 of 14 Xiong et al. [4], as well as primary investigation, a comprehensive map of the KRD control mechanisms in HDA was generated (Figure 2c).Social and economic data-including per capita GDP, population density, and population characteristics, were obtained through a combination of local census data and interviews and surveys with local residents.Hemeroby index (HI) data were derived from the geomorphology and land use types of HDA.HI was divided into five levels, including extremely high, high, moderate, low, and extremely low levels of disturbance.

Design and Calculation of KRD Ecosystem Health Index (K)
Ecosystem health index (K) is a criterion to diagnose the health status of karst regions.A scale from 0-5 was chosen as the basis for this ranking, where a K value of zero indicates the worst possible health state and a K value of 5 indicates the best possible health state.K was further divided into five gradients, with ranges of 0-2, 2-2.5, 2.5-3, 3-4, and 4-5, corresponding to the five health gradients, "Ill", "Unhealthy", "Middle Health", "Healthy", and "Excellent Health", respectively (Table 1).The division of the Ecosystem Health Index was built according to the karst ecological environment quality index and the results of existing research [36].
Ecosystem health index (K) was calculated according to the formula: where K is the comprehensive ecosystem health index, n is the number of diagnostic indicators, W i is the weight of the i-th diagnostic indicator, and E i is the health classification value.

Diagnosis Method in Grid GIS
Diagnoses of the health of the KRD-affected ecosystem was achieved using the following steps (Figure 3).(5) Gridded the layers by first using the "projections" tool in ArcGIS to convert the multi-source data into the same projection coordinates.Then, using the "polygon to raster" tool, all of the indicator data were converted from vector data into raster data.Third, a resolution grid of 5 m × 5 m was chosen as the basic unit of calculation; (6) Stacked the layers and calculated the layers using the "calculator raster" tool to calculate the comprehensive index of ecosystem health (K); (7) Used the "reclassify" tool to classify the maps of Ecosystem Health Index (K) according to classification table; and (8) Analyzed the data.

Diagnosis Indicator System of KRD Ecosystem Health
Karst ecological systems are formed by flows of energy, materials, and information among and between humans and the natural environment [37].The ecosystem health indicator system used in this study was based on 15 indices spanning geological, environmental, and socio-economic factors (Table 2).A healthy ecosystem can be defined in terms of three main features: vigor, organization and resilience [11,38].This definition reflects the comprehensive characteristics of the system, and provides a framework for assessing ecosystem health.Vigor is the ecological metabolism and nutrient recycling capacity of an ecosystem.The vigor of the KRD-affected area is based on the vegetation coverage and coefficient of the cultivated land.The coefficient of the cultivated land is the proportion of the cultivated land area to the total land area.There is a positive linear correlation between NDVI and vegetation coverage.
Fcover is the vegetation coverage value, ranging from 0 to 1.The higher the value, the greater the vegetation coverage.
Organization is the diversity of the ecosystem's composition.This is based on the landscape diversity index, soil erosion, and rocky desertification grades.Soil erosion is a process of external forces (water, wind, gravity) and human activities that lead to soil destruction or soil loss.Soil erosion and rocky desertification are the main problems of ecosystem degradation in HDA.Soil erosion decreases soil quantity and compromises ecosystem functioning, due to the slow rate of soil regeneration.Consequently, the distribution of soil in HDA is not continuous.Rocky desertification in HDA is dominated by moderate and strong rocky desertification.According to the 2010 rocky desertification area statistical data by Xiong and Zhou [35], the rocky desertification area is 27.8 km 2 , >50% of the total area of the HDA.Soil erosion and rocky desertification leads to the decline of soil fertility, the decline of production capacity, and damage of the ecosystem structure and functioning.
Landscape diversity is the foundation for ecological functioning, adaptation, and resilience.Landscape diversity is used to characterize the structure of the ecosystem.In this paper, Shannon's diversity index (H) was selected to quantify landscape diversity.
H is Shannon's diversity index; Pi is the proportion of the total landscape area of each patch type; and m is the total number of patch types.
Resilience refers to the ability of the system to recover gradually once external pressure is eliminated.Comprehensive control measures implanted in HDA that help maintain resilience

Diagnosis Indicator System of KRD Ecosystem Health
Karst ecological systems are formed by flows of energy, materials, and information among and between humans and the natural environment [37].The ecosystem health indicator system used in this study was based on 15 indices spanning geological, environmental, and socio-economic factors (Table 2).A healthy ecosystem can be defined in terms of three main features: vigor, organization and resilience [11,38].This definition reflects the comprehensive characteristics of the system, and provides a framework for assessing ecosystem health.Vigor is the ecological metabolism and nutrient recycling capacity of an ecosystem.The vigor of the KRD-affected area is based on the vegetation coverage and coefficient of the cultivated land.The coefficient of the cultivated land is the proportion of the cultivated land area to the total land area.There is a positive linear correlation between NDVI and vegetation coverage.
F cover is the vegetation coverage value, ranging from 0 to 1.The higher the value, the greater the vegetation coverage.
Organization is the diversity of the ecosystem's composition.This is based on the landscape diversity index, soil erosion, and rocky desertification grades.Soil erosion is a process of external forces (water, wind, gravity) and human activities that lead to soil destruction or soil loss.Soil erosion and rocky desertification are the main problems of ecosystem degradation in HDA.Soil erosion decreases soil quantity and compromises ecosystem functioning, due to the slow rate of soil regeneration.Consequently, the distribution of soil in HDA is not continuous.Rocky desertification in HDA is dominated by moderate and strong rocky desertification.According to the 2010 rocky desertification area statistical data by Xiong and Zhou [35], the rocky desertification area is 27.8 km 2 , >50% of the total area of the HDA.Soil erosion and rocky desertification leads to the decline of soil fertility, the decline of production capacity, and damage of the ecosystem structure and functioning.
Landscape diversity is the foundation for ecological functioning, adaptation, and resilience.Landscape diversity is used to characterize the structure of the ecosystem.In this paper, Shannon's diversity index (H) was selected to quantify landscape diversity.
H is Shannon's diversity index; P i is the proportion of the total landscape area of each patch type; and m is the total number of patch types.
Resilience refers to the ability of the system to recover gradually once external pressure is eliminated.Comprehensive control measures implanted in HDA that help maintain resilience include water resource management initiatives, reforestation projects, clean energy development, community building, and the reallocation of tourism revenues for ecological restoration programs.
Humans have been an important component of the ecological change of ecosystems for millennia, and they have always radically altered the systems of which they have been components [39][40][41].Social goals for sustainable ecosystem management are thus centered on maintaining the "ecological health" of the system [1].We chose to include per capita GDP, population density, and population characteristics.These three indicators were used to characterize the impact of human activities on the ecosystem health of karst areas.Hemeroby index is a measure of the human influence on ecosystems [42].

Assignment of Diagnosis Indicator
In order to normalize the data, an assignment method was chosen to convert each diagnosis indicator to values 1, 3, 5, 7, and 9 corresponding to five health states, "Ill", "Unhealthy", "Middle health", "Healthy", and "Excellent health" (Table 3).The assignment standard of each individual diagnosis indicator was based on national, industrial and local standards.The assignment standard of the lithology, altitude, and slope was built according to the "Comprehensive control of soil and water conservation: General rule of planning (GB/T15772-1995)".The altitude range was assigned from 500-1400 m, based on the elevation of HDA.The annual precipitation factor was built according to the assignment method of precipitation by Xuan and Luo [43].The land use factor was built according to the "Current land use classification in the People's Republic of China (GB/T 201010-2007)".The vegetation coverage and landscape diversity index factors were built according to the study by Fei et al. [44] and Abson et al. [45].The soil erosion factor was built according to the "Techniques standard for comprehensive control of soil erosion and water loss in karst region (SL461-2009)".The rocky desertification degree, coefficient of cultivated land and comprehensive control of KRD was built according to the "Rocky desertification comprehensive treatment planning outline in Karst area (2006-2015)".Social and economic indicators were built according to local construction standards.The assignment standard of the Hemeroby index was built according to the method by Chen et al. [46].

Determination of Weights
Weight is used to measure the relative importance of each factor in the diagnosis system.At present, the analytic hierarchy process (AHP) is a decision making method to analyze complex decision problems with multiple criteria [47,48].AHP is most widely used, as it is relatively mature.Therefore, AHP was used to determine the index weight in the health diagnosis of the ecological system (Table 2).To measure the consistency of the judgments made using AHP about the importance of each factor relative to large samples of purely random judgments, the consistency ratio (CR) was calculated [49].The CR was <0.1, which indicates the reliability of the judgments generated through AHP.

Temporal Dynamic Analysis of KRD Ecosystem Health
Using the health diagnosis index system and diagnosis model, the following conclusions regarding the ecological system health of HDA between 2000 and 2010 were drawn (Table 4).In Figure 4, the horizontal axis represents the health category of the ecological system and the vertical axis represents the proportion of the total area in this health range.

Spatial Distribution of KRD Ecosystem Health
Based on GIS analysis of ecosystem health in the HDA, there were notable temporal and spatial changes in ecological health between 2000 and 2010 (Figure 5).Over those 10 years, excellent health areas did not exist, but the area of healthy ecosystem increased.The most noticeably improved areas were distributed in peak cluster depressions, peak cluster valleys, and villages below 800 m altitude, including Yundongwan, Shuiyanba, Chaeryan, Sanjiazhai, and Kongluoqing villages.The distribution of middle health areas was uniform, located in different types of geomorphology and altitude, however, the area of middle health land of the hill peak platform along the Beipan River increased relative to other regions; unhealthy areas were concentrated in the hill peak platform and erosion platform, and desertification was serious in the area.Unhealthy and middle health were the most common land health types of HDA.Unhealthy areas were mainly distributed north of the Beipan River, located in the northeast part of the demonstration area, where there was erosion of

Spatial Distribution of KRD Ecosystem Health
Based on GIS analysis of ecosystem health in the HDA, there were notable temporal and spatial changes in ecological health between 2000 and 2010 (Figure 5).Over those 10 years, excellent health areas did not exist, but the area of healthy ecosystem increased.The most noticeably improved areas were distributed in peak cluster depressions, peak cluster valleys, and villages below 800 m altitude, including Yundongwan, Shuiyanba, Chaeryan, Sanjiazhai, and Kongluoqing villages.The distribution of middle health areas was uniform, located in different types of geomorphology and altitude, however, the area of middle health land of the hill peak platform along the Beipan River increased relative to other regions; unhealthy areas were concentrated in the hill peak platform and erosion platform, and desertification was serious in the area.Unhealthy and middle health were the most common land health types of HDA.Unhealthy areas were mainly distributed north of the Beipan River, located in the northeast part of the demonstration area, where there was erosion of steep slopes at altitudes 1000-1200 m.The land use types were grassland and bare rocks, characterized by steep slopes, soil erosion, serious rock desertification, and low vegetation coverage.

Discussion
The establishment of a diagnostic indicator system was useful in accurately diagnosing the ecosystem health of the KRD-affected area.Compared with other health indicator systems previously used in China [29,30], the diagnostic method based on grid GIS is more accurate at expressing spatial differences in ecosystem health.The selection of the diagnostic indicators and weight of the ecosystem health varies from area to area [50], and as such, the weighted indicator system applied in this study should only be applied for diagnosing the health of KRD-affected areas.
By analyzing the indicators and spatial characteristics of karst desertification, results show that human activity is the key factor affecting KRD ecosystem health.Appropriate human intervention, such as ecological restoration, bionomic control, and comprehensive control of KRD, can help alleviate the deterioration of karst ecological systems.In order to evaluate the effectiveness of comprehensive control measures on KRD on ecosystem health, ecosystem health must be tracked over time.A comparison of the 2000 map with the 2010 map shows that 78.3% of the total HDA area experienced no change in ecological health, 16.7% of total area experienced an improvement in ecological health, and 5.0% of total area experienced a deterioration in ecological health.The results indicate that the improvement of ecological health in areas with KRD control mechanisms was 36.9%,notably higher than the 9.6% increase in areas without control mechanisms.The results show that the comprehensive control mechanisms implemented in HDA have been effective in improving

Discussion
The establishment of a diagnostic indicator system was useful in accurately diagnosing the ecosystem health of the KRD-affected area.Compared with other health indicator systems previously used in China [29,30], the diagnostic method based on grid GIS is more accurate at expressing spatial differences in ecosystem health.The selection of the diagnostic indicators and weight of the ecosystem health varies from area to area [50], and as such, the weighted indicator system applied in this study should only be applied for diagnosing the health of KRD-affected areas.
By analyzing the indicators and spatial characteristics of karst desertification, results show that human activity is the key factor affecting KRD ecosystem health.Appropriate human intervention, such as ecological restoration, bionomic control, and comprehensive control of KRD, can help alleviate the deterioration of karst ecological systems.In order to evaluate the effectiveness of comprehensive control measures on KRD on ecosystem health, ecosystem health must be tracked over time.A comparison of the 2000 map with the 2010 map shows that 78.3% of the total HDA area experienced no change in ecological health, 16.7% of total area experienced an improvement in ecological health, and 5.0% of total area experienced a deterioration in ecological health.The results indicate that the improvement of ecological health in areas with KRD control mechanisms was 36.9%,notably higher than the 9.6% increase in areas without control mechanisms.The results show that the comprehensive control mechanisms implemented in HDA have been effective in improving ecosystem health.Similar conclusions about the effectiveness of human intervention in curbing KRD have also been made [51,52].

Conclusions
This study has proposed a framework to diagnose the ecological health of KRD-affected areas using grid GIS, with attention to temporal and spatial variations.The ecosystem health status of the Huajiang Demonstration Area is mainly in unhealthy and middle health, which is a result of desertification, soil erosion, and excessive human exploitation.However, the ecosystem health of the Huajiang Demonstration Area has gradually improved over time, especially in areas where comprehensive treatment mechanisms have been applied.The results indicate that rocky desertification protection programs can effectively control the deterioration of the ecosystem in karst areas and improve the health of local ecosystems.Our GIS-based health diagnosis of the HDA area was consistent with field survey results.As such, the comprehensive treatment projects reported in this study may provide an effective model for ecosystem restoration programs in karst areas in South China.
and poor drought tolerance.With the exception of steep peaks with a few shrubs, desertified karst areas occupy ˃ 70% of the area[4].

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

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

Figure 2 .
Figure 2. (a) Geomorphologic map; (b) lithologic map; and (c) map of the comprehensive KRD control mechanisms of HAD.

Figure 2 .
Figure 2. (a) Geomorphologic map; (b) lithologic map; and (c) map of the comprehensive KRD control mechanisms of HAD.

( 1 )
Built a diagnosis indicator system of KRD ecosystem health; (2) Prepared a thematic map of HDA for each of the diagnosis indicators.Each diagnosis indicator represents a layer and each thematic map (layer) corresponds to an indicator in the attribute table; (3) Assigned all diagnosis indicators with corresponding value according to the assignment table; (4) Determined weight of each diagnostic indicator;

Figure 3 .
Figure 3.The procedure for diagnosing the health of KRD-affected areas using Grid GIS.

Figure 3 .
Figure 3.The procedure for diagnosing the health of KRD-affected areas using Grid GIS.
Comparing the comprehensive health diagnosis results of HDA of 2000, 2005, and 2010 shows that an excellent health zone did not exist.The proportion of healthy regions increased over time, from 3.7% in 2000 to 8.2% in 2010; and middle health areas increased to 33.2%, 37.1%, and 43.9% in 2000, 2005, and 2010, respectively.The proportion of unhealthy land was at its peak of 40.4% in 2000.The proportion of unhealthy areas dropped from 39.1% in 2005 to 34.8% in 2010.The proportion of ill areas decreased over time, from 22.7% in 2000, to 17.0% in 2005, and to 13.1% in 2010.Ultimately, these results show that the ecological system health of HDA has improved over time, but the overall ecological system remains in the unhealthy and middle health levels, accounting for 78.7% of the total area in 2010.The rate of ecological health improved faster in 2005-2010 than in 2000-2005.This was because there were no KRD restoration measures used between the years 2000-2005.However, such measures were adapted in 2008.By 2010, it was still too early to measure the full effects of these initiatives.Sustainability 2016, 8, x 10 of 14

Figure 4 .
Figure 4. Ecosystem health diagnoses of the Huajiang Demonstration Area.

Figure 4 .
Figure 4. Ecosystem health diagnoses of the Huajiang Demonstration Area.

Figure 5 .
Figure 5.The spatial distribution map of ecosystem health in the Huajiang Demonstration Area.

Figure 5 .
Figure 5.The spatial distribution map of ecosystem health in the Huajiang Demonstration Area.

Table 1 .
Classification of KRD ecosystem health.
High level of structural integrity, stability, and sustainability; perfect ecological function; none or nearly no negative impact from external human and/or non-human factors; no abnormalities; high resilience.II (3, 4] HealthyModerate level of structural integrity, stability, and sustainability; high ecological function; slight negative impact from external human and/or non-human factors; no abnormalities; moderate resilience.

Table 2 .
Indicators for assessing the health of KRD-affected areas.

Table 3 .
Health diagnosis classification system for KRD-affected areas.

Table 4 .
Diagnosis of ecosystem health in the HDA.
* The proportion of change.