Remote Sensing-Based LULP Change and Its Effect on Ecological Quality in the Context of the Hainan Free Trade Port Plan

: The study of Land Use and Landscape Patterns (LULPs) changes and their ecological quality effects in Haikou city under the background of the Hainan Free Trade Port Plan (HFTPP) helps to promote coordinated development between cities and the environment. To date, most research on ecological quality has focused on areas with extremely fragile ecology and/or is related to LULP analysis. There are few studies in the literature focusing on the impact of high-intensity human activities caused by relevant policies on urban LULPs. The purpose of this research was to design a framework that monitors urban ecological security by considering the effect of the developing free trade port. The proposed framework was constructed by integrating multi-temporal Sentinel-2 remote sensing images, night light remote sensing data, digital elevation model (DEM) data, and spectral index features such as the normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), bare soil index (BSI), and normalized intertidal mangrove index (NIMI), as well as analytical approaches such as the land use transfer matrix, land use dynamic degree, land use degree and transfer matrix, land use gravity center measurement, and landscape pattern index. The framework takes advantage of the Google Earth Engine (GEE) cloud platform and was applied to a highly developed Haikou city, the capital of Hainan province. Maps of brightness (SBI), greenness (GVI), and humidity (WET) were created annually from 2016 to 2021, enabling detailed ecological environment quality evaluation and analysis. The advantages of this study are (1) reliable land cover results obtained automatically and quickly; (2) the strong objectivity of the quantitative research on landscape patterns and land use; and (3) deep integration with free trade port policies. Through the research on the ecological quality problems caused by the change in LULP in the study area, the research results show that, from 2016 to 2021, the spatial distribution of land use and landscape pattern in Haikou city had been constantly changing; the area of construction land has decreased, with most of it having been converted into forest land and farmland; the gravity center of the building land has moved to the northwest; the degree of landscape fragmentation has decreased and the heterogeneity of landscape distribution has increased; the free trade port policies have promoted Haikou’s economic development and ecological civilization construction; and finally, Haikou’s ecological environmental quality has improved significantly.


