Simulation and Assessment of the Capabilities of Orbita Hyperspectral (OHS) Imagery for Remotely Monitoring Chlorophyll-a in Eutrophic Plateau Lakes

: The chlorophyll-a (Chl-a) concentration of eutrophic lakes ﬂuctuates signiﬁcantly due to the disturbance of wind and anthropogenic activities on the water body. Consequently, estimation of the Chl-a concentration has become an immense challenge. Due to urgent demand and rapid development in high-resolution earth observation systems, it has become crucial to assess hyperspectral satellite imagery capabilities on inland water monitoring. The Orbita hyperspectral (OHS) satellite is the latest hyperspectral sensor with both high spectral and spatial resolution (2.5 nm and 10 m, respectively), which could provide great potential for remotely estimating the concentration of Chl-a for inland waters. However, there are still some deﬁciencies that are mainly manifested in the Chl-a concentration remote sensing retrieval model assessment and accuracy validation, as well as signal-to-noise ratio (SNR) estimation of OHS imagery for inland waters. Therefore, the radiometric performance of OHS imagery for water quality monitoring is evaluated in this study by comparing different atmospheric correction models and the SNR with several remote sensing images. Several crucial ﬁndings can be drawn: (1) the three-band model ((1/B15-1/B17)B19) developed by OHS imagery is most suitable for estimating the Chl-a concentration in Dianchi Lake, with the root-mean-square error (RMSE) and the mean absolute percentage error (MAPE) of 15.55 µ g/L and 16.31%, respectively; (2) the applicability of the FLAASH (Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes) atmospheric correction model for OHS imagery in a eutrophic plateau lake (Dianchi Lake) was better than the 6S (Second Simulation of Satellite Signal in the Solar Spectrum) model, and QUAC (Quick Atmospheric Correction) model, as well as the dark pixel method; (3) the SNR of the OHS imagery was similar to that of Hyperion imagery and was signiﬁcantly higher than SNR of the HSI imagery; (4) the spatial resolution showed slight inﬂuence on the SNR of the OHS imagery. The results show that OHS imagery could be applied to remote sensing retrieval of Chl-a in eutrophic plateau lakes and presents a new tool for dynamic hyperspectral monitoring of water quality.


