Retrieval of Chlorophyll-a Concentrations of Class II Water Bodies of Inland Lakes and Reservoirs Based on ZY1-02D Satellite Hyperspectral Data

: Chlorophyll-a is an important parameter that characterizes the eutrophication of water bodies. The advantage of ZY1-02D hyperspectral satellite subdivision in the visible light and near-infrared bands is that it highlights the unique characteristics of water bodies in the spectral dimension, and it helps to assess the Class II water bodies of inland lakes and reservoirs, making it an important tool for reﬁned remote sensing detection of the environment. In this study, the Baiyangdian Nature Reserve in northern China, which contains a typical inland lake and wetland, was chosen as the study area. Using ZY1-02D hyperspectral synchronization transit images and in situ measured chlorophyll-a concentration as the data source, remote sensing of the chlorophyll-a concentration of inland lakes was conducted. By analyzing the correlation between the spectral reﬂectance of the ZY1-02D hyperspectral image and the chlorophyll-a concentration and using algorithms such as the single band, band ratio, and three bands to compare and ﬁlter characteristic wavelengths, a quantitative hyperspectral model of the chlorophyll-a concentration was established to determine the chlorophyll-a concentration of Baiyangdian Lake. The dynamic monitoring of the water body and the assessment of the nutritional status of the water body were determined. The results revealed that the estimation of the chlorophyll-a concentration of Baiyangdian Lake based on the hyperspectral Fluorescence Line Height (FLH) model was ideal, with an R 2 value of 0.78. The FLH model not only comprehensively considers the effects of suspended solids, yellow substances, and backscattering of the water body on the estimation of the chlorophyll-a concentration, but also considers the inﬂuence of the elastic scattering efﬁciency of the chlorophyll. Based on the ZY1-02D hyperspectral data, a spatial distribution map of the chlorophyll-a concentration of Baiyangdian Lake was created to provide new ideas and technical support for monitoring inland water environments.


Introduction
Chlorophyll-a is a pigment contained in all algae. It is an important indicator of the degree of eutrophication of water bodies and it is also monitored in Chinese aquatic environments [1]. It is difficult for traditional water quality monitoring based on on-site sampling and laboratory measurements to capture the temporal and spatial distributions of chlorophyll-a. In comparison, remote sensing technologies have many advantages, including the following: rapidity, ease of access, and large-scale, exhaustive coverage [2]. It can be an effective method of achieving long-term and large-scale monitoring of water eutrophication [3]. In consequence, a series of chlorophyll-a concentration inversion algorithms have been developed [4].
The quantitative inversion of the chlorophyll-a concentration is an important remote sensing method used to assess the nutritional status of water bodies. For example, with ocean-As the largest freshwater lake in the North China Plain, Baiyangdian Lake is characterized by a vast amount of water, complex ground cover types, and strong human influence [19]. The water environment of a lake can reflect the health of the regional ecology to a certain extent. The chlorophyll-a concentration is representative of the primary productivity and richness of the water body, as well as an important indicator of the degree of nutrition. In this study, a combination of Baiyangdian ZY1-02D image data and measured chlorophyll-a data were used to compare and analyze the commonly used chlorophyll-a concentration inversion algorithm and to construct a chlorophyll-a inversion algorithm for Baiyangdian Lake. The objective was to be able to provide suggestions for the eutrophication management of Baiyangdian Lake and guidance for the Xiong'an district

Study Area
Baiyangdian Lake is located in Anxin County, Xiong'an New District, Hebei Province ( Figure 1). It is known as the "Pearl of North China" and is the largest shallow lake wetland in the North China Plain. The water in Baiyangdian Lake is supplied by nine rivers in the Daqing River system and is also called the "nine rivers under the tip" [20]. In modern geography, these nine rivers are the Zhulong River, Xiaoyi River, Tang River, Qingshui River, Baohe River, Fuhe River, Ping River, Caohe River, and Baigou River. In contrast to other Chinese large lakes and reservoirs such as Taihu Lake, Chaohu Lake, and Guanting Reservoir, Baiyangdian Lake is characterized by a complex topography, with many connected rivers and ravines. More than 3700 channels and waterways connect 143 lakes of different sizes to form reed fields and garden terraces, creating a unique natural landscape consisting of a mosaic of depressions and villages. The total area of Baiyangdian Lake is 366 km 2 . There are more than 6 million square meters of lakes, including Baiyangdian Lake, Shaoche Lake, Mapeng Lake, Yangjiao Lake, Chiyu Lake, Houtang Lake, and Xiaobeidian Lake. There are 24 lakes with areas of 10,000 m 2 , and 112 lakes with areas of less than 600,000 m 2 . Among all of these lakes, Baiyangdian Lake has the largest area, so this conglomeration of lakes is referred to as Baiyangdian Lake. With the establishment of the Xiong'an New District, industries have developed rapidly in the upper reaches of Baiyangdian, causing serious water pollution. The main problems include reduced water inflow, insufficient water storage, a low water retention capacity, and poor water self-purification [21].

