An Improved Model for Chlorophyll-a Concentration Retrieval in Coastal Waters Based on UAV-Borne Hyperspectral Imagery: A Case Study in Qingdao, China

: Chlorophyll-a (Chl-a) is an objective biological indicator, which reﬂects the nutritional status of coastal waters. However, the turbid coastal waters pose challenges to the application of existing Chl-a remote sensing models of case II waters. Based on the bio-optical models, we analyzed the suppression of coastal total suspended matter (TSM) on the Chl-a optical characteristics and developed an improved model using the imagery from a hyper-spectrometer mounted on an unmanned aerial vehicle (UAV). The new model was applied to estimate the spatiotemporal distribution of Chl-a concentration in coastal waters of Qingdao on 17 December 2018, 22 March 2019, and 20 July 2019. Compared with the previous models, the correlation coe ﬃ cients (R 2 ) of Chl-a concentrations retrieved by the new model and in situ measurements were greatly improved, proving that the new model shows a better performance in retrieving coastal Chl-a concentration. On this basis, the spatiotemporal variations of Chl-a in Qingdao coastal waters were analyzed, showing that the spatial variation is mainly related to the TSM concentration, wind waves, and aquaculture, and the temporal variation is mainly inﬂuenced by the sea surface temperature (SST), sea surface salinity (SSS), and human activities.


Introduction
The ocean accounts for 71% of the Earth's total area and has 97% of the world's water resources. Ocean waters, especially coastal waters, are important ecosystems and provide a lot of essential ecosystem services for people. However, with the rapid development of the coastal economy, the coastal water quality has gradually deteriorated in the past decades. Especially after 1950, the development of aquaculture has led to significant degradation of the ecological environment, which has threatened biological resources and human health [1][2][3]. Since 2007, green tides have occurred in the Yellow Sea every year due to the water eutrophication caused by ocean pollution. A large amount of floating In this paper, we demonstrated the process of using a UAV equipped with a hyperspectral imager to image coastal waters in Qingdao and successfully estimated the Chl-a concentration. The objective was three-fold: (1) To demonstrate the ability of UAVs equipped with hyperspectral imagers for imaging coastal waters and monitoring water quality; (2) to develop a new model to quantitatively and accurately estimate Chl-a concentration in coastal waters from UAV-borne hyperspectral radiation measurements; and (3) to analyze the spatiotemporal dynamics of Chl-a concentration according to the retrieval results on the sampling dates.

Study Area
The coastal waters of Qingdao are located in the southwest of the Yellow Sea; this area has a warm temperate monsoon climate. The coast is mostly composed of mountain bedrocks, and the main sediment types are silt and sandy silt. The coastal waters are mainly clean sea areas with good hydrodynamic conditions and strong self-purification ability. The deepest location is at the mouth of Jiaozhou Bay, with a water depth of about 64 m [45].
The research area is located in the coastal waters of Huangdao District in Qingdao. There are two specific areas, as shown in Figure