Introduction
Chlorophyll-a (Chl-a) is a crucial parameter which impacts the watercolor of inland lake water, and is a vital indicator measuring the eutrophication degree of lake water [1,2]. The nutrients transported by Chl-a reduce the spread of light through a water column and affect the natural operation of the entire aquatic ecosystem. Therefore, accurate prediction of Chl-a concentration and its spatiotemporal distribution pattern is of great significance to the protection of lake ecosystems and the improvement of water quality monitoring capabilities, and is the basis of water environment remote sensing monitoring and evaluation. Precise prediction of Chl-a concentration is also the focus and challenge of watercolor remote sensing [3][4][5][6].
Dianchi Lake is the largest freshwater lake in Southwest China and the sixth-largest freshwater lake in China. It is an important drinking water source for residents, one of many important habitats for migratory birds, and an essential part of the wetland ecosystem. Due to continuous urbanization, Dianchi Lake has become a "sewage bucket" for domestic sewage, industrial wastewater, and agricultural sewage. The eutrophication of water bodies has become increasingly serious, and this lake has matured into a typical representative of eutrophic lakes in the country. The pollution of the Dianchi Lake has caused an undesirable evolution of the biological population structure and aquatic ecosystems. The resulting cyanobacterial blooms have greatly impacted the production and life of local inhabitants and are a major threat to the ecological environment and biological diversity [7][8][9][10]. Hence, it is particularly important to judiciously monitor the eutrophication level of Dianchi Lake and consider Dianchi Lake as a chief target for eutrophication control.
This eutrophic plateau lake is greatly affected by the disturbance of wind and anthropogenic activities, making it difficult to retrieve the Chl-a concentration. Moreover, most previous studies focused on multi-spectral remote estimation [11][12][13][14][15][16], with few studies for hyperspectral inversion models in eutrophic plateau lake existing [17]. O'Reilly et al. [18] used the blue-green band ratio model to estimate the concentration of Chl-a in the ocean based on Sea-Viewing Wide Field-of-View Sensor (SeaWiFS) images. Gurlin et al. [19] proposed the near-infrared and red band ratio model and applied it to Moderate Resolution Imaging Spectroradiometer (MODIS) and Medium Resolution Imaging Spectrometer (MERIS) images to estimate the Chl-a concentration. Matthews et al. [20] utilized the maximum peak height algorithm to estimate the Chl-a concentration of inland waters. Liu et al. [21] executed semi-analytical algorithms to retrieve the Chl-a concentration of various lakes based on Ocean and Land Color Instrument (OLCI) images. Ioannou et al. [22] implemented neural network algorithms to predict the Chl-a concentration and used them on MODIS images. Due to the complexity of the optical properties of water bodies and the diverse performance of each image, the Chl-a concentration retrieval accuracy was significantly different. Hence, the current dilemmas mainly include the following: (1) the accuracy of multi-spectral inversion models for eutrophic plateau lakes needs improvement; (2) while there is less research on hyperspectral data, its inversion models results may be more accurate compared to multi-spectral data; (3) there is no research on the applicability of OHS (Orbita hyperspectral) data for eutrophic plateau lakes; (4) the problems of atmospheric correction and SNR (signal-to-noise ratio) of OHS imagery have not been resolved.
Retrieving the concentration of Chl-a in eutrophic plateau lakes using the OHS image was a major challenge. This challenge was mainly addressed because OHS imagery has rarely been used in the field of remote sensing and the uncertainty of OHS imagery for Chl-a concentration estimation is still unknown. OHS imagery is an imaging product integrating high spectral resolution, high spatial resolution, and high temporal resolution, and is a hyperspectral satellite with great development potential. Most existing hyperspectral imagers can no longer be used. For instance, the Hyperion hyperspectral sensor mounted on the EO-1 (Earth Observing-1) satellite platform and the AVIRIS (Airborne Visible Infrared Imaging Spectrometer) hyperspectral imager positioned on the ER-2 (European Remote sensing-2) aircraft platform has been terminated; the HJ-1 HSI (Hyperspectral Imager) and the images produced by the CHRIS (Compact High-Resolution Imaging Spectrometer) sensor mounted on the PROBA (Project for On-Board Autonomy) satellite cannot be used due to stripe noise. In addition, the Hyperspectral Precursor of the Application Mission (PRISMA) satellite launched in March 2019 offers great potential for hyperspectral monitoring of water quality in inland lakes, but its imaging products have certain uncertainties [23]; the DESIS (DLR Earth Sensing Imaging Spectrometer) hyperspectral sensor installed on the International Space Station can play a crucial role in monitoring coastal and inland water bodies, but has difficulty monitoring small-scale lake water bodies due to the low spatial resolution [24]; the HISUI (Hyperspectral Imager Suite) hyperspectral instrument onboard the ALOS-3 (Advanced Land Observing Satellite-3) satellite was successfully launched in December 2019, but the quality of its imaging products and its ability to monitor inland water bodies are still unknown. Therefore, the emergence of the OHS satellite provides opportunities for hyperspectral remote sensing monitoring of small-scale inland lakes. This study assesses the performance of OHS imagery for remote sensing inversion of water quality parameters from the perspective of inland lake water quality monitoring.
This study aims to quantitatively reveal and analyze the spatial pattern of Chl-a concentration in eutrophic plateau lakes using OHS imagery in Dianchi Lake, with the following main objectives: (1) Develop a model based on the novel OHS imagery for Chl-a remote sensing retrieval in eutrophic plateau lakes; (2) Explore the applicability and feasibility of OHS imagery for monitoring Chl-a concentration in eutrophic plateau lakes.

Study Area
Located on the Yunnan-Guizhou Plateau, Dianchi Lake (102 Figure 1) is the largest freshwater lake in Yunnan Province and the largest plateau lake in the southwest region of China, with an altitude of 1886 m, a lake area of 330 km 2 , and an average water depth of 5 m. It is a typical plateau lake with eutrophication and is significantly affected by wind and human activities. Dianchi Lake is a vital source of drinking water for local inhabitants. It also functions as a reservoir for industry and agriculture, flood control, shipping, tourism, aquaculture, hydroelectric power generation, and meteorological regulation, and plays a critical role in the development of Kunming's social economy. Over recent decades, the rise of industry, rapid social and economic development, and constant population growth around Dianchi Lake have caused more and more pollutants access to Dianchi Lake, resulting in increasingly severe eutrophication. In current years, although the water quality of Dianchi Lake has improved, it is still in Level-IV (extremely eutrophic water body according to Dianchi Lake Administration Bureau). The management and ecological restoration of Dianchi Lake is threatened.

Remote Sensing Images and Images Pre-Processing
Five types of imagery data were used in this study, including (1) the OHS imagery on April 2, 2019, with both high spectral and spatial resolution (2.5 nm and 10 m, respectively), a band range of 400 nm-1000 nm, and a total of 32 spectral bands [25] (Table A1)

Radiometric Calibration
Radiometric calibration is the process of converting the digital number value (DN) of remote sensing images into top-of-atmosphere (TOA) radiance or apparent reflectance. In this study, radiometric calibration was performed using Environment for Visualizing Images (ENVI) software to convert its DN value into apparent radiance via [25][26]. First, we performed radiometric calibration on OHS imagery in this study.
where, Le is the apparent radiance, gain is the gain coefficient, offset is the offset coefficient, and TDIStage is the integration series (Table A1).

