Fluorescence-based Approach to Estimate the Chlorophyll-a Concentration of a Phytoplankton Bloom in Ardley Cove (antarctica)

A phytoplankton bloom occurred in Ardley Cove, King George Island in January 2016, during which maximum chlorophyll-a reached 9.87 mg/m 3. Records show that blooms have previously not occurred in this area prior to 2010 and the average chlorophyll-a concentration between 1991 and 2009 was less than 2 mg/m 3. Given the lack of in situ measurements and the poor performance of satellite algorithms in the Southern Ocean and Antarctic waters, we validate and assess several chlorophyll-a algorithms and apply an improved baseline fluorescence approach to examine this bloom event. In situ water properties including in vivo fluorescence, water leaving radiance, and solar irradiance were collected to evaluate satellite algorithms and characterize chlorophyll-a concentration, as well as dominant phytoplankton groups. The results validated the nFLH fluorescence baseline approach, resulting in a good agreement at this high latitude, high chlorophyll-a region with correlation at 59.46%. The dominant phytoplankton group within the bloom was micro-phytoplankton, occupying 79.58% of the total phytoplankton community. Increasing sea ice coverage and sea ice concentration are likely responsible for increasing phytoplankton blooms in the recent decade. Given the profound influence of climate change on sea-ice and phytoplankton dynamics in the region, it is imperative to develop accurate methods of estimating the spatial distribution and concentrations of the increasing occurrence of bloom events.


Introduction
Due to the extreme climate and the difficulties of conducting field research above 60 • south, the Southern Ocean (SO), especially around Antarctica, lacks a systematic in situ sampling program of its peculiar bio-optics and micro-organism community structure [1][2][3].As a result, satellite measurements from space still have large errors in estimating phytoplankton biomass [4,5] and global chlorophyll-a satellite algorithms typically underestimate chlorophyll-a in the Southern Ocean [6][7][8][9][10][11].
There are three main reasons for global algorithm underestimation.First is the difference in optical properties between global and SO waters.Compared with the global ocean, the SO has a narrower water leaving radiance in the green band [12].This narrow gap leads to underestimation of chlorophyll-a in blue-green band ratio algorithms such as OC3 or OC4 [13,14].Second, seasonal sea ice causes contaminated pixels, which underestimate chlorophyll-a in > 1.5 mg/m 3 and overestimates in < 1.5 mg/m 3 waters [15,16].These are obvious patterns in the Arctic and are probably similar in Antarctica.Third, inappropriate atmospheric correction in the SO introduces error when converting top of atmosphere radiance into the water leaving radiance.The SO has unique aerosols and cloud coverage [17], and lacks a proper vertical atmospheric simulation to correct for aerosol influences [18].
In addition to band-ratio chlorophyll-a algorithms, state-of-the-art algorithms for global chlorophyll-a estimation include IOP (inherent optical property) bio-optical models (e.g., GSM01 [19] and QAA (quasi-analytical algorithm) [20]) and baseline algorithms, such as nFLH (normalized fluorescence line height).Bands setting at 667, 678, and 748 nm help MODIS become the only satellite to achieve the baseline nFLH approach, which estimates chlorophyll-a through fluorescence intensity from photosynthesis products [21].
where, F(λ em ) is the fluorescence intensity (mol quanta m −3 •s −1 ), E(λ ex ) is the incident intensity (mol quanta m −2 •s −1 ), a * (λ ex ) is are chlorophyll-a absorption coefficients (m 2 •mg•Chla −1 ) per chlorophyll unit (* average per [Chla]), [Chla] is the chlorophyll-a concentration (mg•m −3 ), ϕ f is the fluorescence yield (mol quanta), and Q * a (λ em ) is the re-absorption coefficients in cells.From this equation, a * (λ ex )Q * a (λ em ) reflects the phytoplankton composition features and E(λ ex )ϕ f is from the light acclimation mechanism (photosynthetic adjustment in response to light availability).Under stable phytoplankton community composition and light acclimation, chlorophyll-a will have a linear relationship with fluorescence intensity.The nFLH approach helps to avoid contamination of the chlorophyll-a signal by suspended sediments, detritus, and CDOM (colored dissolved organic matter) and typically produces more accurate results in case 2 coastal waters [22].However, photosynthetic mechanisms are subject to NPQ (non-photochemical quenching) when phytoplankton encounters intense light [23].Non-linearity then occurs between chlorophyll-a and fluorescence intensity, limiting the applications of nFLH in high-intensity light areas, such as midday direct solar radiation.
A phytoplankton bloom occurred in Ardley Cove near King George Island (KGI) in January 2016.Historical records showed that phytoplankton blooms had not previously occurred in this area prior to 2010 and the average of chlorophyll-a between 1991 and 2009 was less than 2 mg/m 3 [24].An obvious phytoplankton bloom was reported until 2010, during which maximum chlorophyll-a reached 20 mg/m 3 .In situ records show an increasing trend of phytoplankton biomass in this area, likely caused by increasing SST (sea surface temperature) related to global warming [25].The aim of this study was to document an algal bloom in Ardley Cove through the validation and regionalization of the MODIS nFLH algorithm.