Image Data and Preprocessing
Images were acquired from a push-broom airborne hyperspectral imager named Pika L (Resonon, Inc., Bozeman, MT, USA), which covers a spectral range of 400-1000 nm with a spectral resolution of 2.1 nm. Pika L was mounted on an M600 PRO UAV (DJI, Co., Ltd., Shenzhen, Guangdong, China), and an Ellipse2-D inertial sensor (SBG Systems, S.A.S., Carrières-sur-Seine, France) was integrated to acquire orientation and position in real time during the flight, as shown in Figure 2. Before image acquisition, the flight path was planned, and flight speed and height were set Images were acquired in coastal waters of northern Golden Beach (35 • 57 N, 120 • 15 E) in Huangdao District on 17 December 2018, i.e., research area A. Golden Beach is located in the south of Shandong peninsula, facing the Yellow Sea in the south and extending towards the southeast and is crescent-shaped. It has a total length of about 3500 m and a width of about 300 m. The water here is clear, and the sand is fine. It is a barrier-free coast formed by weathering and transformation of bedrock. There are two high and low tides every day. The annual average high tide is 3.5 m, and the average low tide is 0.9 m [46]. The primary sediments in coastal waters of northern Golden Beach are medium sands and coarse sands, so it has high visibility. There is a lot of marine life here, as well as marine farms, where one can see many fishing boats and fishers working. The long dike extending into the sea is convenient for ships to dock. There is also a broad platform on the shore for taking off and landing of UAVs equipped with hyper-spectrometers.
Images were acquired in coastal waters of northern Zhangtunzui (35 •  Tidal flats are mainly formed by shellfish deposits, and weathered and eroded bedrock gravel are deposited on lagoon beaches. Compared with the coast of Golden Beach, the hydrodynamic conditions are generally better, but there are many reefs exposed. South of the study area is close to Zhangtunzui that extends along Fenghuang mountains to the ocean, and hydrodynamic force where reef exposure is weakened, so that Enteromorpha gather easily in summer. There is also a broad platform on the shore, which is convenient for unmanned aerial vehicles (UAVs) equipped with airborne hyper-spectrometers to take off and land.

Image Data and Preprocessing
Images were acquired from a push-broom airborne hyperspectral imager named Pika L (Resonon, Inc., Bozeman, MT, USA), which covers a spectral range of 400-1000 nm with a spectral resolution of 2.1 nm. Pika L was mounted on an M600 PRO UAV (DJI, Co., Ltd., Shenzhen, Guangdong, China), and an Ellipse2-D inertial sensor (SBG Systems, S.A.S., Carrières-sur-Seine, France) was integrated to acquire orientation and position in real time during the flight, as shown in Figure 2. Before image acquisition, the flight path was planned, and flight speed and height were set using the flight control software. The flight height, speed, and sampling frame rate of hyper-spectrometer jointly determine the ground spatial resolution of the pixels on the image.

Image Data and Preprocessing
Images were acquired from a push-broom airborne hyperspectral imager named Pika L (Resonon, Inc., Bozeman, MT, USA), which covers a spectral range of 400-1000 nm with a spectral resolution of 2.1 nm. Pika L was mounted on an M600 PRO UAV (DJI, Co., Ltd., Shenzhen, Guangdong, China), and an Ellipse2-D inertial sensor (SBG Systems, S.A.S., Carrières-sur-Seine, France) was integrated to acquire orientation and position in real time during the flight, as shown in Figure 2. Before image acquisition, the flight path was planned, and flight speed and height were set using the flight control software. The flight height, speed, and sampling frame rate of hyperspectrometer jointly determine the ground spatial resolution of the pixels on the image. In Qingdao, hyperspectral images of coastal waters were acquired three times. The flight days are all clear and cloudless, with plenty of light and gentle wind on the shore and offshore. The flight In Qingdao, hyperspectral images of coastal waters were acquired three times. The flight days are all clear and cloudless, with plenty of light and gentle wind on the shore and offshore. The flight parameters are shown in Table 1. Each of the flight tasks needed to be divided into multiple flights to cover the whole study area because of battery limitations.
Image preprocessing includes radiation calibration, atmospheric correction, image filtering, geometric correction, and flight belts mosaicking. The hyperspectral images were radiometrically calibrated to covert the digital number to radiance. Atmospheric correction was implemented on the images after radiation calibration to eliminate the influence of the atmosphere and light using the Fast Line-of-Sight Atmospheric Analysis of Spectral Hypercubes (FLAASH) module via ENVI 5.4 software (Exelis Visual Information Solutions, Inc., Boulder, CO, USA). Savitzky-Golay filtering based on local polynomial least square fitting in the time domain was adopted in this study because it can not only smooth jagged noise but also keep the basic spectral characteristics of the curves unaffected [47,48]. The flight belts were then geometrically corrected based on the principle of homography [49] and configured with geographic information according to the coordinates of imaging centers provided by the Global Positioning System (GPS) receiver connected with the Ellipse2-D sensor. All geometrically corrected flight belts were mosaicked to an image that covers the whole study area. The true color synthesis images of sea surface after preprocessing are displayed in Figure 3.

Shipborne Spectrum Sampling and Spectral Correction
To improve the quality of the sea surface spectra collected by Pika L, a shipborne marine hyperspectrometer named QE Pro (Ocean Optics, Inc., Dunedin, FL, USA) was used to automatically collect the spectra every two minutes while collecting water samples in the study area on 22 March

Water Sampling and Measurement
Seawater sampling and image acquisition were carried out at the same time. A total of 20, 24, and 20 water samples were collected on 17 December 2018, 22 March 2019, and 20 July 2019, respectively. The location of the sampling sites is presented in Figure 3. Because the research area is close to shore and water is shallow, sediment was easily stirred up when the ship passed. Therefore, we needed to wait for a few minutes after reaching the planned sampling sites. Surface water was collected using a 5000-mL water sampler. Water samples were injected from the water sampler into bottles and quickly frozen and stored for measurements of concentration in the laboratory.
The TSM concentration and absorption coefficient of CDOM at 440 nm were also necessary in this study. Therefore, three portions of each sample were used for collecting Chl-a, TSM, and CDOM. For Chl-a, Whatman GF/F glass microfiber filters (0.7 µm pore size) were used to filter water samples with a capacity of 2000 mL, and acetone with a volume fraction of 90% was added to extract the pigments. Chl-a concentration was determined by spectrophotometry using a Cary 5000 UV-Vis-NIR spectrophotometer (Agilent, Inc., Santa Clara, CA, USA). For TSM, water samples with a capacity of 1000 mL were filtered through Whatman cellulose acetate filters (0.45 µm pore size) and the filters were weighted gravimetrically to determine the TSM concentration. For CDOM, after the samples were filtered through Isopore membrane filters (0.2 µm pore size), the absorption coefficient of CDOM was measured in a 10-cm optical cuvette against distilled water with a Cary 5000 UV-Vis-NIR spectrophotometer. During the laboratory measurements, the test equipment was thoroughly cleaned and dried before each new sample was measured to ensure that no external pollution was introduced.

Shipborne Spectrum Sampling and Spectral Correction
To improve the quality of the sea surface spectra collected by Pika L, a shipborne marine hyperspectrometer named QE Pro (Ocean Optics, Inc., Dunedin, FL, USA) was used to automatically collect the spectra every two minutes while collecting water samples in the study area on 22 March 2019. QE Pro has a sensitivity range of 350-975 nm at an increment of 0.7 nm. It was fixed on the shipboard, measuring the radiance of skylight (L sky ), upwelling radiance from water (L sw ), and reference panel (L p ) according to the above-water method. The observation geometry was automatically determined by the direction of sun incidence and sea surface normal. A total of 34 spectra were obtained, as highlighted in Figure 4. In situ measurement processing included two parts: spectral preprocessing and remote sensing reflectance calculation. First, some preprocessing work (e.g., removing abnormal spectra affected by solar flares and shadows, correcting dark noise, and averaging spectra collected multiple times) was conducted to get more realistic spectra. Then, remote sensing reflectance (R rs ) was calculated according to Equation (1): where ]ρ p is the irradiance reflectance of the reference panel, and δ is a proportionality coefficient that relates L sky to the reflected sky radiance determined when the detector viewed the water surface. Images were acquired approximately 150-200 m above the water surface using Pika L. The impact of the atmosphere between the sensor and sea surface cannot be ignored. Even if the hyperspectral images are atmospheric corrected, it still cannot represent the real sea surface reflectance. As mentioned before, a total of 34 spectra were collected by QE Pro on 22 March 2019. One of them was removed because it was located in the specular reflection area of the image, so there were 33 valid sites remaining in the dataset. Figure 5a,b display the spectral reflectance curves of 33 sites collected by Pika L and QE Pro, respectively. It can be seen that the trends of two spectral curves at the same site are similar, but the reflectance and spectral characteristics between them are different. QE Pro was closer to the sea surface than Pika L, and the observation geometry was automatically optimized. Accordingly, spectra from QE Pro can be considered more real. The spectra acquired by Pika L need to be corrected using those collected by QE Pro. Based on the coordinates of spectrum sampling sites, the differences of the reflectance from hyperspectral images and QE Pro were calculated. For each band, the least square fitting was performed on the average differences to obtain the spectral correction parameters, which were then used to correct the reflectance of images. It can not only correct the low radiance of airborne spectra due to the long observation distance but also maintain the small differences of reflectance among different sampling sites as much as possible. Figure 5c,d show the results of the least square fitting and the spectral curves of the sampling sites after correction. Considering that the differences between airborne and shipborne spectra at the same point were mainly due to factors, including the observation distance and atmospheric environment, and the distance and weather condition for each observation were basically similar, the correction parameters calculated from the where is the irradiance reflectance of the reference panel, and δ is a proportionality coefficient that relates to the reflected sky radiance determined when the detector viewed the water surface. Images were acquired approximately 150-200 m above the water surface using Pika L. The impact of the atmosphere between the sensor and sea surface cannot be ignored. Even if the hyperspectral images are atmospheric corrected, it still cannot represent the real sea surface reflectance. As mentioned before, a total of 34 spectra were collected by QE Pro on 22 March 2019. One of them was removed because it was located in the specular reflection area of the image, so there were 33 valid sites remaining in the dataset. Figure 5a,b display the spectral reflectance curves of 33 sites collected by Pika L and QE Pro, respectively. It can be seen that the trends of two spectral curves at the same site are similar, but the reflectance and spectral characteristics between them are different. QE Pro was closer to the sea surface than Pika L, and the observation geometry was automatically optimized. Accordingly, spectra from QE Pro can be considered more real. The spectra acquired by Pika L need to be corrected using those collected by QE Pro. Based on the coordinates of spectrum sampling sites, the differences of the reflectance from hyperspectral images and QE Pro were calculated. For each band, the least square fitting was performed on the average differences to obtain the spectral correction parameters, which were then used to correct the reflectance of images. It can not only correct the low radiance of airborne spectra due to the long observation distance but also maintain the small differences of reflectance among different sampling sites as much as possible. Figure 5c,d show the results of the least square fitting and the spectral curves of the sampling sites after correction. Considering that the differences between airborne and shipborne spectra at the same point were mainly due to factors, including the observation distance and atmospheric environment,

