remote Assessment of Urban Ecological Quality and Spatial Heterogeneity Based on Remote Sensing: A Case Study of the Rapid Urbanization of Wuhan City

: Rapid urbanization signiﬁcantly affects the productivity of the terrestrial ecosystem and the foundation of regional ecosystem services, thereby detrimentally inﬂuencing the ecological environment and urban ecological security. The United Nations’ Sustainable Development Goals (SDGs) also require accurate and timely assessments of where people live in order to develop, implement and monitor sustainable development policies. Sustainable development also emphasizes the process of protecting the ecological environment for future generations while maintaining the current needs of mankind. We propose a comprehensive evaluation method for urban ecological quality (UEQ) using Landsat TM/ETM+/OLI/TIRS images to extract remote sensing information representing four ecological elements, namely humidity, greenness, heat and dryness. An improved comprehensive remote sensing ecological index (IRSEI) evaluation model is constructed by combining the entropy weight method and principal component analysis. This modeling is applied to the city of Wuhan, China, from 1995 to 2020. Spatial autocorrelation analysis was conducted on the geographic clusters of the IRSEI. The results show that (1) from 1995 to 2015, the mean IRSEI of Wuhan city decreased from 0.60 to 0.47, indicating that environmental deterioration overwhelmed improvements; (2) the global Moran’s I for IRSEI ranged from 0.535 to 0.592 from 1995 to 2020, indicating signiﬁcant heterogeneity in its spatial distribution, highlighting that high and low clusters gradually developed at the edge of the city and at the city center, respectively; (3) the high clusters are mainly distributed in the Huangpi and Jiangxia districts, and the low clusters at the city center, which exhibits a dense population and intense human activity. This paper uses remote sensing index methods to evaluate UEQ as a scientiﬁc theoretical basis for the improvement of UEQ, the control of UEQ and the formulation of urban sustainable development strategies in the future. Our results show that the UEQ method is a low-cost, feasible and simple technique that can be used for territorial spatial control and spatiotemporal urban sustainable development.


Introduction
Urban ecological quality (UEQ) evaluation is an important field of urban ecology research and the basis of urban planning and ecological management. With the continuous expansion of urbanization, China's cities have achieved medium-high quality development. However, social problems, such as resource exhaustion, an imbalance of economic structure and environmental pollution, do appear frequently. It is urgent to improve the capacity to implement urban sustainable development. In 2015, United Nations (UN) member states unanimously committed to achieving the Sustainable Development Goals (SDGs) by thereby compensating for the deficiencies of the EI. This paper uses existing research from a new perspective to more accurately study urban socio-economic activity intensity and its relationship with the regional ecological environment. Using the RSEI will help in studying the interactions between human activities and natural ecology, and the resulting knowledge of theory, concepts and methods is expected to benefit local governments [20]. In recent years, the RSEI has been applied in ecological quality monitoring in 35 cities of China [19,21,22], Eurasia [23] and America [21,23]. The RSEI and the results of principal component analysis (PCA) have been combined to develop an ecological index [19,24]. However, using the PCA results in insufficient information utilization, as the adaptive nature of PCA algorithms inevitably limits the full use of the available information. For example, the RSEI obtained in two studies using only the first component for normalization ranged from 60% to 90%, which cannot guarantee adequate contribution rates.
Accordingly, the aim of the current study is to improve the RSEI calculation method by proposing an improved-comprehensive remote sensing ecological index (IRSEI) constructed by employing PCA and equal weights (EW). Our study overcomes the shortcomings of previous studies, which only considered the application of PCA in ecological quality assessment, and the resolved knowledge gaps are reflected in the comprehensive consideration of EW and the PCA method to determine the UEQ. The contribution rates of the eigenvalues of PCA and EW are taken as the weights. This method enables the full use of the available data and ensures that the value of the calculated IRSEI is ecologically optimal. In addition, more indicators could be integrated and the IRSEI reduces noise interference and makes optimal use of practical image information. These factors facilitate the reliable and quantitative monitoring of the regional ecological environment.
A comparison and evaluation of the differences in quality in large cities can improve the cognitive ability of the internal mechanism of the reciprocal feed-back relationship between the construction of megacities and regional ecological balance, and can provide a scientific reference for controlling the scale of urban sustainable development and ecological planning and regulation. Wuhan is one of the fastest growing cities in central China, but few studies have been conducted on quantitative UEQ monitoring based on remote sensing data. Therefore, we used a series of parameters obtained from remote sensing imagery to construct the IRSEI for the evaluation of the UEQ of Wuhan city from 1995 to 2020. In addition to the UEQ, we determined the temporal and spatial changes in the city. We present a discussion of the ecological changes caused by economic and social developments and natural conditions. Finally, we provide theoretical guidance and a scientific basis for ecological construction in Wuhan city.
The objectives of this study are to: (1) Use GIS and RS technology to construct the IRSEI efficiently by integrating multiple sensors, including the Landsat Thematic Mapper (TM), Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS); (2) Monitor spatial and temporal changes in UEQ in Wuhan from 1995 to 2020; (3) Explore the spatial differentiation characteristics of the IRSEI in Wuhan.

