Assessment of VIIRS on the Identiﬁcation of Harmful Algal Bloom Types in the Coasts of the East China Sea

: Visible Infrared Imaging Radiometer Suite (VIIRS) data were systematically evaluated and used to detect harmful algal bloom (HAB) and classify algal bloom types in coasts of the East China Sea covered by optically complex and sediment-rich waters. First, the accuracy and spectral characteristics of VIIRS retrieved normalized water-leaving radiance or the equivalent remote sensing reﬂectance from September 2019 to October 2020 that were validated by the long-term observation data acquired from an offshore platform and underway measurements from a cruise in the Changjiang Estuary and adjacent East China Sea. These data were evaluated by comparing them with data from the Moderate-Resolution Imaging Spectroradiometer. The bands of 486, 551, and 671 nm provided much higher quality than those of 410 and 443 nm and were more suitable for HAB detection. Secondly, the performance of four HAB detection algorithms were compared. The Ratio of Algal Bloom (RAB) algorithm is probably more suitable for HAB detection in the study area. Importantly, although RAB was also veriﬁed to be applicable for the detection of different kinds of HAB ( Prorocentrum donghaiense , diatoms, Ceratium furca, and Akashiwo sanguinea ), the capability of VIIRS in the classiﬁcation of those algal species was limited by the lack of the critical band near 531 nm.


Introduction
The Visible Infrared Imaging Radiometer Suite (VIIRS) is the new generation of medium resolution image radiometer of the United States, which combines most features of its predecessor, such as the Moderate-resolution Imaging Spectroradiometer (MODIS). It has nine visible/infrared bands with an average bandwidth of 20 nm, twelve mid-far infrared bands, and one day/night band, among which only five bands have a spatial resolution of 375 m and the spatial resolution of the remaining bands is 750 m [1][2][3][4]. VIIRS is mainly used to monitor the radiation changes in ocean, atmosphere, ice, and land in visible and infrared bands, and provide data for monitoring ocean color, ocean surface temperature, and other surface changes [4][5][6][7][8]. The first VIIRS onboard Suomi National Polar-orbiting Partnership (SNPP) was successfully launched into an 824-km sun synchronous polar orbit on 28 October 2011; subsequently, the first Joint Polar Satellite System (JPSS-1) satellite launched into space on 18 November 2017 with VIIRS. In December remains unclear. Furthermore, whether or not the absence of these bands limits the VIIRS application in discriminating bloom types also needs to be investigated. Therefore, in this study, we evaluated VIIRS data in bloom type detection and selected the ECS coasts as the study area. This assessment involves three main objectives. The first is to validate R rs data products from VIIRS (both MSL12 and L2GEN processing) when compared to Moderate Resolution Imaging Spectroradiometer of Aqua (MODISA; L2GEN processing only) and in situ datasets, since the reliability of operational monitoring of algal blooms primarily depends on the degree of the quantification of R rs products. Second is to compare different bloom detection methods using VIIRS data in coasts of the ECS. The last but the most important is to evaluate the capabilities of VIIRS data in the discrimination of various bloom types, such as P. donghaiense and diatoms, along with other blooms.