Atmospheric Correction
Atmospheric correction eliminates and corrects the absorption or scattering of sunlight caused by the atmosphere, the reflection from the target or absorption and scattering

Remote Sensing Images and Images Pre-Processing
Five types of imagery data were used in this study, including (1) the OHS imagery on 2 April 2019, with both high spectral and spatial resolution (2.5 nm and 10 m, respectively), a band range of 400 nm-1000 nm, and a total of 32 spectral bands [25] (Table A1) Among them, OHS imagery was used for all operations in this study, and other imageries were used to estimate the SNR.

Radiometric Calibration
Radiometric calibration is the process of converting the digital number value (DN) of remote sensing images into top-of-atmosphere (TOA) radiance or apparent reflectance. In this study, radiometric calibration was performed using Environment for Visualizing Images (ENVI) software to convert its DN value into apparent radiance via [25,26]. First, we performed radiometric calibration on OHS imagery in this study.
where, L e is the apparent radiance, gain is the gain coefficient, offset is the offset coefficient, and TDIStage is the integration series (Table A1).

Atmospheric Correction
Atmospheric correction eliminates and corrects the absorption or scattering of sunlight caused by the atmosphere, the reflection from the target or absorption and scattering of scattered light, and the incident light to the sensor in addition to the light from the target. Its basic algorithm is given in Equation (2) [26]. We performed atmospheric correction on OHS imagery in this study.
where, L is the total radiance received by the pixel at the sensor, ρ is the real surface reflectance of the pixel, ρ e is the average surface reflectance around the pixel, S is the atmosphere spherical albedo, and L a is the atmospheric backscattered radiance (atmospheric range radiation); A and B are two coefficients depending on atmospheric conditions and geometric conditions. In this study, the FLAASH (Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes) model, QUAC (Quick Atmospheric Correction) model, 6S (Second Simulation of Satellite Signal in the Solar Spectrum) model, and dark pixel method were used to perform atmospheric correction on the OHS image. Among them, FLAASH is an atmospheric correction module based on the MODTRAN5 (Moderate Resolution Atmospheric Transmission) radiation transfer model, which can be applied to hyperspectral and multispectral images; QUAC model is a quick atmospheric correction method incorporated in ENVI software; 6S model is used to simulate the radiation transfer process of 0.25~4.00 µm solar reflectance wave band and is an improved version of the 5S (Simulation of Satellite Signal in the Solar Spectrum) model; dark pixel method assumes that the radiance of water body dark pixel in the infrared band is approximately 0 and the radiation received by the sensor in this band comes from atmospheric effects from which the atmospheric parameters and aerosol factors are derived. The visible bands are extrapolated and atmospherically corrected [26].

Field Measurements
Three cruise field surveys conducted on 19 to 20 September 2009, 13 to 16 April 2017, and 13 to 15 November 2017 over Dianchi Lake measured and collected remote sensing reflectance (R rs ) and water sample data. We collected 93 samples in 3 field surveys ( Figure 1). At each station, the longitude and latitude coordinates were recorded using a Global Positioning System (GPS) receiver and roughly 1000-1500 mL surface water samples from a 50 cm water depth were collected for laboratory measurement and analysis [27,28]. In addition, we obtained the measured Chl-a concentration data collected by Dianchi Administration Bureau on 2 April 2019. Among them, 93 observation stations datasets from three field measurements (19 to 20 September 2009, 13 to 16 April 2017, and 13 to 15 November 2017) were used to develop and validate the Chl-a retrieval models; the in situ Chl-a concentration data (N = 10) of Dianchi Lake synchronized with OHS imagery (2 April 2019) were used to validate the performance of the OHS imagery retrieval results (Table 1).

. Radiometric Measurements
Spectral reflectance was measured using an ASD FieldSpec spectroradiometer. The instrument has 512 spectral channels with a spectral resolution of 1.5 nm across a spectral range of 350-1050 nm. The radiance spectra of the reference panel, water, and sky were measured using the above-water measurement method followed by the Ocean Optics Protocols [29,30]. The observation geometry with an azimuth of 135 • and a nadir angle of 40 • was used to avoid direct sunlight reflection on the water [30]. Ten spectra were obtained for each sampling site, from which we eliminated abnormal spectra due to occasional factors and averaged valid ones to calculate the remote sensing reflectance (R rs (λ)) via the following equation [30,31] (Figure 2): where, L t is the measured total radiance from the water surface; ρ is the skylight reflectance received from the air-water surface, which is considered as 2.2% for calm weather, 2.5% for wind speed reaching 5 m·s −1 , and 2.6%-2.8% for wind speed of 10 m·s −1 ; L sky is the determined radiance of the sky; L g is the radiance measured with the reference panel; R g is the reflectance from the reference panel, and its value is 30%.

