Remote Sensing and Argo Float Observations Reveal Physical Processes Initiating a Winter-Spring Phytoplankton Bloom South of the Kuroshio Current Near Shikoku

: BIO-Argo ﬂoat (chlorophyll a (Chl-a), temperature, and salinity proﬁles) and remote sensing data (Chl-a, photosynthetic available radiation (PAR), and wind) located south of the Kuroshio current near Shikoku from September 2018 to May 2019 were used to study phytoplankton bloom and their mechanisms of development in open oceans. Results show that higher (lower) Chl-a concentrations are correlated with a deeper (shallower) mixed layer (R Pearson = 0.77, R crit = 0.12 (alpha = 0.05, n = 263)) compared to the average of Chl-a and mixed layer depth (0.13 mg / m 3 and 105 m). The average net accumulation rates (r) of phytoplankton were close to 0.08 d − 1 . An increasing r corresponds to a gradually increasing surface Chl-a (S (Chl-a): 0–20 m average Chl-a) and integrated Chl-a inventory (I (Chl-a): integrated Chl-a from surface to euphotic depth). These phenomena indicate that the mechanism of winter-spring phytoplankton blooms is consistent with the dilution-recoupling hypotheses (DRH). During the bloom formation, winter deep mixing and eddy-wind Ekman pumping are enhanced by a strong winter monsoon. The enhancement may disturb predator–prey interactions and dilute zooplankton in deep mixed layers. Moreover, winter deep mixing and eddy-wind Ekman pumping can cause the nutrients to be transported into the euphotic layer, which can promote the growth of phytoplankton and increase grazing. During the bloom extinction, the stratiﬁcation strengthens and the intensity of light increases; this increases grazing and nutrient consumption, and decreases the phytoplankton bloom signiﬁcantly (S (Chl-a) and I (Chl-a) increase by 0.3 mg / m 3 and 27 mg / m 2 , respectively). The output from a biogeochemistry model shows that nutrients are consistent with the temporal distribution of S (Chl-a) and I (Chl-a). Our results suggest that physical processes (deep winter mixing and eddy-wind Ekman pumping) under the DHR framework are critical factors for winter-spring blooms in open oceans with an anticyclone eddy. methodology, T.W.; software, T.W.; validation, T.W., F.C., W.Z. and H.N.; formal analysis, T.W., F.C., and S.Z.; investigation, T.W., F.C., W.Z. and H.N.; resources, T.W. and F.C.; data curation, T.W. and S.Z.; writing—original draft preparation, T.W.; writing—review and editing, S.Z., F.C., J.P. and A.T.D.; visualization, T.W., H.N. and W.Z.; supervision, S.Z. and J.P.; project administration, S.Z. and J.P.; funding acquisition, S.Z. and J.P.


Introduction
Chlorophyll a (Chl-a) concentration is affected by the changes in intracellular pigment levels driven by light and nutrients. Recently, satellite-estimated phytoplankton carbon has been used as an

BIO-Argo Float
BIO-Argo float (2902750) data observed south of the Kuroshio current near Shikoku (28.6 • N-31.5 • N, 132.7 • E-135.9 • E) from the China Array for Real-time Geostrophic Oceanography (ARGO) project were used. The study area is located in the subtropical region far from land (Figure 1a), and can be influenced by the Kuroshio recirculation [17]. The BIO-Argo float observations used in this study were taken from the period of 9 September 2018 16:45:00 to 31 May 2019 23:33:12, a duration of 265 d, using SBE (Sea-Bird Electronics), CTD (Conductivity, Temperature and Depth), and wet labs multi-parameter (measuring Chl-a fluorescence and BBP700 backscattering coefficient) sensor data and sampled profiles of Chl-a concentrations, taken in trajectories of "ring" shapes ( Figure 1b), trapped in an anticyclone eddy (Supplementary Materials Figure S1). A complete vertical profile was sampled once per day before 20 March 2019 and every two days afterwards. Chl-a concentrations were determined following the method of Welschmeyer [18]: Chl-a = (Fluorescence Chl-a _ − Dark Chl_a ) * Scale Chl-a ; Scale Chl-a = 0.0073; Dark Chl-a = 48 (1)