Study Area and In Situ Measurements
The study area covers most areas of the ECS (26)(27)(28)(29)(30)(31)(32)(33) • N, 119-126 • E), including the Yangtze River Estuary (YRE) and the coasts of Shanghai Municipality and Zhejiang Province ( Figure 1). Most of the coastal ocean has a depth of less than 50 m and is influenced by the Yangtze River plume, the Taiwan Strait Warm Current, the Kuroshio Current, and the coastal current along the Zhejiang coast [29,30]. Some of algal bloom records for the time of occurrence, location, and causative species were obtained from the Wenzhou Marine Environmental Monitoring Center Station (WMEMC) of State Oceanic Administration. These records cover the period from April to May 2020. Other records come from field measurements carried out during the cruise of Project of Long-term Observation and Research Plan in the Changjiang Estuary and Adjacent East China Sea (LORCE) in July 2013, August 2020, and August 2021. At each bloom station, 1 L or 500 mL water samples collected from Niskin bottles on the surface (2-3 m depth) were fixed with formalin to a final concentration of 2%. After at least 24 h of sedimentation, preserved samples were concentrated to 10-50 mL by slowly siphoning off the supernatant [31,32]. Subsequently, the concentrated samples were transferred into a 0.1 mL scaled slide using a pipette. The phytoplankton taxa were identified and counted on the entire slide using light microscopes (Leica DM2500 and DM6B) at 200×, 400×, or 630× magnification. At least 300 units (individual cells or colonies) were counted for each sample. All the collection, sedimentation, species identification, and abundance calculation of phytoplankton samples were performed complying with the Chinese National Standard [33].
Moreover, a fixed oceanographic platform called Dongou was used to collect a longterm series of radiometric measurements (such as R rs and normalized water-leaving radiance (L wn )). This platform was located in the southern part of the ECS (27.675 • N and 120.358 • E) at approximately 80 km southeast of Wenzhou in an area that minimized the adjacency effects in satellite data. In addition, the water depth around the platform was around 50 m, so that the bottom effects in L wn were negligible. This oceanographic platform, owned and operated by WMEMC, has served as an offshore laboratory since 2017 and is an excellent site supporting ocean color validation activities through a comprehensive collection of bio-optical data. On this offshore platform, a Sea-Viewing Wide Field-of-View Sensor Photometer Revision for Incident Surface Measurements (SeaPRISM) autonomous radiometer system was deployed; its sea-viewing measurement sequence was executed every 30 min from 08:00 to 16:00 local time. The system collected water surface and sky radiation data; after processing to obtain L wn data, it divided by the mean extraterrestrial solar irradiance (F 0 ) of the corresponding waveband to obtain R rs data. The SeaPRISM system configuration performed ocean color measurements at 11 wavelengths from 400-1020 nm. The SeaPRISM data used in this study were collected during September 2019 to October 2020 and were manually checked to make sure that no corrupted spectra were present in this dataset. Only the L wn and R rs data at 412, 442, 490, 560, and 667 nm center wavelengths were used for the match-up comparison analysis with VIIRS. The overall uncertainty of SeaPRISM L wn data was estimated to be around~5% for all wavelengths shorter than 668 nm [34], which meets the requirement of VIIRS data evaluation. By taking advantage of high frequency measurements of SeaPRISM, all in situ data used in the quantitative match-up comparison analysis were selected from the measurements made within a short (±0.5 h) time window of the satellite overpass time of the locations of the Dongou site, in order to avoid unexpected effects induced by the highly dynamic changes in coastal waters. Another part of R rs data comes from the underway observation of the LORCE cruise in the ECS in August 2021; the instrument used a shipboard fully automated radiometric measurement system called CrusieAOP-T (Figure 2a), collecting various types of data such as the measured total radiance leaving the sea surface L sfc (λ), sky radiance L sky (λ), surface incident irradiance E s (λ), and R rs (λ), which was calculated from above-water measurements of the above three radiometric parameters following the NASA Ocean Optics Protocols [35]. The information about HAB events obtained by the above equipment is shown in Table 1 Moreover, a fixed oceanographic platform called Dongou was used to collect a longterm series of radiometric measurements (such as Rrs and normalized water-leaving radiance (Lwn)). This platform was located in the southern part of the ECS (27.675°N and 120.358°E) at approximately 80 km southeast of Wenzhou in an area that minimized the adjacency effects in satellite data. In addition, the water depth around the platform was around 50 m, so that the bottom effects in Lwn were negligible. This oceanographic platform, owned and operated by WMEMC, has served as an offshore laboratory since 2017 and is an excellent site supporting ocean color validation activities through a comprehensive collection of bio-optical data. On this offshore platform, a Sea-Viewing Wide Fieldof-View Sensor Photometer Revision for Incident Surface Measurements (SeaPRISM) autonomous radiometer system was deployed; its sea-viewing measurement sequence was executed every 30 min from 08:00 to 16:00 local time. The system collected water surface and sky radiation data; after processing to obtain Lwn data, it divided by the mean extraterrestrial solar irradiance (F0) of the corresponding waveband to obtain Rrs data. The Sea-PRISM system configuration performed ocean color measurements at 11 wavelengths from 400-1020 nm. The SeaPRISM data used in this study were collected during September 2019 to October 2020 and were manually checked to make sure that no corrupted spectra were present in this dataset. Only the Lwn and Rrs data at 412, 442, 490, 560, and 667 nm  of the LORCE cruise in the ECS in August 2021; the instrument used a shipboard fully automated radiometric measurement system called CrusieAOP-T (Figure 2a), collecting various types of data such as the measured total radiance leaving the sea surface Lsfc(λ), sky radiance Lsky(λ), surface incident irradiance Es(λ), and Rrs(λ), which was calculated from above-water measurements of the above three radiometric parameters following the NASA Ocean Optics Protocols [35]. The information about HAB events obtained by the above equipment is shown in Table 1. The four stations numbered S1-S4 were selected to represent turbid water, medium turbid water, algal bloom water, and clear water, respectively.  The four stations numbered S1-S4 were selected to represent turbid water, medium turbid water, algal bloom water, and clear water, respectively.

