Monitoring Sea Ice in Liaodong Bay of Bohai Sea during the Freezing Period of 2017/2018 Using Sentinel-2 Remote Sensing Data

It is of great significance to monitor sea ice for relieving and preventing sea ice disasters. In this paper, the growth and development of sea ice in Liaodong Bay of Bohai Sea in China were monitored using Sentinel-2 remote sensing data during the freezing period from January toMarch in 2018. Based on the comprehensive analysis of the spectral characteristics of seawater and sea ice in visible bands, supplemented by the Normalized Difference Snow Index (NDSI) and the Normalized Difference Vegetation Index (NDVI), we proposed a newmethod based on decision tree classification for extracting sea ice types in Liaodong Bay of Bohai Sea. Using the remote sensing data of eight satellite overpasses acquired from Sentinel-2A/B satellites, the distribution and area of the different sea ice types in Liaodong Bay during the freezing period of 2017/2018 were obtained. Compared with the maximum likelihood (ML) classification method and the support vector machine (SVM) classification method, the proposedmethod has higher accuracy when discriminating the sea ice types, which proved the newmethod proposed in this paper is suitable for extracting sea ice types from Sentinel-2 optical remote sensing data in Liaodong Bay. And its classification accuracy reaches 88.05%.)ewhole process of evolution such as the growth and development of sea ice in Liaodong Bay during the freezing period from January to March in 2018 was monitored. )e maximum area of sea ice was detected on 27 January 2018, about 10,187 km. At last, the quantitative relationship model between the sea ice area and themean near-surface temperature derived by MODIS data in Liaodong Bay was established. )rough research, we found that the mean near-surface temperature was the most important factor for affecting the formation and melt of sea ice in Liaodong Bay.