Study Area and Data Preprocessing
We select the rapidly urbanizing city of Wuhan as study area for ecological monitoring and assessment. Wuhan is the capital city of Hubei Province. Its geographical location is 29 • 58 -31 • 22 N and 113 • 41 -115 • 05 E (Figure 1). From the perspective of Wuhan's geographical location and the location of its basin, the development of Wuhan has had great impact on the environment of the whole Yangtze River basin, and even the whole country. Therefore, ecological assessment and policy-based restoration and protection in Wuhan are vital for the ecological restoration of the Yangtze River basin. The city has jurisdiction over six central urban areas and seven distant urban areas. The land area comprises 8494.41 km 2 . The permanent population was 10.91 million in 2018. The Yangtze and Han rivers meet there, forming a geographical pattern referred to as "two rivers and Remote Sens. 2021, 13, 4440 4 of 23 three towns". Wuhan has a subtropical humid monsoon climate, with abundant rainfall and sufficient heat throughout the year. The average annual temperature is 15.8 to 17.5 • C. The area is rich in ecological resources, with nearly 40% green coverage and more than 10 m 2 of green space per capita. These ecological resources are crucial for Wuhan to build an ecological civilization city and, therefore, are critical factors in the protection of the ecological environment. The area is rich in ecological resources, with nearly 40% green coverage and more than 10 m 2 of green space per capita. These ecological resources are crucial for Wuhan to build an ecological civilization city and, therefore, are critical factors in the protection of the ecological environment. In order to consider the quality of the remote sensing data, such as cloud cover and vegetation condition, we use data from Landsat 5 TM in 1995, Landsat 5 ETM in 2005 and Landsat 8 OLI in 2015 and 2020 as the main remote sensing data. RS data are particularly useful because they can be used for temporal and spatial monitoring [25]. Details on the satellite images used in this study are provided in Table 1. The source dates of the images are relatively close; therefore, differences caused by different seasons and vegetation growth states can be ignored. Owing to topographic differences in images at different times and the influence of illumination and atmospheric factors on surface reflectance, the selected images required preprocessing with radiometric calibration and atmospheric and geometric correction prior to the calculation of the IRSEI. The corrections were performed using the Environment for Visualizing Images (ENVI) software. The Fast Line-of-sight Atmospheric Analysis of Hypercubes (FLAASH) model was used for atmospheric correction to eliminate the radiation error caused by atmospheric absorption and scattering. The accuracy of radiation calibration was more than 95%, and that of atmospheric correction exceeded 85%. Further, the error in geometric correction was controlled to less than 1 pixel. The quadratic polynomial and the nearest neighbor methods were used to correct the geometry of the images and the preprocessed images of the study area were clipped In order to consider the quality of the remote sensing data, such as cloud cover and vegetation condition, we use data from Landsat 5 TM in 1995, Landsat 5 ETM in 2005 and Landsat 8 OLI in 2015 and 2020 as the main remote sensing data. RS data are particularly useful because they can be used for temporal and spatial monitoring [25]. Details on the satellite images used in this study are provided in Table 1. The source dates of the images are relatively close; therefore, differences caused by different seasons and vegetation growth states can be ignored. Owing to topographic differences in images at different times and the influence of illumination and atmospheric factors on surface reflectance, the selected images required preprocessing with radiometric calibration and atmospheric and geometric correction prior to the calculation of the IRSEI. The corrections were performed using the Environment for Visualizing Images (ENVI) software. The Fast Line-of-sight Atmospheric Analysis of Hypercubes (FLAASH) model was used for atmospheric correction to eliminate the radiation error caused by atmospheric absorption and scattering. The accuracy of radiation calibration was more than 95%, and that of atmospheric correction exceeded 85%. Further, the error in geometric correction was controlled to less than 1 pixel. The quadratic polynomial and the nearest neighbor methods were used to correct the geometry of the images and the preprocessed images of the study area were clipped using the vector data of the administrative districts of Wuhan. Other data sources included the administrative zoning map of Wuhan, digital elevation model data of Wuhan from the geospatial data cloud (http://www.gscloud.cn/sources/accessdata/310?pid=302 (accessed on 15 May 2021)) and the cloud platform of geographical national condition monitoring of China We combine principal component analysis (PCA) and the entropy value method to design synthetic indicators that facilitate quick and quantitative assessment of UEQ, based on humidity, greenness, dryness and the heat index. Using this method enables prioritizing the natural factors of the ecological evaluation system. The overall framework of IRSEI modeling, as shown in Figure 2, includes four main steps. First, we obtain Landsat Enhanced Thematic Mapper Plus (ETM+)/OLI/TIRS images and perform preprocessing, including atmospheric correction, radiometric calibration and image mosaic (see Section 2.1). Second, we derive four remote sensing indicators: humidity, greenness, dryness and heat. Third, we calculate the PCA components, obtain PC1 and use the entropy method to calculate the results that are used for the construction of the IRSEI. Finally, the characteristics of the spatial and temporal changes in the ecological quality of Wuhan over the past 25 years are determined, and the spatial heterogeneity of the city is analyzed.