Introduction
The quality of the ecological environment is an important basis for the sustainable development of human society and an important part of national security [1].It is also an important guarantee for economic security, political security, and social stability.The evaluation of the ecological environment quality combined with remote sensing technology and geographic information technology shows the unique characteristics and advantages of interdisciplinary synthesis.Ecological quality evaluation is an important basis for China to promote the construction of ecological civilization and achieve sustainable development, and it is also a hot issue that has been widely considered by academic and political circles in recent years.Haikou, also known as "Coconut City", is the core city of Hainan Free Trade Port.On 13 April 2018, General Secretary Xi Jinping delivered an important speech at the conference to celebrate the 30th anniversary of Hainan's construction of a special economic zone.He supported Hainan in gradually exploring and steadily promoting the construction of a free trade port with Chinese characteristics, and he called for the establishment of a free trade port policy and institutional system step by step [2].
Google Earth Engine (GEE) offers several advantages over other platforms for geospatial analysis and remote sensing.The advantages of GEE can be summarized as follows: (1) GEE hosts a vast collection of satellite imagery and geospatial datasets, including Landsat, Sentinel, MODIS, and more.This extensive archive allows users to access petabytes of data without worrying about storage or data management (https: //developers.google.com/earth-engine/datasets(accessed on 3 May 2021)).(2) The Earth Engine is built on Google's cloud infrastructure, which enables scalable and parallelized processing of large-scale geospatial data.Users can analyze terabytes of data in a fraction of the time it would take using traditional computing resources.Earth Engine also provides a wide range of analysis tools and functions for image processing, time series analysis, machine learning, and statistical modeling.These tools are integrated into the platform, making it easier for users to perform complex analyses without switching between different software packages (https://developers.google.com/earth-engine/(accessed on 3 July 2021)).( 3) Earth Engine offers a Code Editor with a JavaScript-based integrated development environment (IDE) that allows users to write, test, and execute code directly in the browser.Additionally, Earth Engine provides APIs for Python and JavaScript, enabling seamless integration with existing workflows and programming languages (https://developers.google.com/earth-engine/playground(accessed on 3 July 2021)).( 4) Earth Engine facilitates collaboration and sharing of geospatial data, algorithms, and analysis results through its cloud-based platform.Users can share their code, scripts, and visualizations with others, fostering a collaborative and open-source approach to geospatial research (https://groups.google.com/g/google-earth-engine-users(accessed on 14 November 2021)).(5) In my opinion, the most import advantage is that GEE provides free access to its platform for non-commercial and educational use.This lowers the barrier to entry for researchers, educators, and students who want to explore and analyze geospatial data for various applications, including environmental monitoring, land cover mapping, disaster response, and more (https://earthengine.google.com/terms/(accessed on 17 November 2021)).These advantages make GEE a powerful platform for geospatial analysis and remote sensing applications, suitable for a wide range of users from researchers and scientists to educators and students [3].Users can use Java Script and Python programming languages to call upon various images of the platform and develop classification algorithms for remote sensing applications.GEE can support portable and interactive applications, providing great convenience for remote sensing workers [4].Therefore, in recent years, the GEE cloud platform has been widely used in the field of ecology, based on different scales of the world [5], country [6,7], and region [8,9], involving cities [10], cultivated land [11], watersheds [12,13], and other aspects of ecology environmental studies.
The change in land use and landscape pattern is the result of the interaction between man, nature, and other factors, and it is also an important carrier of economic and social activities [14].The study of LULP evolution plays a vital role in regional development.At present, research on the change in LULPs is mainly concentrated in ecologically fragile areas [2,15], focusing on long time-series studies [16,17].However, few researchers have considered the impact of high-intensity human activities and related policies on LULPs in combination with the conditions of the study area and the general environment of the study area, especially for Haikou City, where the urbanization process has developed rapidly in respect to the free trade port construction in recent years.
In this context, it is of great significance to clarify the impact of the free trade port construction on LULP changes and ecological quality evaluation in Haikou.Therefore, this study took Haikou City as an example, used the GEE cloud platform, selected Sentinel-2 data from 2016 to 2021, and then constructed a feature set to obtain the land cover results of the study area.Then, this study used the land use transfer matrix, land use dynamic degree, land use degree and transfer matrix, land use gravity center measurement, and landscape pattern index to carry out land use and landscape pattern change analysis.Finally, we used indicators such as "brightness (SBI)", "greenness (GVI)", and "humidity (WET)" to evaluate the ecological health of Haikou City and analyze the ecological impact of the free trade port construction on Haikou City.Research outcomes can significantly enhance decision-making for governmental actors striving for more sustainable cities.By utilizing historical satellite imagery and geospatial data from 2016 to 2021 in this research, governmental actors can understand urban growth patterns, identify areas of rapid development, and predict future expansion.In addition, the research overcomes the spatial and temporal limitations of traditional methods and can better reflect the landscape pattern information of surface cover, provide important reference for ecological environmental protection and restoration, and further promote the application and development of remote sensing technology in urban planning and environmental management.

Study Area
Haikou City (N 19 ) is located in the north of Hainan province, and has a tropical marine monsoon climate.The topography of Haikou is slightly long and heart-shaped, and it is dominated by plains and terraced fields, with rich wetlands and numerous rivers.Haikou has jurisdiction over the districts Xiuying, Longhua, Qiongshan, and Meilan.Haikou is the political, economic, technological, and cultural center of Hainan Province, and it is also the largest transportation hub [18].The geographical location of the study area is shown in Figure 1.
present, research on the change in LULPs is mainly concentrated in ecologically fragile areas [2,15], focusing on long time-series studies [16,17].However, few researchers have considered the impact of high-intensity human activities and related policies on LULPs in combination with the conditions of the study area and the general environment of the study area, especially for Haikou City, where the urbanization process has developed rapidly in respect to the free trade port construction in recent years.
In this context, it is of great significance to clarify the impact of the free trade port construction on LULP changes and ecological quality evaluation in Haikou.Therefore, this study took Haikou City as an example, used the GEE cloud platform, selected Sentinel-2 data from 2016 to 2021, and then constructed a feature set to obtain the land cover results of the study area.Then, this study used the land use transfer matrix, land use dynamic degree, land use degree and transfer matrix, land use gravity center measurement, and landscape pa ern index to carry out land use and landscape pa ern change analysis.Finally, we used indicators such as "brightness (SBI)", "greenness (GVI)", and "humidity (WET)" to evaluate the ecological health of Haikou City and analyze the ecological impact of the free trade port construction on Haikou City.Research outcomes can significantly enhance decision-making for governmental actors striving for more sustainable cities.By utilizing historical satellite imagery and geospatial data from 2016 to 2021 in this research, governmental actors can understand urban growth pa erns, identify areas of rapid development, and predict future expansion.In addition, the research overcomes the spatial and temporal limitations of traditional methods and can be er reflect the landscape pa ern information of surface cover, provide important reference for ecological environmental protection and restoration, and further promote the application and development of remote sensing technology in urban planning and environmental management.

Study Area
Haikou City (N 19°31′32″-20°04′52″, E 110°07′22″-110°42′32″) is located in the north of Hainan province, and has a tropical marine monsoon climate.The topography of Haikou is slightly long and heart-shaped, and it is dominated by plains and terraced fields, with rich wetlands and numerous rivers.Haikou has jurisdiction over the districts Xiuying, Longhua, Qiongshan, and Meilan.Haikou is the political, economic, technological, and cultural center of Hainan Province, and it is also the largest transportation hub [18].The geographical location of the study area is shown in Figure 1.

Data Sources
The remote sensing data used in this study mainly include Sentinel-2 satellite data and night light data (VIIRS).These data were acquired through the GEE platform combined with the research area (Haikou City) and time information (2016)(2017)(2018)(2019)(2020)(2021).The digital elevation model (DEM) data and slope data of Haikou City come from the Resource and Environment Science and Data Center (https://www.resdc.cn/(accessed on 18 January 2022); the administrative division data come from the National Geomatics Center of China (http: //www.ngcc.cn/ngcc/(accessed on 18 January 2022); and Haikou Statistical Yearbooks come from The People's Government of Hainan Province (https://www.hainan.gov.cn/(accessed on 27 January 2022).

Data Preprocessing
The remote sensing data were Sentinel-2 high-resolution multi-spectral images, which are generated by two complementary satellites, with a revisit period of 5 days and an imaging resolution of 10-60 m [19].The Sentinel-2 images used in this study are products stored in the GEE platform after atmospheric correction and orthorectification.Because the study area has a weather feature with a lot of cloud cover, in order to obtain high-quality image data, we carried out cloud processing.Based on the official method provided by GEE, the study used the QA band to quickly perform cloud remove processing on the Sentinel-2 images.When acquiring Sentinel-2 images of the study area, we limited the percentage of cloud cover (<20%) and chose the median method to synthesize the image data with the minimum cloud cover.The preprocessing of the remote sensing images was completed through steps such as mosaic and cropping.Finally, the basic remote sensing image dataset was obtained.
According to the classification standards in the "Classification of Land Use Status" (GB 21010-2017) [20] and the actual land use/cover classification status of Haikou City, the classification system was divided into eight land use types: grassland, wetland, water, forest land, building land, farmland, bare land, and mangrove.For the high-resolution remote sensing image data in the GEE platform, the collection strategy of "small area continuous collection, global uniform distribution" was used to visually select samples.

Research Route
Firstly, with the help of the GEE cloud platform, Sentinel-2 images from 2016 to 2021 in the study area were screened out and clouds were removed, feature sets were constructed, and random forest classification was performed; accuracy evaluation and filtering processing were performed on the classification results.Then, with the help of ArcGIS 10.7 software, based on the land cover classification results of Haikou City, land use change analysis was carried out, with methods including transfer matrix, land use dynamic degree, land use degree and transfer matrix, and land use gravity center measurement.Using Fragstats 4.2 software, a series of landscape pattern indices were selected and analyzed.Finally, the greenness, brightness, and humidity indicators were extracted on the GEE cloud platform.Normalized with the help of ENVI 5.3 software, the ecological quality of Haikou City was analyzed.The research results provide a reference for the construction of the Haikou Free Trade Port. Figure 2 shows the overall methodology for this study.

Construction of Feature Set
The selection of feature variables is an important step in remote sensing classification.Using the powerful cloud computing capabilities of the Google Earth Engine platform, for the selected Sentinel-2 remote sensing image dataset, comprehensively considering the characteristics of each band, this study selected the reflectance data of the blue band, green band, red band, near-infrared band, and short-wave infrared band to construct spectral band features, which were used to analyze the spectral difference characteristics of land cover types in Haikou City (Table 1).Using a variety of feature variables and feature variable combinations can effectively improve the accuracy of remote sensing classification.We highlighted ground features by adding spectral index features such as the normalized difference vegetation index (NDVI) [21], enhanced vegetation index (EVI) [22], and bare soil index (BSI) [23]; importantly, we introduced the normalized intertidal mangrove index (NIMI) [24] to extract mangroves.The introduction of night light data can effectively improve the determination of construction land.Topographic features are important features in remote sensing image classification.The elevation information extracted based on DEM data was used as topographic features, which can effectively assist in judging land cover type information at different altitudes.

Land Use Information Extraction 4.2.1. Construction of Feature Set
The selection of feature variables is an important step in remote sensing classification.Using the powerful cloud computing capabilities of the Google Earth Engine platform, for the selected Sentinel-2 remote sensing image dataset, comprehensively considering the characteristics of each band, this study selected the reflectance data of the blue band, green band, red band, near-infrared band, and short-wave infrared band to construct spectral band features, which were used to analyze the spectral difference characteristics of land cover types in Haikou City (Table 1).Using a variety of feature variables and feature variable combinations can effectively improve the accuracy of remote sensing classification.We highlighted ground features by adding spectral index features such as the normalized difference vegetation index (NDVI) [21], enhanced vegetation index (EVI) [22], and bare soil index (BSI) [23]; importantly, we introduced the normalized intertidal mangrove index (NIMI) [24] to extract mangroves.The introduction of night light data can effectively improve the determination of construction land.Topographic features are important features in remote sensing image classification.The elevation information extracted based on DEM data was used as topographic features, which can effectively assist in judging land cover type information at different altitudes.

Classification and Filtering
In addition to the data source, the choice in classification algorithm also directly affects the classification results.In recent years, with the development of machine learning classification algorithms, the random forest algorithm (RF) has high calculation accuracy and less model training time.RF is able to determine the relative importance of variables in the model with low sensitivity to the number and quality of training samples [25].Therefore, RF is widely used in land cover identification.In this study, the random forest algorithm was used for land use classification in Haikou City; 70% of the sample points were selected as training samples for training the classifier, and 30% of the sample points were used as validation samples for accuracy verification.
Because the topography of the study area breaks up the land use units, the classification results obtained have a "salt and pepper phenomenon", so the classification results need to be post-processed.In this paper, the method of mode filtering was used to remove, merge, and cluster the fragmented image patches after classification.We used the nearest eight-pixel neighborhood as the filter mode, that is, when at least half of the eight pixels are the same, the corresponding pixel value can be changed [23].This ensures that the spatial structure unit in the classification result is relatively complete.

Accuracy Evaluation
The accuracy of land use classification in Haikou City was verified by the confusion matrix and sample points collected from high-resolution images.The evaluation incorporated several precision metrics, including producer's accuracy (PA), user's accuracy (UA), overall accuracy (OA), and the Kappa coefficient (Kappa).The confusion matrix provides a clear visual representation of the classification results, but it does not readily reveal the accuracy of individual category classifications.With the help of these multiple accuracy metrics, we can gain a comprehensive understanding of the strengths and limitations of our land use classification model in Haikou City.Among these, OA and the Kappa coefficient (Kappa) are particularly noteworthy.Overall accuracy (OA) refers to the ratio of the number of pixels, which is the correctly classified category to the total category, and the Kappa coefficient (Kappa) is an indicator used to measure the consistency of the classification result of the verification sample [26].

Overall accuracy =
Number o f correct plots(value) Total number o f plots(value) where i is the class number, n is the total number of classified pixels that are being compared to actual data, nii is the number of pixels belonging to the actual data class i which were classified with a class i, C i is the total number of classified pixels belonging to class i, and G i is the total number of actual data pixels belonging to class i.

Land Use Change 4.3.1. Land Use Transition Matrix
The land use transfer matrix can quantitatively reflect the structure and transfer changes in land types in different periods.Its mathematical expression is shown in (3).Based on ArcGIS 10.7 software, using the Dissolve Tool and the Intersect Tool, combined with Excel 2016 software, a total of six periods of land use transfer data from 2016 to 2021 were generated.These data were used to analyze the characteristics of land use quantity and direction in Haikou City.
In this formula, S ij is the area of land type i converted into land type j, and n is the number of land use types.

Land Use Dynamic Degree
(1) The Single Land Use Dynamic Degree.
The dynamic degree of a single land use type expresses the quantity changes in a certain land use type within a certain time range in a certain research area.In this paper, it is expressed as the annual change rate [27], and its expression formula is (4).
In the formula, K is the annual change rate of a certain type of land use during the study period; U a and U b are the areas of a certain type of land use at the beginning and end of the study period, respectively; and T is the length of the study period.
(2) The Integrated Land Use Dynamic Degree The integrated dynamic index of land use is the transfer rate among land categories during the research period.The calculation formula of the integrated land use dynamic degree in Haikou [28] is shown in (5).
In this equation, U bi and U ai are the areas of a certain land use category at the end and beginning moment of the research, respectively, T is the duration of the study, and n is the quantity of land use categories and is equal to eight in this study.When T is year(s), k total is the annual integrated change rate of area change in all land use categories.

Land Use Degree
The land use degree reflects the breadth and depth of land use.It also reflects the comprehensive effects of human activities and the natural environment on land use [29].According to the comprehensive analysis method of land use degree proposed by Zhuang Dafang et al. [30], this paper divides the land use degree into four grades according to the land natural complex in a state of natural balance.The assigned grading indices are shown in Table 2. (1) Land Use Degree Index The land use degree index model [31] is (6): In the formula, L j is the index of land use degree; A i is the grading index of the i-level land use degree; C i is the area percentage of the i-level land use degree; and n is the grade of the land use degree.The change interval of the land use degree index is [100, 400]; the higher the land use degree, the larger the index.
(2) Rate of land use degree The rate of the land use degree model [32] is formula (7): In the formula, R is the change rate of the land use degree; A i is the grading index of the i-level land use degree; and C ib and C ia are the percentages of the area occupied by the i-level land use degree at time b and time a.If R is positive, the land use in this area is in the development period, and if R is negative, it means that it is in the adjustment period or the decline period.

Land Use Gravity Center Measurement
Land use gravity center, that is, the geometric balance point of a land use type, is an important indicator to measure the overall distribution of a certain land use type in a region [33].By comparing the gravity center distribution of land use types at the beginning and end of the research period, we can obtain the characteristics of land use change and distribution in a certain period of time.This index mainly reflects the concentrated distribution of land use types in space.The calculation formulas are shown in Formula (8) and Formula (9): In the formula, X t and Y t are the gravity center coordinates of a certain land use type in year t; n is the number of patches; C n is the area of a certain land use type; X i and Y i are the barycentric longitude and latitude of the i-th patch, respectively; and M it is the area of the i-th patch in year t.

Landscape Pattern Index
The landscape pattern index is a quantitative index that reflects the spatial layout and structural composition characteristics of landscape elements.It can highly reflect the landscape pattern information [34].The landscape pattern index is divided into three levels, which are patch level index, class level index, and landscape level index [35].Referring to the previous research, at the patch-type level, this paper selected the number of patches index (NP), the mean patch area index (AREA_MN), the largest patch index (LPI), and the landscape shape index (LSI) to analyze the changes in landscape pattern.At the landscape level, indices such as the mean patch area index (AREA_MN), the largest patch index (LPI), the landscape shape index (LSI), the patch density index (PD), Shannon's diversity index (SHDI), and Shannon's evenness index (SHEI) were selected to analyze the landscape pattern changes.Fragstats 4.2 software provides a large number of spatial metrics, and a specific subset of them were specifically selected for this study; they are described in Table 3 [36,37].

AREA_MN Mean patch area
This indicates the degree of landscape fragmentation.The smaller the value, the more fragmented the landscape is.

LPI
Largest patch index 0 < LPI ≤ 100.LPI approaches 0 when the largest patch in the landscape is increasingly small.LPI = 100 when the entire landscape consists of a single patch; that is, when the largest patch contains 100% of the landscape.Its change can reflect the change in the intensity of human disturbance activities.

Ecological Quality Assessment
This study conducted rapid monitoring and evaluation of the regional ecological environment based on remote sensing information.On the basis of the obtained tasseled cap transformation results, the research used the GEE platform to synthesize the humidity, brightness, and greenness to construct and describe the spatial distribution of the ecological environment quality.
(1) Tasseled Cap Transformation Tasseled cap transformation is an orthogonal transformation of images based on the distribution structure of soil, vegetation, and other information in multidimensional spectral space [38].The transformed main components are "Soil Brightness Index (SBI)", "Vegetation Greenness Index (GVI)", and "Wetness (WET)", which, respectively, reflect information such as soil rocks, vegetation, and moisture in soil and vegetation.The information represented by these components can be applied to ecological environment monitoring.The tasseled cap transformation depends on the band settings and characteristics of the sensor itself, and the transformation coefficients of different sensors are not universal.In this study, referring to previous research results [39], the corresponding transformation coefficients were set to calculate the SBI, GVI, WET, and other components of Sentinel-2 data.The transformation matrix formula is shown in (10).
In the formula, B 2 , B 3 , B 4 , B 8 , B 11 , and B 12 are the band values corresponding to Sentinel-2.
(2) Normalization Processing According to the correlation between each index value and the ecological environment, the corresponding method was used to normalize the initial calculation results [40].For positive correlation indicators, the normalization formula is (11).For negative correlation indicators, the normalization formula is (12).
In the formula, y i is the index score after normalization, and x i , x min , and x max are the initial value, minimum value, and maximum value of the indicator, respectively.

The Spatiotemporal Changes in Land Use in Haikou City
The spatiotemporal distribution of land use in Haikou City was evaluated utilizing the GEE platform and the random forest classification algorithm, followed by spatial connectivity processing with the filter mode.The resulting map, after the removal of small patches, is presented in Figure 3.The accuracy evaluation of the land use data of Haikou City in each period was carried out, and the results are shown in Table 4.We found that the overall accuracy (OA) and the Kappa coefficient of land use in each period were greater than 0.85, meeting the relevant accuracy requirements.By showing the classified statistics of the data on land use types in Haikou City from 2016 to 2021 (see Table 5), this table can succinctly express the area changes in different years.The results show that there are obvious regional differences in the distribution of land use types in Haikou City.The user accuracy (UA) of building land from 2016 to 2021 was 0.92, 1, 1, 0.85, 1, 0.77, and 0.97, respectively, while the production accuracy (PA) was 0.79, 0.95, 0.94, 1, 1, and 1, respectively, all maintained at a high level.These results show that our classification algorithm can accurately distinguish building land from other land types in most cases.On the whole, forest land, building land, and farmland are the main types, among which forest land occupies the largest area.From the perspective of spatial distribution, forest land is widely distributed, mainly in the central, western, and southeastern regions of Haikou, except for buildings; mangroves are mainly distributed in the northeastern coastal area; and building land is mainly distributed in the northern part of Haikou.The water network of rivers and lakes in Haikou City is densely distributed, and the water area is mainly located in the Nandu River flow area, which runs from south to north.The proportion of bare land and wetland is relatively small, indicating that Haikou has a high degree of land development and utilization, and the reserve land resources are not sufficient.

Analysis of Land Use Transfer Matrix
We used the ARCGIS software to export the attribute table, and then obtained the land use transfer matrix table in Excel.The table clearly shows the characteristics and rules of the internal conversion among land use types in Haikou City.In the conversion of land use types, grassland is the most converted, and it is mainly converted into forest land; the converted grassland area amounts to 52.755 km 2 .The second most converted is building land, whose converted area is 43.885 km 2 .Wetland is mainly converted into farmland; water, bare land, and mangroves are mainly converted into building land; building land is mainly converted into farmland and forest land; and farmland is mainly converted into forest land.The land use transfer matrix from 2016 to 2021 is shown in Table 5.

Analysis of Land Use Dynamic Degree
It can be seen from the land use dynamic degree results (Table 6) that the integrated land use dynamic degree in the five periods shows a trend of decreasing first and then increasing.From 2016 to 2017, the integrated land use dynamic degree was the highest, with a value of 13.789%, mainly due to the increase in farmland.At this time, the single land use dynamic degree was 90.482%.The integrated land use dynamic degree during 2018-2019 was 5.727%, which was the lowest value among the five periods.This change was mainly caused by the increase in farmland area.During 2017-2018, the area of farmland decreased the most, and its single land use dynamic degree was −18.215%.During the period of 2019-2020, the mangroves changed greatly, and their single land use dynamic degree was 290.424%.The type of surface cover with the largest increase in area during 2020-2021 was forest land, and its single dynamic degree was 27.633%.The research divided the land use types into four grades, and Figure 4 shows a map of the land use degree changes in the six phases in this area.From the analysis of the figure, we can conclude that the degree of land use presents a trend of first decreasing, then increasing, and then decreasing.The change interval of the land use degree index is [100, 400]; the higher the land use degree, the greater the index value.It can be seen from the land use degree change map of Haikou City that the degree dropped from 280.652 in 2016 to 271.682 in 2019, then increased to 279.050 in 2020, and then slightly decreased to 273.646 in 2021.On the whole, the degree of land use showed a decline, indicating that land use has been in a stage of constant adjustment during this period.Haikou is gradually figuring out its own development mode, which requires us to pay attention to the construction of ecological civilization, constantly adjust the way the land is used, and finally, improve the efficiency of the land use.The decrease in the land use degree index also indirectly proves that Haikou's ecological environment protection has achieved certain results.
The calculated results of the land use degree change rate in Haikou City in the five stages are shown in Table 6.From this table, we can see that the land use degree change rate in the five periods has positive and negative values.The land use degree change rate in 2016-2020 and 2020-2021 are both negative, and in these two periods, Haikou belonged to the adjustment stage, indicating that the degree of land use by human activities was low, the land load in these periods was small, and the natural state was good.From 2019 to 2020, the land use degree change rate in Haikou City was 2.712%, which is a positive value, indicating that Haikou City was in the development period during this time, meaning that the degree of land use by human activities was increasing.At this stage, the output value of land was constantly changing, and Haikou was gradually entering into a stage of adjustment and development.
to the adjustment stage, indicating that the degree of land use by human activities was low, the land load in these periods was small, and the natural state was good.From 2019 to 2020, the land use degree change rate in Haikou City was 2.712%, which is a positive value, indicating that Haikou City was in the development period during this time, meaning that the degree of land use by human activities was increasing.At this stage, the output value of land was constantly changing, and Haikou was gradually entering into a stage of adjustment and development.

Analysis of the Gravity Center Shift of Building Land in Haikou
In ArcGIS, the average center of building land distribution was calculated for the layers of the past years, so as to represent the gravity center of the building land area.This made it easier to further understand the spatial distribution of the building land.It can be seen from the gravity center transfer map (Figure 5) that the gravity center of the building land during the research period was entirely in Qiongshan District.The gravity center of the building land moved from the southeast to the northwest of Haikou as a whole.During the three stages of 2016-2017, 2018-2019, and 2020-2021, the gravity center of building land moved to the northwest, while the gravity center of building land during 2017-2018 and 2019-2020 moved to the southeast.In ArcGIS, the average center of building land distribution was calculated for the layers of the past years, so as to represent the gravity center of the building land area.This made it easier to further understand the spatial distribution of the building land.It can be seen from the gravity center transfer map (Figure 5) that the gravity center of the building land during the research period was entirely in Qiongshan District.The gravity center of the building land moved from the southeast to the northwest of Haikou as a whole.During the three stages of 2016-2017, 2018-2019, and 2020-2021, the gravity center of building land moved to the northwest, while the gravity center of building land during 2017-2018 and 2019-2020 moved to the southeast.

Analysis of Landscape Pa ern Changes at Landscape Level
The results of the landscape pa ern index of Haikou City at the landscape level from 2016 to 2021 are shown in Figure 6.It can be seen from the figure that the landscape pa ern of Haikou City from 2016 to 2021 generally shows the characteristics of gradually increasing landscape fragmentation, simplifying the shape of the landscape, and increasing landscape heterogeneity.From 2016 to 2021, the mean patch area and patch density in Haikou City have alternating changes of increase and decrease.The mean patch area reached its maximum in 2019, with an area of 0.93 hm 2 .Overall, the mean patch area increased from 0.612 hm 2 in 2016 to 0.698 hm 2 in 2021.The patch density reached 1556.857pieces/100 hm 2 in 2020.These show that under the influence of the external environment, Haikou's land-

Analysis of Landscape Pa ern Changes at Class Level
From 2016 to 2021, the landscape pa ern of Haikou at a class level underwent great changes, and the results are shown in Figure 7. From Figure 7a, it can be seen that the number of patches of grassland, bare land, building land, mangroves, and water showed a decreasing trend from 2016 to 2021, while the number of patches of wetland, forest land, and farmland showed an increasing trend.Among them, the number of forest land patches decreased sharply during 2019-2020; the number of wetland patches continued to increase from 2016 to 2018, which may be due to the issuance of the "Notice of Haikou City on the Implementation of the Key Points for the Construction of Ecological Civilization in Hainan Province in 2016".We must continue to strengthen the construction of ecological civilization, carry out the construction of ecological restoration in a planned way, and vigorously protect wetland resources.It can be seen from Figure 7b that the mean patch area of wetland, bare land, mangrove, and grassland during the study period is relatively small, indicating that their patches were relatively fragmented.In 2017, the mean patch area of wetland was the smallest, which was 0.0913 hm 2 .The mean patch area

Analysis of Landscape Pattern Changes at Class Level
From 2016 to 2021, the landscape pattern of Haikou at a class level underwent great changes, and the results are shown in Figure 7. From Figure 7a, it can be seen that the number of patches of grassland, bare land, building land, mangroves, and water showed a decreasing trend from 2016 to 2021, while the number of patches of wetland, forest land, and farmland showed an increasing trend.Among them, the number of forest land patches decreased sharply during 2019-2020; the number of wetland patches continued to increase from 2016 to 2018, which may be due to the issuance of the "Notice of Haikou City on the Implementation of the Key Points for the Construction of Ecological Civilization in Hainan Province in 2016".We must continue to strengthen the construction of ecological civilization, carry out the construction of ecological restoration in a planned way, and vigorously protect wetland resources.It can be seen from Figure 7b that the mean patch area of wetland, bare land, mangrove, and grassland during the study period is relatively small, indicating that their patches were relatively fragmented.In 2017, the mean patch area of wetland was the smallest, which was 0.0913 hm 2 .The mean patch area index of building land, water, and farmland showed alternating changes of increasing and decreasing, but generally increased slightly.Forest land fragmentation decreased during 2020-2021.This indicates that the degree of fragmentation of these three landscape types decreased.From the largest patch index in Figure 7c, it can be seen that the building landscape has increased rapidly overall, indicating that its dominant position has continuously improved.The largest patch index of forest land fluctuated greatly, with a maximum of 19.7107 in 2016 and a minimum of 3.4054 in 2020.The reason may be the rapid economic development and the acceleration of urbanization in recent years.It can be seen from Figure 7d that the landscape shape of the bare land is the most complex, and the landscape shape of the mangrove is relatively simple.Except for bare land, the landscape shape index of all landscape types increased during 2019-2020, indicating that the landscape shape became more complex during this period.7d that the landscape shape of the bare land is the most complex, and the landscape shape of the mangrove is relatively simple.Except for bare land, the landscape shape index of all landscape types increased during 2019-2020, indicating that the landscape shape became more complex during this period.

Analysis of Landscape Pattern Changes
The normalized results of the three indicators (brightness, greenness, and humidity) of Haikou City from 2016 to 2021 are shown in Figure 8. Combined with the statistical results (Table 7) of the three indicators (SBI, GVI, and WET) from 2016 to 2021, we can determine that the mean brightness of Haikou in 2016 was 0.5997.From then, the standard deviation of brightness value first decreased and then increased until it reached 0.5848 in 2021.Comparing the brightness distribution maps from 2016 to 2021, it can be clearly seen that the brightness values in the northeast and northwest corners of Haikou City decreased significantly.As for the greenness distribution maps, it is easy to observe that the southern part of Haikou City has a higher greenness value.One reason for this is that forest land and farmland are mostly distributed here.Another reason is that the building land is gathered in the northern part of Haikou, and more human activities have resulted in the destruction of the ecological environment.The mean value of greenness first decreased and then increased, and the overall greenness value increased, showing that the green area of Haikou City has generally increased and the ecological environment has gradually improved.The standard deviation of humidity in Haikou City in 2016 was 0.0386.The mean value of WET first decreased and then increased during 2016-2021, with small-scale changes showing a downward trend overall.From the spatial distribution maps of humidity, it can be seen that the areas with high humidity in Haikou City are distributed along the northern coast of Haikou.The changes in humidity are related to the population density and regional economic development.The decrease in humidity value may be caused by the increase in population and the acceleration of regional economic development.It can be clearly observed from the ecological environment quality distribution maps of Haikou City in 2016 and 2021 (Figure 9) that the ecology of Haikou City in 2021 improved significantly compared to 2016, indicating that the overall quality of the ecological environment improved.The areas with better ecology are mainly distributed in the central and southern vegetation coverage areas, while the ecology of rivers, urban centers, and coastal areas is poor, which is closely related to the scope of human activities.

Conclusions
This paper analyzed the land use classification results of Haikou City based on the GEE platform.In the context of the Hainan Free Trade Port policy, Haikou City has undergone tremendous changes in LULP.The results of this study will play an important role in maintaining the sustainable development of Haikou City in order to provide a basis and guidance for the ecological restoration and free trade port construction of Haikou City.
The main conclusions are as follows:

Conclusions
This paper analyzed the land use classification results of Haikou City based on the GEE platform.In the context of the Hainan Free Trade Port policy, Haikou City has undergone tremendous changes in LULP.The results of this study will play an important role in maintaining the sustainable development of Haikou City in order to provide a basis and guidance for the ecological restoration and free trade port construction of Haikou City.
The main conclusions are as follows: (1) This paper uses methods such as the land use transfer matrix, land use degree, land use dynamic degree, and land use gravity center measurement to quantitatively analyze the land use changes in Haikou during the period from 2016 to 2021.The building land in Haikou City is mainly distributed in the northern coastal area; the central and southern parts of Haikou are rich in forest land and farmland; and the mangroves are mostly distributed in the northeastern coastal area of Haikou.In the past six years, the area of building land in Haikou City has decreased overall, while the area of forest land, farmland, and wetland has increased.In recent years, Hainan Province's idle land disposal policy has been continuously improved, and the province's farmland, permanent basic farmland, forest land, etc., are being approved and incorporated into the national land space planning, which is conducive to promoting land use changes in Haikou City.The overall shift in the building's gravity center to the northwest of Haikou is also inseparable from the rapid economic development.
(2) Regarding the landscape, in the research stage, the degree of landscape fragmentation in Haikou City has decreased, the landscape shape has become increasingly simple, and the heterogeneity of landscape distribution has increased.During the development of the free trade port plan, in terms of the operation mode, Haikou has explored a more flexible policy system in terms of domestic and foreign trade and investment.Based on the economic development of various industries such as tourism, the building land has increased rapidly, and its dominant position has been continuously improved, making Haikou's island-based economy and society develop with high quality.While developing the economy, people are advocating for green development and insisting on thick planting and openness, so that the gap between various landscape types is gradually reduced, which gradually reduces the gap between various landscape types and reduces the degree of landscape fragmentation.(3) After comparing the normalized ecological quality maps of the two periods, 2016 and 2021, it is easy to conclude that the ecological quality of Haikou City improved significantly during the study period.Since the construction of the free trade port requires the protection of the ecological environment, the six iconic projects such as the construction of national parks and the construction of clean energy islands have also achieved certain results.Most of the areas where the ecological environment has improved are forest parks, wetland parks, woodlands, and grasslands; most of the areas with poor ecological environment quality are vegetation-free areas and new construction areas.This shows that the ecological construction promoted by the

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

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

Figure 2 .
Figure 2. The flow chart of the methodology.

Figure 2 .
Figure 2. The flow chart of the methodology.

Figure 4 .
Figure 4. Line chart of land use degree changes in Haikou.

Figure 4 .
Figure 4. Line chart of land use degree changes in Haikou.

Sustainability 2024 ,
16, x FOR PEER REVIEW 14 of 22 5.1.4.Analysis of the Gravity Center Shift of Building Land in Haikou

Figure 5 .
Figure 5.The gravity center shift of building land in Haikou City.

Figure 5 .
Figure 5.The gravity center shift of building land in Haikou City.

5. 2 .
Analysis of Landscape Pattern Changes5.2.1.Analysis of Landscape Pattern Changes at Landscape LevelThe results of the landscape pattern index of Haikou City at the landscape level from 2016 to 2021 are shown in Figure6.It can be seen from the figure that the landscape pattern of Haikou City from 2016 to 2021 generally shows the characteristics of gradually increasing landscape fragmentation, simplifying the shape of the landscape, and increasing landscape heterogeneity.From 2016 to 2021, the mean patch area and patch density in Haikou City have alternating changes of increase and decrease.The mean patch area reached its maximum in 2019, with an area of 0.93 hm 2 .Overall, the mean patch area increased from 0.612 hm 2 in 2016 to 0.698 hm 2 in 2021.The patch density reached 1556.857pieces/100 hm 2 in 2020.These show that under the influence of the external environment, Haikou's landscape distribution is becoming increasingly fragmented, and the degree of interlacing is deepening.As we can see from Figure6c, the landscape shape index dropped from 380.431 in 2016 to 315.534.The overall decline in the landscape shape index indicates that the distribution of various types of patches in Haikou is becoming more and more compact, and the shape of the patches is more regular and simplified.The constant changes in the largest patch index, Shannon's diversity index, and Shannon's evenness index indicate that the landscape planning of Haikou City has been constantly adjusted.The largest patch index decreased from 19.711 in 2016 to 19.106 in 2021, while Shannon's diversity index and Shannon's evenness index decreased overall.This indicates that the proportional gap between various types of landscapes was slightly increasing.Some dominant landscapes are more prominent, the diversity and balance of the landscapes are relatively reduced, and the heterogeneity of landscape distribution is increasing.Sustainability 2024, 16, x FOR PEER REVIEW 15 of 22prominent, the diversity and balance of the landscapes are relatively reduced, and the heterogeneity of landscape distribution is increasing.

Figure 6 .
Figure 6.Changes in landscape pa ern index at the landscape level: (a) mean patch area index; (b) largest patch index; (c) landscape shape index; (d) patch density; (e) Shannon's diversity index; and (f) Shannon's evenness index.

Figure 6 .
Figure 6.Changes in landscape pattern index at the landscape level: (a) mean patch area index; (b) largest patch index; (c) landscape shape index; (d) patch density; (e) Shannon's diversity index; and (f) Shannon's evenness index.

Sustainability 2024 ,
16, x FOR PEER REVIEW 16 of 22 development and the acceleration of urbanization in recent years.It can be seen from Figure

Figure 7 .
Figure 7. Changes in landscape pa ern index at the class level: (a) number of patches; (b) mean patch area; (c) largest patch index; and (d) landscape shape index.

Figure 7 .
Figure 7. Changes in landscape pattern index at the class level: (a) number of patches; (b) mean patch area; (c) largest patch index; and (d) landscape shape index.

Figure 8 .
Figure 8.The SBI, GVI, and WET maps of Haikou in 2016 and 2021.Figure 8.The SBI, GVI, and WET maps of Haikou in 2016 and 2021.

Figure 9 .
Figure 9. Distribution maps of Haikou's ecological environment quality in 2016 and 2021: (a) in 2016 and (b) in 2021.

Table 1 .
Classification indicators used in the study.

Table 1 .
Classification indicators used in the study.

Table 2 .
Classification values of land use degree classification index.

Table 3 .
Spatial metrics selected in this study.

Table 4 .
Accuracy evaluation of classification results.

Table 5 .
The land use transfer matrix of Haikou from 2016 to 2021 (km 2 ).

Table 4 .
Accuracy evaluation of classification results.

Table 6 .
The land use dynamic degree and change rate of Haikou from 2016 to 2021.

Table 7 .
The results of SBI, GVI, and WET of Haikou from 2016 to 2021.