Radiometric Measurements
Spectral reflectance was measured using an ASD FieldSpec spectroradiometer. T instrument has 512 spectral channels with a spectral resolution of 1.5 nm across a spectr range of 350-1050 nm. The radiance spectra of the reference panel, water, and sky we measured using the above-water measurement method followed by the Ocean Optics Pr tocols [29,30]. The observation geometry with an azimuth of 135° and a nadir angle of 4 was used to avoid direct sunlight reflection on the water [30]. Ten spectra were obtain for each sampling site, from which we eliminated abnormal spectra due to occasional fa tors and averaged valid ones to calculate the remote sensing reflectance (Rrs(λ)) via t following equation [30,31] (Figure 2): where, Lt is the measured total radiance from the water surface; ρ is the skylight reflectan received from the air-water surface, which is considered as 2.2% for calm weather, 2.5 for wind speed reaching 5 m•s −1 , and 2.6%-2.8% for wind speed of 10 m•s −1 ; Lsky is the d termined radiance of the sky; Lg is the radiance measured with the reference panel; Rg the reflectance from the reference panel, and its value is 30%.

Water Sample Analysis
Water samples were filtered using Whatman GF/F fiberglass filters, and phytoplan ton pigments were extracted in 90% hot ethanol at 80 °C. Chl-a concentration was mea ured using hot ethanol methods [32], that is, the concentration of Chl-a was extracted wi 90% hot ethanol at 80 °C, the resulting extract was acidified with 1% dilute hydrochlor acid, and the absorbance at 665 nm and 750 nm was measured by a UV2550 spectroph tometer. The Chl-a concentration was determined using the absorbance at 665 nm and 7

Water Sample Analysis
Water samples were filtered using Whatman GF/F fiberglass filters, and phytoplankton pigments were extracted in 90% hot ethanol at 80 • C. Chl-a concentration was measured using hot ethanol methods [32], that is, the concentration of Chl-a was extracted with 90% hot ethanol at 80 • C, the resulting extract was acidified with 1% dilute hydrochloric acid, and the absorbance at 665 nm and 750 nm was measured by a UV2550 spectrophotometer. The Chl-a concentration was determined using the absorbance at 665 nm and 750 nm. The concentrations of total suspended matter (TSM), organic suspended matter (OSM), and inorganic suspended matter (ISM) were measured using detailed descriptions from American Public Health Association [33]. The water sample was filtered on a Whatman GF/F glass fiber filter and weighed according to the method of Huang [34]. The filter was then re-heated at 550 • C for 4 h to remove the organic ingredients and weighed again to determine ISM. The OSM concentration can be obtained by subtracting ISM from TSM.

Model Calibration Based on Simulated OHS Imagery
To estimate the Chl-a concentration using the OHS image, the spectral reflectance of the field measurements and the spectral response function of the OHS image were used by the atmosphere-surface scene simulation method to simulate the spectral characteristics of the in situ remote sensing reflectance to the spectral characteristics of OHS image. R rs (λ) was used to represent the remote sensing reflectance received by the water-surface reflecting sunlight and reaching the ASD sensor detector. After considering the spectral response function f (λ), the simulation result can be derived by the following formula [35,36]: where, R rs (B i ) denotes the simulated field-measured R rs (λ) for the i-th band of OHS image, with integration from λ min to λ max for the i-th band.

Signal-to-Noise Ratio Estimation
SNR is a measure of the relationship between signal and noise and is a particularly important indicator to evaluate the performance of a sensor. According to Tobler's first law of geography, adjacent features of the same kind have similar spectral characteristics, corresponding to a similar gray value in the image of the same band. The SNR of the OHS image was determined using the variance method [37,38]. However, the premise of the method must ensure that the pixels involved in the operation were homogeneous features. The central area of Dianchi Lake selected in this study was rather homogeneous, and it was assumed that their bio-optical characteristics were spatially stable.
The specific process was: (1) extracted a pure window, that is, an area with relatively homogeneous optical properties in the center of Dianchi Lake; (2) determined the average DN value and standard deviation of all pixels in the pure window; (3) the ratio of the mean value and the standard deviation was SNR of remote sensing imagery. The formula of SNR is as follows: where, MDN is the mean DN value of the pixels, DN i is the DN value of each pixel in the image of the selected area, N is the number of pixels, S is the standard deviation of the DN value of the pixels, and SNR is the signal-to-noise ratio of the OHS image.

Statistical Analysis
Statistical analyses including calculations of the mean, maximum, and minimum values, correlation analysis, and linear and non-linear regressions were performed using the Microsoft Excel software. Additionally, several metrics were used to assess the algorithm's performance and uncertainties in this study, which include (1) the Pearson correlation coefficient (r), (2) the determination coefficient (R 2 ), (3) the root-mean-square error (RMSE), (4) the mean absolute percentage error (MAPE), and (5) p-value that represents the significance level of the association between two variables, and was reported to be significant (p < 0.05) or not significant (p > 0.05) with a t-test. Among them, the r and R 2 indicate the degree of correlation between the coupled data; the RMSE represents the difference between the predicted value and the actual value, smaller RMSE denotes higher prediction accuracies; and the MAPE is a measure of the prediction accuracy of a forecasting method, smaller MAPE indicates better modeling results. These accuracy criteria are defined as [39,40]: where, n refers to the number of samples, and Chl-a i and Chl-a i ' are measured and predicted values of Chl-a concentration, respectively.