Calculation of Component Indices
• Humidity index ( ) The Kauth-Thomas transform (also called the tasseled hat transform) is a linear transformation method based on multispectral imaging [26,27]. This method is widely used in ecological monitoring for data compression and removal of redundancy. The moisture The Kauth-Thomas transform (also called the tasseled hat transform) is a linear transformation method based on multispectral imaging [26,27]. This method is widely used in ecological monitoring for data compression and removal of redundancy. The moisture component obtained by this transform reflects moisture information in the soil and vegetation. A low humidity value indicates severe land degradation, low vegetation cover and a poor ecological environment. A high humidity value indicates sufficient soil moisture, rich surface vegetation cover and a good ecological environment.
• Greenness index (I ndvi ) The normalized difference vegetation index (NDVI) is often used to monitor vegetation growth [30] and directly reflects the quality of the regional ecological environment. This index is used in the classification of regional land cover, environmental change and vegetation. The NDVI greenness index is computed as follows [31]: where ρ4 represents the reflectance of the near-infrared band and ρ3 represents the reflectance of the red band.
• Heat index (I heat ) Land surface temperature (LST) refers to heat, which is related closely to vegetation growth, crop yield, surface water circulation, urbanization, other natural phenomena and processes and human activities [32]. LST can be used as a heat index to reflect the surface ecological environment. Several algorithms use thermal infrared technology to retrieve LST, including the atmospheric correction, single-window and single-channel algorithms. Comparison between LST retrieval results obtained using the atmospheric correction method and the actual measurement of LST indicates that the error is within 1 • C, thereby meeting research accuracy requirements. LST is generated using Equations (5)-(9) [33,34]: where DN is the pixel gray value, gain and bias are thermal infrared band excursions and L is the radiation brightness value. Equation (7) is a simplified form of the inverse function of Planck's formula, with K1 and K2 being the calibration parameters. All of the parameter values are available from the metadata file (MTL) of the satellite data. ε is the specific infrared emissivity and Remote Sens. 2021, 13, 4440 7 of 23 is calculated with the method proposed by Min [35]. λ is the central wavelength of the thermal infrared band and ρ = s 1.438 10 −2 mK.
Vegetation coverage (F veg ) refers to the ratio (%) of the vertical projection area of vegetation on the ground to the total statistical area. P veg is based on Landsat NDVI and adopts the dichotomy model of mixed pixels [36]. The calculation formula is as follows [37]: where NDVI is the normalized vegetation index, NDV I soil is the normalized vegetation index value of bare land and NDV I veg is the normalized vegetation index value of complete vegetation coverage. NDV I soil and NDV I veg were selected as NDV I max and NDV I min with a confidence level of more than 95%.
• Dryness index (I dry ) The dryness index refers to the quantification of soil desiccation, which is a condition detrimental to the ecological environment. As most urban construction land is located in our study area, the dryness index can be represented by combining the bare soil index (SI) and the built-up index (IBI) into a normalized building-bare-soil index (NDBSI) [29]. We proposed extracting the bare soil and building area by setting an appropriate threshold and, subsequently, calculating the NDBSI as a weighted average and employing the area ratio as the weight. NDBSI = (SI + IBI)/2 (10) (12) where ρ1, ρ2, ρ3, ρ4 and ρ5 have been defined earlier in the context of the humidity index.

Water Mask and Standardization
The humidity index reflects the moisture of the vegetation and soil. The area covered by water in the study area occupies a large proportion of the I wet , which reduces the advantage of vegetation and soil in I wet . Therefore, the calculated I wet is not a true reflection of the vegetation and soil moisture, and it is necessary to mask the water bodies present in the study area. We use a modified normalized difference water index (MNDWI) to mask these water bodies. The formula is: where ρ Green represents the reflectance of the near-infrared band and ρ MIR represents the reflectance of the red band.

Construction of the Improved Remote Sensing Ecological Index (IRSEI) Evaluation Model
First, we obtain the primary remote sensing ecological index based on PCA. The four indices are standardized to the range [0-1] and PCA is used to combine these indices. PCA1 is obtained from the four RSEIs to build a preliminary assessment model. Generally, the first PCA collects most of the information on the four indicators, and PC1 can be used to represent the characteristics of the regional ecological environment. Therefore, we use only one PC in further analyses. To facilitate index measurement and comparison, the initial RSEI is standardized, as follows: where I ndvi represents the green component; I wet represents the humidity component; I heat represents heat; I dry represents dryness; and PC1 is the first principal component. The obtained RSEI value is within the [0-1] range. e i is the characteristic value contribution rate of the index corresponding to PC1. The closer RSEI is to 1, the better the UEQ of the region. The first principal component analysis index values are listed in Table 2. A detailed description of the calculation steps is available in the relevant literature [11,22,24,29]. Second, we introduce the entropy value method, which determines the weight of each index according to the information provided by the observed values of each index [38,39]. The evaluation index system includes N indices (NDVI, WET, NDBSI and LST). This is a problem that consists of m samples (cell) and uses N indicators for comprehensive evaluation. The initial data matrix A of the evaluation system is formed and X ij is the value in i cell of the j remote sensing ecological indicator. The detailed procedures of the entropy method are described as follows [22,38,40]: Proportion of the value in i cell of the indicator j. 2.
Entropy value of the j th index.