Introduction
e Bohai Sea in China is the southernmost frozen sea in the Northern Hemisphere. Affected by the Siberian cold current and the stable cold high-pressure in Mongolia, the sea ice occurs mostly in winter, especially in the north of Bohai Sea, the Liaodong Bay [1,2]. Sea ice disaster is a serious natural disaster. e freezing and drifting of sea ice on a large scale are a threat to coastal aquaculture, fishery production, navigation, offshore production operations, and oil-gas exploration, and the resulting economic losses are also serious [3][4][5]. It is of great significance to effectively monitor the types, the distribution, and the spatiotemporal evolution of sea ice and to analyze the mechanism of its growth, development, and melt for relieving and preventing the sea ice disaster in the Bohai Sea.
Satellite remote sensing technology has the characteristics of economy, timeliness, and large-area simultaneous observation, so it is an important technical means for monitoring the sea ice [2][3][4][5][6]. e satellite remote sensing technology for monitoring sea ice mainly includes microwave radiometer, microwave scatterometer, synthetic aperture radar (SAR), microwave altimeter, hyperspectral remote sensing, and optical remote sensing [2][3][4][5][6][7][8][9][10][11][12][13][14][15]. ese technologies have been applied to monitor the sea ice. At present, one of the important parameters of monitoring the sea ice is the sea ice types. It is of great significance to accurately discriminate the sea ice types for calculating the area of sea ice, assessing the conditions of sea ice, navigation, marine production operations, and so on.
Many experts and scholars have done a lot of research work in monitoring the sea ice types by satellite remote sensing technology [6][7][8][9][10][11][12][13][14][15]. Most of the studies have been carried out on sea ice in polar and high latitudes, and a series of research results have been achieved [7][8][9][10]. Dabboor et al. used four-polarized RADARSAT-2 data to classify sea ice of the Arctic region by using a classification method that quantitatively selects the best subset of CP (compact polarimetry) parameters and achieved high classification accuracy [8]. Angelika et al. proposed a semiautomatic classification algorithm and used EM-bird to measure the thickness of sea ice, which improved the accuracy of aerial photographic sea ice classification [9]. Liu et al. used the support vector machine (SVM) method based on texture features and sea ice concentration to classify by using RADARSAT-2 dual-polarized ScanSAR data [10]. A review of the main approaches developed for sea ice classification using satellite imagery was presented [11]. e fully automatic design of a numerically optimized decision tree algorithm and demonstration of its application to sea ice classification from SAR data was introduced [12]. ere are relatively few studies on monitoring the sea ice in Bohai Sea. Zhang et al. extracted the polarization scattering characteristics of sea ice using full-polarimetric SAR images and then carried out the research in the classification algorithm of sea ice types by polarimetric SAR [13]. Combining the advantages of support vector machine (SVM) and Markov random field (MRF), Shen et al. proposed a classification system of MRF-v-SVC for sea ice classification of SAR images [14]. Karvonen et al. estimated sea ice parameters including sea ice saturation and sea ice thickness in Bohai Sea, based on thermodynamic sea ice model and SAR, radiometer and other Earth observation data [15]. Microwave sensors, such as SAR, radar altimeter and microwave radiometer, are the main sensors used for monitoring sea ice in polar regions [6][7][8][9][10][11][12]. And they are used to monitor the sea ice in the Bohai Sea [13][14][15], but these sensors have limited capabilities in discriminating sea ice types.
Optical remote sensing has obvious advantages in discriminating the types of terrestrial objects due to its high spatial resolution and abundant spectral information, and it has been applied to monitor the sea ice [16][17][18][19][20][21][22]. For example, Shi et al. used MODIS data to obtain the distribution information of sea ice in Bohai Sea [16,17]. Su et al. used the method of gray level cooccurrence matrix texture analysis to obtain the distribution of sea ice in Bohai Sea based on the MODIS data [18]. Chen et al. used HJ-1 CCD data and support vector machine classification method to extract sea ice in Liaodong Bay of Bohai Sea [19]. Yan et al. used MODIS data to estimate the amount of sea ice resources in the winter of 2012-2013 [20]. Yan et al. used Geostationary Ocean Color Imager (GOCI) data to monitor the characteristics of the Bohai Sea ice during the winter of 2012-2013, including area and thickness [21]. rough analysis, it can be found that the current research on the sea ice in the Bohai Sea mainly focuses on the distribution, concentration, ice edge, and ice thickness.
ere were fewer studies on monitoring the sea ice types and the situation of sea ice growth and development in the Bohai Sea. Bohai Sea ice appears every winter in the past 60 years, and this area is an important economic development zone in China and undertaking busy maritime traffic tasks. So it is very important to monitor the sea ice in the Bohai Sea. One of the tasks of monitoring sea ice is to extract the spatial distribution information of sea ice types, which is the sea ice classification mentioned in this paper. Most of the existing studies were based on MODIS data. While the spatial resolution of MODIS data is low, it cannot meet the needs of discriminating sea ice types, and it cannot well be used to analyze the whole process of sea ice growth, development, and melt through the evolution of sea ice types in the Bohai Sea. High-resolution optical satellites (such as SPOT, QuickBird, GeoEye, Worldview, and geostationary satellites) have a long repetition period, which cannot meet the needs of the high temporal resolution in monitoring the growth and development of sea ice. To monitor the growth and development of sea ice in Liaodong Bay of Bohai Sea in detail, it is necessary to find some satellite data with short repetition period and high spatial resolution. With the development of satellite network observation modes, especially in the form of constellation established by multiple satellites, such as the Sentine1-2A/2B optical remote sensing satellite launched by the European Space Agency (ESA) [23][24][25], its repetition period is greatly shortened, which is conducive to analyze the whole process of sea ice growth and development. In addition, there are many kinds of sea ice, including new ice, white ice, gray-white ice, fast ice, and so on. e difference of spectral reflectance characteristics of different sea ice is much smaller than that of terrestrial objects. It is another difficult problem to be solved to accurately discriminate various sea ice types from similar spectral reflectance. In addition, the typical reflectivities of optical images may be affected by the effect of snow on sea ice and the issue that the waters in Liaodong Bay are turbid and not clear, containing a lot of sediments, complicating optical remote sensing of sea ice in this bay. erefore, to solve the problems in monitoring growth and development of sea ice based on remote sensing technique in Liaodong Bay, this paper used Sentine1-2A/2B optical remote sensing data to find a method of extracting sea ice types in Liaodong Bay of Bohai Sea that is suitable for medium-high resolution optical remote sensing satellites. On this basis, the multitemporal data were used to monitor the temporal and spatial evolution of sea ice in Liaodong Bay of Bohai Sea during the freezing period in 2017/2018. e main sections of this paper are organized as follows. In Section 2, the study area and data are introduced in detail. In Section 3, based on the analysis of the spectral characteristics of sea ice, we propose a new method of extracting the sea ice types of Liaodong Bay that is suitable for mediumhigh resolution optical images. In Section 4, the results are showed. A total of 32 Sentinel-2 optical images from January to March 2018, were classified by decision tree classification.