Data Descriptive Statistics
The water quality parameters obtained from field sampling sites showed a high dynamic range and considerable variability ( Table 2). The mean Chl-a during the sampling period was 87.35 µg/L, with a maximum of 187.01 µg/L and a minimum of 38.97 µg/L. TSM ranged from 20.98 to 66.60 mg/L (mean = 37.04 mg/L, standard deviation = 8.19 mg/L). OSM and ISM varied from 11.16 to 52.10 mg/L and 0.00 to 28.10 mg/L, with a mean value of 23.54 mg/L and 13.50 mg/L. A total of 93 sample data were randomly divided into two groups: 61 of them were used as training samples for the development of the algorithm, and the remaining 32 samples were used as the validation data set to evaluate the accuracy of the algorithm [41].  This band ratio model can eliminate the atmospheric influence to some extent, and can also reduce the interference of the roughness from the water surface [42]. The blue-green band ratio model is mainly used for Case-I water bodies with phytoplankton as the dominant. The blue-green band ratio model commonly adopted include R rs (490) > R rs (510)/R rs (560), R rs (490)/R rs (560), etc. We referred to the R rs (490)/R rs (560) model in this study, the reflectance ratio factor B2/B7 was used as the independent variable, and the corresponding Chl-a concentration was the dependent variable: where, B2 and B7 represent the simulated OHS-based remote sensing reflectance at spectral channels of 2nd and 7th band, respectively (the following spectral bands expressions are similar to this), and coefficient a = −154.84 and b = 156.71 are determined empirically.

Near-Infrared and Red (NIR-Red) Band Ratio Model
The near-infrared and red band ratio model is commonly used in optically complex Case-II water bodies mainly affected by Chl-a concentration, total suspended solids, and colored dissolved organic matter [19]. The generally used near-infrared and red band ratio models are R rs (700)/R rs (670), R rs (719)/R rs (667), R rs (706)/R rs (682), R rs (708)/R rs (665), etc. We referred to the R rs (708)/R rs (665) model in this study: where, coefficient a = 56.226 and b = 0.2191.