3.
Difference coefficient of the first index.
For the j th index, the more significant the difference is in the index value X ij , the greater the effect on the scheme evaluation and the smaller the entropy value.
The larger the g j value, the more critical the indicator.

5.
Ecological index score based on the entropy method.
The PCA effectively removes redundant information between bands and compresses multiband image information into a few independent bands that are more effective than the original band. The entropy method can effectively remove deficiencies caused by a lack of PCA information. The weights for all of the indicators are listed in Table 3. Finally, the IRSEI integrates humidity, greenness, heat and dryness through PCA and EW, which is calculated according to Equation (19): In this formula, RSEI PCA is the main component, RSEI EW is the weighted result of the entropy method and the final IRSEI is calculated as their arithmetic average. The IRSEI for each year has to be standardized to accurately compare the remote sensing images of different time frames. The closer IRSEI is to 1, the better the UEQ (and vice versa). The IRSEI for the four years is classified into five groups employing the ArcGIS software (Esri, USA). Referring to previous studies [22][23][24]29,41], these groups are labeled "Excellent, Good, Moderate, Fair, and Poor" and they facilitate comparisons across the study area (Table 4). Global spatial autocorrelation (SA) measures the average correlation, spatial distribution pattern and significance of all of the objects in the entire study area. SA visualizes spatial aggregations and exceptions to the IRSEI. The Moran's index is commonly used to calculate SA [42]. The main calculation indices for spatial autocorrelation are the global Moran's index and the local Moran's index. We analyze both the "global" spatial clustering and the "local" spatial clustering of the IRSEI. The formula for calculating the global Moran's index is: where n is the total number of grid cells in the study area (500 m × 500 m); W ij represents the spatial weight of elements i and j; x i and x j are the attribute values of cell i and cell j, respectively; x represents the average value of the attributes across all cells; and S 2 is the sample variance. The value of the global Moran index I varies between −1 and 1, where I > 0 indicates positive SA, i.e., a high value corresponds to high-value clusters, whereas a low value corresponds to low-value clusters. The closer I is to 1, the smaller the overall spatial difference. When I < 0, there is negative SA, i.e., there is significant spatial difference between a cell and its surrounding cells. The closer I is to −1, the greater the overall spatial difference. When I = 0, there is no SA.
Due to the fact that the global Moran's index describes the overall aggregation situation, it cannot accurately determine where the place of aggregation is located and is unable to indicate the hot spots and cold spots of the entire region. Accordingly, we use the local indicator of SA to measure local SA and determine hot and cold spots. The formula for the local Moran's I i for cell i is: When the local Moran index I i > 0, the spatial difference between the cell and its surrounding cells is minor. When the local Moran index I i < 0, the spatial difference between cell i and its surrounding cells is significant. When the local Moran index I i = 0, there is no spatial difference between cell i and its surrounding cells. In this study, we use the software GeoDA to calculate and obtain the global and local Moran's indices.

Attributing Factors
A comparison of the spatial distributions of the four ecological factors in the study area ( Figure 3) shows high levels of land surface moisture close to and alongside the Yangtze River, which extends in the central part of Wuhan from west to east. The NDVI is high on the northeast side, along the Yangtze and Han rivers, the central part of Wuhan and in patches in the south and east. Comparing the NDVI, LST and moisture maps shows that moderate temperature and moisture are the most favorable conditions for vegetation growth, whereas extreme weather conditions can damage plant vitality. The temperature and moisture conditions are moderate in the study area and the NDVI is remarkably high. A high LST is detected in the southern part of Wuhan, with some patches in the north and east. A moderate LST is detected in the central part of Wuhan. The NDBSI does not display much variation, as most of the study area is covered by agricultural land (Figure 3). To test the representativeness of the index IRSEI, we calculate the correlation coefficients among IRSEI, WET, NDVI, NDSI and LST in the same period (Table S1, Supplementary Materials), and test the applicability of the model through average correlations. From 1995 to 2015, the average correlation of IRSEI with the other variables is the highest, ranging from 0.60 to 0.70. The mean correlation of IRSEI over this period was 0.64, which indicates that IRSEI integrates most of the information embodied in all four indicators. It is more representative than any single indicator and can better reflect the ecological situation. To test the representativeness of the index IRSEI, we calculate the correlation coefficients among IRSEI, WET, NDVI, NDSI and LST in the same period (Table S1, Supplementary Materials), and test the applicability of the model through average correlations. From 1995 to 2015, the average correlation of IRSEI with the other variables is the highest, ranging from 0.60 to 0.70. The mean correlation of IRSEI over this period was 0.64, which indicates that IRSEI integrates most of the information embodied in all four indicators. It is more representative than any single indicator and can better reflect the ecological situation.