Satellite Data
VIIRS onboard SNPP (VIIRS/SNPP) and VIIRS onboard NA20 (VIIRS/NA20) granules covering the dates from September 2019 to October 2020 were downloaded from Level-2 from NASA's Goddard Space Flight Center archives (https://oceancolor.nasa.gsfc.gov (accessed on 31 May 2021)). The two sensors are generally identical and have only minor difference in the central wavelength of each band. In addition, data from MODIS onboard Aqua (MODISA) of the same period were downloaded to compare with VIIRS. The specific parameters of the three sensors are given in Table 2. These files conform to calibration updates for which atmospheric correction was performed with the iterative Near Infrared approach [36][37][38][39]. Furthermore, only VIIRS/SNPP "science quality" data for these dates and locations were also acquired from National Oceanic and Atmospheric Administration (NOAA) CoastWatch (https://coastwatch.noaa.gov (accessed on 31 May 2021)). These data correspond to the April 2017 Software-Defined Radio (SDR) and calibration update, with atmospheric correction performed using the Near Infrared-Short Wave Infrared procedure [6,36,40,41]. More importantly, one of major differences between NASA and Remote Sens. 2022, 14, 2089 6 of 23 NOAA processing schemes of VIIRS/SNPP data is the vicarious calibration gains, detail values of which for 410, 443, 486, 551, 671, 745, and 862 nm channels can be found in the Level-2 granules (NetCDF4 attributes). Within this manuscript, VIIRS/SNPP data from these two sources are termed 'VIIRS-SNPP L2GEN' and 'VIIRS-SNPP MSL12 , respectively, as well as VIIRS/NA20 Level-2 data from the NASA website which is termed 'VIIRS-NA20 L2GEN'. The VIIRS/SNPP band centers (410, 443, 486, 551, and 671) are slightly different from associated band centers of MODIS (412, 443, 488, 547, and 667), SeaPRISM, and even VIIRS/NA20 (411, 445, 489, 556, and 667). For convenience, we describe the bands of other sensors based on the central band name of VIIRS/SNPP, for example, 410 nm of VIIRS/NA20 means it is the 411-nm band.

Data Processing
Barnes et al. proposed that in situ-and satellite-derived R rs data should be assessed according to a variety of exclusion criteria, including Level-2 Processing Flags, quality assessment method, and spatial homogeneity [11]. Following this proposal, collocated Level-2 satellite pixel(s) were first masked by Level-2 Processing Flags (ATMFAIL, HILT, and CLDICE (termed "CLOUD" in MSL12 datasets)) in order to remove low-quality data at the time of data extraction. Recently, a quality assessment method for in situ-and satellitederived R rs was provided by Wei et al., which can yield an R rs quality score (hereafter termed 'QA_Wei', meaning the quality assurance system proposed by Wei et al.) and water type for each spectrum [42]. Note that R rs (λ) or L wn data could be converted to each other (L wn (λ) = R rs (λ)×F 0 (λ)) using spectral response integrated F 0 values [43]. Thus, the remaining spectra were assessed according to their QA_Wei scores and water types, and then all pixels with a QA_Wei score ≥ 0.5 were remained in the dataset.
Additionally, for the matchup analysis with in situ measurements from SeaPRISM, data which remained after masking by the L2 Flags and QA_Wei were partitioned according to spatial homogeneity, assessed as the coefficient of variation (CV = standard deviation/mean) of the 3 × 3 pixel boxes with the matchup location in the center [44]. Note that CV calculations were carried out only if ≥5 of the nine pixels in the 3 × 3 pixel boxes were not flagged and the associated R rs at 551 nm or the closest center wavelength in each sensors was used for the spatial homogeneity CV calculation. Matchup analysis was implemented for those pixels with a low CV. In fact, the spatial homogeneity test and QA criteria served as a good effort to minimize the impact of differences between in situ measurements (point sample) and satellite measurements (integrated ≥ 1 km 2 pixel). The coefficient of temporal variation was calculated as follows: (1) where CV i,j is the coefficient of spatial variation in the (i, j) pixel, x is the mean value of the window, and S i,j is the statistical standard deviation of the image data within the window. The match-up comparison analyses were based on linear regression between any two datasets being compared. The average absolute percent differences, denoted as APD, and average relative percent differences, denoted as RPD of N total number of matchups, were the primary measures used to assess satellite uncertainty and bias, respectively, as where x i and y i are the in situ and satellite data, respectively, for matchup i of N. For direct comparison of other published validation results, coefficients of correlation (R 2 ) and Root Mean Squared Error (RMSE) were also computed at each wavelength for the comparisons in order to provide information on how well the data being compared match. Furthermore, the two-standard deviation filtering procedure was employed to filter out some extreme cases where the statistics were out of the range of the majority of cases (only 6 matchup points out of 141 were filtered out) [45].

HAB Detection Algorithm
Currently used algal bloom detection algorithms can be divided into two main types, the band ratio and spectral shape difference methods; the available typical algorithms are limited due to a lack of characteristic terrestrial bands and fluorescence bands in VIIRS. Two band ratio methods, Ratio of Algal Bloom (RAB) [28], RGCI [21], and two spectral shape index methods, Spectral Shape index at 486 nm (SS(486)) [19] and Red tide Index (RI) [46], were selected for subsequent comparison in this paper (Table 3). However, as mentioned above, VIIRS lacks the 531 nm band, so the 486 nm band was used in the RAB algorithm instead of the MODIS 531 nm band. Table 3. Four algorithms related to harmful algal blooms and their formulas for comparison. (Note: RAB, Ratio of Algal Bloom; SS(486), Spectral Shape index at 486 nm; RI, Red tide Index; RGCI, Red-Green-Chlorophyll-a Index.)