Three-Band Model
The three-band model is based on bio-optical theory and is used for remote sensing retrieval of Chl-a concentration in turbid Case-II water bodies, and its manifestation is [R rs −1 (λ 1 ) − R rs −1 (λ 2 )] R rs (λ 3 ), the model is based on three assumptions [43]: (1) the backscattering coefficients of the three bands are approximately equal, that is, represent the backscattering coefficients at λ 1 , λ 2 , λ 3 , respectively; (2) λ 3 is dominated by the absorption coefficient of pure water, and the absorption coefficients in-water constituents could be ignored, i.e., a w (λ 3 ) >> a nap (λ 3 ) + a COOM (λ 3 ) + b b (λ 3 ), where a w (λ 3 ), a nap (λ 3 ), a COOM (λ 3 ) are the absorption coefficient of pure water, non-pigmented particulate matter and colored dissolved organic matter at λ 3 , respectively; (3) the absorption coefficients of 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, namely, a nap (λ 1 ) ≈ a nap (λ 2 ), a COOM (λ 1 ) ≈ a COOM (λ 2 ) and a ph (λ 1 ) >> a ph (λ 2 ) where a ph (λ 1 ) and a ph (λ 2 ) are the absorption coefficient of phytoplankton pigment particles at λ 1 and λ 2 , respectively. Regularly used band combination forms of three-band model are [ (708), etc. By analyzing the optical properties of the Dianchi Lake waters, we found that the above conditions were satisfied when λ 1 = 680 nm, λ 2 = 710 nm, and λ 3 = 745 nm, and the three-band model used in this study is: where, a = 137.35 and b = 59.741.

Four-Band Model
Due to the optically complex properties of inland eutrophic lakes, the assumptions in the three-band model cannot be established simultaneously. For this reason, a new nearinfrared band was introduced in the four-band model to reduce the impact of the absorption of pure water and non-pigmented particles on estimation of the Chl-a concentration [44]. We referred to the method described by Le et al. [44] for the four-band model, namely, where, a = −0.0002 and b = 89.498.

Fluorescence Line Height (FLH) Model
The fluorescence line height model was first proposed by Neville & Gower [45] to be applied to the estimation of Chl-a concentration in various types of water. This method is based on the relationship between the peak distance from the baseline and the Chl-a concentration. The fluorescence line height model is as follows: Chl-a_FLH = a × FLH + b where, a = −11510 and b = 75.573; λ B14 , λ B15 , and λ B17 are central wavelength of 14th, 15th, and 17th band of OHS imagery, respectively.

Model Validation Based on Simulated-OHS Imagery
Model validation was used to assess the effect of the Chl-a concentration estimation model and the existing uncertainty and can be used to evaluate whether the retrieval algorithm has practical application value. In this study, the root mean square error (RMSE) and the mean absolute percentage error (MAPE) were used to assess the accuracy of the Chl-a concentration estimation model via Formulas (8) and (9), and to analyze the accuracy difference and uncertainty among different algorithms ( Figure 3 and Table 3). Table 3. Estimation model of Chl-a concentration based on simulated OHS imagery. R 2 , RMSE, and MAPE represent the determination coefficient, root-mean-square error, and mean absolute percentage error, respectively. As shown in Figure 3 and Table 3, the effect of the four-band model was the worst (RMSE = 24.44 µg/L, MAPE = 28.17%); the three-band model had the highest accuracy (RMSE = 15.55 µg/L, MAPE = 16.31%) with the error within the acceptable range, which was suitable for the remote estimation of the Chl-a concentration in Dianchi Lake. The results of the estimation model were validated based on the simulated OHS imagery.  . Scatter-plots of Chl-a model development (a-i) between measured in situ Chl-a data and simulated OHS-based derived R rs , and Chl-a model validation (b-j) between measured and simulated OHS-based derived Chl-a through Equations (10)- (15). Chl-a is chlorophyll-a concentration and OHS is OHS imagery.

Validation of Derived Chl-a for OHS Imagery
To further assess the effectiveness of the three-band algorithm for the OHS image, the measured Chl-a concentration data synchronized with OHS imagery from 10 observation stations evenly distributed in Dianchi Lake on 2 April 2019 (OHS imaging date) were used for model evaluation. In order to determine the influence of atmospheric correction on Chl-a concentration inversion, we used satellite-ground synchronization data to validate the Chl-a concentration retrieval results obtained by applying different atmospheric correction algorithms such as the FLAASH model, 6S model, QUAC model, and dark pixel method (Figure 4). Among them, the FLAASH model had the best validation effect (m MAPE = 8.26%) and could be applied to the atmospheric correction of the OHS imagery of Dianchi Lake, while the Chl-a concentration obtained by the 6S atmospheric correction model was significantly higher than the in situ Chl-a concentration (RMSE = 31.49 µg/L, MAPE = 56.86%), and was not suitable for the OHS imagery of Dianchi Lake. The Chl-a concentration estimated based on the three-band model had a significant correlation with the measured Chl-a concentration (R 2 = 0.858, p < 0.05) with RMSE of 5.71 µg/L and MAPE of 8.26% of the relative difference ( Figure 4). Comparisons between the in situ Chl-a measured and the Chl-a estimated using the proposed three-band model in this study showed that these values were acceptable. These results indicate that the three-band model could be used with satisfactory performance to retrieve Chl-a concentration in this inland water body.

Spatial Patterns of Chl-a in Dianchi Lake
The three-band model was applied to the OHS imagery to obtain the spatial distribution pattern of Chl-a concentration in Dianchi Lake (Figure 5a), which revealed that the concentration of Chl-a in Dianchi Lake showed an increasing trend from the center of the lake to the coast of the lake. To quantitatively reveal the spatial distribution characteristics of Chl-a concentration in Dianchi Lake, three cross-section lines were randomly selected to analyze the variation of Chl-a concentration with the distance from the lakeshore. Fig-Figure 4. Comparison of OHS-derived and measured Chl-a for OHS imagery. (a-d) respectively represent the validation of Chl-a concentration inversion results after using different atmospheric correction algorithms such as the FLAASH model, 6S model, QUAC model, and dark pixel method.

Spatial Patterns of Chl-a in Dianchi Lake
The three-band model was applied to the OHS imagery to obtain the spatial distribution pattern of Chl-a concentration in Dianchi Lake (Figure 5a), which revealed that the concentration of Chl-a in Dianchi Lake showed an increasing trend from the center of the lake to the coast of the lake. To quantitatively reveal the spatial distribution characteristics of Chl-a concentration in Dianchi Lake, three cross-section lines were randomly selected to analyze the variation of Chl-a concentration with the distance from the lakeshore. Figure 5b intuitively describes the overall decreasing trend of Chl-a concentration as the distance from the lake shore increases. The center of the lake had the lowest concentration of Chl-a, being 20~40 µg/L; the Chl-a concentration along the lake was the highest, reaching 100~180 µg/L; and the mean Chl-a concentration of the whole Dianchi Lake on 2 April 2019, was estimated to be 66.43 ± 12.98 µg/L.

Signal-to-Noise Ratio of OHS Imagery
In order to compare the SNR of the OHS image with other multispectral and hyperspectral images, Landsat-8 OLI (Operational Land Imager) (spatial resolution of 30 m), Sentinel-2 MSI (Multispectral Instrument, spatial resolution of 10 m, 20 m and 60 m), EO-1 Hyperion (Earth Observing-1, spatial resolution of 30 m), and HJ-1 HSI (Hyperspectral Imager, spatial resolution of 100 m) were selected to perform comparable analysis with OHS imagery. A comparatively large area with relatively stable bio-optical properties in Dianchi Lake was selected as the SNR study area (Figure 6a). Simultaneously, the OHS imagery was resampled to 30 m and 100 m to explore the effect of spatial resolution on the SNR of the OHS imagery. The SNR of the Landsat-8 OLI imagery was significantly higher than that of the OHS imagery, Hyperion imagery, HSI imagery, and Sentinel-2 imagery, but it gradually decreased with wavelength increase. The SNR of multispectral images was generally higher than that of hyperspectral images. The SNR of the HSI image in all bands was relatively low, indicating that as the spectral resolution advances, the noise decreases accordingly. The SNR of the Hyperion imagery gradually declined with the increase of the wavelength, which was similar to that of the OHS imagery. Although the SNR of OHS imagery after resampling has improved, the change was not obvious, reflecting that the spatial resolution had no significant effect on the OHS imagery ( Figure  6b,c, Table 4). This result is different from the change of multi-spectral image SNR with

Signal-to-Noise Ratio of OHS Imagery
In order to compare the SNR of the OHS image with other multispectral and hyperspectral images, Landsat-8 OLI (Operational Land Imager) (spatial resolution of 30 m), Sentinel-2 MSI (Multispectral Instrument, spatial resolution of 10 m, 20 m and 60 m), EO-1 Hyperion (Earth Observing-1, spatial resolution of 30 m), and HJ-1 HSI (Hyperspectral Imager, spatial resolution of 100 m) were selected to perform comparable analysis with OHS imagery. A comparatively large area with relatively stable bio-optical properties in Dianchi Lake was selected as the SNR study area (Figure 6a). Simultaneously, the OHS imagery was resampled to 30 m and 100 m to explore the effect of spatial resolution on the SNR of the OHS imagery. The SNR of the Landsat-8 OLI imagery was significantly higher than that of the OHS imagery, Hyperion imagery, HSI imagery, and Sentinel-2 imagery, but it gradually decreased with wavelength increase. The SNR of multispectral images was generally higher than that of hyperspectral images. The SNR of the HSI image in all bands was relatively low, indicating that as the spectral resolution advances, the noise decreases accordingly. The SNR of the Hyperion imagery gradually declined with the increase of the wavelength, which was similar to that of the OHS imagery. Although the SNR of OHS imagery after resampling has improved, the change was not obvious, reflecting that the spatial resolution had no significant effect on the OHS imagery (Figure 6b,c, Table 4). This result is different from the change of multi-spectral image SNR with spatial resolution [46]. The SNR of the OHS image generally decreased with increasing wavelength, and there was an obvious minimum near 760 nm ( Figure 6). As a result of the weak water signal in watercolor and the low SNR of high spatial resolution sensors, it was difficult to distinguish the subtle differences among signals, limiting its application in watercolor remote sensing to some extent and reducing the accuracy of Chl-a concentration estimation.

Limitation
The three-band model based on simulated OHS imagery proposed in this study has successfully retrieved the Chl-a concentration in eutrophic plateau lakes but still has the following limitations: (1) Limitations of the application area. Compared with the existing Chl-a concentration retrieval algorithms [19,[47][48][49][50][51][52] (Table A2 and A3, Figure A1-A4), the three-band model developed in this study is most suitable for remote sensing retrieval of Chl-a concentration in eutrophic plateau lakes. However, due to the optically complex properties of the water body, whether the three-band model proposed in this study is applicable to other optical water types such as Dongting Lake and Poyang Lake is still unknown. In addition, the three-band model developed in this study is an empirical algorithm, which appears to rely on measured data. Therefore, for different in situ data sets, the model coefficients may vary, and the model needs to be re-calibrated. (2) Limitations  The SNR of the OHS image generally decreased with increasing wavelength, and there was an obvious minimum near 760 nm ( Figure 6). As a result of the weak water signal in watercolor and the low SNR of high spatial resolution sensors, it was difficult to distinguish the subtle differences among signals, limiting its application in watercolor remote sensing to some extent and reducing the accuracy of Chl-a concentration estimation.

Limitation
The three-band model based on simulated OHS imagery proposed in this study has successfully retrieved the Chl-a concentration in eutrophic plateau lakes but still has the following limitations: (1) Limitations of the application area. Compared with the existing Chl-a concentration retrieval algorithms [19,[47][48][49][50][51][52] (Tables A2 and A3, Figures A1-A4), the three-band model developed in this study is most suitable for remote sensing retrieval of Chl-a concentration in eutrophic plateau lakes. However, due to the optically complex properties of the water body, whether the three-band model proposed in this study is applicable to other optical water types such as Dongting Lake and Poyang Lake is still unknown. In addition, the three-band model developed in this study is an empirical algorithm, which appears to rely on measured data. Therefore, for different in situ data sets, the model coefficients may vary, and the model needs to be re-calibrated. (2) Limitations of OHS imagery. This study only discusses the capability of OHS imagery to monitor the water quality in eutrophic plateau lakes; further research needs to be undertaken on whether OHS imagery could be applied to other types of water bodies or other application fields such as vegetation remote sensing.

Conclusions
A new and customized Chl-a concentration retrieval algorithm based on measured spectral data and OHS imagery is proposed in this study. The model improved Chl-a estimation accuracy using the OHS imagery in extremely eutrophic plateau water bodies compared to the existing models proposed for clear open sea and turbid coastal waters. This was then used to document the temporal and spatial patterns of Chl-a distributions in Dianchi Lake, leading to several key findings.
First, the three-band model based on OHS imagery is most suitable for estimating the concentration of Chl-a in Dianchi Lake, with RMSE and MAPE of 15.55 µg/L and 16.31%, respectively.
Second, the patterns of Chl-a retrieval from OHS imagery exhibit a significant spatial heterogeneity in Dianchi Lake, which is significantly higher in the surrounding region than the central waters.
Third, the FLAASH model is appropriate for estimating the concentration of Chl-a in Dianchi Lake, while the 6S model, QUAC model, and dark pixel method are not suitable for eutrophic Dianchi waters.
Moreover, the SNR of the OHS image (mean value of 19.43) is similar to that of the Hyperion image (mean value of 27.87) but nearly 2.5 times SNR of the HSI image (mean value of 8.47), and the variations in spatial resolution do not show a significant effect on the SNR of the OHS image. The findings indicated that the performance of SNR is much improved from HSI to OHS imagery; however, further improvement is still needed since there is still a gap compared with advanced hyperspectral imageries, such as Hyperion, etc.
Lastly, OHS imagery can be applied to remote sensing retrieval of Chl-a concentration in eutrophic plateau lakes.

Data Availability Statement:
The data presented in this study are available upon request from the corresponding author.
Remote Sensing Research Group of the Key Laboratory of Virtual Geography Environment of the Ministry of Education of Nanjing Normal University. We sincerely thank Yunmei Li for her help and the hard work of the experimental staff. We would like to appreciate Zhuhai Orbita Aerospace Technology Co., Ltd. and Dianchi Administration Bureau for providing remote sensing data and in situ chlorophyll-a concentration data respectively. We are deeply grateful to the three anonymous reviewers and academic editors for their constructive comments and suggestions to improve the quality of this work.

Conflicts of Interest:
The authors declare no conflict of interest.
Appendix A Table A1. OHS image bands setting and radiometric calibration coefficient. Wavelength represents the central wavelength of each band, FWHM (Full-Width Half-Maximum) is the wave width of each band, Radiance Gains and Radiance Offsets are the radiometric calibration gain coefficient and offset coefficient respectively, TDIStage is the integration series. The parts in bold are the bands used in this study.

Bands
Wavelength ( . Figure A1. Comparison of existing Chl-a concentration retrieval algorithms based on measured data. (a-n) are the models reconstructed by applying the algorithms in Table A2 to Dianchi Lake through calibration parameters. Figure A1. Comparison of existing Chl-a concentration retrieval algorithms based on measured data. (a-n) are the models reconstructed by applying the algorithms in Table A2 to Dianchi Lake through calibration parameters. Figure A2. Comparison of the validation results of existing Chl-a concentration retrieval algorithms based on measured data. (a-n) are the validation of the reconstruction models in Figure A1 in turn. Figure A2. Comparison of the validation results of existing Chl-a concentration retrieval algorithms based on measured data. (a-n) are the validation of the reconstruction models in Figure A1 in turn. Figure A3. Comparison of existing Chl-a concentration retrieval algorithms based on OHS imagery. (a-n) are the models reconstructed by applying the algorithms in Table A3 to Dianchi Lake through calibration parameters. Figure A3. Comparison of existing Chl-a concentration retrieval algorithms based on OHS imagery. (a-n) are the models reconstructed by applying the algorithms in Table A3 to Dianchi Lake through calibration parameters. Figure A4. Comparison of the results of existing Chl-a concentration retrieval algorithms based on OHS imagery. (a-n) are the validation of the reconstruction models in Figure A3 in turn. Figure A4. Comparison of the results of existing Chl-a concentration retrieval algorithms based on OHS imagery. (a-n) are the validation of the reconstruction models in Figure A3 in turn.