Materials and Methods
All the samples were collected on a Zodiac (an inflatable boat) in the Great Wall Cove and Ardley Cove near the China Great Wall Station on KGI between 6 and 27 January 2016 (Figure 1a,b).Clear sky, high solar evaluation (40 • -50 • in 10 am-2 pm), and light breeze (wind speed < 5m/s) were chosen as the threshold conditions during sampling to reduce the impact on above-water optical property retrieval.

Apparent Optical Properties
We applied (45°, 135°) angles for above-water measurements recommended by NASA Ocean color protocol [26], to normalize water radiance from simultaneous sky radiance and to remove interference from environmental light conditions.A hand-held VNIR spectroradiometer (HH2, ASD Inc., Boulder, Colorado, USA), with a high spectral resolution of 1 nm, was used to measure water leaving radiance.When measuring, the Zodiac stopped its engine to prevent white capping.However, without power, the boat did not drift far from its original location.The black hull of the boat decreased the probability of reflected light from contaminating the spectrometer readings.Furthermore, being close to the water allowed measurements to be made from 20 cm above the surface of the water, qualifying the measurement as an 'at surface' reading.Every station had duplicate measurements and each duplicate had 15 samples for water, sky and standard plaque.Averaging of the 15 samples improved the signal-noise ratio.Upwelling radiance was converted into normalized water leaving reflectance with Equation (2), where,   ,   ,   are signals for water, sky, and standard plaque respectively,   is the reflectance rate for the standard plaque on whole bands (%, which is provided by factory calibration), and  is the sea-air interface reflectance rate (%).As a response to wind speed, r varies from 2.5%-2.7%(see details in Tang, et.al [27]).Solar irradiance has significant impacts on water upwelling radiance, while lesser influences on water reflectance (see Figures 4 and 5 in Mobley [28]).Since reflectance takes into account environmental factors from solar irradiance after normalization and all measurements were conducted at the same angles, they can be directly compared with each other.Water reflectance then went through a baseline correction to shift the infrared band to 0 (Rrs(λ) − Rrs(763)), and was smoothed with a five-point median filter.We also simultaneous collected sky iPAR (instant photosynthetically available radiation) from a cosine receptor setup on an ASD HH2 at every station.

Apparent Optical Properties
We applied (45  ) angles for above-water measurements recommended by NASA Ocean color protocol [26], to normalize water radiance from simultaneous sky radiance and to remove interference from environmental light conditions.A hand-held VNIR spectroradiometer (HH2, ASD Inc., Boulder, Colorado, USA), with a high spectral resolution of 1 nm, was used to measure water leaving radiance.When measuring, the Zodiac stopped its engine to prevent white capping.However, without power, the boat did not drift far from its original location.The black hull of the boat decreased the probability of reflected light from contaminating the spectrometer readings.Furthermore, being close to the water allowed measurements to be made from 20 cm above the surface of the water, qualifying the measurement as an 'at surface' reading.Every station had duplicate measurements and each duplicate had 15 samples for water, sky and standard plaque.Averaging of the 15 samples improved the signal-noise ratio.Upwelling radiance was converted into normalized water leaving reflectance with Equation (2), where, S water , S sky , S p are signals for water, sky, and standard plaque respectively, ρ p is the reflectance rate for the standard plaque on whole bands (%, which is provided by factory calibration), and r is the sea-air interface reflectance rate (%).As a response to wind speed, r varies from 2.5%-2.7%(see details in Tang, et al. [27]).Solar irradiance has significant impacts on water upwelling radiance, while lesser influences on water reflectance (see Figures 4 and 5 in Mobley [28]).Since reflectance takes into account environmental factors from solar irradiance after normalization and all measurements were conducted at the same angles, they can be directly compared with each other.Water reflectance then went through a baseline correction to shift the infrared band to 0 (Rrs(λ) − Rrs(763)), and was smoothed with a five-point median filter.We also simultaneous collected sky iPAR (instant photosynthetically available radiation) from a cosine receptor setup on an ASD HH2 at every station.

In Situ Chlorophyll-a Concentration
Duplicate surface water samples (500 ml for each) were collected at each station (Figure 1b).Water samples were filtered under 50 kpa onto a GF/F filter.The filter was immersed into a flask with 10 ml 90% acetone and wrapped with aluminum foil.All steps were conducted under low light conditions to prevent chlorophyll-a decomposition.Extracted chlorophyll-a was stored for 24 h in a freezer.The flask was then placed in a centrifuge for 10 min at 4000 revolutions/min (TDL-60B, Anke Ins., Ninbo, China).Total chlorophyll-a was measured from its supernatant with three readings using a fluorometer (AquaFluor®Handheld Fluorometer and Turbidimeter, Turner Design, San Jose, CA, USA), with an excitation band of 430 nm and an emission band of 660 nm.Pheophytin concentration was measured again with the fluorometer after the supernatant had a chemical reaction with 10% HCl for 1 min.The chlorophyll-a concentration responsible for photosynthesis is total chlorophyll-a concentration minus the pheophytin concentration after converting the fluorescence intensity to the chlorophyll-a value from a chlorophyll-a: fluorescence curve.All duplicates show a mean deviation lower than 5.3%.Prior to chlorophyll-a evaluation, a spectrophotometer (HITACHI F-2700, Hitachi High-Technologies Corporation, Tokyo, Japan) was used to calibrate the various chlorophyll-a concentration values on the fluorometer using standard stock chlorophyll-a diluted gradient solutions.