Index
Algorithm MODISA VIIRS-SNPP VIIRS-NA20 After correct detection of the HAB area, this paper used VIIRS data for the taxonomic identification of P. donghaiense, diatoms, and other algal species in the ECS based on the PDI and DI indices of MODIS [28]. The PDI index and DI index were calculated as follows: Since the NASA L2GEN products of VIIRS were missing a band near 640 nm, the VIIRS data used in the subsequent algal bloom types identification were level 2 data of VIIRS-SNPP MSL12 from the NOAA CoastWatch website; in addition, since the 531 nm band was missing in VIIRS, the PDI index of MODIS was replaced by SS(486). The MODIS based PDI and DI indices were also calculated for the comparison with VIIRS results in the same day.

Assessment of Normalized Water-Leaving Radiance
Since the reliability of remote sensing monitoring of algal bloom mainly depends on the degree of the quantification of satellite derived R rs , the R rs at bands with lowest uncertainty and best degree of quantification in the study area should be preferred for constructing the algorithm. Validating the accuracy and spectral characteristics of VIIRS retrieved L wn (λ) or the equivalent R rs (λ) is, thus, an important prerequisite for carrying out algal bloom monitoring. A first step when evaluating the quality of the satellite data is to examine the spectral consistency of the water-leaving radiances from the VIIRS and MODIS missions based on the in situ data collected from the Dongou site. The Dongou site is frequently characterized by a large variability in bio-optical quantities because of its position in a transitional region between open sea and coastal waters (Figure 3a,e,i). Thus, the L wn spectra from VIIRS and MODIS (Figure 3b,f,j) also show high variability. It can be found that the spectra exhibiting large variations show minima at 410 and 671 nm, with the majority of values below 20 W/m 2 /nm/sr and unique maxima at 551 nm with values in the range of 4.0-50 W/m 2 /nm/sr, and suggest the presence of seawater that is moderately dominated by sediments, as well as the derived water type provides a further confirmation with an occurrence of roughly 70% Case 2 water (of which MODISA, VIIRS-SNPP, and VIIRS-NA20 spectral data accounted for 75%, 83%, and 79% of this type, respectively) [42]. The comparison of L wn spectra indicates that the spectral variation ranges of MODISA, VIIRS/SNPP, and VIIRS/NA20 are very consistent to each other. Furthermore, the matchup L wn spectra from SeaPRISM was selected (Figure 3c,g,k); matchups of the overall average L wn spectra were also calculated from all available spectral data of satellites and SeaPRISM, as seen in the fourth column of Figure 3d,h,l. These matchup comparisons indicate that a general spectral concordance exists between satellite and in situ data. However, in the average L wn spectra match-up comparisons, a spectral discrepancy also can be observed at the shorter wavelength bands (410 and 443 nm).
Thus, additional insight from the qualitative matchup analysis of satellite L wn (λ) data against in situ SeaPRISM for each individual spectral band is presented in Figure 4 and Table 4. After the quality control methods were applied, a total of 135 in situ L wn spectra matched up with at least one satellite dataset, in which there are 34, 60, and 41 matchups available for MODIS/Aqua, VIIRS/SNPP MSL12, and VIIRS/NA20, respectively. It is noticeable that only 44 matchups were available for VIIRS-SNPP L2GEN, less than that for VIIRS-SNPP MSL12, because much tighter restrictions in the flag conditions applied in the filtering cloud and bright pixels were implemented in the L2GEN procedure. Additionally, the number of matchups of MODSA was less than the other three VIIRS sources, probably because the coastal water of the ECS is dominated by sediments (as seen in Figure 3a,e,i), Remote Sens. 2022, 14, 2089 9 of 23 which will more easily result in saturation of MODIS 862 nm and thus makes fewer matchups with L wn spectra available.
confirmation with an occurrence of roughly 70% Case 2 water (of which MODISA, VIIRS-SNPP, and VIIRS-NA20 spectral data accounted for 75%, 83%, and 79% of this type, respectively) [42]. The comparison of Lwn spectra indicates that the spectral variation ranges of MODISA, VIIRS/SNPP, and VIIRS/NA20 are very consistent to each other. Furthermore, the matchup Lwn spectra from SeaPRISM was selected (Figure 3c,g,k); matchups of the overall average Lwn spectra were also calculated from all available spectral data of satellites and SeaPRISM, as seen in the fourth column of Figure 3d,h,l. These matchup comparisons indicate that a general spectral concordance exists between satellite and in situ data. However, in the average Lwn spectra match-up comparisons, a spectral discrepancy also can be observed at the shorter wavelength bands (410 and 443 nm). Thus, additional insight from the qualitative matchup analysis of satellite Lwn(λ) data against in situ SeaPRISM for each individual spectral band is presented in Figure 4 and Table 4. After the quality control methods were applied, a total of 135 in situ Lwn spectra matched up with at least one satellite dataset, in which there are 34, 60, and 41 matchups available for MODIS/Aqua, VIIRS/SNPP MSL12, and VIIRS/NA20, respectively. It is noticeable that only 44 matchups were available for VIIRS-SNPP L2GEN, less than that for VIIRS-SNPP MSL12, because much tighter restrictions in the flag conditions applied in    In addition, percent differences (APD and RPD) suggest good qualities of Lwn retrievals from satellites at 488 nm, 551 nm, and 671 nm, although high uncertainties and bias still exist at 410 nm and 443 nm ( Table 4). The highest APD values were observed at 443 nm (with 39.48%, 61.55%, 59.89%, and 48.88% for MODISA, VIIRS-NPP L2GEN, VIIRS-NPP MSL12, and VIIRS-N20 L2GEN, respectively), followed by 410 nm with a minimum Overall, the correlation coefficients (R 2 ) between the SeaPRISM and satellite datasets at each wavelength were relatively high, and all comparisons, particularly at 488 nm, 551 nm, and 671 nm, were also very close to a 1:1 line (Figure 4). Notably, the MODISA data performed quite well even though this MODIS sensor has been in operation during orbit for more than 17 years and has long passed its expected lifetime. At individual wavelengths, correlations for the four satellites were very close to each other, which indicate good consistency between MODIS and the two VIIRS sensors; in addition, the variations in the water-leaving radiance data for the Dongou location were well captured by all satellite sensors. Nevertheless, the resulting R 2 values showed strong spectral dependencies exhibiting the tendency of displaying larger differences in the shorter end of the spectrum. For example, VIIRS-SNPP MSL12 achieved stronger R 2 values of 0.91, 0.91, and 0.95 at 486 nm, 551 nm, and 671 nm, respectively, while at 410 nm and 443 nm only moderate correlations were attained with R 2 values equal to 0.61 and 0.82, respectively. Similar spectral behavior and values of R 2 can be found in MODISA, VIIRS-SNPP L2GEN, and VIIRS-NA20 L2GEN comparisons. This degradation in the correlation for 410 nm and 443 nm agrees with more recent findings from the Long Island Sound coasts, USA, which was considered to have originated from the data processing procedure [10].
In addition, percent differences (APD and RPD) suggest good qualities of L wn retrievals from satellites at 488 nm, 551 nm, and 671 nm, although high uncertainties and bias still exist at 410 nm and 443 nm ( Table 4). The highest APD values were observed at 443 nm (with 39.48%, 61.55%, 59.89%, and 48.88% for MODISA, VIIRS-NPP L2GEN, VIIRS-NPP MSL12, and VIIRS-N20 L2GEN, respectively), followed by 410 nm with a minimum of 28.14%. The larger uncertainties further confirmed the moderate correlations at the two deep blue bands. In the individual band match-up comparisons (Figure 4), it can be found in all the four datasets that had satellite-derived L wn at 410 nm and 443 nm as a whole were larger than that from SeaPRISM; a similar trend was shown in comparisons of the average L wn spectra ( Figure 3). Although these overestimates could be partially attributed to the uncertainties in the SeaPRISM datasets, limited retrieval accuracies were achieved at the two blue bands in both MODIS and VIIRS, indicating the challenge of using current remote sensing L wn at 410 nm and 443 nm to derive bio-optical products for the ECS. Despite these problems, as compared to shorter wavelengths, matchup statistics were much better for the 486 nm, 551 nm, and 671 nm wavelengths.
In order to fully access the VIIRS R rs products, the supplementary underway R rs measurements in the ECS were used. Figure 5 shows that the section of underway R rs measurements crossed the coastal area from turbid to clean water bodies; it fortunately covered turbid, moderately turbid, clean, and HAB water bodies at the same time. Although it does not cover extremely turbid water bodies, the HAB events were very infrequent in such turbid waters due to light limitations; the satellite R rs data were usually not available. Figure 5a-f shows the validation results of the section of underway R rs observation against the VIIRS MSL12 R rs data. Similar to the validation results at the Dongou platform (Figure 4), the VIIRS R rs data of 410 nm and 443 nm showed a relatively high bias, while those at the 486 nm, 551 nm, and 671 nm bands were in very good agreement with the underway data. The spectral comparison in the four sites associated with turbid, moderately turbid, HAB, and clean water bodies is also presented in Figure 5g-j. In terms of magnitude and spectral shape, R rs results obtained from VIIRS were very consistent with field measurements.
In the aggregate, these analyses show that as attributing to the most recent calibration efforts (and associated atmospheric correction routines) and reprocessing efforts of both NOAA and NASA, the MODISA, VIIRS/SNPP, and VIIRS/NA20 datasets from either NASA or NOAA yield reliable performance and consistent environmental data records, showing the band data of 486 nm, 551 nm, and 671 nm had better potential for the long-term monitoring in coasts of the ECS.  In the aggregate, these analyses show that as attributing to the most recent calibration efforts (and associated atmospheric correction routines) and reprocessing efforts of both NOAA and NASA, the MODISA, VIIRS/SNPP, and VIIRS/NA20 datasets from either (g-j) Spectra of S1-S4, which represent turbid water, medium turbid water, algal bloom water, and clear water, respectively.