2
Journal of Spectroscopy And our results were compared with that obtained by the maximum likelihood classification method and support vector machine classification method. In Section 5, the distribution of sea ice types in Liaodong Bay during the freezing period of 2017/2018 is obtained. en the spatialtemporal evolution of sea ice growth and development in Liaodong Bay was analyzed. e mechanism of growth and development of sea ice in Liaodong Bay was given. Finally, some valuable conclusions drawn from this study are given in Section 6.

Study Area.
e Bohai Sea, located in the northeast of China, is the only frozen sea in China. It includes Liaodong Bay, Bohai Bay, and Laizhou Bay. Among them, Liaodong Bay is located in the north of the Bohai Sea. As shown in Figure 1, Liaodong Bay is located at 39-41°N and 120.5-122.5°E. It has the highest latitude and the most serious sea ice disaster area in China. Liaodong Bay freezes every winter in the past 60 years, and the main sea ice types are new ice, white ice, gray-white ice, fast ice.

Experimental Data.
In order to monitor the distribution of sea ice and the whole process of sea ice growth and development in Liaodong Bay of Bohai Sea during the freezing period of 2017/2018, a total of 32 satellite remote sensing images with medium-high resolution acquired from Sentinel-2A/2B were used as the main data source.
e Sentinel-2A satellite was launched on 23 June 2015. It was mainly equipped with the optical sensor called a Multispectral Imager (MSI), including 13 spectral bands, ranging from optical, near infrared to short-wave infrared, with several spatial resolutions (e.g., 10 meters, 20 meters, and 60 meters) [23][24][25]. e spatial resolution of the panchromatic band is 10 meters, and the resolution of the multispectral band is 20 meters. e central wavelength and resolution of each band of Sentinel-2A/2B is shown in Table 1. e swath width of the Sentinel-2 satellite image is about 290 km, so it is especially suitable for monitoring the sea ice with a wide range. On 7 March 2017, ESA launched Sentine1-2B, a sister satellite of Sentinel-2A. Its satellite parameters are exactly the same as those of Sentinel-2A. e repetition period of a single satellite is 10 days. In the form of constellation established by two satellites, the period of image acquisition is greatly shortened, and the two satellites can revisit every 5 days in the equatorial region of the Earth. In the case of high latitude regions, the repetition period is only 3 days. e short repetition period of 5 days in Bohai Sea also can provide more data with medium-high resolution when studying the growth and development of sea ice in Liaodong Bay of Bohai Sea. Currently, Sentinel-2A/2B data are officially available for users around the world for free.
We collected 8 satellite overpasses of the Sentine1-2A/2B optical remote sensing data in the study area during the freezing period of 2017/2018, as shown in Table 2. And the cloud covers of all data were less than 10%. ere are 3 satellite overpasses data of Sentinel-2A and 5 satellite overpasses data of Sentinel-2B. e swath-width of MSI images from Sentinel-2 is 290 km, but ESA cuts the complete data of each strip into 100 km-sized images and stores them in different files, each of which is less than 1 GB, in order to reduce the amount of data and facilitate data management. So the data of each overpass needs 4 scenes to completely cover Liaodong Bay in Bohai Sea. e details of the experimental data are as follows: It is not enough to only use the data in the paper to study the growth and evolution of sea ice in Liaodong Bay. However, the focus of this paper is to study an automatic classification method of sea ice based on decision tree, so as to monitor the growth and evolution of sea ice in Liaodong Bay in the winter of 2018.
To make full use of the spectral characteristics of Sentinel-2A/B for discriminating sea ice types, six bands were mainly used, including band B2 (blue), B3 (green), B4 (red), B5 (vegetation red edge), B8A (vegetation red edge), and B11 (SWIR, short-wave infrared). e data were uniformly sampled to the pixel size of 20 meters for discriminating the sea ice types.
In addition, in order to derive the near-surface temperature, we used the MODIS L1B (MOD02 1 km) data on the same day as Sentinel-2 data. MODIS data can be obtained at NASA's official website: http://modis.gsfc.nasa. gov/. e imaging time of the MODIS data was all later less than 0.5 hours from that of the Sentinel-2 data. e SRTM (Shuttle Radar Topography Mission) DEM (Digital Elevation Model) and the meteorological data also were used to analyze the growth and development of sea ice in Liaodong Bay. And the SRTM data were downloaded from http://srtm.csi.cgiar.org. e topographic map of study area is shown in Figure 2.

Methods
At present, there are many classification methods for optical remote sensing images, such as the maximum likelihood classification method, support vector machine (SVM) classification method, and decision tree classification method [4,22,26].