Model Assessment
In this paper, a semi-analytical and semi-empirical algorithm was used to establish the retrieval model of Chl-a concentration. Accuracy evaluation indicators for polynomial regression were used to evaluate the algorithm, including the correlation coefficient (R 2 ), root mean square error (RMSE), mean absolute error (MAE), and average percentage difference (APD). They were calculated by Equations (2)-(5), respectively:

Model Assessment
In this paper, a semi-analytical and semi-empirical algorithm was used to establish the retrieval model of Chl-a concentration. Accuracy evaluation indicators for polynomial regression were used to evaluate the algorithm, including the correlation coefficient ( ), root mean square error (RMSE), mean absolute error (MAE), and average percentage difference (APD). They were calculated by Equations (2), (3), (4), and (5), respectively:

Analysis of Water Samples and Spectra
Image acquisition was carried out on cloudless and sunny days. Sea surface fluctuations easily cause direct reflection of sun rays, also known as specular reflection or solar flare, appearing as a piece or scattered bright area on hyperspectral images. The water spectra of sampling sites located within the specular reflection area of images cannot represent the real reflectance and need to be eliminated.
In conclusion, a valid sampling site that can be used to establish and verify the retrieval model must meet the following requirements: 1.
It is located within the study area covered by UAV-borne hyperspectral images after geometric correction;  The spectral characteristics of seawater are jointly affected by pure water, pigment particles, non-pigment particles, and CDOM [50]. In this paper, spectra were collected from the coastal waters in Qingdao with shallow water and a sandy beach. The spectral characteristics of seawater are relatively complex because of not only the water condition but also some external factors (e.g., observation geometry and illumination). According to the coordinates of the water sampling sites, the spectral curves of valid sites for three samplings were obtained from the spectral-corrected images, as shown in Figure 6.
It can be seen that the trends of the spectral curves at the wavelength of 400-800 nm are similar in different dates. The spectral characteristic is consistent with that of case II waters, but it has its uniqueness to some extent. Due to the absorption of CDOM and Chl-a, the reflectance is relatively low at the wavelength of 400-500 nm. Particularly, a valley caused by the absorption of blue-violet light by Chl-a is located at the wavelength of 420-440 nm. The absorption valley is more evident on 20 July, which is also associated with the increase of algae at this time. There is a significant reflection peak near 550 nm due to the scattering of TSM, and an unclear absorption valley of Chl-a at around 670 nm. The reflection peak appears again at the wavelength of 680-720 nm, which is generally regarded as the fluorescence peak of Chl-a. Its position is shifted to the long-wave direction as the concentration of Chl-a increases [51]. However, the high TSM concentration of coastal waters weakens the response of the fluorescence peak to the Chl-a concentration and the fluorescence peak is not significant. The peak at 750-800 nm is another reflection peak of TSM. When the concentration of TSM increases, the reflection peak will also move to the long-wave direction.
Reflectance between 400 and 800 nm is normalized to reduce the impact of external factors on the spectrum collection by the UAV-borne hyper-spectrometer and determine the sensitive bands for retrieval of the Chl-a concentration. The Pearson correlations between the normalized reflectance and Chl-a concentration or TSM concentration are demonstrated in Figure 7. It can be seen that the trends of the spectral curves at the wavelength of 400-800 nm are similar in different dates. The spectral characteristic is consistent with that of case II waters, but it has its uniqueness to some extent. Due to the absorption of CDOM and Chl-a, the reflectance is relatively low at the wavelength of 400-500 nm. Particularly, a valley caused by the absorption of blue-violet light by Chl-a is located at the wavelength of 420-440 nm. The absorption valley is more evident on 20 July, which is also associated with the increase of algae at this time. There is a significant reflection peak near 550 nm due to the scattering of TSM, and an unclear absorption valley of Chl-a at around 670 nm. The reflection peak appears again at the wavelength of 680-720 nm, which is generally regarded as the fluorescence peak of Chl-a. Its position is shifted to the long-wave direction as the concentration of Chl-a increases [51]. However, the high TSM concentration of coastal waters weakens the response of the fluorescence peak to the Chl-a concentration and the fluorescence peak is not significant. The peak at 750-800 nm is another reflection peak of TSM. When the concentration of TSM increases, the reflection peak will also move to the long-wave direction.
Reflectance between 400 and 800 nm is normalized to reduce the impact of external factors on the spectrum collection by the UAV-borne hyper-spectrometer and determine the sensitive bands for retrieval of the Chl-a concentration. The Pearson correlations between the normalized reflectance and Chl-a concentration or TSM concentration are demonstrated in Figure 7. From Figure 7a, it can be concluded that there is a weak correlation between the normalized reflectance and Chl-a concentration at some bands, the most noticeable of which are the maximum positive correlation at 520-560 nm caused by the strong reflection of Chl-a at the green band and the minimum correlation at approximately 420 nm caused by the strong absorption of Chl-a at the blueviolet band. However, the correlations between the normalized reflectance and TSM concentration in different dates are not only consistent but also strong around 560 nm, as shown in Figure 7b. Overall, the correlation between Chl-a and the normalized reflectance is weaker than that between TSM and the normalized reflectance. The spectral information of seawater is dominated by a high level of TSM concentration in the coastal area of Qingdao.