HAB Algorithm Development and Validation
Comparative analysis of various VIIRS algorithms was carried out to assess advantages and limitations in the application of these techniques for the detection of algal blooms in the sediment-dominated coasts of the ECS. Recently, different forms of HAB algorithms have become available, which specially include reflectance band-ratios and spectral band difference models [47]. As seen in Table 2, the two band-ratio algorithms, Algal Bloom Ratio (RAB) [28] and RGCI, and the two spectral band difference ones, Red Tide index (RI) [46] and normalized spectral shape at 488 nm (SS(488)), were selected for the comparison. Due to the similar band configuration of VIIRS with MODIS, most of these MODIS based algorithms can be directly applied to VIIRS. However, as mentioned above, VIIRS lacks the 531 nm band so that the 488 nm and 551 nm bands were used instead of MODIS 531 and 555 nm bands in the RAB algorithm. Figure 6 shows scatterplots of in situ/satellite matchups for the four HAB indices, in which the algal bloom observations confirmed by field measurements at the Dongou site are plotted as red circles, as well as green circles indicating turbid waters with R rs (551) greater than 0.014 sr −1 , and blue circles were associated with clear ones with R rs (551) lower than 0.014 sr −1 . General consistency between satellite data and the SeaPRISM data can be found for all the four indices, but a difference exists in the capacity of the separation of bloom waters from normal ones. However, for sediment-dominated waters, the optical properties in the red bands were determined by both phytoplankton and sediment, and thus the separation of bloom areas from turbid waters by RGCI was not satisfactory. Although the performance of RI and SS(486) seems to be better, the separation in the two methods makes both not clear enough, probably because the blue band of 443 nm with relative high uncertainties was used in both of them. Thus, it can be clearly found that the RAB method probably yielded the best discrimination and the threshold value of 1.6 seems suitable for both MODIS and VIIRS for the identification of algal blooms. This success can be attributed not only to good quality of 488 nm and 551 nm data (Figure 4), but also to the much shearer slope between 488 nm and 551 nm of phytoplankton absorption [28]. To further access the performance of the simple VIIRS RAB method in identifying HAB in the ECS, image series and near-concurrent field surveys for large-scale P. donghaiense blooms in April to May 2020 were used for method validation (Figure 7). These blooms were massive, extending over thousands of square kilometers, and persisted for To further access the performance of the simple VIIRS RAB method in identifying HAB in the ECS, image series and near-concurrent field surveys for large-scale P. donghaiense blooms in April to May 2020 were used for method validation (Figure 7). These blooms were massive, extending over thousands of square kilometers, and persisted for nearly one month. The field measurements of cell counts from WMEMC confirmed that the RAB method can successfully differentiate bloom and non-bloom waters. Particularly, even in a very near-shore region, as seen in Figure 7e,f, small-area bloom events can also be captured in the RAB images. Similar success was achieved in the diatom bloom event occurring near the YRE in August 2020 (Figure 8a-c), indicating that RAB is also suitable for diatom bloom detection. In addition, RAB can also identify blooms dominated by other algal species, such as C. furca and A. sanguinea, validation results of which are shown in Figure 8d,e,g-i, respectively. Although there is an hourly error in the matchups of VIIRS derived RAB data and field surveys, these results are still sufficient to confirm the ability of VIIRS to detect HABs in the ECS.  Figure A1 in the Appendix A. In the algal bloom ratio images, only the areas positively flagged as bloom waters are shown in yellow and red. The black and green circles indicate bloom and non-bloom sites (5 × 10 5 cells/L is taken as the reference concentration of HAB [33]), which have been confirmed by field surveys; the cell count data are plotted close to the circle in each image. (Satellite data from the square area on the map were used for the subsequent classification of algal bloom types.).