Technical Method.
To meet the needs of discriminating the sea ice types in Liaodong Bay of Bohai Sea, on the basis of analyzing the features of medium-high resolution of Sentinel-2 optical remote sensing image and the spectral characteristics of sea ice, we proposed a method for extracting the sea ice types based on decision tree. en the whole process of sea ice growth and development in Liaodong Bay of Bohai Sea was monitored by using multitemporal Sentinel-2 optical remote sensing data and the proposed method based on the decision tree. e specific technical method is shown in Figure 3. Firstly, the original Sentinel-2 images should be preprocessed, including radiometric calibration, atmospheric correction, image fusion, and mosaic. en the training samples of different sea surface objects were selected to analysis spectral characteristics of different sea ice types and to derive the parameters of NDVI and NDSI. At last, the classification rules suitable for extracting the sea ice types were established by the reflectance differences of various types of sea ice.

Preprocessing.
Firstly, the original Sentinel-2 remote sensing images should be preprocessed, including radiometric calibration, atmospheric correction, image fusion, and mosaic [24]. Radiometric calibration is to convert the grayscale brightness value of an image into absolute radiance brightness value. Atmospheric correction is mainly to eliminate the radiation error caused by atmospheric affect. We used the Sen2Cor plug-in, developed by ESA specifically, to carry out the radiometric calibration and atmospheric correction for the Sentinel-2 MSI images. Due to the wide range of Liaodong Bay in the Bohai Sea, at least 4 scenes of Sentinel-2 remote sensing images were needed to completely cover the whole study area. erefore, image mosaic was necessary. e color differences between different images need to be dealt with during the mosaic process. e true color composite images of Sentinel-2 data after preprocessing are displayed in Figure 4.