Influence of TSM on Chl-a Retrieval
In order to prove that a high TSM concentration weakens the contribution of Chl-a on the spectra in the retrieval of coastal Chl-a concentration, the bio-optical models were used to simulate the reflectance of coastal waters. The influence of pure water, CDOM, and the external environment on the reflectance was excluded, and only the interference of TSM on the retrieval of Chl-a was discussed. From Figure 7a, it can be concluded that there is a weak correlation between the normalized reflectance and Chl-a concentration at some bands, the most noticeable of which are the maximum positive correlation at 520-560 nm caused by the strong reflection of Chl-a at the green band and the minimum correlation at approximately 420 nm caused by the strong absorption of Chl-a at the blue-violet band. However, the correlations between the normalized reflectance and TSM concentration in different dates are not only consistent but also strong around 560 nm, as shown in Figure 7b. Overall, the correlation between Chl-a and the normalized reflectance is weaker than that between TSM and the normalized reflectance. The spectral information of seawater is dominated by a high level of TSM concentration in the coastal area of Qingdao.

Influence of TSM on Chl-a Retrieval
In order to prove that a high TSM concentration weakens the contribution of Chl-a on the spectra in the retrieval of coastal Chl-a concentration, the bio-optical models were used to simulate the reflectance of coastal waters. The influence of pure water, CDOM, and the external environment on the reflectance was excluded, and only the interference of TSM on the retrieval of Chl-a was discussed.
The total absorption coefficient a(λ) is expressed as a sum of the absorption coefficients for pure water a w (λ), Chl-a a φ (λ), and CDOM a g (λ), whereas the total backscattering coefficient b b (λ) is expressed as a sum of the backscattering coefficients for pure water b bw (λ), Chl-a b bp (λ), and TSM b bs (λ), as Equation (6) describes: The absorption and backscattering coefficients for pure water come from the literature [52,53]. The coefficients for other components can be simulated using the empirical models as follows [54][55][56][57]: where [Chla] and [TSM] are the concentrations of Chl-a and TSM, a g (440) is the absorption coefficient of CDOM at 440 nm, a 0 (λ) and a 1 (λ) are empirical values and taken from Lee, and λ is the wavelength. Above-surface remote-sensing reflectance R rs (λ) can be calculated as follows [58,59]: where u(λ) is the ratio of the backscattering coefficient to the sum of the absorption and backscattering coefficients and r rs (λ) is the subsurface remote-sensing reflectance.
The in situ measurements were substituted into the above models to simulate water reflectance, as shown in Figure 8a. The magnitude and trend of reflectance from the simulation are not completely consistent with those from the UAV-borne images because of the empirical parameters and the exclusion of the observation environment. However, the influence of changes in various parameters on remote sensing reflectance can be studied by simulation.
According to the measurements on 17 December 2018, when the TSM concentrations were set to be the reference concentration (2 mg/L), the measured minimum (~8 mg/L) and the measured maximum (~12 mg/L), the change rates of reflectance caused by the increase of Chl-a from the measured minimum (0.8409 mg/m 3 ) to the measured maximum (1.1372 mg/m 3 ) were simulated, as shown in Figure 8b. Similarly, Figure 8c,d show the simulations based on the measurements on 22 March 2019, and 20 July 2019, respectively. From the figures, we can see that the increase of the Chl-a concentration disturbs the water reflectance in varying degrees at varying bands. A positive change rate indicates that the water reflectance increases while the concentration of Chl-a increases, and the scattering effect of Chl-a is greater than the absorption effect. On the contrary, a negative value indicates that the absorption effect is greater than the scattering effect.
It can be seen from Figure 8b-d that when the TSM concentration is 2 mg/L and the Chl-a concentration rises (red lines), the water reflectance shows a noticeable negative rate of change between 430 and 520 nm. At this time, the absorption of Chl-a is significant. Meanwhile, the water reflectance shows a positive change after 540 nm and the scattering of Chl-a is dominant now. In particular, the change rate of reflectance is positive at 550 nm, because the reflection of Chl-a for green light is enhanced as the Chl-a concentration increases. The optical properties of water are mainly affected by Chl-a when the TSM concentration is maintained at a low level. When the TSM concentration increases to the actual level (blue and green lines), the variations of reflectance caused by the change of the Chl-a concentration become larger at the blue-green bands and smaller at the red bands, which proves that the blue and green bands are more sensitive to Chl-a concentration changes in this area. In addition, the change rate of reflectance becomes negative at 550 nm, because the contribution of Chl-a is no longer dominant due to a high level of TSM concentration. At this time, the scattering caused by the increase of Chl-a is much weaker than the TSM scattering, i.e., u(λ) decreases, and more consideration should be given to the TSM in the retrieval. backscattering coefficients and is the subsurface remote-sensing reflectance. The in situ measurements were substituted into the above models to simulate water reflectance, as shown in Figure 8a. The magnitude and trend of reflectance from the simulation are not completely consistent with those from the UAV-borne images because of the empirical parameters and the exclusion of the observation environment. However, the influence of changes in various parameters on remote sensing reflectance can be studied by simulation.
According to the measurements on 17 December 2018, when the TSM concentrations were set to be the reference concentration (2 mg/L), the measured minimum (~8 mg/L) and the measured maximum (~12 mg/L), the change rates of reflectance caused by the increase of Chl-a from the measured minimum (0.8409 mg/m 3 ) to the measured maximum (1.1372 mg/m 3 ) were simulated, as shown in Figure 8b. Similarly, Figure 8c,d show the simulations based on the measurements on 22 March 2019, and 20 July 2019, respectively. From the figures, we can see that the increase of the Chla concentration disturbs the water reflectance in varying degrees at varying bands. A positive change rate indicates that the water reflectance increases while the concentration of Chl-a increases, and the scattering effect of Chl-a is greater than the absorption effect. On the contrary, a negative value indicates that the absorption effect is greater than the scattering effect.