In Situ Chlorophyll-a Data
In this study, 21 sets of water quality sample data were collected simultaneously in Baiyangdian on 7 October 2020 according to ZY1-02D transit time, and the sampling points are shown in Figure 1. A Hach Hydrolab MS5 multi-probe water quality measuring Remote Sens. 2022, 14, 1842 4 of 16 instrument was used to conduct the on-site water quality measurements. The MS5 multiprobe can simultaneously measure the water temperature, Total Dissolved Solids (TDSs), chlorophyll-a (Chl-a), turbidity, and other parameters. The water quality was measured for 10 s at each sampling point, three replicate sets of measures were obtained at each sampling point. The average value was used in subsequent analysis.

In Situ Remote Sensing Reflectance Data
The research team carried out satellite-ground synchronous observation experiments in the study area according to the transit time of the satellite to obtain a water-color matching dataset. An ASD FieldSpec3 portable field hyperspectral ground spectrometer was used. The spectrum was collected under dry, windless, clear, and cloudless or less cloudy conditions, and the standard board (25% gray board) was corrected in time according to the weather conditions. The measurement time was taken during the period (10:00-14:00) when the Sun's altitude angle was greater than 45 • . The measurement of the water spectrum was conducted following the international standard by measuring above the water surface [22].
where R rs is the spectral reflectance value of the water body; S sw , S sky , and S p are the measured signal values when the spectrometer was facing the water, the sky, and the standard plate, respectively; and r is the reflectance of the air-water interface. Based on previous studies, r = 2.2% for calm water. ρ p is the reflectivity of the standard plate.
To obtain a ZY1-02D satellite-ground synchronously measured reflectance dataset for Baiyangdian, the following two principles were followed to match the reflectivity R rs of the sampling point and the corresponding ZY1-02D satellite image pixel. Spatially, it was ensured that the position of the pixel and the ground point were consistent. In addition, it was necessary to ensure that the pixel had a good uniformity, and the ratio of the standard deviation of the 3 × 3 neighborhood pixel value around the pixel to the mean value did not exceed 0.4. Temporally, it was ensured that the time interval between the pixel and the measured point was no more than 5 days. Therefore, the research team measured 8 sets of spectral data of water bodies in Baiyangdian on 7 October 2020. Each set of data was measured with 10 spectral curves, with curve processing including spectral smoothing, denoising, and averaging. The measured spectral curves were obtained as shown in Figure 2. Comparison of the in situ measured water spectrum with the spectrum of the water bodies before and after atmospheric correction of the same pixel.

ZY1-02D Image Data
The ZY1-02D satellite is China's first civilian hyperspectral service satellite. ZY1-02D is equipped with a 9-band multi-spectral camera and a 166-band hyperspectral camera, including 2.5 m panchromatic, 10 m multi-spectral, and 30 m hyperspectral cameras. The