In Vivo Fluorescence
A fluorometer (AquaFluor®Handheld Fluorometer and Turbidimeter, Turner Design, San Jose, CA, USA.) was also used to measure in vivo fluorescence at all stations.Surface water was directly placed in the measuring window to retain phytoplankton light acclimation and physiology information.All samples were measured 3 times.The in vivo fluorometer has an excitation band at 430 nm and emission band at 660 nm.This instrument directly transfers water fluorescence intensity into chlorophyll-a concentration from a preset linear equation based on factory calibration.

Satellite/In Situ Match-Ups
We limited satellite/in situ match-ups with three criteria.First, we used a mean 3 × 3 spatial windows on the satellite image.Second, the gap between in situ and satellite image was limited to less than four hours.Lastly, an AOT (aerosol optical thickness) index lower than 0.15 was chosen.All satellite data were L2 products from the Ocean Color Website collected by MODIS Aqua and Terra.Before converting to water leaving reflectance, the L2 products were subject to pixel-by-pixel atmospheric correction with the 6S model and kept at a spatial resolution of 1.1 km without any resampling.Matchups between in situ chlorophyll-a concentration and water optical properties (Rrs) were processed with global algorithms (Table 1) to estimate chlorophyll-a concentration accuracy and error around Ardley and Great Wall Cove.
To obtain regional coefficients for the nFLH algorithm, we applied the 'Leave One Out' cross-validation method [29].For each iteration, a matching pair was left out and residual error and curve fitting coefficients were estimated.The 27 coefficients and residual errors were then averaged to obtain the unbiased estimation.
where, RMSE is the residual error, nFLH in situ is the fluorescence value from the field and nFLH est is the value estimated from the 26 matching pairs.Equation ( 4) is the fitting function and c 1 , c 2 , and c 3 are the fitting coefficients.The leave-one out approach helps to limit the problem of overfitting [29].

Phytoplankton Absorption Coefficent Spectrum
We applied bio-optical models to determine phytoplankton absorption coefficients from water leaving reflectance.There are two commonly used approaches in bio-optics, QAA [20] and GSM01 [19].GSM01 retrieves simultaneous estimates for chlorophyll-a concentration, the absorption coefficients for dissolved, and detrital materials and particulate backscatter.Model parameters are then tuned through simulated annealing.The GSM01 model decreases residual errors from satellite and estimated Rrs through a multiple iteration approach (Maritorena,et al. [19]).The QAA (Lee, et al. [19]) analytically calculates coefficient values of total absorption and backscattering from remote sensing reflectance.In comparison to GSM01, the QAA approach has the benefit of not requiring any prior information about the spectral shape of a Φ (λ) and thereby reduces potential errors and uncertainties with spectral models or inappropriate spectral shapes [19].
Considering the uniqueness of SO water optics [13,14], we intended to keep more original Rrs signals in phytoplankton absorption coefficients.Therefore, we applied QAA approach to invert the in situ Rrs and expanded the single band absorption into whole visible bands by introducing some of GSM01 equations.The QAA approach we applied here followed its V5 coefficients [30] and pure water absorption [31].Details are shown in Appendix A.
Ciotti et al. [32] developed a micro-pico cell-size composition estimation approach from absorption coefficient spectrum.They used least-square fitting to gain abundance (S f ) of micro and pico cells using whole bands.
where a ph (λ) are phytoplankton absorption coefficients, S f is abundance of pico-phytoplankton, a pico (λ) and a micro (λ) are standard pico-and micro-absorption coefficients per unit from Table 3 in Ciotti, et al. [32].This algorithm was used to estimate micro-phytoplankton fractions.* nlw is water leaving radiance; WAP(West Antarctic Peninsula)

Spatial Distribution of Water Optical Properties
Chlorophyll-a varied significantly between Ardley Cove and Great Wall Cove (Figure 1c).Ardley Cove had a large phytoplankton bloom, with a maximum of chlorophyll-a at 9.87 mg/m 3 .Relatively low chlorophyll-a concentrations occurred in the Great Wall Cove, reaching only 1.37 mg/m 3 .Simultaneous in vivo fluorescence did not follow the same pattern, with the maximum occurring on the east side of Ardley Island instead of the north side (Figure 1d).
Remote Sens. 2017, 9, 210 6 of 15 KGI Rrs spectrum produced 4 different shapes, with varied slopes between 550-560 nm (Figure 2a).These shapes were driven by chlorophyll-a concentration, with chlorophyll-a concentration increasing the overall Rrs value across all wavelengths less than 600 nm.For red bands higher than 600 nm, 4 Rrs spectra decreased significantly due to pure water absorption [32].
Despite the different Rrs spectral shapes for the various chlorophyll-a concentrations, the phytoplankton absorption coefficients showed similar curves after removing non-phytoplankton absorptions and particle backscattering coefficients (Figure 2b).Two obvious peaks appear in the blue and red bands of phytoplankton absorption spectra.Higher chlorophyll-a concentration has higher absorption coefficients, and the 660 nm absorption band has a good agreement with chlorophyll-a concentration (52.6%, p < 0.05) (Figure 2c).
Remote Sens. 2017, 9, 210 6 of 15 increasing the overall Rrs value across all wavelengths less than 600 nm.For red bands higher than 600 nm, 4 Rrs spectra decreased significantly due to pure water absorption [32].Despite the different Rrs spectral shapes for the various chlorophyll-a concentrations, the phytoplankton absorption coefficients showed similar curves after removing non-phytoplankton absorptions and particle backscattering coefficients (Figure 2b).Two obvious peaks appear in the blue and red bands of phytoplankton absorption spectra.Higher chlorophyll-a concentration has higher absorption coefficients, and the 660 nm absorption band has a good agreement with chlorophyll-a concentration (52.6%, p < 0.05) (Figure 2c).