Retrieval Model and Comparison
As we all know in the Tassan model [42], Chl-a concentration is retrieved by: where C is the retrieved concentration of Chl-a; c 0 , c 1 , and c 2 are fitting coefficients, x = R rs(443) /R rs(555) R rs(412) /R rs (490) a ; a is a constant related to the study area; and R rs(i) is the remote sensing reflectance at band i. The selection of bands depends on the sensitive bands of different components. The first term of x is used to express the Chl-a concentration by the ratio of the absorption and reflection peak, and the second one is used to eliminate the influence of CDOM and TSM on Chl-a concentration retrieval of case II waters.
Considering that Chl-a still contributes to the spectra at the wavelength of 490 nm in the Yellow and East China Seas, the impact of TSM was removed using 510 nm instead of 490 nm, i.e., x = R rs(443) /R rs(555) R rs(412) /R rs(510) a . Therefore, the NSOAS model was proposed, which was more suitable for concentration retrieval in the Yellow and East China Seas [43].
Compared with the case II waters with a certain distance from the coast, the water is shallower and the TSM concentration is higher in coastal waters of Qingdao. The previous analysis proved that the optical properties of Chl-a are suppressed by the TSM, which causes a decline in the predictive accuracy of the Tassan and NSOAS models. Therefore, in order to improve the retrieval accuracy and further reduce the impact of TSM, two improvements were made on the NSOAS model: multiple logarithmic power of the second ratio is added to the right side of Equation (9), and lg 2 x is decomposed to better reflect the concentration of Chl-a. After verification, a new retrieval model is obtained as follows: where c 0 , c 1 , . . . , c 7 are fitting coefficients, x 1 = R rs(443) /R rs (555) , and x 2 = R rs(412) /R rs(510) . Subsequently, the accuracy evaluation was performed to compare the Tassan, NSOAS, and new models, as illustrated in Figure 9. Figure 9a shows the comparison of the three models on 17 December 2018, and R 2 increased from less than 0.20 to 0.65; Figure 9b shows the comparison on 22 March 2019, and R 2 increased from 0.21 to 0.66; and Figure 9c shows the comparison on 20 July 2019, and R 2 increased from about 0.70 to 0.79. Table 3 summarizes the accuracy analysis of different models. It can be concluded that for the retrieval of the Chl-a concentration in the study area on 17 December 2018, 22 March 2019, and 20 July 2019, R 2 of the new model was higher, and RMSE, MAE, and APD of the new model were lower than those of Tassan and NSOAS models. It is proved that the new model has higher accuracy and more advantages when applied to retrieve Chl-a concentration in coastal waters of Qingdao.