Satellite Observations
We used a daily multisensory composite (surface Chl-a concentration (mg/m 3 ), photosynthetic available radiation (PAR) (I 0 (einstein/m 2 /day)), and diffused attenuation coefficient at 490 nm of the downwelling irradiance (K 490 (m −1 )) in GLOBAL 4 km pixel size resolution. Multisensory imagery merging was undertaken following three merging techniques (simple averaging, weighted averaging, GSM (Garver-Siegel-Maritorena) mode) [19,20], downloaded from https://resources.marine.copernicus.eu/ and http://hermes.acri.fr/. The relative performance of merging techniques depends on the conditions (water types, region, and glint/aerosol conditions). The reprocessing version was NASA R2018.0 for MODIS Aqua, VIIRS NPP, and VIIRS JPSS-1 sensors; and ESA PB 2.16 to 2.55 for the OLCI-A sensor. The Chl-a multi-sensor product aims to merge chlorophyll-a fields that have been rebuilt with a similar approach for each sensor. For oligotrophic, mesotrophic, and coastal water, the product relies on algorithms of Hu et al. [21], OC3-OC4 [22,23], and OC5 [24], respectively. A multi-sensor product is obtained by combining the processed data from each sensor. K 490 and I 0 rely on the Morel algorithm [25] and Frouin et al. algorithm [26], respectively. The number of match-ups with in situ data is significant (R Pearson = 0.84, Rcrit = 0.01 (alpha = 0.05, n = 27,917)) [27].
We used the following equation to convert attenuation coefficient at 490 nm wavelength (K 490 ) to the attenuation coefficient at the I 0 wavelength (K PAR ): Euphotic depth (Z eu (m)), defined here as the depth of 1% of surface irradiance, was estimated from PAR profiles assuming Lambert-Beer's relationship [28]: The daily 10 m wind speed (m/s) data were obtained from ftp://ftp.remss.com/ccmp/v02.0/, with a spatial resolution of 25 km. The daily surface heat flux (W/m 2 ) was obtained from https://www.esrl.noaa.gov/, with a spatial resolution of 200 km. Net surface heat flux was calculated as the sum of net longwave radiation, net shortwave radiation, sensible heat flux, and latent heat flux [29].

Model Design and Numerical Products
The Global Analysis Forecast Bio 001_028 product is based on NEMO version (v3.6_STABLE), provided by Copernicus Marine Environment Monitoring Service (CMEMS) by the Global Analysis Forecast Phys 001_024 product at daily frequency. Initial conditions were taken from the World Ocean Atlas 2018 (https://www.nodc.noaa.gov/OC5/woa18/) for nitrate (N (mmol/m 3 )), phosphate (P (mmol/m 3 )), silicate (S (mmol/m 3 )) and dissolved oxygen (DO (mmol/m 3 )). The Operational Mercator Ocean biogeochemical global ocean analysis and forecast system can provide up to ten days of 3D global ocean forecasts, with a spatial resolution of 25 km, updated weekly. The time series are aggregated to obtain a two-year time series sliding window. This product includes daily and monthly means of biogeochemical parameters (Chl-a, N, P, S, DO, dissolved iron, primary production, phytoplankton, pH, and surface partial pressure of carbon dioxide) over the global oceans. The spatial and temporal resolutions of outputs are 25 km and daily, respectively. Twenty vertical levels were used, ranging from 0 to 400 m.

Mixed Layer Depth
Mixed layer depths (Z mixed (m)) were calculated using temperature and salinity data from BIO-Argo floats profiles. Z mixed were based on two criteria: the temperature change is 0.2 • C relative to the value at 10 m depth, and the density change is 0.03 kg/m 3 relative to the value at 10 m depth [30,31]. For each profile, the mean of the two values was chosen for Z mixed .

Turbulence Regimes, Mixing Depth and Mixing Time Scales
We calculated and distinguished turbulence regimes, mixing depth, and mixing time scales according to the method mentioned by Rumyantseva et al. [31], and turbulence of convective boundary layers under the forcing of wind stress, τ, can be quantified by the Monin-Obukhov length scale [32] given by: where u * = τ ρ 0 is the friction velocity, ρ 0 is the seawater reference density, k = 0.4 is von Karman's constant, and B 0 = Qλg C p ρ 0 is the surface buoyancy flux, where λ is the thermal expansion coefficient, Q is the surface heat flux, C p is the heat capacity, and g is the gravity parameter [32].
The effects of convection and wind stress on mixing of the ocean surface layer can be classified into two different regimes [33]: Regime 1: Wind mixing: Regime 2: Convective mixing: Z mixed > C 1 Z MO and Q < 0.
For Regime 1, Q < 0, Z mixed can be deeper than active mixing depth (Z mixing ) [34]: where C 2 = 2 and f is the Coriolis parameter [31]. When it is in surface warming (Q > 0): where C 3 = 1 and C 4 = 0.57 are determined by the output of the large-eddy simulation (LES) model [35]. For Regime 2: Z mixing = Z mixed (7) The time scale of vertical mixing affected by wind and convection can be estimated based on Deardorff [36], respectively: where C 5 = 1. C i (i = 1,2,3,4,5) represents scaling constants for Z mixed and Z mixing , adopted from Thorpe [33], Enriquez and Taylor [35], and Deardorff [36]. The choice of constants (C i ) in these equations does not affect the main conclusions [31].

Phytoplankton Specific Growth Rate
Generally, the growth rate of phytoplankton is mainly controlled by the nutrient supply, light conditions, and temperature [37]. The Kuroshio, a western boundary current, is known as a nutrient-poor warm salty northward stream [38]. For the Kuroshio current region south of Shikoku, the nitrate concentration is smaller than 1 mmol/m 3 within 100 m of the upper ocean ( Figure S4), which can be regarded as an oligotrophic region, as can be seen in WOA2018. The maximum specific growth rate can be calculated with following equation: where u max is the maximum phytoplankton growth rate and T is the averaged mixed layer temperature [37]. The phytoplankton growth rate [39] as a function of light conditions is determined by the following equations: where α chl is the Chl-a specific slope of the phytoplankton-irradiance curve, α chl = 11 (mol/m 2 ) −1 , and θ = 0.03 (g chl/g C) is the cellular chlorophyll-to-carbon mass ratio [31,39]. The phytoplankton growth time scale (t g ) is used to compare with the mixing time scale:

Net Accumulation Rate of Phytoplankton
We obtained the phytoplankton net accumulation rate (r) according to the method described by Behrenfeld [11]: (15) where S (Chl-a) is the 0-20 m average Chl-a, and I (Chl-a) is the integrated Chl-a from the surface to euphotic depth.
The critical phytoplankton depth (Z CR ) [9,40] is derived from PAR (I 0 ), K PAR , and the compensation irradiance (I C ), as follows:

Distribution of Chl-a Concentration
The obtained values of surface Chl-a based on BIO-Argo correlated well (R Pearson = 0.65, R crit = 0.12 (alpha = 0.05, n = 263)) with the satellite-derived daily Chl-a over the sampling site ( Figure 2a). However, the values were somewhat different in magnitude. The BIO-Argo observations may be overestimated and may represent the same change trend with satellite-derived daily Chl-a. Chl-a was further corrected based on match-up satellite-derived daily Chl-a using an exponential function (Chl-a = 0.2534 * Chl-a argo 0.0049 ) [41]; the result (R Pearson = 0.81, R crit = 0.12 (alpha = 0.05, n = 263)) is shown in Figure 2b. In the following content, our research is mainly based on BIO-Argo float observations (temperature, salinity, corrected Chl-a). Chl-a concentrations' seasonal cycles of variable amplitude were found; in Figure 3a, Chl-a peaks typically range from 0.15 to 0.35 mg/m 3 in winter-spring. This order of magnitude was similar to that of Rumyantseva et al. [31]. When the Chl-a concentration first reaches 0.15 mg/m 3 , the bloom starts each year and can last for up to 24 days. The variations of Chl-a south of the Kuroshio current to the south of Shikoku (Figures 1 and 3), at a water depth of 800 m, clearly display seasonal cycles of blooms (>0.15 mg/m 3 ) from September 2018 to May 2019, with high Chl-a in winter and spring. The peak (i.e., blooming in winter and spring) of Chl-a (Figure 3c) typically appears in December-May. There were seasonal variations evident in the magnitude of Chl-a concentration. The phytoplankton bloom (Chl-a > 0.15 mg/m 3 ) first appeared at the end of December; that is, one month after the initial increase in surface Chl-a (caused by the erosion of the deep chlorophyll maximum) and the bloom decayed in April-May. During the period before the coming of winter, S (Chl-a) and I (Chl-a) were approximately 0.05 mg/m 3 and 9 mg/m 2 , respectively (Figure 3c), and then increased gradually to 0.3 mg/m 3 and 27 mg/m 2 , respectively, with the mixed layer deepening in winter-spring. Compared with the average of Chl-a and Z mixed (0.13 mg/m 3 and 105 m), higher (lower) Chl-a generally correspond to a deeper (shallower) mixed layer. Notably, Chl-a in May was lower than before the winter-spring blooms.

Distribution of Ocean Conditions during Winter-Spring Blooms
The mixed layer gradually deepened from the end of September until mid-February, and then became deeper than Z eu and remained consistent, until a rapid transition to the shallow seasonal stratification was observed at the end of March (Figure 3a). During this blooming period, a lower SST (Sea Surface Temperature) coincided with the bloom, Z eu , and PAR reduced slightly, and the surface net heat flux was mainly negative (Figure 3b,e,f). Surface wind stress was stronger in winter and early spring, and generally positive trends of S (Chl-a) and I (Chl-a) were present at that time, which is consistent with gradually increasing surface heat flux and PAR (Figure 3c,e,f) [31].

The Comparison of Different Spatial Scale Data
Satellite observations (Chl-a, PAR, K 490 ), wind components, heat flux, and model outputs with different spatial scales were extracted for the grid point centered on the location of the BIO-Argo float, the closest pixel to the sampling site. This approach has been applied in the work of Rumyantseva et al. [31], and may also introduce some errors and uncertainties. However, this study investigates relative changes (seasonal variation) in Chl-a, rather than its absolute magnitude. To further evaluate the error, we calculated the standard deviation and average of all grid points in the study area, and compared it with the closest pixel to the sampling site ( Figures S2 and S3). It was found that there is a good correlation between the closest pixel to the sampling site value and the average value of this area, and the corresponding standard deviation is small. The closest pixel to the sampling site value can roughly indicate surface Chl-a (R Pearson = 0.91, R crit = 0.12 (alpha = 0.05, n = 263)), wind speed (R Pearson = 0.44, R crit = 0.12 (alpha = 0.05, n = 263)), PAR (R Pearson = 0.71, R crit = 0.12 (alpha = 0.05, n = 263)), and heat flux (R Pearson = 0.98, R crit = 0.12 (alpha = 0.05, n = 263)) at the sampling site ( Figure S2). The vertical profile of the model outputs (Chl-a and nitrate) at the sampling site is shown in Figure S3; there was little difference between the vertical distribution of the closest pixel to the sampling site and the average of this area.
Some parameters (t m,wind , euphotic depth, Monin-Obukhov length scale, critical phytoplankton depth and Z mixing ) were derived from remote sensing data. Because calculation of other parameters (t m,convection , Z mixed , maximum specific growth rate, phytoplankton growth time scales, and phytoplankton net accumulation rate) require subsurface data, these parameters were derived from BIO-Argo float observations (temperature, salinity, Chl-a). The work of Rumyantseva et al. [31] also used the same calculation method and similar data to investigate the impact of atmospheric forcing and light conditions on phytoplankton blooms in the temperate North Atlantic.

The Exclusion of the CDH and CTH for Bloom Initiation
As can be seen in Figure 3, the phytoplankton bloom occurred in weakly-stratified conditions before spring, and compared with the average of Chl-a and Z mixed (0.13 mg/m 3 and 105 m), higher (lower) Chl-a generally corresponded to a deeper (shallower) mixed layer, which demonstrated a good correlation (R Pearson = 0.77, R crit = 0.12 (alpha = 0.05, n = 263)). This suggests that physical processes may control phytoplankton blooms in this area.
Here, three hypotheses (i.e., CDH, CTH, and DRH) were used to analyze the mechanism of the blooms. According to Sverdrup's theory [9], the standard for the beginning of a spring bloom is that the depth of the mixed layer becomes shallower than the critical depth. At this critical depth, the integrated production is equal to the integrated respiration [31]. As surface PAR gradually increases (decreases) (Figure 3e), the critical depth becomes deep (shallow) enough, even though low PAR and shallow critical depth may result in a net accumulation of phytoplankton for several months before seasonal re-stratification. As shown in Figure 4a, the estimated ranges of Z CR and Z mixing were compared according to Sverdrup's theory of critical depth; the result showed Z CR ≤ Z mixing from 23 December to 8 March. During this period, in which the phytoplankton bloom appeared, the corresponding r was mostly positive (Figure 4b). At the arrival of spring, Z CR > Z mixing , the bloom began to decay. Our results suggest that the period of the decreasing cumulative sum of r corresponded to a shallower mixed layer and improving light conditions for phytoplankton during spring. However, the critical depth criterion did not match well with the timing of the bloom, and the cumulative sum of r showed a rapid decay (Figure 4c). Thus, the CDH hypothesis can be ruled out for phytoplankton blooms in our study.
For the CTH hypothesis, the starting time of the bloom is defined as the time when phytoplankton growth time scales are shorter than mixing time scales [10]. We thereby compared the estimated mixing (t m ) and growth (t g ) time scales. As shown in Figure 4e, the growth time scales were approximately two orders of magnitude greater than turbulent mixing time scales for wind and convective mixing regimes, suggesting that the critical turbulence criterion was insufficient to support the bloom initiation. In particular, the growth and mixing time scales have a large difference, namely, the difference is as high as one order of magnitude. Therefore, the CTH hypothesis does not appear valid for the bloom in our study.
The BIO-Argo observations also do not support the hypothesis as predicted by CDH [9], that bloom initiation can be explained by the depth of the mixed layer (Z mixing ) being shallower than the critical depth as light conditions improve (Figure 4a). According to the estimated mixing and growth time scales, t m are approximately two orders of magnitude smaller than t g (Figure 4e), therefore, CTH also does not support the bloom initiation.

The Winter-Spring Bloom Development Revealed by DRH
Because CDH and CTH cannot account for the phytoplankton blooms in our study, we now consider DRH. According to DRH, zooplankton grazing can exert a huge impact on the regulation of the seasonal phytoplankton cycle [8]. DRH regarding phytoplankton blooms focuses on physical and ecological processes influencing the balance between phytoplankton division and loss rates, in particular the losses to grazing [11]. The blooms are initiated by physical processes (deep winter mixing, freshwater input, upwelling, or polar night) that disrupt the balance between phytoplankton division and grazer consumption [12]. Some studies have found that the basis of net phytoplankton population growth in winter is diluted by winter convective mixing in deep mixed layers because intense mixing can reduce encounter rates of phytoplankton with grazers, thus resulting in a positive net accumulation rate of phytoplankton [11,13,42]. Unfortunately, there are no observational data or literature about zooplankton quantity and composition in the study region. Investigations of the role of such events in the long-term changes of this region's zooplankton quantity and composition should be the subject of future studies.
In our study, more positive r was detected in winter, associated with the development of upper-ocean stratification when µ mean of phytoplankton was low due to light availability (Figure 4). Of the three hypotheses examined here, only DRH can account for the changes observed in Chl-a (Figure 4d). By early winter, phytoplankton cells were homogeneously dispersed by vertical mixing, and the impact from grazers exceeds light-driven decreases in cell division, thus helping to induce phytoplankton bloom [11,42,43]. From the end of November to the beginning of March, the convective mixing regime was the dominant process, inducing the deepening of the mixed layer (Figure 4a). During this period, mean net accumulation rates were close to 0.08 d −1 , and the cumulative sum of r continued to increase (Figure 4b,c), corresponding to a gradually increasing S (Chl-a) and I (Chl-a), even though the µ mean of phytoplankton was low (Figure 4d). The strength of positive net accumulation rates were mostly less than negative net accumulation rates after the beginning of March, and corresponded to a decreasing cumulative sum of r, gradually improving light conditions and an increasing µ mean (Figure 4c,d), but a gradually decreasing S (Chl-a) and I (Chl-a). Most of the time (90%) convective regimes dominated (Figure 4e). With an increase in winter-convective mixing, the bloom was gradually aggravated even though light conditions did not begin to improve (before the beginning of March). The subsequent bloom period lasted for three months, and intense wind stress and positive eddy-wind Ekman pumping velocity events became more frequent, further contributing to deepening the mixed layer and convective mixing (Figures 3d and 5) to sustain the bloom. With the onset of spring, the mixed layer stopped deepening. Thus, the spring bloom slightly increased before the end of March and lasted for only 20 d. However, a significant decrease (0.2 mg/m 3 and 18 mg/m 2 ) in S (Chl-a) and I (Chl-a) was observed from the end of March to the end of May. In the meantime, r < 0 was detected, seasonal stratification developed, and the accumulation of biomass decreased. Behrenfeld and Boss [12] also found that the phytoplankton in the mixed layer increased, but concentrations remained low due to the effect of dilution by convective mixing. Therefore, it is clear evidence that the winter bloom initiation occurred due to the dilution effect, which could sustain the spring bloom (from March to April). However, the improvement of light conditions did not result in a more pronounced phytoplankton bloom, as discussed below. As can be seen from Figure 3a,d, and Figure 4d, the calculated average growth rate is generally low for the wind and convective mixing regimes, but corresponds to an increasing phytoplankton bloom, and shows that sustained winter convective mixing and prevailing winter monsoon (more strong wind stress) can induce decreasing grazer pressure and vertical nutrient transport. The "eddy pumping" effect can result in downwelling in an anti-cyclone eddy associated with nutrient-depleted and reduced Chl-a, whereas eddy Ekman pumping can inhibit this process through eddy-wind interaction [44][45][46]. The positive Ekman pumping velocity resulting from the interaction of the northward monsoon wind and anti-cyclone eddy further contributed to providing nutrients for phytoplankton (Figure 5a,b). The spring bloom only lasted for one month despite improving light conditions until nutrients were depleted, and then began to decay. A comparison of the decay initiation of S (Chl-a) and I (Chl-a) indicated that S (Chl-a) before I (Chl-a) was closely related to vertical nutrient distribution. In their examination of the DRH framework, Behrenfeld and Boss [12] suggested that spring blooms were initiated by physical processes that disrupted the balance between phytoplankton division and grazer consumption, such as deep winter mixing and upwelling. In the recoupling stage, the change of phytoplankton population was determined by the balance between the specific growth rate of light-dependent phytoplankton and the loss caused by grazing.
The outputs from the biogeochemistry model are shown in Figure 6 in winter and spring, which assimilate the satellite observation data and reveal that the temporal evolution of Chl-a was consistent with the observed result of the BIO-Argo float. Therefore, we infer that the r calculated based on the model results and BIO-Argo float observations were also consistent, and it can be said that the model results also provide support for the dilution-recoupling hypothesis. In addition, based on the consistency of Chl-a concentrations, the corresponding modeled nutrient changes can further explain the phytoplankton bloom process. With deep winter-convective mixing, nitrate concentration was gradually improving before light conditions began to improve (before the beginning of March) in the upper ocean, which was consistent with the winter phytoplankton bloom. The remaining nutrients in the mixed layer can allow the spring phytoplankton bloom to last for one month, without the need for intense convective mixing. This was true even though light conditions improved after March, and S (Chl-a) and I (Chl-a) continued to decrease owing to increased grazing and nutrient depletion in the stably-stratified upper ocean though phytoplankton photosynthesis. Overall, physical processes (deep winter mixing and eddy-wind Ekman pumping) were the critical factors for the sustaining of the winter-spring bloom considered here for the open ocean with an anticyclone eddy. However, these influences will decrease under the background of global warming, which appears likely to cause more stable water stratification (and thereby relatively mild winter deep mixing and eddy-wind Ekman pumping) in the future. This phenomenon will reduce the available nutrients in seawater and decrease the primary productivity, which in turn may affect the whole ecosystem. Therefore, future investigations of the long-term changes in phytoplankton dynamics are critical.

Conclusions
We studied the winter-spring bloom initiation and development using biogeochemical observations from BIO-Argo floats (2902750) south of the Kuroshio current near Japan and tested three commonly discussed hypotheses (CDH, CTH, or DRH) to determine which can best explain the observed results.
According to Sverdrup' s theory of critical depth and growth (mixing) time scales, we ruled out two hypotheses (CDH, CTH). Our results confirmed that the DRH hypothesis can best explain the blooms, which were initiated when winter-convective mixing increased, and were associated with dilution effects due to winter deep mixing, which can disturb predator-prey interactions associated with the seasonal stratification in the upper ocean. The development of the bloom from winter to spring was markedly influenced by sustained winter convective mixing and meteorological conditions (winter monsoon) through their effect on nutrients necessary for phytoplankton. More frequent strong wind stress associated with upwelling induced by eddy-wind Ekman pumping can contribute to decreasing grazing pressure and injecting nutrients in the euphotic layer. A biogeochemistry model product suggested that the development of bloom was consistent with the observed result of BIO-Argo floats, and that nutrient availability was a key factor for the development of the winter-spring bloom. This model supported the conclusions of the earlier work by Behrenfeld and Boss [12]. The bloom began to decay when increased grazing and nutrient depletion accompanied by strong stratification occurred after April, even under enhanced light conditions. In addition, physical conditions within the anticyclone eddy would have included significant downwelling, which would not have been representative of open waters. At present, we are using the existing data to conduct a preliminary exploration of individual cases. The low temporal and spatial resolution of buoy data limits the further exploration of the universality of DHR theory, and more Argo floats and high-resolution model data may be needed to increase knowledge about the temporal variability of biogeochemical processes. Investigations of the role of such events on the long-term changes of this region's phytoplankton dynamics should be the subject of future studies.