Establishing the Decision Trees.
e decision tree classification algorithm is an inductive learning algorithm, which is based on samples. e algorithm infers classification rules from a set of unordered and irregular samples and establishes a classifier or a prediction model and then classifies the data of unknown categories. e physical basis of discriminating targets is to distinguish them according to the differences of spectral characteristics between targets. Firstly, sea ice monitoring is actually to discriminate between seawater and sea ice. It is required to analyze the spectral characteristics of seawater and various sea ice types statistically. e spectral curves can be used to analyze the reflectance differences of seawater and various sea ice types in different spectral regions, and then   Analyzing the temporal and spatial evolution of sea ice over multiple overpasses   the classification rules of decision tree classification method can be established to achieve the purpose of sea ice classification [26]. According to the spectral characteristics of sea ice in optical images and the existing research in the Liaodong Bay of Bohai Sea [1,4,11,[18][19][20][21], the feature types are mainly divided into land area (LA), seawater (SW), fast ice (FI), new ice (NI), white ice (WI), and gray-white ice (GWI).
In order to analyze the spectral characteristics of seawater and various sea ice types, we selected the image acquired on 1 February 2018, in Liaodong Bay as the training sample data set. rough our visual interpretation, multiple ROIs (region of interest) of each type of sea surface target are selected for the samples, and each ROI contains 20-200 pixels. e samples have been indicated in Figure 4(  We should build the rules for the classification of decision tree firstly. As shown in Figure 5, the reflectivity of seawater in each band is low, while the difference between sea ice and seawater in band B2 and band B5 is obvious. So the band B2 and band B5 were selected to discriminate between seawater and sea ice. We used the Normalized Difference Snow Index (NDSI) to extract snow information [27]. It mainly showed that snow has higher reflectivity in the visible band and lower reflectivity in the short-wave infrared band. We found that characteristics of sea ice are similar to snow. erefore, NDSI is regarded as important auxiliary information to discriminate the sea ice types. e formula of NDSI is as follows: where B3 represents the green band and B11 represents the short-wave infrared band of Sentinel-2A/B. Taking band B2 and band B5 as the characteristic bands and using NDSI as the characteristic parameter, the appropriate classification rules were established to monitor the sea ice. rough a lot of experimental analysis, we found that the threshold of the reflectivity of B2 band for discriminating seawater and sea ice is 0.065, that is, less than 0.065 for seawater and greater than 0.065 for sea ice. e reflectivity of band B5 has a threshold of 0.08, which is less than 0.08 for seawater and greater than 0.08 for sea ice. e two thresholds are very close to the typical albedo value of open sea water in existing literature [28]. e albedo threshold of band B5 is larger than that of band B2, which may be due to sediment. e Normalized Difference Snow Index (NDSI) has a threshold of 0.65, which is less than 0.65 for seawater and greater than 0.65 for sea ice. On the basis of sea ice monitoring, the spectral reflectance characteristics of various sea ice types are further explored to achieve sea ice classification.
As can be seen from Figure 6, the reflectivity of fast ice and white ice in each band is high, and the spectral characteristics are not different, so it is difficult to discriminate them only by spectral differences. It is found that the Normalized Difference Vegetation Index (NDVI) has better effect on separating the fast ice from others [29]. e formula for NDVI is as follows: where B4 represents the red band and B8A represents the near infrared band of Sentinel-2A/B. If the value of NDVI is greater than −0.09, it is classified into fast ice. After repeated classification experiments and spectral analysis, the optimal threshold of decision tree classification model can be selected. e reflectivity of new ice in band B8A is lower than 0.07, so a classification rule of decision tree is established to discriminate new ice by band B8A. It can be seen from Figure 6 that the reflectivity of white ice is the most different from that of gray-white ice in the band B5. erefore, the band B5 is selected as the basis for discriminating white ice from gray-white ice. e threshold value is 0.31; that is, the reflectivity of gray-white ice is less than 0.31, and the reflectivity of white ice is greater than 0.31.
According to the above analysis of reflectance characteristics, the classification rules suitable for extracting the sea ice types were established by the reflectance differences of various sea ices in each band. Finally, the decision tree classification model for discriminating the sea ice types was generated. e final decision tree of the classification model was established as shown in Figure 7.

Classification Results.
Based on the proposed method for discriminating the sea ice types, the multiple satellite overpasses data during the freezing period of 2017/2018 in Liaodong Bay were used for classification. e classification results are displayed in Figure 8.

Comparison of Different Methods.
In order to evaluate the effectiveness of the proposed method based on decision tree for sea ice classification, our result of the proposed method was compared with those of maximum likelihood (ML) classification method and support vector machine (SVM) classification method. e images acquired from 11 February 2019, were selected to evaluate the classification accuracy. e ML classification method and SVM classification method are carried out in the same samples. e results of maximum likelihood classification method and support vector machine classification method are shown in Figures 9(a) and 9(b), respectively. e result of the proposed method is shown in Figure 9(c).
By comparing the classification results with the original optical images, we found the distribution of sea ice edge from the proposed method based on decision tree in this paper was clear and basically consistent with the actual distribution. At the same time, it could be seen that the degree with which different ice types are discriminated from each other was also relatively large, and the outline of fast ice around the island and the coastline was substantially the same as that of the original image. In addition, the outline of the white ice in the middle part of the ice region was also better displayed. In terms of efficiency, the proposed method based on decision tree was far superior to the other two algorithms, and it was an effective way to achieve the sea ice classification.
For the images acquired from 11 February 2018, a total of 500 validation samples were generated randomly using the sea ice type maps obtained by different classification methods.
en 500 samples were interpreted by manual visual interpretation. Combining the original optical image, the classification accuracies of various sea ice types were calculated which were determined by the proposed method, the maximum likelihood classification method, and the support vector machine classification method, as shown in Table 3. 8 Journal of Spectroscopy From Table 3, we can see that the classification accuracies of these methods are between 74.30% and 88.05%. e classification accuracy of the proposed method is higher than the other two methods. And its classification accuracy reaches 88.05%, which indicates that the proposed method is suitable for Sentinel-2 optical remote sensing data to discriminate the sea ice types in Liaodong Bay of Bohai Sea.

Spatiotemporal Evolution of Sea Ice of Liaodong Bay.
e spatiotemporal evolution analysis of sea ice growth, development, and melt was carried out based on the classification results from 2 January to 8 March 2018 (Figure 8). e area of sea ice and its proportion in seawater was calculated, as shown in Table 4.
According to Table 4 and Figure 10, the evolution trend of sea ice in Liaodong Bay was as follows: with the decrease of near-surface temperature, the thickness of sea ice gradually increased. e types of sea ice gradually developed from new ice to gray-white ice and white ice. And the distribution area of sea ice gradually expanded. On 27 January 2018, the sea ice area reached the maximum, which was about 10,187 km 2 , covering almost the whole Liaodong Bay. en, as the near-surface temperature rose, the sea ice slowly melted. And the sea ice area gradually decreased until it finally disappeared completely.

e Mechanisms of Growth and Development of Sea Ice in Liaodong Bay.
e near-surface temperature was derived using the method introduced in reference [30,31] through the MODIS L1B (MOD02 1 km) data of different time. Since the interval of imaging time between MODIS and Sentinel-2 did not exceed half an hour, it could approximately be considered that the mean near-surface temperature derived from the MODIS data represented that of Sentinel-2 imaging. In addition, according to the weather forecast information issued by China Meteorological Administration, the maximum and minimum temperatures in the study area were counted, as shown in Figure 11. It showed the mean near-surface temperature, minimum temperature, and maximum temperature at the eight satellite overpasses in Liaodong Bay during the freezing period in 2017/2018.
According to the distribution of sea ice of Liaodong Bay during the freezing period in 2017/2018, the growth and development of sea ice in Liaodong Bay was mainly affected by the mean near-surface temperature from the MODIS data.          When the mean near-surface temperature decreases to 0°C, it begins to freeze in Liaodong Bay. e sea ice type is mainly new ice, and the thickness of sea ice is thin. As per the further decrease of mean near-surface temperature, the frozen area is expanding, and the sea ice types gradually develop from new ice to gray-white ice and white ice, and the thickness of sea ice gradually increases. If the mean near-surface temperature continues to increase above 0°C, sea ice will melt and sea ice area will be greatly reduced. rough linear regression analysis, a quantitative relationship model between the sea ice area and the mean nearsurface temperature derived from the MODIS data was established in Liaodong Bay (as shown in Figure 12). e specific relationship model was as follows:

Journal of Spectroscopy
where y is the sea ice area at a certain time (unit: km 2 ) and x is the mean near-surface temperature at this time (unit:°C).

Conclusions
In this paper, we used a total of 32 Sentinel-2 optical remote sensing images from January to March 2018 to monitor the sea ice in Liaodong Bay of the Bohai Sea in China. According to the spectral characteristics of sea ice and seawater, the classification rules based on decision tree were established to discriminate the sea ice types at different time. A new method based on decision tree classification for discriminating the sea ice types was proposed. e validity of the new method proposed in this paper was verified. At the same time, combined with the near-surface temperature derived from the MODIS data and meteorological data, the evolution of sea ice distribution and area in Liaodong Bay was analyzed. rough the study, some valuable conclusions can be drawn as follows: (1) e new method makes full use of the spectral characteristics, NDSI, and NDVI to discriminate the sea ice types. e classification accuracy of the proposed method was higher than those of the maximum likelihood classification method and the support vector machine classification method. Its classification accuracy reached 88.05%. e test proved that the new method is suitable for Sentinel-2 optical images with medium-high resolution to discriminate the sea ice types.
(2) e growth and development of sea ice in Liaodong Bay during the freezing period in 2017/2018 was analyzed through 8 Sentinel-2 satellite overpasses. As the mean near-surface temperature decreases, the sea ice in Liaodong Bay generally begins to freeze. e maximum area of sea ice was detected on 27 January 2018, about 10,187 km 2 . e sea types include fast ice, new ice, white ice, and gray-white ice in Liaodong Bay. (3) e mean near-surface temperature was the most important factor for affecting the formation and melt of sea ice in Liaodong Bay. e quantitative relationship model between sea ice area and mean nearsurface temperature in Liaodong Bay was established, which can provide a theoretical basis for analyzing the growth and development of sea ice in Liaodong Bay. e results of sea ice classification were affected by image quality and preprocessing. Because it needed 4 scenes of Sentinel-2 images to completely cover the study area, the uniform color processing during image mosaic has a certain influence on the classification accuracy of sea ice. In addition, the selection of training samples would also affect the accuracy of sea ice classification. e sample of the gray-white ice was relatively few, so the classification accuracy of gray-white ice was not very high. e evolution of sea ice was affected by many factors and there were existing randomness and variability.
Although the repetition period of Sentinel-2A/2B satellite network has been reduced to 5 days, the satellite remote sensing data of some satellite overpasses are greatly contaminated by clouds and fogs because of weather reasons, so they cannot be used. is made the data time interval larger in some satellite overpasses, which is not conductive to analyze the growth and development of sea ice in Liaodong Bay in detail. However, using the eight satellite overpasses of Sentinel-2 remote sensing data, the type change of the sea ice growth and development in Liaodong Bay can be basically ascertained. In the next, we will use Sentinel-1 radar satellite data that is not affected by clouds and fogs [32] to analyze the growth and development of sea ice in Liaodong Bay.

Data Availability
e Sentinel-2 data, MODIS data, and SRTM DEM data can be available for free.