Spatial Variation
The new model was applied to estimate Chl-a concentration based on the UAV-borne hyperspectral images, and the spatial distributions of Chl-a concentration on the sampling dates were obtained, as shown in Figure 10. Chl-a concentration was displayed in chromaticity. In order to demonstrate the spatial variation of Chl-a concentration in more detail, the land and retrieved outliers were excluded from the images, as well as the highest and lowest concentrations, which accounted for 5% of the total number of pixels in the study area. All excluded pixels were colored gray or white, and the remaining pixels were colored and displayed according to the concentration.    Figure 10a shows the estimated spatial distribution of Chl-a concentration in study area A on 17 December 2018. In this area, the range of Chl-a concentration was approximately 0.82-1.41 mg/m 3 on this day, and the average concentration was 1.09 mg/m 3 . The average concentration of Chl-a is the highest of these three days. The average concentration at latitudes > 35.97 • (1.13 mg/m 3 ) is slightly higher than that at latitudes < 35.97 • (1.07 mg/m 3 ), which is mainly caused by human activities. Figure 10b shows the estimated spatial distribution of Chl-a concentration in study area B on 22 March 2019. Affected by the marine climate, the seawater temperature reaches the lowest level of these three days, as well as the Chl-a concentration, i.e., 0.05-0.32 mg/m 3 , with an average concentration of 0.30 mg/m 3 . The lowest concentrations of Chl-a are distributed in the south of the study area, where it is close to the coast of Zhangtunzui. Figure 10c shows the estimated spatial distribution of Chl-a concentration in study area B on 20 July 2019. On this day, the spatial changes of Chl-a are significant, ranging from 0.13 to 1.58 mg/m 3 . The average Chl-a concentration is about 0.72 mg/m 3 . The higher concentrations are distributed in the areas where algae aggregation occurs. obtained, as shown in Figure 10. Chl-a concentration was displayed in chromaticity. In order to demonstrate the spatial variation of Chl-a concentration in more detail, the land and retrieved outliers were excluded from the images, as well as the highest and lowest concentrations, which accounted for 5% of the total number of pixels in the study area. All excluded pixels were colored gray or white, and the remaining pixels were colored and displayed according to the concentration.   The distribution of Chl-a showed clear spatial variations in coastal waters of Qingdao. For study area A, a large number of fish and shrimp farms are distributed on the coast, resulting in a higher Chl-a concentration near the coast than in the distant sea on 17 December 2018. Aquaculture farms provide organic matter and nutrients to the surrounding waters, promoting the development of phytoplankton biomass, which shows the impact of human activities on the concentration of Chl-a [19,60]. For study area B, the concentration of Chl-a in the south is lower than that in the north on 22 March 2019, because the southern area is close to the coast of Zhangtunzui, and seabed sediments are deposited. A high level of TSM concentration becomes a key factor that restricts the growth of phytoplankton [24,61,62]. On 20 July 2019, the spatial changes of Chl-a concentration in this area are significant, which is related to the increase of seawater temperature and the accumulation of phytoplankton caused by wind waves [63].