Algorithim Chlorophyll-a Estimation and In Situ/Satellite Match-Ups
Chlorophyll-a estimation from various algorithms were compared to in situ Rrs and residual errors were calculated.The results showed poor relationships with almost all previously developed algorithms (Table 1).All correlations were lower than 30%.An exception was the algorithm of Mitchell et.al [33].Relative errors from estimated and in situ chlorophyll-a were greater than 60%.Although all these algorithms were built from the SO dataset, scarce samples in high latitude and high chlorophyll-a coastal water resulted in poor performance of chlorophyll-a estimation in KGI waters.
In addition, satellite-derived chlorophyll-a and fluorescence intensity from global empirical algorithms were evaluated.The comparison between in situ and satellite data showed a relatively good correlation with fluorescence intensity (55.35%) (Figure 3a).Direct band-ratio chlorophyll-a estimation showed a non-linear relationship with in situ chlorophyll-a with an 11.43% correlation (Figure 3b).Poor estimation was also apparent in the QAA bio-optical approach at a correlation of 1.62% (Figure 3b).

Algorithim Chlorophyll-a Estimation and In Situ/Satellite Match-Ups
Chlorophyll-a estimation from various algorithms were compared to in situ Rrs and residual errors were calculated.The results showed poor relationships with almost all previously developed algorithms (Table 1).All correlations were lower than 30%.An exception was the algorithm of Mitchell et.al [33].Relative errors from estimated and in situ chlorophyll-a were greater than 60%.Although all these algorithms were built from the SO dataset, scarce samples in high latitude and high chlorophyll-a coastal water resulted in poor performance of chlorophyll-a estimation in KGI waters.
In addition, satellite-derived chlorophyll-a and fluorescence intensity from global empirical algorithms were evaluated.The comparison between in situ and satellite data showed a relatively good correlation with fluorescence intensity (55.35%) (Figure 3a).Direct band-ratio chlorophyll-a estimation showed a non-linear relationship with in situ chlorophyll-a with an 11.43% correlation (Figure 3b).Poor estimation was also apparent in the QAA bio-optical approach at a correlation of 1.62% (Figure 3b).