Spatial and Temporal Distribution of UEQ in Wuhan
Generally, higher IRSEI values are associated with higher levels of greenness and moisture, whereas lower IRSEI values are directly proportional to dryness and temperature. This implies that high IRSEI values represent positive ecological conditions. As shown in Figure 4, the IRSEI increases from 0.79 to 0.98 from 2010 to 2015, indicating improved ecological conditions. However, from the second half of 2015 up to 2020, its value drops to 0.82, indicating deterioration. Comparing the values from 2010 to 2020 indicates overall improved conditions, as the IRSEI increases from 0.79 to 0.82. However, the maximum values (1.09, 1.03 and 0.96) decline continuously, indicating that high-quality IRSEI conditions are declining continuously. Further, low-quality IRSEI conditions improve in the first half of the study period (1995 to 2005); however, in the second half (2005 to 2020), these conditions decline and reach their previous stage. Our findings also show maximal variation in the median IRSEI values, i.e., indicating the recovery of favorable conditions (moderate to high temperature, moderate to low moisture and higher vegetation) for all factors during the study period.
Generally, higher IRSEI values are associated with higher levels of greenness moisture, whereas lower IRSEI values are directly proportional to dryness and tem ture. This implies that high IRSEI values represent positive ecological conditions. As shown in Figure 4, the IRSEI increases from 0.79 to 0.98 from 2010 to 2015, cating improved ecological conditions. However, from the second half of 2015 up to its value drops to 0.82, indicating deterioration. Comparing the values from 2010 to indicates overall improved conditions, as the IRSEI increases from 0.79 to 0.82. How the maximum values (1.09, 1.03 and 0.96) decline continuously, indicating that highity IRSEI conditions are declining continuously. Further, low-quality IRSEI condition prove in the first half of the study period (1995 to 2005); however, in the second half to 2020), these conditions decline and reach their previous stage. Our findings also maximal variation in the median IRSEI values, i.e., indicating the recovery of favo conditions (moderate to high temperature, moderate to low moisture and higher ve tion) for all factors during the study period.
The mean IRSEI value and the area and percentage of each evaluation grade in han from 1995 to 2020 are displayed on Figures 4 and 5. Overall, the proportion of with average and good IRSEI ratings is the highest during the study period (>57%) proportions of average and above average regions are 82.33%, 87.14%, 74.21% and 57 indicating that the ecological environment of Wuhan was unstable from 1995 to 2020, ecological conditions first improving and subsequently deteriorating. The UEQ of the zhou, Hanyang, Qiaokou, Huangpi and Caidian districts show the most obvious de with reduction rates of 32.32%, 30.18%, 27.84%, 27.67% and 27.24%, respectively. From the perspective of a single year (see Table 5), the area share of good ecolo environment in 1995 was the highest, reaching 43.67% of the total area. The area sha poor ecological environment was the lowest in 1995, comprising an area of only 371 or less than 5% of the total area. The share of poor ecological environment was app mately 12% of the total area. The area with a good ecological environment rating in was larger than that of 1995 and accounted for the highest proportion (45.57%), com ing an area of 3494 km 2 . The percentage of area rated excellent was the smallest (6  From the perspective of a single year (see Table 5), the area share of good ecological environment in 1995 was the highest, reaching 43.67% of the total area. The area share of poor ecological environment was the lowest in 1995, comprising an area of only 371 km 2 , or less than 5% of the total area. The share of poor ecological environment was approximately 12% of the total area. The area with a good ecological environment rating in 2005 was larger than that of 1995 and accounted for the highest proportion (45.57%), comprising an area of 3494 km 2 . The percentage of area rated excellent was the smallest (6.23%) after 1995. In 2020, the poor ecological environment generally accounted for the highest   As governments and social organizations have become increasingly aware of en ronmental protection, Wuhan has strengthened its enforcement of ecologically relev laws and regulations, effectively halting the trend of environmental deterioration. Th reflected in the varying ecological evaluation grades. The differences in rating reflec increase in area from 643 km 2 in 2005 to 1684 km 2 in 2015 (area expansion of 14%) to 2 km 2 in 2020 (area expansion of 7%).
The spatial distribution (Table 6 and Figure 6) shows that areas with a good ecol cal environment are distributed mainly in the surrounding urban areas of Wuhan. Th areas have a relatively weak economy and the land-use types are mainly cultivated l and woodland, with rich vegetation and high biodiversity levels. The areas with poor e logical environments are concentrated in Hongshan, Hanyang, Wuchang and Qingsh   Ecological problems ascribed to human activities, such as vegetation damage and soil pollution, have become increasingly prominent.
As governments and social organizations have become increasingly aware of environmental protection, Wuhan has strengthened its enforcement of ecologically relevant laws and regulations, effectively halting the trend of environmental deterioration. This is reflected in the varying ecological evaluation grades. The differences in rating reflect an increase in area from 643 km 2 in 2005 to 1684 km 2 in 2015 (area expansion of 14%) to 2221 km 2 in 2020 (area expansion of 7%).
The spatial distribution (Table 6 and Figure 6) shows that areas with a good ecological environment are distributed mainly in the surrounding urban areas of Wuhan. These areas have a relatively weak economy and the land-use types are mainly cultivated land and woodland, with rich vegetation and high biodiversity levels. The areas with poor ecological environments are concentrated in Hongshan, Hanyang, Wuchang and Qingshan. According to the different functions of each administrative region of Wuhan, Hongshan is based mainly on the education industry. Several colleges and universities are located in the area, and it is densely populated. Qingshan, Hanyang and Wuchang are primarily industrial areas. Heavy industrial companies, such as Wuhan Iron & Steel Co., Ltd., Wushi Chemical Co., Ltd. and Dongfeng Motor Co., Ltd., are located in these areas. Industrial production and human economic activities have a direct detrimental effect on the environment of these areas.

Dynamic Monitoring of UEQ in Wuhan
Based on the IRSEI grade classification, the detected changes were divided further into nine levels and seven classes. The range for the levels of detected changes was −4 to +4, with a positive value indicating that the UEQ had improved, 0 indicating no change and a negative value indicating deterioration. For the classes with no detected changes, level 0 was classified as unchanged, level −4 as significantly worse and levels −2 and −3 as worse; level −1 as slightly worse; level 1 as slightly better; levels 2 and 3 as better; and level 4 as significantly better (Table 7).

Dynamic Monitoring of UEQ in Wuhan
Based on the IRSEI grade classification, the detected changes were divided further into nine levels and seven classes. The range for the levels of detected changes was −4 to +4, with a positive value indicating that the UEQ had improved, 0 indicating no change Remote Sens. 2021, 13, 4440 16 of 23 and a negative value indicating deterioration. For the classes with no detected changes, level 0 was classified as unchanged, level −4 as significantly worse and levels −2 and −3 as worse; level −1 as slightly worse; level 1 as slightly better; levels 2 and 3 as better; and level 4 as significantly better (Table 7).  Table 8 presents the ecological changes in Wuhan from 1995 to 2020. The size of the area representing both UEQ and ecological deterioration (obviously worse and slightly worse) is 3636 km 2 , accounting for the highest proportion (39.44%) over 2015-2020. The size of the area with the same UEQ (no change) is 2984 km 2 , accounting for 35.56% of the total area. Among the areas with deteriorating UEQ, most (69.51%) deteriorated by one grade. Deterioration in UEQ accounted for 25.60%. Most of the areas showing improved environmental conditions improved by one grade, accounting for 79.41% of the entire improved area. The areas improving by two grades account for 18%. The areas representing levels 3 or 4 are relatively small, indicating gradual changes. The areas with significant changes are related to direct economic activities, such as the transformation of cultivated land and woodland into construction and industrial land. The spatial distribution of UEQ (Figure 7) shows that the deteriorating areas are located mainly around cities and most water bodies. The deterioration of the ecological environment around water bodies is related to a leakage of urban domestic sewage and enterprise wastewater and a rise in aquaculture in recent years. Moreover, the areas with a deteriorating ecological environment are expanding along both sides of the Yangtze and Han rivers. Except for the water area, the UEQ in the central metropolitan area remains mainly unchanged and several areas show signs of improvement. This result indicates that environmental governance in the main urban area of Wuhan has played a positive role in recent years.

Spatial Autocorrelation Analysis
We explore the spatial autocorrelation (SA) of the IRSEI at a grid cell scale of 500 m × 500 m and our results indicate the existence of SA. The Moran's I was 0.568 in 1995, and 0.535 in 2020. All four IRSEI maps (1995,2005,2015 and 2020) display an extremely low

Spatial Autocorrelation Analysis
We explore the spatial autocorrelation (SA) of the IRSEI at a grid cell scale of 500 m × 500 m and our results indicate the existence of SA. The Moran's I was 0.568 in 1995, and 0.535 in 2020. All four IRSEI maps (1995,2005,2015 and 2020) display an extremely low probability (p-value < 0.01) of completely random spatial distribution. Therefore, the statistical significance test shows that SA exists for all of the ecological factors. The IRSEI increased in places where spatial distribution was favorable to the UEQ. In 1995, high-value clustering of the IRSEI in Wuhan was distributed mainly in the south and north of the study area, whereas low-value clustering was concentrated in the middle of the study area. In 2005, high IRSEI values started gathering gradually in the southern region, and low IRSEI values became more concentrated in the clustering distribution. By 2015, the high/high clustering and low/low clustering of the IRSEI in the study area became more dispersed and tended to spread in every direction. In 2020, low/low clusters had spread from the middle to the east and west, whereas high/high clusters were concentrated mainly in the south and north of Wuhan City.
The Moran's I scatter graph is divided into four quadrants, corresponding to four different spatial distribution types ( Figure 8). The first quadrant represents high/high clustering, the second quadrant low value and high-value aggregation, the third quadrant low/low aggregation and the fourth quadrant high-value and low-value aggregation. The IRSEI of Wuhan is concentrated mainly in the first and third quadrants. This result indicates that the IRSEI spatial distribution in Wuhan represents positive spatial autocorrelation, and high IRSEI agglomeration zones are mainly distributed in outer suburban areas, mainly in the north and southeast Wuhan.

Literature, Policy and Practice
We have reviewed previous studies and demonstrated that it is feasible to evaluate the quality of the urban ecological environment through remote sensing. This research proposes a feasible method. Other remote sensing images could also have been used as data in this research, such as Tiangong-2 WIS images [11]. In terms of method improvement, we mainly improved the integration of quantitative factors. A related similar index,

Literature, Policy and Practice
We have reviewed previous studies and demonstrated that it is feasible to evaluate the quality of the urban ecological environment through remote sensing. This research proposes a feasible method. Other remote sensing images could also have been used as data in this research, such as Tiangong-2 WIS images [11]. In terms of method improvement, we mainly improved the integration of quantitative factors. A related similar index, RSUSEI, has primarily increased remote sensing ecological factors by adding the impervious surface cover (ISC) [15]. ISC is also one of the most important factors that distinguish different types of land use/land cover characteristics in urban environments, and has a strong impact on UEQ. However, in our study we also consider the dryness index (NDBSI), the bare soil index (SI), the building index (IBI) and the normalized buildings-bare-soil index. However, there are strong correlations between the impervious surface, bare soil and building indices. Previous studies have found that the relationship between ISC and LST has the form of an exponential function, rather than a simple linear function, as commonly believed [43]. This exponential relationship has been confirmed by many subsequent studies [44,45]. Our IRSEI index takes into account the bare soil, building index and surface temperature. We suggest that the correlations of remote sensing ecological indicators affecting the regional ecological environment should be introduced into comprehensive indicators, or different indicators should be set according to the characteristics of the study region.
There are few high-quality ecological environment patches in Wuhan (IRSEI > 0.8), with close to zero over the past five years, and most of the patches are in the center of the ecological environment. Therefore, we propose a policy whereby Wuhan would focus on protecting forest land and gardens, build high-quality ecological corridors and coordinate the management of rivers in the future, so as to guide sustainable urban development and achieve sustainable development goals (such as SDG 11, sustainable cities and communities). Lake and wetland protection and ecological restoration and management will optimize the pattern of ecological security. Further analyses of the results indicate that there was a negative correlation between LSI, NDBSI and urban ecological quality. The ecological environment in areas with a high surface temperature, such as the Wuhan downtown area and coastal area around the Yangtze River, has tended to deteriorate; however, the humidity indices in these areas were also relatively high, which is conducive to ecological protection. Low vegetation index values in the central urban area also affect the quality of the ecological environment of Wuhan to a certain extent. The IRSEI can macro-evaluate the quality of the regional ecological environment, which is more convenient and efficient. In the future, higher precision can be introduced at the block level. Data, such as Google Street View data, could be used with machine learning algorithms to further identify the proportion of regional urban green space, trees, etc., and improve the accuracy of ecological environment assessment. The index has a high ability to distinguish between different land cover uses. The framework can also be easily extended to a global scale or to map other gridded socio-economic variables (such as GDP and population) to monitor and assess progress towards the SDGs [25]. The assessment and modelling of uses is critical to supporting sustainability assessment in achieving Sustainable Development Goals (SDGs), such as sustainable cities and communities. Therefore, IRSEI can be used to assess the spatial and temporal sustainability of cities.

Analysis of the Factors Affecting the UEQ
The regression least squares method (OLS) can be used to quantitatively describe the relationship between the ecological index and natural, economic and social factors in Wuhan. The data include temperature, precipitation, elevation, slope and DMSP as explanatory variables. The night light variable reflects the human footprint and fundamentally affects the urban ecological environment. Impervious surfaces and roads and a high population density are not conducive to UEQ. The regression coefficients represent the contribution of six independent variables to the dependent variable. The regression coefficients of precipitation and elevation are equal to 0.522 and 0.441, respectively, indicating that precipitation and elevation positively contribute to the IRSEI.
In contrast, the regression coefficients of night light and slope are negative, indicating that these variables contribute negatively to the IRSEI. The night light variable has a regression coefficient of -0.619, indicating a negative effect. The R 2 is 0.901 and p < 0.05, indicating that climate, precipitation, elevation, slope and night light data account for 90% of the variations of the IRSEI.

Method Framework and Validation Analysis
Weighting is an important process in the development of aggregated ecological indices that help promote sustainability. Different weighting methods have different characteristics, and the method employed could reflect the subjectivity of the decision makers. However, such methods combined with remote sensing index data can facilitate decisions and reduce the calculations required.
PCA is widely used in the evaluation of the RSEI. In several studies, the three principal components obtained after dimensionality reduction did not show any obvious effects (contribution was below 15%). However, including all of the pixels in extensive data calculations is a time-consuming process. The RSEI employs a covariance-based (unstandardized) PCA to determine the importance of each indicator involved. The weight of each indicator can be assigned objectively and automatically based on the load (contribution) of each indicator to PC1. In this study, we used PCA and EW to comprehensively calculate the IRSEI. After improvement, the combined method was able to reflect the degree of change in the index, and the calculation was quick and uncomplicated. The spatial distribution of the UEQ over the study period (1995-2020) is consistent with the information in the bulletin on the eco-environmental situation in China in that year. The current, more popular assessment method is based on habitat quality (HQ) [46][47][48][49][50]. In further research, we intend to include HQ in this quantitative assessment.

Limitations and Future Prospects
The proposed UEQ evaluation model is feasible and straightforward, providing a new idea for ecological protection and comprehensively reflecting the changes in UEQ in Wuhan. From 1995 to 2020, the UEQ of Wuhan declined overall, probably owing to a combination of natural factors and human activities. However, the ecological level in the eastern and southeastern mountainous areas has increased because of the influence of forest resource protection, desertification land management and the warm and humid climate. In contrast, the regional ecological level has declined, owing to the overexploitation and overgrazing of lake resources in the northwest and southwest of Wuhan. The constantly rising levels of urbanization and construction over nearly 20 years have resulted in a downward trend in the UEQ. Overall, the ecology of Wuhan is in a fragile state. In 2020, the proportion of areas with poor ecological environment grades remained high, accounting for 42.68% of the total area.
In future social and economic development, we should follow the laws of nature, prioritize protection and rationally develop and utilize natural resources. The IRSEI effectively revealed the spatial distribution of and change in the UEQ in Wuhan, based on remote sensing images. Although four types of ecological factors closely related to the ecological environment were selected in the calculation process, the ecological environment is a complex and comprehensive variable. Areas with a deteriorating ecological environment tend to be spread along the Yangtze and Han rivers and around the central urban area. Urban expansion has damaged the ecological environment, and urban planning should integrate more ecological concepts to promote a harmonious coexistence and sustainable development for humans, nature and society.
Comprehensive quantitative evaluation requires selecting several impact factors that reflect the actual situation in the study area. We aimed to conduct the UEQ evaluation by employing a scientific, objective and feasible method. Nevertheless, choosing the UEQ evaluation index remains exploratory work. Determining the index weight affects the accuracy of the evaluation results. Accordingly, expanding research to a more scientific multifunction performance index system and determining the index weights require further work. Furthermore, the limited availability of data and a lack of longitudinal comparison of urban data have affected the scientific nature of our research results. In addition, our next step will be exploring how the IRSEI changes at different spatial scales. The rapid development of cities will inevitably lead to a series of ecological and environmental problems, and the deterioration of the ecological environment may further affect the surrounding environment, forming a cycle and harming urban sustainability. This study also demonstrates that IRSEI is characterized by spatial heterogeneity; that is, the poor UEQ patches will focus on areas where the ecological environment is poor and the urbanization is also highest.

Conclusions
In this study, the IRSEI model was used to evaluate and monitor the ecological environment in Wuhan from 1995 to 2020. The IRSEI is an ecological environmental quality assessment method based on remote sensing technology. The method has many advantages, such as the ease of obtaining parameters, a large time sequence span and a wide evaluation range. The UEQ method employing remote sensing technology is feasible and simple, and provides a new tool for territorial spatial control and spatiotemporal urban sustainable development. Our proposed UEQ assessment framework can also help to develop potentially relevant additional sub-indicators, which could help to address one of the current challenges in SDG monitoring, namely how to implement SDG indicators. We have implemented the proposed workflow in this study based on an open-source platform and free satellite data, making it an appealing option that is applicable in almost all countries.
The main conclusions from the results of this study are: • The mean IRSEI value in Wuhan decreased annually from 2005 to 2015. The UEQ continued to decline, mainly because of the rapid economic development of Wuhan, reduction in vegetation coverage caused by human activities, gradual decrease in lake area and transformation of the land-use structure caused by urban expansion; • From the perspective of spatial patterns, the UEQ of the central urban areas, such as Qingshan, Hanyang, Hongshan and Wuchang, was lower than that of the surrounding metropolitan areas, such as Huangpi, Jiangxia and Caidian. In terms of time series, the UEQ in the central city of Wuhan has been mainly unchanged or improved, indicating that the management of the ecological environment in the central city had achieved specific results; • The global Moran's I value range from 0.535 to 0.592 from 1995 to 2020, respectively, indicating that the IRSEI spatial distribution displays significant spatial heterogeneity. This finding indicates that high clustering gradually developed to the edge of the city, whereas low clustering gradually developed to the center of the city. The spatial correlation and local index cluster diagram of the IRSEI show that the high points are located mainly in the Huangpi and Jiangxia districts; • The UEQ evaluation model constructed in this study is feasible and simple and could be implemented at no or negligible cost, making it applicable to most regional areas. Our model, therefore, could be considered a new tool for ecological management and protection, and for assessing progress toward urban sustainable development.

Conflicts of Interest:
The authors declare no conflict of interest.