Temporal Variation
Qingdao is located in the north extratropical monsoon climatic region, and its climatic conditions are favorable for many kinds of aquatics. Due to the influence of the southeast monsoon, ocean currents, and water masses, Qingdao has significant features of a maritime climate. It is humid, hot, and rainy in summer; the low temperature lasts for a long time in winter; and the temperature rises slowly in spring, about one month later than inland. Sea surface temperature (SST) has obvious seasonal variation characteristics, with high temperatures in summer and low temperatures in winter and spring [8]. Sea surface salinity (SSS) is mainly affected by the rainfall, and it gradually decreases from the maximum in spring with an increase of the rainfall [64,65]. Many studies have proven that the concentration of Chl-a can be influenced by various factors, such as SST, SSS, and human activities [65][66][67][68].
The Chl-a concentration in coastal waters of Qingdao has changed greatly with time. Figure 11 shows the temporal variations of the average Chl-a concentration, SST, and SSS in the study area. SST and SSS data are from the in situ measurements of the nearby Xiaomaidao ocean observation station released by the National Marine Science Data Center. As we can see from the figure, when SST decreased from 17 December to 22 March, the average Chl-a concentration also decreased. However, when SST increased drastically from 22 March to 20 July, the average Chl-a concentration also increased slightly, but it was still lower than that on December 17. It proves that SST is a critical factor that affects the concentration of Chl-a. Within a certain temperature range, the positive correlation between Chl-a concentration and SST is significant, i.e., decreasing SST will inhibit the growth of phytoplankton, and Chl-a concentration will also decrease; on the contrary, as SST increases, phytoplankton grows vigorously, and the concentration of Chl-a also increases. However, excessively high temperatures will inhibit the rise of the Chl-a concentration, which is the reason why the concentration of Chl-a on 20 July is lower than that on December 17; this finding is consistent with previous research results [69][70][71][72][73]. The change of SSS is almost the opposite of that of the Chl-a concentration. From 17 December to 22 March, SSS increased and Chl-a concentration decreased, but from 22 March to 20 July, SSS decreased and Chl-a concentration increased, showing that high salinity is a limitation for phytoplankton growth [61,65].