Capability in Bloom Type Discrimination
Although the reduction in the number of VIIRS bands has little impact on the identification of HABs in the ECS, it was found that it has a significant impact on the discrimination of bloom types or dominant species. Figure 9a shows the distribution of HAB events on 24 May 2019; note that the dominant algal species in region A was identified as P. donghaiense, while the one of region B in the other end of the stretch was caused by diatoms. In the above two regions, the typical Rrs spectra of P. donghaiense and diatom blooms from both MODIS and VIIRS were normalized by the maximum value of each spectrum in the green bands (Figure 9b,c). Other Rrs spectra of` P. donghaiense blooms that occurred on 28 April 2020 ( Figure 7a) and a diatom bloom that occurred on 18 August 2020 (Figure 8b) are also presented. According to the finding of Tao et al. [28], the MODIS derived spectra of P. donghaiense and diatoms show a large difference mainly near the 531 nm and 645 nm bands; the Rrs of diatoms has a high shoulder peak at the 645 nm while that of P. donghaiense has a more obvious trough at 531 nm (Figure 9b). Therefore, the P. donghaiense bloom can be concluded as having extremely high PDI value and relatively low DI value compared to those of diatoms; then, the two type blooms can be clearly classified in the scatter plot of DI against PDI (Figure 9d). For VIIRS, due to a lack of a band

Capability in Bloom Type Discrimination
Although the reduction in the number of VIIRS bands has little impact on the identification of HABs in the ECS, it was found that it has a significant impact on the discrimination of bloom types or dominant species. Figure 9a shows the distribution of HAB events on 24 May 2019; note that the dominant algal species in region A was identified as P. donghaiense, while the one of region B in the other end of the stretch was caused by diatoms. In the above two regions, the typical R rs spectra of P. donghaiense and diatom blooms from both MODIS and VIIRS were normalized by the maximum value of each spectrum in the green bands (Figure 9b,c). Other R rs spectra of' P. donghaiense blooms that occurred on 28 April 2020 ( Figure 7a) and a diatom bloom that occurred on 18 August 2020 (Figure 8b) are also presented. According to the finding of Tao et al. [28], the MODIS derived spectra of P. donghaiense and diatoms show a large difference mainly near the 531 nm and 645 nm bands; the R rs of diatoms has a high shoulder peak at the 645 nm while that of P. donghaiense has a more obvious trough at 531 nm (Figure 9b). Therefore, the P. donghaiense bloom can be concluded as having extremely high PDI value and relatively low DI value compared to those of diatoms; then, the two type blooms can be clearly classified in the scatter plot of DI against PDI (Figure 9d). For VIIRS, due to a lack of a band near 531 nm, the difference between P. donghaiense and diatoms in their VIIRS R rs spectra was only observed at 638 nm (Figure 9c). After the replacement of the PDI index by SS(486) for VIIRS, the separation of P. donghaiense and diatoms in the scatter plot of SS(486) against DI is not clear enough when compared with that in MODIS (Figure 9e). In contrast to PDI, SS(486) did not have obvious data that could be used to distinguish between P. donghaiense and diatoms; nevertheless, owing to the presence of the DI index, a clear trend still exists that can be used to separate the two types of blooms. Although VIIRS still has the ability to distinguish between P. donghaiense and diatoms in the ECS to some extent, the lack of a 531 nm band greatly reduces its ability to discriminate other types of algal blooms. near 531 nm, the difference between P. donghaiense and diatoms in their VIIRS Rrs spectra was only observed at 638 nm (Figure 9c). After the replacement of the PDI index by SS(486) for VIIRS, the separation of P. donghaiense and diatoms in the scatter plot of SS(486) against DI is not clear enough when compared with that in MODIS (Figure 9e). In contrast to PDI, SS(486) did not have obvious data that could be used to distinguish between P. donghaiense and diatoms; nevertheless, owing to the presence of the DI index, a clear trend still exists that can be used to separate the two types of blooms. Although VIIRS still has the ability to distinguish between P. donghaiense and diatoms in the ECS to some extent, the lack of a 531 nm band greatly reduces its ability to discriminate other types of algal blooms. Based on the observation of the C. furca blooms in July 2013, its Rrs spectra derived from MODIS showed some distinguishing spectral features that can be used to separate them from P. donghaiense and diatom blooms (Figure 10a). Compared with diatoms, a C. furca bloom does not have a prominent shoulder peak at 645 nm, while compared with P. donghaiense, it has no reflection trough at 531 nm. These features can be directly reflected in the scatter plot of PDI to DI (Figure 10c). The clusters of C. furca blooms are roughly located in the lower left quarter of Figure 10c so that the separation between C. furca and diatom blooms is very clear, although some points of a C. furca bloom are not well separated from clusters associated with P. donghaiense. Nevertheless, the classification of C. furca, P. donghaiense, and diatom blooms cannot be achieved based on the VIIRS data, since the VIIRS derived Rrs spectra between 488 nm and 551 nm of a C. furca bloom do not yield different features from the other two types of blooms when the 531 nm band is not available (Figure 10b,d). Based on the observation of the C. furca blooms in July 2013, its R rs spectra derived from MODIS showed some distinguishing spectral features that can be used to separate them from P. donghaiense and diatom blooms (Figure 10a). Compared with diatoms, a C. furca bloom does not have a prominent shoulder peak at 645 nm, while compared with P. donghaiense, it has no reflection trough at 531 nm. These features can be directly reflected in the scatter plot of PDI to DI (Figure 10c). The clusters of C. furca blooms are roughly located in the lower left quarter of Figure 10c so that the separation between C. furca and diatom blooms is very clear, although some points of a C. furca bloom are not well separated from clusters associated with P. donghaiense. Nevertheless, the classification of C. furca, P. donghaiense, and diatom blooms cannot be achieved based on the VIIRS data, since the VIIRS derived R rs spectra between 488 nm and 551 nm of a C. furca bloom do not yield different features from the other two types of blooms when the 531 nm band is not available (Figure 10b,d). For A. sanguinea blooms, there is another interesting feature found near 488 nm in Rrs spectra of both MODIS and VIIRS (Figure 11a,b), which makes both PDI and SS(486) of A. sanguinea become significantly lower than that of P. donghaiense. Similar to the classification results of C. furca, the distribution of A. sanguinea points in the PDI-DI scatter plot of MODIS (Figure 11c) is at the lower left of P. donghaiense and diatom clusters, and more importantly, the distinction between the A. sanguinea and other two type blooms is more evident. Comparatively, the distribution of A. sanguinea points in the SS(486)-DI scatter plot of VIIRS partially overlaps with the P. donghaiense and diatom points (Figure 11d), but the separations of A. sanguinea from the P. donghaiense and diatom blooms is also better than that of C. furca (Figure 10d). Although the classification results of VIIRS showed a certain trend of separation, but was still not enough to separate A. sanguinea from P. donghaiense and diatoms. Based on the above results, one can conclude that the 531 nm band is essential for the bloom type discrimination, a lack of which significantly limits the capability of VIIRS. For A. sanguinea blooms, there is another interesting feature found near 488 nm in R rs spectra of both MODIS and VIIRS (Figure 11a,b), which makes both PDI and SS(486) of A. sanguinea become significantly lower than that of P. donghaiense. Similar to the classification results of C. furca, the distribution of A. sanguinea points in the PDI-DI scatter plot of MODIS (Figure 11c) is at the lower left of P. donghaiense and diatom clusters, and more importantly, the distinction between the A. sanguinea and other two type blooms is more evident. Comparatively, the distribution of A. sanguinea points in the SS(486)-DI scatter plot of VIIRS partially overlaps with the P. donghaiense and diatom points (Figure 11d), but the separations of A. sanguinea from the P. donghaiense and diatom blooms is also better than that of C. furca (Figure 10d). Although the classification results of VIIRS showed a certain trend of separation, but was still not enough to separate A. sanguinea from P. donghaiense and diatoms. Based on the above results, one can conclude that the 531 nm band is essential for the bloom type discrimination, a lack of which significantly limits the capability of VIIRS. Remote Sens. 2022, 14, x FOR PEER REVIEW 20 of 25

Conclusions
This paper involved a systematic assessment of VIIRS for its ability to detect HAB and classify bloom types in coastal waters of the ECS. Although the data in the 412 and 443 nm bands had relatively high uncertainties, we found that the data quality at the 486, 551, and 671 nm was reliable. The above three bands show strong stability and consistency in the comparison of data accuracy and spectral characteristics of VIIRS. More importantly, we determined the accuracy of VIIRS data in different water bodies through a comparison with Rrs data collected during LORCE underway observation, which showed that the three bands were fully competent for construction of the algorithm for long-term HAB monitoring due to their high data quality. The applicability of different HAB algorithms in the ECS also varied. The two spectral shape algorithms, RI and SS(486), were less effective because they used the 443 nm with poor stability, and because of the red light absorption characteristics of suspended sediment, RGCI was not suitable for the detection of HABs in sediment-rich coasts of the ECS. Only RAB, thanks to the high quality of 486 nm, 551 nm, and the shearer slope of the two bands, has the best recognition effect and can accurately identify different types of algae. It is noteworthy that, due to the absence of the band near 531 nm, VIIRS was able to separate P. donghaiense from diatoms, but the ability to discriminate a third type of algae was significantly limited more than that of MODIS. We could hardly discriminate C. furca or A. sanguinea from P. donghaiense and diatoms using VIIRS data, and although A. sanguinea showed some separation trend, Figure 11. (a,b) Normalized spectra of MODIS and VIIRS derived R rs (λ) associated with P. donghaiense, diatoms, and A. sanguinea, respectively. (c,d) Scatter plot of DI against PDI or SS(486) derived from MODIS or VIIRS data, respectively. The satellite data of P. donghaiense and diatom are the same as Figure 9, while the data of A. sanguinea are taken from the rectangle regions of Figure 8g,h on 28 September 2021 and 2 October 2021, respectively.

Conclusions
This paper involved a systematic assessment of VIIRS for its ability to detect HAB and classify bloom types in coastal waters of the ECS. Although the data in the 412 and 443 nm bands had relatively high uncertainties, we found that the data quality at the 486, 551, and 671 nm was reliable. The above three bands show strong stability and consistency in the comparison of data accuracy and spectral characteristics of VIIRS. More importantly, we determined the accuracy of VIIRS data in different water bodies through a comparison with R rs data collected during LORCE underway observation, which showed that the three bands were fully competent for construction of the algorithm for long-term HAB monitoring due to their high data quality. The applicability of different HAB algorithms in the ECS also varied. The two spectral shape algorithms, RI and SS(486), were less effective because they used the 443 nm with poor stability, and because of the red light absorption characteristics of suspended sediment, RGCI was not suitable for the detection of HABs in sediment-rich coasts of the ECS. Only RAB, thanks to the high quality of 486 nm, 551 nm, and the shearer slope of the two bands, has the best recognition effect and can accurately identify different types of algae. It is noteworthy that, due to the absence of the band near 531 nm, VIIRS was able to separate P. donghaiense from diatoms, but the ability to discriminate a third type of algae was significantly limited more than that of MODIS. We could hardly discriminate C. furca or A. sanguinea from P. donghaiense and diatoms using VIIRS data, and although A. sanguinea showed some separation trend, the effect was still insufficient to distinguish them. Therefore, it is suggested that subsequent ocean color   Appendix B