ZY1-02D Image Data
The ZY1-02D satellite is China's first civilian hyperspectral service satellite. ZY1-02D is equipped with a 9-band multi-spectral camera and a 166-band hyperspectral camera, including 2.5 m panchromatic, 10 m multi-spectral, and 30 m hyperspectral cameras. The detailed parameters are shown in Table 1. Hyperspectral data can provide more accurate retrieval of standard aquatic products, and this sensor can potentially provide more detailed information about optically active components. The remote sensing reflectance is the basis of the application of hyperspectral remote sensing data. The hyperspectral reflectance of the ZY1-02D satellite is directly related to the efficiency and quality of applications involving hyperspectral data. The preprocessing of ZY1-02D hyperspectral data mainly includes coincident band elimination, radiometric calibration, and atmospheric correction.
According to the ZY1-02D satellite spectrum settings and the characteristics of the water spectrum curve, the last five bands of the visible near-infrared (NIR) band overlap with the first three bands of the NIR band, so that the redundant bands (77-79) are eliminated, and the total number of bands is 163.
Radiometric calibration is used to establish the quantitative relationship between the radiances of the different wavebands and the original pixel value of the sensor. The equation for the radiometric calibration of remote sensing images is ( where L is the apparent radiation measurement value; DN is the original value of the remote sensing image; Gain is the gain value of the absolute radiation calibration coefficient; and offset is the offset value. The signal detected by the remote sensor is affected by the absorption and scattering caused by atmospheric components (by approximately 90%). These phenomena are one of the uncertain factors in the remote monitoring of the chlorophyll-a concentration [23]. The purpose of atmospheric correction is to eliminate atmospheric molecular reflection. The apparent radiance or apparent reflectance after the radiation calibration is inverted into the true reflectance close to the ground object. There are many atmospheric correction methods, which can be divided into absolute atmospheric correction methods and relative atmospheric correction methods, according to their corrected results. The atmospheric radiation transfer model is the most accurate correction method among the many models of absolute atmospheric correction methods, including fast line-of-sight atmospheric analysis of hypercubes (FLAASH), Second Simulation of the Satellite Signal in the Solar Spectrum (6S), and Quick Atmospheric Correction (QUAC) models. The FLAASH model is suitable for multi-spectral and hyperspectral images and can accurately compensate for atmospheric effects. Its applicable wavelength range is from visible light to NIR and short-wave infrared, with a maximum wavelength of 3 µm. In the moderate resolution atmospheric transmission version 4 (MODTRAN4) model, the atmospheric correction parameters include the satellite transit time, aerosol type, and atmospheric model. The FLAASH model provides three wavebands for water vapor inversion: 1050-1210 nm (1135 nm), 870-1020 nm (940 nm), and 770-870 nm (820 nm). The aerosol optical thickness is dark and is obtained using the pixel inversion method. The use of the MODTRAN standard atmospheric model and the aerosol model to replace real-time atmospheric parameters broadens the application range of the FLAASH model and improves its stability [24].
The 6S model is used to simulate the influence of the atmosphere during the transmission of sunlight from the Sun to the ground target to the sensor. Through simulations, the degree of influence of the atmosphere on the solar radiation can be obtained to carry out atmospheric correction. The latest 6S atmospheric correction code requires auxiliary data, including the moisture content, ozone content, aerosol optical thickness, and terrain height.
The QUAC model is an algorithm based on dark pixels, and it considers three hypotheses. (1) The image must contain more than 10 spectrally different pixels. (2) The standard deviation of the end-member pixel reflectance is spectrally independent and can be used to calculate the transmittance. (3) The existence of a related number of dark pixels is used to calculate the constant baseline, which is assumed to be the measurement of the attenuation (scattering and absorption) and remote sensing. The calculation speed of the atmospheric correction method is significantly faster than that of the method based on the radiation transmission model because the QUAC model does not involve first-principle radiation transmission calculations, and it only requires an approximate description of the sensor band position (the central wavelength) and its radiation calibration.
In this study, the FLAASH module of the ENVI 5.5 software was chosen to perform the atmospheric correction of the water bodies, and atmospheric correction of the ZY1-02D images was carried out after radiation calibration because the aerosols in the study area were not affected by urban or industrial inputs. Due to the strong influence of the source, the rural aerosol model in the FLAASH module and the 2-band (KT) aerosol inversion method were chosen. Figure 2 shows a comparison of the water reflectance data before and after the atmospheric correction for the same pixel and a comparison between the water reflectance after the atmospheric correction and the measured water reflectance. After the FLAASH atmospheric correction, the spectrum of the satellite image retains the major features of the measured spectrum.

Water Extraction
For small and medium-sized inland lakes, the lake area is greatly affected by the proximity of the land and vegetation [25]. The water body extraction area directly affects the accuracy of the chlorophyll-a inversion of the water body, which leads to distortion of the chlorophyll-a inversion model. Most lake water body methods use masks to solve the neighbor effect. A small mask range will increase the number of land pixels and the uncertainty of the remote sensing reflectance. If the mask range is too large, it will reduce the amount of effective data. The ZY1-02D data were collected using an AHSI sensor with a 30 m spatial resolution and a VNIC sensor with a 10 m spatial resolution.
Based on multi-spectral Landsat data, Mcfeeters proposed the use of the NIR band and the green light band to construct the Normalized Difference Water Index (NDWI) based on a large number of experiments [26]. The definition of the NDWI is where λ 1 is the reflectance in the green light band and λ 2 is the reflectance in the NIR band. The NDWI can quantitatively describe the performance characteristics of the water body in the image, so it can effectively realize the automatic extraction of the water body information. The underlying surface of Baiyangdian Lake is complicated by human activities, and the water area is small and fragmented. Based on the narrow-band water spectrum information about Baiyangdian Lake ZY1-02D, the band combination for the NDWI was determined. Baiyangdian Lake is an optically shallow water area. The water reflectance curve is shown in Figure 3. In the range of 400-500 nm, due to the strong absorption of blue light by chlorophyll-a and yellow substances, the water reflectance in this wavelength range is generally low. In the range of 510-620 nm, due to the weak absorption of chlorophyll-a and carotene and the scattering of cells and suspended particles, a reflection peak is formed. The size of the reflection peak is related to the pigment composition, and the higher the chlorophyll-a concentration of the water, the larger the peak. Thus, it can be used as a quantitative indicator of the chlorophyll-a concentration. Another reflection peak appears in the range of 685-730 nm. This peak is generally considered to be the fluorescence peak of chlorophyll-a. Its position and reflectivity increase as the chlorophyll-a concentration increases. Displacement occurs in the direction of the longer wavelengths, and this peak is the most prominent feature of algae-bearing water bodies. Its presence or absence is usually used as the basis for judging whether the water contains algae-related chlorophyll-a [27]. The position and value of the reflection peak are important indicators of the chlorophyll-a concentration. Beyond 730 nm, the reflectivity of the water body is relatively low due to the absorption by the water body.  The NDWI band was determined based on the spectral curve characteristics of the features of Baiyangdian Lake. The reflection peak and absorption valley in Figure 2 correspond to = 576 nm and = 730 nm in the ZY1-02D image. Therefore, the water extraction formula for the ZY1-02D hyperspectral data is ZY1-02D is equipped with a multi-spectral sensor with a spatial resolution of 10 nm. The central wavelength of B2 is 555 nm, and the central wavelength of B7 is 725 nm. The water extraction formula for the ZY1-02D multi-spectral data is The multi-spectral data bright surface water extraction results and the hyperspectral bright surface water extraction results were superimposed and analyzed, which effec- The NDWI band was determined based on the spectral curve characteristics of the features of Baiyangdian Lake. The reflection peak and absorption valley in Figure 2 correspond to λ 1 = 576 nm and λ 2 = 730 nm in the ZY1-02D image. Therefore, the water extraction formula for the ZY1-02D hyperspectral data is ZY1-02D is equipped with a multi-spectral sensor with a spatial resolution of 10 nm. The central wavelength of B2 is 555 nm, and the central wavelength of B7 is 725 nm. The water extraction formula for the ZY1-02D multi-spectral data is The multi-spectral data bright surface water extraction results and the hyperspectral bright surface water extraction results were superimposed and analyzed, which effectively eliminated the information about the banks, improved the accuracy of the water extraction, and obtained the water surface distribution of Baiyangdian Lake (Figures 4 and 5).
The NDWI band was determined based on the spectral curve characteristics of the features of Baiyangdian Lake. The reflection peak and absorption valley in Figure 2 correspond to = 576 nm and = 730 nm in the ZY1-02D image. Therefore, the water extraction formula for the ZY1-02D hyperspectral data is = rs (567) − rs (730) rs (567) + rs (730) ZY1-02D is equipped with a multi-spectral sensor with a spatial resolution of 10 nm. The central wavelength of B2 is 555 nm, and the central wavelength of B7 is 725 nm. The water extraction formula for the ZY1-02D multi-spectral data is The multi-spectral data bright surface water extraction results and the hyperspectral bright surface water extraction results were superimposed and analyzed, which effectively eliminated the information about the banks, improved the accuracy of the water extraction, and obtained the water surface distribution of Baiyangdian Lake (Figures 4  and 5).

Screening of Sensitive Bands
Single Band The correlation coefficients between the chlorophyll-a concentration and each band of the ZY1-02D pixel remote sensing reflectivity were calculated and the best band for the chlorophyll-a concentration inversion was identified. Figure 6 shows that the chlorophylla concentration at 400-680 nm is negatively correlated with the reflectance, while the chlorophyll-a concentration at 680-900 nm is positively correlated with the reflectance. The maximum correlation coefficient (0.6518) occurred at 713.79 nm, so the reflection at

Screening of Sensitive Bands Single Band
The correlation coefficients between the chlorophyll-a concentration and each band of the ZY1-02D pixel remote sensing reflectivity were calculated and the best band for the Remote Sens. 2022, 14, 1842 9 of 16 chlorophyll-a concentration inversion was identified. Figure 6 shows that the chlorophylla concentration at 400-680 nm is negatively correlated with the reflectance, while the chlorophyll-a concentration at 680-900 nm is positively correlated with the reflectance. The maximum correlation coefficient (0.6518) occurred at 713.79 nm, so the reflection at 713.79 nm was selected. The rate and chlorophyll-a concentration were used for the regression modeling.

Single Band
The correlation coefficients between the chlorophyll-a concentration and each band of the ZY1-02D pixel remote sensing reflectivity were calculated and the best band for the chlorophyll-a concentration inversion was identified. Figure 6 shows that the chlorophylla concentration at 400-680 nm is negatively correlated with the reflectance, while the chlorophyll-a concentration at 680-900 nm is positively correlated with the reflectance. The maximum correlation coefficient (0.6518) occurred at 713.79 nm, so the reflection at 713.79 nm was selected. The rate and chlorophyll-a concentration were used for the regression modeling.

Band Calculations
The band ratio model factors in the chlorophyll-a inversion mainly include the bluegreen band ratio and the NIR band ratio. Among them, the NIR and red-band ratio models are usually used for optically complex Case-II water bodies, which are mainly affected by the chlorophyll-a concentration, the total suspended solids, and the influence of colored dissolved organic matter. The commonly used NIR and red-band ratio models are R rs (700)/R rs (670), R rs (719)/R rs (667), R rs (706)/R rs (682), and R rs (708)/R rs (665). According to the band characteristics of the ZY1-02D data and the existing band combinations, the SPSS Statistics software (IBM) was used to perform the correlation analysis between the different band combinations and the measured chlorophyll-a concentrations to obtain the best band combination using Equation (6).

Three-Band Model
The three-band model is based on the bio-optical theory for the remote sensing inversion of the Chl-a concentrations in turbid Case-II water bodies [28].
The model is based on three assumptions. (1) The backscattering coefficients of the three bands are approximately equal. That is, b , and b b (λ 3 ) are the backscattering coefficients at λ 1 , λ 2 , and λ 3 , respectively.
(2) λ 3 is mainly determined by the absorption coefficient of pure water, and the absorption coefficient of the water components can be ignored. That is, a w (λ 3 ) >> a d (λ 3 ) + a CDOM (λ 3 ) , and a CDOM (λ 3 ) are the absorption coefficients of pure water, non-pigmented particulate matter, and colored dissolved organic matter at λ 3 , respectively. (3) The absorption coefficients of the non-pigmented particulate matter and colored dissolved organic matter at λ 1 and λ 2 are approximately equal, and λ 1 is located near the absorption peak of Chl-a. That is, a ph (λ 1 ) ≈ a ph (λ 2 ), a CDOM (λ 1 ) ≈ a CDOM (λ 2 ), and a ph (λ 1 ) >> a ph (λ 2 ), where a ph (λ 1 ) and a ph (λ 2 ) are the absorption coefficients of phytoplankton pigment particles at λ 1 and λ 2 , respectively. λ 1 should be located at the peak of a ph according to the principle of the standard three-band model. λ 1 is close to λ 2 , λ 2 ≥ 690 nm, λ 3 is located in the NIR band, and λ 3 ≥ 730 nm.
By analyzing the spectral characteristics of the water in Baiyangdian Lake and the commonly used bands of the three-band model, it was found that when λ 1 = 679 nm, λ 2 = 705 nm, and λ 3 = 756 nm, the measured Chl-a concentration in Baiyangdian had the highest correlation with the model, and the correlation coefficient was 0.7.

Four-Band Model
Due to the complex optical properties of inland eutrophic lakes, the assumptions of the three-band model could not simultaneously hold true. A new NIR band was introduced, i.e., the four-band model, to reduce the effect of the absorption by pure water and nonpigmented particles on the estimation of the Chl-a concentration [27].
Based on the characteristics of the spectral curve for the water in Baiyangdian Lake, the selection of the bands in the commonly used four-band model, and the spectral channels of ZY1-02D used, in the study, the best bands were determined to be λ 1 = 679 nm, λ 2 = 696 nm, λ 3 = 739 nm, and λ 4 = 705 nm through repeated experiments and comparisons.

Normalized Difference Chlorophyll Index (NDCI)
Sachidananda Maisha et al. proposed a novel index for estimating the chlorophyll-a concentration, the Normalized Difference Chlorophyll Index (NDCI) [29]. The NDCI index is constructed by applying the remotely sensed reflectance in both the λ 1 and λ 2 bands (Equation (9)) in a normalized form, which can partially eliminate the effects of the solar altitude angle, atmospheric radiation, and other factors. Among them, λ 1 and λ 2 are chosen at the characteristic spectra and their optical properties are dominated by chlorophyll-a. This also reduces the influence of the other water components to a certain extent.
Based on the characteristics of the spectral curve for the water in Baiyangdian Lake and the spectral channels of ZY1-02D used, the best bands were determined to be λ 1 = 662 nm and λ 2 = 705 nm through repeated tests and comparisons.

Fluorescence Line Height
The Fluorescence Line Height (FLH) model proposed by Neville and Gower is used to estimate the chlorophyll-a concentrations of various types of water bodies [30], and the basic principle of the FLH found using the model is to use the line connecting the off-water irradiance of the left and right channels of the fluorescence peak to the baseline, to calculate the distance between the off-water irradiance of the fluorescence channel and the baseline (Figure 7), and to establish the correlation between the chlorophyll-a mass concentration and the FLH. Using this method, the correlation between the chlorophyll-a mass concentration and the FLH was established. The FLH model is described by Equation (10): water irradiance of the left and right channels of the fluorescence peak to the baseline, to calculate the distance between the off-water irradiance of the fluorescence channel and the baseline (Figure 7), and to establish the correlation between the chlorophyll-a mass concentration and the FLH. Using this method, the correlation between the chlorophyll-a mass concentration and the FLH was established. The FLH model is described by Equation (10):  Based on the spectral curve characteristics of the water in Baiyangdian Lake and the spectral channel of ZY1-02D, λ max = 705 nm, λ 1 = 679 nm, and λ 2 = 748 nm.

Model Construction
Using the stepwise regression method in the SPSS22 software (gradually eliminating independent variables and increasing the correlation coefficient), a single-band model, ratio model, three-band model, and four-band model were constructed based on the measured Chl-a concentration data and the narrow band of the ZY1-02D AHSI images. The normalized chlorophyll-a difference index and the FLH model were used to establish linear, polynomial, and exponential models. Finally, the fitting equation and R 2 values were obtained. Figure 8 shows the estimated results of the different models. Among them, the fluorescent line height model has the best estimation accuracy, with y = 0.028x + 8.023 and R 2 = 0.84.

Accuracy Assessment
The Leave-One-Out Crossover (LOOCV) method was used to verify the accuracy of the inversion. This verification method has been proven to be effective in evaluating the inductive characteristics of statistical models, and the result is almost unbiased [8]. This method involves taking out one of the N samples s, the remaining N−1 samples are used to design the prediction model, and then, the samples taken are used for the testing. This process is repeated N times, the tests are conducted N times, and the average forecast deviation is counted.
The SPSS Statistics (IBM) software was used for all statistical analysis, including calculating the average, maximum, and minimum values, correlation analysis, and linear and non-linear regression. In addition, several indicators were used to evaluate the performance and uncertainty of the algorithm, including (1) the coefficient of determination (R 2 ), (2) the Root Mean Square Error (RMSE), and (3) the Mean Absolute Percentage Error (MAPE). Among them, the R 2 value represents the degree of correlation between the coupled data. The RMSE represents the difference between the predicted value and the actual value, and a smaller RMSE indicates a higher prediction accuracy. The MAPE is used to measure the accuracy of the prediction, and a smaller MAPE value indicates a better modeling result. These accuracy standards are defined in Equations (11) and (12).
where n is the number of samples, and Chl_a i and Chl_a i are the measured and predicted Chl-a concentrations.

Model Validation
Based on the ground-truth (i.e., measured) Chl-a concentrations and the satellite spectral profiles of the corresponding image elements, the inverse model construction and accuracy validation were carried out using the five models introduced in Section 3.2. The obtained model formulas, R 2 values, and error indexes of the validation set are presented in Table 3. The RMSE of the fluorescence height model was 15.55 µg/L, and the MAPE was 16.31%. Figure 9 shows that the scatter points of the measured and inverse Chl-a concentrations were evenly distributed on both sides of the 1:1 line, indicating that the constructed FLH model can be used to invert the Chl-a concentration of Baiyangdian Lake.

Inverted Chlorophyll-a Results
Based on the FLH model, the chlorophyll-a concentrations of the Baiyangdian water bodies were inverted, and the distribution of the chlorophyll-a concentration in Baiyangdian Lake was obtained. The maximum value of the inversion results was 15.40 μg/L, and the minimum value was 2.05 μg/L. Figure 10a-c shows the spatial distribution of the chlorophyll-a concentration in Baiyangdian. The difference in the chlorophyll-a concentrations in different areas was obvious, as well as the establishment of a wet-land in the inlet of Fu River. The natural reserve at the mouth of the Fuhe River was established to purify the water, the chlorophyll-a concentration was low, and the water quality was good. The spatial distribution of chlorophyll-a concentrations in Shaoche Lake and Laowang Lake shows a low center and high level on both sides. Jinlong Lake has fishponds and perimeter dikes. Figure 10 shows that the chlorophyll-a concentrations of the fishponds inside the perimeter dike were significantly different from those outside the perimeter dike, and the chlorophyll-a concentrations of the fishponds were lower after treatment than before treatment, indicating at the spatial distribution of the chlorophyll-a concentration in the study area based on the inversion of ZY1-02D hyperspectral data is relatively reasonable.

Inverted Chlorophyll-a Results
Based on the FLH model, the chlorophyll-a concentrations of the Baiyangdian water bodies were inverted, and the distribution of the chlorophyll-a concentration in Baiyangdian Lake was obtained. The maximum value of the inversion results was 15.40 µg/L, and the minimum value was 2.05 µg/L. Figure 10a-c shows the spatial distribution of the chlorophyll-a concentration in Baiyangdian. The difference in the chlorophyll-a concentrations in different areas was obvious, as well as the establishment of a wet-land in the inlet of Fu River. The natural reserve at the mouth of the Fuhe River was established to purify the water, the chlorophyll-a concentration was low, and the water quality was good. The spatial distribution of chlorophyll-a concentrations in Shaoche Lake and Laowang Lake shows a low center and high level on both sides. Jinlong Lake has fishponds and perimeter dikes. Figure 10 shows that the chlorophyll-a concentrations of the fishponds inside the perimeter dike were significantly different from those outside the perimeter dike, and the chlorophyll-a concentrations of the fishponds were lower after treatment than before treatment, indicating at the spatial distribution of the chlorophyll-a concentration in the study area based on the inversion of ZY1-02D hyperspectral data is relatively reasonable.

ZY1-02 -Based and Sentinel-Based Chl-a Products
Sentinel-2 carries a high spatial (10-20 m) and temporal (2-3 days) resolution of MSI and contains a sensitive "red-edge band", which is a sensitive band for inversion of chlorophyll-a concentration by fluorescence line height model. In this study, the sentinel-2 data on 4 October 2020 were selected, and the hyperspectral inversion model (FLH) was applied to the sentinel-2 data to obtain the sentinel-2 multispectral inversion results. The inversion results were simulated and verified with the inversion results of the measured chlorophyll-a concentration, and the verification R 2 = 0.66 was obtained ( Figure 11). This further illustrates the advantage of hyperspectral inversion of chlorophyll-a concentration.

Uncertainties and Limitations
In this study, there are two main factors to improve the accuracy of chlorophyll-a concentration assessment in Baiyangdian: the advantage of the spectral resolution of ZY1- Figure 10. Results of the chlorophyll-a inversion for Baiyangdian Lake. (a) Shaochedian Lake; (b) Laowangdian Lake; (c) Jinlongdian Lake.

ZY1-02 -Based and Sentinel-Based Chl-a Products
Sentinel-2 carries a high spatial (10-20 m) and temporal (2-3 days) resolution of MSI and contains a sensitive "red-edge band", which is a sensitive band for inversion of chlorophyll-a concentration by fluorescence line height model. In this study, the sentinel-2 data on 4 October 2020 were selected, and the hyperspectral inversion model (FLH) was applied to the sentinel-2 data to obtain the sentinel-2 multispectral inversion results. The inversion results were simulated and verified with the inversion results of the measured chlorophyll-a concentration, and the verification R 2 = 0.66 was obtained ( Figure 11). This further illustrates the advantage of hyperspectral inversion of chlorophyll-a concentration.  Figure 10. Results of the chlorophyll-a inversion for Baiyangdian Lake. (a) Shaochedian Lake; (b) Laowangdian Lake; (c) Jinlongdian Lake.

ZY1-02 -Based and Sentinel-Based Chl-a Products
Sentinel-2 carries a high spatial (10-20 m) and temporal (2-3 days) resolution of MSI and contains a sensitive "red-edge band", which is a sensitive band for inversion of chlorophyll-a concentration by fluorescence line height model. In this study, the sentinel-2 data on 4 October 2020 were selected, and the hyperspectral inversion model (FLH) was applied to the sentinel-2 data to obtain the sentinel-2 multispectral inversion results. The inversion results were simulated and verified with the inversion results of the measured chlorophyll-a concentration, and the verification R 2 = 0.66 was obtained ( Figure 11). This further illustrates the advantage of hyperspectral inversion of chlorophyll-a concentration. Figure 11. Scatter plot of in situ versus sentinel-2-derived Chl-a concentrations.

Uncertainties and Limitations
In this study, there are two main factors to improve the accuracy of chlorophyll-a concentration assessment in Baiyangdian: the advantage of the spectral resolution of ZY1-

Uncertainties and Limitations
In this study, there are two main factors to improve the accuracy of chlorophyll-a concentration assessment in Baiyangdian: the advantage of the spectral resolution of ZY1-02D and the chlorophyll-a concentration dataset of simultaneous measurements. However, there are limitations in the model assessment. First, the acquisition of the chlorophyll-a dataset is mainly focused on autumn, the next step needed is to acquire synchronous chlorophyll-a concentration datasets in spring and summer by star field to validate the model. Second, it is found that the chlorophyll-a inversion results obtained using ZY1-02D data are affected by the image noise and remote sensing reflectance accuracy. The image-based remote sensing reflectance estimation method for water surface in this study produced good results, but more accurate remote sensing reflectance acquisition is needed, which will further improve the accuracy of chlorophyll-a concentration inversion. Further development of the high-precision atmospheric correction method for ZY1-02D water body images is also needed. Finally, due to the long revisit period (approximately 55 days) and low temporal resolution of the ZY1-02D single star, it is difficult to meet the demand for multi-frequency water monitoring, the next work can increase the time series with the help of multispectral image data.

Conclusions
The spectral resolution of ZY1-02D is 10 nm within 400-900 nm. Its high spectral resolution is more likely to reflect the sensitive spectrum of the water-color absorption and reflection characteristics and to achieve high-precision chlorophyll-a concentration inversion. In this study, ZY1-02D hyperspectral data and multi-spectral data were used to clarify the surface water area to obtain the clear water surface area in Baiyangdian and to monitor and evaluate the chlorophyll-a concentration of the open water in Baiyangdian at the optimal wavelength band. Then, based on the semi-empirical method for determining the Baiyangdian chlorophyll-a concentration, multi-model regression fitting was performed. Finally, the FLH model was determined as having the highest inversion accuracy. Applying the inversion model to Baiyangdian Lake and its evaluation improved the inversion accuracy of the chlorophyll-a concentration of the water bodies, and becomes of great significance for the monitoring and evaluation of the chlorophyll-a concentrations of the Baiyangdian water bodies.