Conclusions
Considering the low predictive accuracy of the existing Chl-a concentration retrieval algorithms for case II waters, an improved retrieval algorithm for coastal waters in Qingdao was established based on UAV-borne hyperspectral images. The particularity of coastal water constituents was considered in the new algorithm, i.e., a high level of TSM concentration will suppress the optical properties of Chl-a and result in a low correlation between Chl-a concentration and spectral reflectance. Similar to the principle of the existing algorithms, the new algorithm was improved to further reduce the impact of TSM on the retrieval of Chl-a concentration. Compared with the previous algorithms, the correlation between retrieved and in situ Chl-a concentrations was higher, and the errors were lower using the new algorithm, which proved that the new algorithm was more suitable for Chl-a concentration retrieval of coastal waters.
The Chl-a concentration in the research area was estimated using the new algorithm on the sampling dates. The spatial and temporal variation characteristics of the Chl-a concentration and the

Conclusions
Considering the low predictive accuracy of the existing Chl-a concentration retrieval algorithms for case II waters, an improved retrieval algorithm for coastal waters in Qingdao was established based on UAV-borne hyperspectral images. The particularity of coastal water constituents was considered in the new algorithm, i.e., a high level of TSM concentration will suppress the optical properties of Chl-a and result in a low correlation between Chl-a concentration and spectral reflectance. Similar to the principle of the existing algorithms, the new algorithm was improved to further reduce the impact of TSM on the retrieval of Chl-a concentration. Compared with the previous algorithms, the correlation between retrieved and in situ Chl-a concentrations was higher, and the errors were lower using the new algorithm, which proved that the new algorithm was more suitable for Chl-a concentration retrieval of coastal waters.
The Chl-a concentration in the research area was estimated using the new algorithm on the sampling dates. The spatial and temporal variation characteristics of the Chl-a concentration and the relationships with various factors (e.g., SST, SSS, wind waves, and human activities) were analyzed. It is concluded that the spatial variation of Chl-a concentration in coastal waters is mainly related to the growth and reproduction of algae influenced by TSM concentration, wind waves, and aquaculture. A high level of TSM concentration will suppress the increase of Chl-a concentration, while wind waves and aquaculture will promote the rise of the Chl-a concentration of the surrounding waters. The temporal variation of the Chl-a concentration in the study area is the result of combined effects of SST, SSS, and human activities, which was the highest on December 17, followed by 20 July, and lowest on 22 March. As SST rises, phytoplankton grows vigorously and Chl-a concentration increases, but excessively high SST inhibits the rise of Chl-a concentration. Increased SSS inhibits the growth of phytoplankton, resulting in a decrease in the Chl-a concentration.