Fluorescence Approach Estimation for Chlorophyll-a
Based on the relatively good linear relationship between in situ and satellite fluorescence measurements, we explored the feasibility of using the fluorescence approach for estimating chlorophyll-a in KGI.This issue of non-linearity between in vivo fluorescence and chlorophyll-a, caused by NPQ under intense solar irradiance, was considered.Our measurements also showed a decreasing exponential phase of fluorescence per chlorophyll-a: iPAR (Figure 4a).Increasing solar irradiance decreases the fluorescence of chlorophyll-a yield, illustrating the impact of NPQ.
Previous research [23] reported a turning point of NPQ when iPAR reaches 100 μ mol/sec.After that point, NPQ will grow slowly.Since our measurements were collected in the southern summer with long daytime and high solar irradiance, the in vivo fluorescence influence on chlorophyll-a decreased (Figure 4a).Their correlation was 40.64%.High chlorophyll-a concentration samples had larger underestimation and deviation from 1:1 line than low chlorophyll-a concentration measurements (Figure 4b).
For the relationship between global empirical fluorescence from satellite and in situ chlorophylla (Figure 4c), the global coefficients do not work well in estimating chlorophyll-a in Ardley Cove since the fluorescence reacts a lot on phytoplankton physiology and iron stress [34].Therefore, we optimized the three coefficients in the nFLH algorithm from the 30 in vivo fluorescence samples for the KGI region to produce the following algorithm: where at 678, 667, and 748 nm, respectively.This optimization improved fluorescence estimation by about 10% from global empirical coefficients according to Figure 3a.Then, NPQ (FLH/iPAR) and phytoplankton absorption package correction (FLH/  ℎ ̅̅̅̅̅, where  ℎ ̅̅̅̅̅ is absorption mean between 300-700 nm per chlorophyll-a unit following Babin et al. [35]) were applied to the fluorescence estimation.The final fluorescence: chlorophyll-a relationship increased from original 40.64% to 59.46% (Figure 4d).Those two corrections referenced Behrenfeld et al. [23] correction in their Figure 2b,d.We do not introduce a global satellite average for iPAR and  ℎ to non-dimensionlization as they The satellite fluorescence algorithm here is based on the global fluorescence coefficients from NASA ocean color group (https://oceancolor.gsfc.nasa.gov/atbd/nflh/).(b) Correlation between satellite and in situ chlorophyll-a concentrations.The satellite algorithm applied here is the OC3 global band-ratio (http://oceancolor.gsfc.nasa.gov/cms/atbd/chlor_a) and NASA IOP model derived absorption coefficient for phytoplankton (http://oceancolor.gsfc.nasa.gov/cms/atbd/giop).

Fluorescence Approach Estimation for Chlorophyll-a
Based on the relatively good linear relationship between in situ and satellite fluorescence measurements, we explored the feasibility of using the fluorescence approach for estimating chlorophyll-a in KGI.This issue of non-linearity between in vivo fluorescence and chlorophyll-a, caused by NPQ under intense solar irradiance, was considered.Our measurements also showed a decreasing exponential phase of fluorescence per chlorophyll-a: iPAR (Figure 4a).Increasing solar irradiance decreases the fluorescence of chlorophyll-a yield, illustrating the impact of NPQ.
Previous research [23] reported a turning point of NPQ when iPAR reaches 100 µ mol/sec.After that point, NPQ will grow slowly.Since our measurements were collected in the southern summer with long daytime and high solar irradiance, the in vivo fluorescence influence on chlorophyll-a decreased (Figure 4a).Their correlation was 40.64%.High chlorophyll-a concentration samples had larger underestimation and deviation from 1:1 line than low chlorophyll-a concentration measurements (Figure 4b).
For the relationship between global empirical fluorescence from satellite and in situ chlorophyll-a (Figure 4c), the global coefficients do not work well in estimating chlorophyll-a in Ardley Cove since the fluorescence reacts a lot on phytoplankton physiology and iron stress [34].Therefore, we optimized the three coefficients in the nFLH algorithm from the 30 in vivo fluorescence samples for the KGI region to produce the following algorithm: nFLH = where at 678, 667, and 748 nm, respectively.This optimization improved fluorescence estimation by about 10% from global empirical coefficients according to Figure 3a.Then, NPQ (FLH/iPAR) and phytoplankton absorption package correction (FLH/a ph , where a ph is absorption mean between 300-700 nm per chlorophyll-a unit following Babin et al. [35]) were applied to the fluorescence estimation.The final fluorescence: chlorophyll-a relationship increased from original 40.64% to 59.46% (Figure 4d).Those two corrections referenced Behrenfeld et al. [23] correction in their Figure 2b,d.We do not introduce a global satellite average for iPAR and a ph to non-dimensionlization as they did in their correction.In our equation, the magnitude of in situ iPAR (10 3 ) and a ph (10 −3 ) balanced out and produced agreeable results.

Nutrient, Light, and Phytoplankton Composition
Behrenfeld et al's.[23] algorithm was useful in providing phytoplankton physiological status for our samples.According to Behrenfeld et al. [23], the remaining fluorescence signal divided by chlorophyll-a after NPQ and phytoplankton absorption package corrections are defined as fluorescence quantum yield and reflects the factors limiting phytoplankton growth (e.g., iron, light, or nitrogen/phosphorus). NPQ correction removes the photo-protection impacts on fluorescence intensity, and phytoplankton absorption package correction decreases the gaps between various phytoplankton communities.Fluorescence quantum yield is then quantified to global phytoplankton growth limitation factors where > 1.4% stands for iron limitation and < 1.4% means light, nitrogen/phosphorus limitation (in Figure 4 from their publication).The quantum yield in our samples was < 0.7% indicating that light and nitrogen/phosphorus limitation is greater in KGI Great Wall Cove than Ardley Cove (Figure 5a).
The spatial pattern of iPAR (Figure 5b) is similar to fluorescence quantum yield and chlorophyll-a in Figure 1c., which showed a stronger positive link between solar irradiance and phytoplankton biomass.

Nutrient, Light, and Phytoplankton Composition
Behrenfeld et al's.[23] algorithm was useful in providing phytoplankton physiological status for our samples.According to Behrenfeld et al. [23], the remaining fluorescence signal divided by chlorophyll-a after NPQ and phytoplankton absorption package corrections are defined as fluorescence quantum yield and reflects the factors limiting phytoplankton growth (e.g., iron, light, or nitrogen/phosphorus). NPQ correction removes the photo-protection impacts on fluorescence intensity, and phytoplankton absorption package correction decreases the gaps between various phytoplankton communities.Fluorescence quantum yield is then quantified to global phytoplankton growth limitation factors where > 1.4% stands for iron limitation and < 1.4% means light, nitrogen/phosphorus limitation (in Figure 4 from their publication).The quantum yield in our samples was < 0.7% indicating that light and nitrogen/phosphorus limitation is greater in KGI Great Wall Cove than Ardley Cove (Figure 5a).after 2002 (Figure 6b,c).When we focus on the year with extended sea ice coverage and high sea ice concentration (2009, 2011, and 2015), fluorescence intensity also increased in the following southern summer (2010, 2012, and 2016).Those years with high phytoplankton concentrations included 2010 [7].It is noted that the fluorescence data shown here did not go through NPQ and phytoplankton absorption package corrections, so Figure 6b illustrates the relative trend with in the most recent decade.

Fluorescence Approach Estimation for Chlorophyll-a
Light energy is an important factor regulating phytoplankton growth in KGI though photoinhibition.The spatial pattern of IPAR (Figure 5b), which showed good agreement with chlorophylla concentration (Figure 1c), clearly illustrated the influence of light limitation on phytoplankton growth.This physical condition prevents optimal light for phytoplankton growth and this along with intense wind and terrigenous particles, typically leads to low primary productivity KGI [25].From the perspective of this study, low light conditions ensured the reliability of the fluorescence approach and the nFLH algorithm.
Studies have shown improved application of the baseline fluorescence approach in case II waters over the global blue-green band-ratio algorithms [36].The fluorescence approach also has higher effectiveness than the semi empirical algorithms like QAA (this study).Studies have also shown that inefficiencies in the SO bio-optics algorithms is generated from green band reflectance, leading to an The spatial pattern of iPAR (Figure 5b) is similar to fluorescence quantum yield and chlorophyll-a in Figure 1c., which showed a stronger positive link between solar irradiance and phytoplankton biomass.
Phytoplankton size, deduced from the phytoplankton absorption coefficients, showed that the dominant phytoplankton in the bloom area are micro-size algae (> 20µm) (Figure 6a).The abundance of micro-size algae occupied 79.58% of the entire bloom phytoplankton community in Ardley Cove during our study.The mean residual error for Sf (pico-phytoplankton fraction estimation) is 16.54%.Phytoplankton size, deduced from the phytoplankton absorption coefficients, showed that the dominant phytoplankton in the bloom area are micro-size algae (> 20μm) (Figure 6a).The abundance of micro-size algae occupied 79.58% of the entire bloom phytoplankton community in Ardley Cove during our study.The mean residual error for Sf (pico-phytoplankton fraction estimation) is 16.54%.
Additionally, satellite-estimated chlorophyll-a and sea ice in this area show an increasing trend after 2002 (Figure 6b,c).When we focus on the year with extended sea ice coverage and high sea ice concentration (2009, 2011, and 2015), fluorescence intensity also increased in the following southern summer (2010, 2012, and 2016).Those years with high phytoplankton concentrations included 2010 [7].It is noted that the fluorescence data shown here did not go through NPQ and phytoplankton absorption package corrections, so Figure 6b illustrates the relative trend with in the most recent decade.

Fluorescence Approach Estimation for Chlorophyll-a
Light energy is an important factor regulating phytoplankton growth in KGI though photoinhibition.The spatial pattern of IPAR (Figure 5b), which showed good agreement with chlorophylla concentration (Figure 1c), clearly illustrated the influence of light limitation on phytoplankton growth.This physical condition prevents optimal light for phytoplankton growth and this along with intense wind and terrigenous particles, typically leads to low primary productivity KGI [25].From the perspective of this study, low light conditions ensured the reliability of the fluorescence approach and the nFLH algorithm.
Studies have shown improved application of the baseline fluorescence approach in case II waters over the global blue-green band-ratio algorithms [36].The fluorescence approach also has higher effectiveness than the semi empirical algorithms like QAA (this study).Studies have also shown that inefficiencies in the SO bio-optics algorithms is generated from green band reflectance, leading to an Additionally, satellite-estimated chlorophyll-a and sea ice in this area show an increasing trend after 2002 (Figure 6b,c).When we focus on the year with extended sea ice coverage and high sea ice concentration (2009, 2011, and 2015), fluorescence intensity also increased in the following southern summer (2010, 2012, and 2016).Those years with high phytoplankton concentrations included 2010 [7].It is noted that the fluorescence data shown here did not go through NPQ and phytoplankton absorption package corrections, so Figure 6b illustrates the relative trend with in the most recent decade.

Fluorescence Approach Estimation for Chlorophyll-a
Light energy is an important factor regulating phytoplankton growth in KGI though photo-inhibition.The spatial pattern of IPAR (Figure 5b), which showed good agreement with chlorophyll-a concentration (Figure 1c), clearly illustrated the influence of light limitation on phytoplankton growth.This physical condition prevents optimal light for phytoplankton growth and this along with intense wind and terrigenous particles, typically leads to low primary productivity KGI [25].From the perspective of this study, low light conditions ensured the reliability of the fluorescence approach and the nFLH algorithm.
Studies have shown improved application of the baseline fluorescence approach in case II waters over the global blue-green band-ratio algorithms [36].The fluorescence approach also has higher effectiveness than the semi empirical algorithms like QAA (this study).Studies have also shown that inefficiencies in the SO bio-optics algorithms is generated from green band reflectance, leading to an obvious underestimation in chlorophyll-a when using global empirical coefficients [9].In addition, the blue band is highly affected by CDOM, whose absorption coefficient was reported to occupy 70% of the non-water components in WAP (west of the Antarctica Peninsula) [37].The absorption spectrum of CDOM has decreasing exponential influence, showing less influence on longer bands.Since the red band has longer wavelengths in the visible band, the selection of red bands will introduce less contamination in CDOM rich waters.
Conversely, this band selection has some issues.The red/NIR bands are more susceptible to atmospheric interference.Two recent programs, SOCTRATES [16] and WCRP [38] focused on the SO atmosphere, intend to improve the understanding of clouds, aerosols, air-sea exchanges over the SO.The improvement of SO atmospheric model will further meet the demands of input parameters for ocean color estimation in the near future.
SST-based corrections of fluorescence intensity have also been applied to SO waters.Browning et al. [39] found that extremely low SST in the SO has a strong connection with NPQ parameters and can be used to correct field chlorophyll-a fluorescence signals.They showed highly variable phytoplankton physiology features and communities under regional irradiance conditions and suggested more regional studies on NPQ capacity impacts.However, this correction needs an empirical relationship between SST and physiology coefficient B and was not applied in this study.
The byproduct of fluorescence quantum yield from the fluorescence approach also has a significant role in estimating primary productivity.State-of-the-art algorithms for primary productivity estimation are mainly based on chlorophyll-a absorption-fluorescence quantum yield equation or a photosynthesis-irradiance equation [40].If we improve the estimation of fluorescence quantum yield from space, we will improve estimation for SO phytoplankton photosynthesis and carbon export.This will lead to a clearer understanding of the role of the SO in global climate change.

Factors on Phytoplankton Bloom in KGI
As an island hosting several international research stations, the ecological trend of the coastal waters surrounding KGI deserves attention.Historical records in KGI have shown that no phytoplankton blooms [6] occurred in these waters until 2010, during which maximum chlorophyll-a reached 20 mg/m 3 .Research has assumed that the increasing phytoplankton growth was exacerbated by global warming [25].This manuscript documented another phytoplankton bloom in 2016 and attempted to determine its potential cause from optics and satellite history records.Iron limitation was widely observed in SO phytoplankton growth, which leads to ferredoxin deficiency in photosystem I and decreases photosynthetic efficiency.Iron depletion is reported as a major restriction for SO primary productivity [41,42].However, previous research in KGI found no obvious nutrient or iron deficiencies in phytoplankton growth [43].The microphytoplankton dominant features in this bloom (Figure 6a) further confirms iron enrichment because micro-size cells favor the iron abundant water.Micro-phytoplankton maintains a high growth rate in iron rich water, and their large size prevents zooplankton grazing [44].
Diatoms found around the WAP are larger in size than those found in other regions [44,45].The phytoplankton bloom reported in 2010 was dominated by micro-algae (e.g., Porosira glacialis, Thalassiosira antarctica, and T. ritscheri [25]).Phytoplankton tends to adapt their cellular physiology to optimize light harvesting and photoreception [44].Generally, small size algae benefit from light limited conditions due to their high absorption-photosynthesis effectiveness and small cellular shading [46].Micro-size algae have more self-protection under high light conditions because their large surface area will cause low absorption-photosynthesis effectiveness.Therefore, for the micro-phytoplankton dominant bloom in KGI, nutrient support and zooplankton grazing has more influence on phytoplankton community composition than light limitation.
Owing to extended sea ice coverage and high sea ice concentration (Figure 6b,c), KGI has more fresh water input following southern summer to form strong stratification and keep phytoplankton in surface waters.This shallow mixed layer increases the probability for phytoplankton to absorb light energy.In addition, sea ice melt brings nutrients for algal growth [46][47][48][49].The increasing trend of sea ice increases nutrients input (e.g., iron or micronutrients), which will lead to phytoplankton blooms in spring and summer.In addition, phytoplankton tends to bloom in marginal sea ice regions in the coastal ocean of the SO [46][47][48][49].Therefore, increasing sea ice extent and coverage in KGI is a likely trigger for increased phytoplankton blooms.

Errors for Fluorescence Approach Estimating Chlorophyll-a
Two factors are responsible for accurate satellite estimation of chlorophyll-a.These include the sea-air model for calculating chlorophyll-a from water leaving Rrs and the atmospheric radiation transfer model to transfer top-of-atmosphere radiance into water leaving Rrs.Atmospheric correction has a large probability of introducing errors, because its signal occupies 90% of the total signal received at the satellite [50].Therefore, atmospheric correction deserves great attention, particularly in improving its accuracy.Until now, there has not been a proper high latitude atmospheric correction model for the southern hemisphere and high cloud coverage in the SO creates a challenge for accurate atmospheric correction [16,17,[51][52][53].
Atmospheric correction requires several steps to remove the influence of the atmosphere, such as single scattering, multiple molecular scattering, aerosol/ozone, and other gas absorption [54].The procedure usually removes cloud contaminated pixels by the threshold of the near infrared band albedo.However, this method does not work well on cloud edges.Current satellite algorithms provide estimates of cloud optical thickness [55], which can be used to quantify the impacts of clouds on ocean color data in the near future.
The ocean color signal reflected from water is generally the integrated light reflected from the light penetration depth [56,57].With increasing depth, deep waters have less contribution to the water modulated reflectance signals.However, we still cannot ignore it in chlorophyll-a inversion algorithms.In the SO, the deep chlorophyll-a maximum increases the contribution from deeper waters.Research seldom applies chlorophyll-a profiles and integrated chlorophyll-a signals to compare in situ and satellite match-ups.Current reports have already shown the obvious underestimation from satellite algorithms [6][7][8][9][10][11].If considering the deep chlorophyll-a maximum and chlorophyll-a profile, there is a larger underestimation in SO chlorophyll-a retrieval.Chlorophyll-a profile estimation should be considered in future research.In fact, due to wind mixing, a deep chlorophyll-a maximum is a common situation in the SO [58,59].Sullivan et al. [60] built a model on polar phytoplankton growth and reduction.They found that there is a high probability of maintaining a subsurface chlorophyll-a maximum during the post-bloom period.

Conclusions
Decades of development in satellite ocean color has produced multiple algorithms to detect, among other things, chlorophyll-a, fluorescence intensity, particulate organic carbon, photosynthetic available radiation.Satellites provide almost every parameter needed for observation and understanding the global ocean.However, the uncertainty of those parameters varies between regions.Our research validated the fluorescence approach to estimate chlorophyll-a concentration in KGI, resulting in a good agreement (59.46%) at this high latitude, high chlorophyll-a region.The phytoplankton bloom in 2016 showed that the dominant phytoplankton is micro-phytoplankton, occupying 79.58% of the total phytoplankton community in Ardley Cove.Increasing sea ice coverage and sea ice concentration are possible reasons for the increasing occurrence phytoplankton blooms in the recent decade.Due to NPQ, the fluorescence approach does not accurately estimate chlorophyll-a in places that have intense solar irradiance.However, the fluorescence approach and red band selection has notable advantages in avoiding CDOM interference from blue bands and decreasing gaps from the peculiar bio-optics of SO green bands.Our future work intends to validate and extend the application of this algorithm to the entire SO.The photosynthesis mechanism revealed by the fluorescence approach will provide more information on SO primary productivity estimation and its role on global carbon cycling.

Figure 1 .
Figure 1.(a) Sampling area location.(b) Locations of in situ samples.(c) Chlorophyll-a spatial pattern interpolated from the in situ samples during January 2016.(d) In vivo fluorescence interpolated from the in situ samples during January 2016.

Figure 1 .
Figure 1.(a) Sampling area location.(b) Locations of in situ samples.(c) Chlorophyll-a spatial pattern interpolated from the in situ samples during January 2016.(d) In vivo fluorescence interpolated from the in situ samples during January 2016.

Figure 2 .
Figure 2. (a) Four selected above water Rrs spectra.(b) Four featured phytoplankton absorption coefficients derived from Figure 2a using QAA v5 presented in Appendix.A. (c) Correlation between phytoplankton absorption coefficients and in situ chlorophyll-a concentration.

Figure 2 .
Figure 2. (a) Four selected above water Rrs spectra.(b) Four featured phytoplankton absorption coefficients derived from Figure 2a using QAA v5 presented in Appendix A. (c) Correlation between phytoplankton absorption coefficients and in situ chlorophyll-a concentration.

Figure 4 .
Figure 4. (a) Relationship between fluorescence per chlorophyll-a unit and iPAR.(b) Relationship between in situ chlorophyll-a and fluorescence intensity.(c) Relationship between in situ chlorophylla and in situ Rrs estimated fluorescence.(d) Relationship between in situ chlorophyll-a and in situ fluorescence estimated chlorophyll-a after regional coefficients improvement, NPQ correction and phytoplankton absorption package correction.

Figure 4 .
Figure 4. (a) Relationship between fluorescence per chlorophyll-a unit and iPAR.(b) Relationship between in situ chlorophyll-a and fluorescence intensity.(c) Relationship between in situ chlorophyll-a and in situ Rrs estimated fluorescence.(d) Relationship between in situ chlorophyll-a and in situ fluorescence estimated chlorophyll-a after regional coefficients improvement, NPQ correction and phytoplankton absorption package correction.

Figure 5 .Figure 6 .
Figure 5. (a) Fluorescence quantum yield interpolated from in situ samples during January 2016, and (b) iPAR spatial pattern interpolated from in situ samples during January 2016

Figure 5 .
Figure 5. (a) Fluorescence quantum yield interpolated from in situ samples during January 2016, and (b) iPAR spatial pattern interpolated from in situ samples during January 2016.

Figure 5 .Figure 6 .
Figure 5. (a) Fluorescence quantum yield interpolated from in situ samples during January 2016, and (b) iPAR spatial pattern interpolated from in situ samples during January 2016

Figure 6 .
Figure 6.(a) Pico-phytoplankton fraction spatial pattern retrieved from the in situ spectrum algorithm (Equation (5)) during January 2016 (b) Hovmoller longitude-average plots of satellite records on fluorescence estimated chlorophyll-a from Giovanni (http://giovanni.gsfc.nasa.gov/giovanni/).(c) Hovmoller longitude-average plots of satellite records on sea ice concentration from Giovanni.

Table 1 .
Various models application on chlorophyll-a estimations in KGI phytoplankton bloom area.