Comparison of Primary Production Using in situ and Satellite-Derived Values at the SEATS Station in the South China Sea

Satellite-based observations of primary production (PP) are broadly used to assess carbon fixation rate of phytoplankton in the global ocean with small spatiotemporal limitations. However, the remote sensing can only reach the ocean surface, the assumption of a PP vertically exponential decrease with increasing depth from the surface to the bottom of euphotic zone may cause a substantial and potential discrepancy between in situ measurements and satellite-based observations of PP. This study compared euphotic zone integrated PP derived from measurements based on ship-based in situ incubation (i.e., PPin situ) and those derived from the satellite-based vertically generalized production model (VGPM; PPVGPM) for the period 2003∼2016 at the South East Asian Time-series Study (SEATS) station. PP values obtained during the NE-monsoon (NEM: Nov∼Mar; PPin situ = 323 ± 134; PPVGPM = 443 ± 142 mg-C m–2 d–1) were ∼2-fold higher than those recorded during the SW-monsoon (SWM: Apr∼Oct; PPin situ = 159 ± 58; PPVGPM = 250 ± 36 mg-C m–2 d–1), regardless of the method used for derivation. The main reason for the higher PP values during the NEM appears to have been a greater abundance of inorganic nutrients were made available by vertical advection. Note that on average, PPin situ estimates were ∼50% lower than PPVGPM estimates, regardless of the monsoon. These discrepancies can be mainly attributed to differences from the euphotic zone depth between satellite-based and in situ measurements. The significantly negative relationship between PP measurements obtained in situ and sea surface temperatures observed throughout this study demonstrates that both methods are effective indicators in estimating PP. Overall, our PPin situ analysis indicates that a warming climate is unfavorable for primary production in low-latitude open ocean ecosystems.


INTRODUCTION
Primary production at the bottom of the marine food web plays a key role in the ocean ecosystem (Pauly and Christensen, 1995) and represents a major pathway for sequestration and/or cycling of atmospheric CO 2 by the oceans. However, this process is highly susceptible to environmental and climatic changes (Buitenhuis et al., 2013;Hung et al., 2013Hung et al., , 2016Liu et al., 2021;Zhong et al., 2021). It is important for oceanographers to gain a comprehensive understanding of the spatiotemporal characteristics of primary production (Field et al., 1998;Campbell et al., 2002;Tang et al., 2008;Tilstone et al., 2015).
The quantification of oceanic primary production is generally based on in situ measurements pertaining to incubation wherein the euphotic zone integrated primary production (PP) is calculated via trapezoidal integration (i.e., the PP inventory is evaluated by integrating PP divided into small trapezoids from the surface to the bottom of the euphotic zone). Measurements of PP by the radiolabeled carbon (C-14) uptake method has been extensively used in marine environments since the introduction of this method to determine the carbon uptake rate of phytoplankton (Nielsen, 1952;Hama et al., 1983;Gong, 1993;Shiah et al., 2000). Despite considerable research in the western North Pacific (WNP), i.e., East China Sea, South China Sea (SCS), and Taiwan Strait, to establish phytoplankton carbon fixation rates, researchers continue to debate whether phytoplankton growth conditions at the surface are indicative of conditions in deeper waters. Obtaining measurements of temperature at depth is straightforward; however, the in situ observation of PP (PP in situ ) requires considerable effort in terms of manpower and time. Specifically, seawater samples collected at discrete depths must be incubated on deck within an incubator using running surface water for cooling (Shiah et al., , 2003(Shiah et al., , 2005Chen, 2005;Lai et al., 2014;Chen et al., 2016). Note also that gaps in PP in situ coverage inevitably lead to discrepancies in corresponding estimates. The first empirical algorithm for PP predictions based on remote sensing was proposed by Balch et al. (1989). Behrenfeld and Falkowski (1997) developed an attractive alternative approach to estimating global PP using a small number of inputs. Their vertically generalized production model (VGPM; PP VGPM ) is widely regarded as the most highly optimized yet usable methods for PP estimation (Kameda and Ishizaka, 2005;Yamada et al., 2005;Ishizaka et al., 2007;Hill and Zimmerman, 2010).
The semi-analytical VGPM uses satellite-based data as an input to calculate PP. Thus, it is conceivable that data measured in situ could be used as an alternative to satellite-based input data, and vice versa (Behrenfeld and Falkowski, 1997;Hill and Zimmerman, 2010). The VGPM can be used to derive PP using satellite-based observations, such as remote passive ocean color, sea surface temperature (SST), surface optimal carbon fixation rate (P B opt ) per unit chlorophyll a (Chl) of the euphotic zone (Z eu ), surface Chl concentration (Chl s ), surface light intensity (E 0 ), and the surface light diffuse attenuation coefficient (K d ) (Yamada et al., 2005;Ishizaka et al., 2007). Nonetheless, it is still problematic whether the assumption based on peak values of PP and Chl s occur at the surface and decrease exponentially with depth until reaching the bottom of the Z eu is true (Hill and Zimmerman, 2010;Buitenhuis et al., 2013). Furthermore, the reliability of VGPM in evaluating PP diminishes when geographic features, vertical hydrographic distributions, regional characteristics, extreme weather events, and climate change are taken into consideration (Dierssen, 2010;Friedland et al., 2012;Hung et al., 2013Hung et al., , 2016Chen et al., 2015;Shih et al., 2015Shih et al., , 2020b. In the absence of a robust estimation method, the results derived from PP VGPM cannot be relied upon to reflect the actual situation throughout the oceans. Empirical models developed by Dunne et al. (2005); Laws et al. (2011), andHenson et al. (2011) are widely used by oceanographers to estimate particulate organic carbon (POC) export flux; however, reliance on PP VGPM also calls into question all corresponding estimates pertaining to global carbon export flux and oceanic carbon sequestration.
The South East Asian Time-series Study (SEATS) conducted in the South China Sea (SCS) was the lowest latitude time-series program implemented during the Joint Global Ocean Flux Study era (Karl et al., 2003;Wong et al., 2007). Numerous studies have characterized the upper ocean of the SCS as stratified and oligotrophic (Chen et al., 2004;Chen, 2005;Wong et al., 2007). Thus, biological activity and regulation of biogeochemical responses depend heavily on dynamic perturbations, including the yearly monsoon, typhoons, storms, internal waves, Kuroshio intrusion, and atmospheric deposition as well as nutrient supply in the form of phytoplankton nitrogen fixation (Liu et al., 2002;Chou et al., 2006;Chen et al., 2008Chen et al., , 2020Du et al., 2013;Yang et al., 2014;Li et al., 2018;Shih et al., 2020a).
Throughout the SCS basin, the annual modeled PP ranges from 280 to 343 and PP VGPM ranges from 308 to 354 mg-C m −2 d −1 (Table 1; Liu et al., 2002;Tan and Shi, 2009;Ma et al., 2014). High PP values are generally associated with the strong NE-monsoon (NEM) system during the cold season, whereas low PP values are associated with a relative weak SWmonsoon (SWM) during the warm season (Table 1). Based on long-term satellite-based SST records, Chen et al. (2020) reported that declining PP can be attributed at least in part to rising SST (0.012 • C y −1 ). Their and several selected researches in our study went a long way toward establishing a connection between monsoons and PP; however, all of the suppositions are based on indirect measurements; i.e., remote sensing (Liu et al., 2002;Hao et al., 2007;Zhao et al., 2008;Tan and Shi, 2009;Pan et al., 2012;Ma et al., 2014;Chen et al., 2020), rather than direct in situ incubation, such as the C 13 and C 14 methods (Liu et al., 2002;Chen et al., 2004Chen et al., , 2007Ning et al., 2004). Note also that even PP research based on in situ oceanographic analysis is limited to short-term observations rather than long-term timeseries (Chen, 2005;Chen et al., 1998Liu et al., 2002;Tseng et al., 2005). Remote sensing has made it possible to conduct large-scale PP monitoring over extended durations; however, most existing sensing technology is limited to the upper ocean and there is little evidence supporting the assumption of an exponential decrease in PP as a function of depth (Buitenhuis et al., 2013;Shih et al., 2020b). In the current study, we employed data obtained in the SEATS time-series study to elucidate long-term variations in PP in situ under the prevailing monsoon system. We also looked  (Chen et al., 1998). for discrepancies between satellite-based and in situ observations, which could potentially influence PP estimates derived using the VGPM algorithm.  (Karl et al., 2003;Wong et al., 2007). On-deck incubation of PP in accordance with 14 C protocols was also conducted during 20 of 25 expeditions between 2003 and 2016 (Figure 1). Conductivity-temperature-depth (CTD) (SBE9/11, SeaBird) and quantum scalar irradiance (QSP-200L, Biospherical) were, respectively, used to obtain vertical profiles of temperature and underwater photosynthetically available radiation (PAR). A Biospherical instrument (QSR-240) was used to obtain daily PAR measurements (i.e., daily surface light intensity). In that study, Z eu was defined as the depth at which underwater PAR reached 1%. The mixed layer depth (MLD) was defined as the deepest depth at which the temperature was 0.8 • C lower than at the surface (Kara et al., 2000;Chou et al., 2006). Seawater samples were collected at discrete depths to determine Chl concentrations and PP values. Briefly, 2L seawater samples were filtered using 47-mm GF/F filters before obtaining Chl concentrations (via acidification) using a Turner 10-AU-005 fluorometer Shiah et al., 2003). PP was measured using the 14 C assimilation method (Parsons et al., 1984) following incubation under an artificial light source (2,000 µE m −2 s −1 with full spectrum from 350 to 2450 nm, similar to sunlight) for ∼3 h in a proprietary isolation container using flowing surface seawater for cooling. Following incubation and acidification (0.5 N HCl), radioactive substances collected in 25-mm GF/F filters were transported to a lab for quantification using a scintillation counter (Packard 2200) Shiah et al., 2003;Lai et al., 2014). Finally, PP was calculated via trapezoidal integration (i.e. integrated from surface to 1% PAR). Note that SST, Chl and PP measurements from the Ocean Data Bank (Ministry of Science and Technology, Taiwan 1 ) were used to compensate for deficiencies in in situ observations. Some of the hydrographic records (i.e., SST) and biogeochemical parameters (i.e., PP, Chl s ) of May-October from 2003 to 2006 have been present by Pan et al. (2012). 1 http://www.odb.ntu.edu.tw/

Remote Sensing Observation
Between 2003 and 2016, the daily data of level-3 products were obtained using a passive ocean color MODerate resolution Imaging Spectroradiometer (Aqua sensor; MODIS-Aqua) with spatial resolution of 4 km for the region covering 17.5-18.5 • N and 115.5-116.5 • E. We also collected PP VGPM , SST, E 0 , K d , and Chl s values from the Environmental Research Division's Data Access Program (ERDDAP), National Oceanic and Atmospheric Administration, Department of Commerce, U. S. 2 . Note that 8-day composite data were applied in situations where daily observations were missing. Note also that PP VGPM values obtained from ERDDAP were estimated using VGPM developed by Behrenfeld and Falkowski (1997).
Plausible deviations between PP VGPM and PP in situ were examined by comparing data obtained from satellites vs. data obtained in situ, including i.e., SST sate , (P B opt ) sate , E 0−sate , Note that this analysis was conducted using the VGPM algorithm proposed by Behrenfeld and Falkowski (1997) (with far fewer input variables), as Equation 2-1: where the D irr refers to the photoperiod. Z eu was estimated from the average K d of the water column (Behrenfeld and Falkowski, 1997 and references therein). P B opt was computed using the empirical equation proposed by Behrenfeld and Falkowski (1997), as Equation 2-2, where T refers to the incubation temperature measured in the deck incubator as SST in situ or SST sate . In cases where SST exceeded 28.5 • C, a constant value of 4.0 mg-C mg-Chl −1 h −1 was taken as P B opt (Behrenfeld and Falkowski, 1997;Hu et al., 2014).
To evaluate the impact of monsoonal force on time-series variations, all data sets were divided into two groups: SWmonsoon (SWM, Apr. to Oct., including inter-monsoons of Apr and Oct) and NE-monsoon (NEM, Nov. to Mar). Variations between satellite-based observations and in situ measurements were examined by averaging all relevant values and reporting them as mean ± standard deviation (SD). Linear regression was used to assess relationships between any two of the variables. The t-tests were used to compare sets of observations with the significance level set at p = 0.10.

RESULTS
Time-Series Distributions: PP, SST, E 0 , K d , and Chl s As shown in Figure 2A, PP in situ was ∼2.0 times higher during the NEM (117-528, average = 323 ± 134 mg-C m −2 d − ) than during the SWM (76-247, average = 159 ± 58 mg-C m −2 d −1 ) ( Table 2). PP VGPM values between 2003 and 2016 were as follows: SWM was 250 ± 36; 154-771 mg-C m −2 d −1 and NEM was 443 ± 142; 190-1,153 mg-C m −2 d −1 (Figure 2A). Note that PP VGPM during the NEM was ∼1.8 times higher than during the SWM (p < 0.01) ( Table 2). The differences in PP in situ and PP VGPM between the NEM and SWM are in line with those reported by Ning et al. (2004); Chen (2005), and Hao et al. (2007), all of which were obtained from the same area of the SCS using different methods ( Table 1). PP VGPM exhibited monsoonal variations resembling the curve derived from PP in situ ; however, the magnitudes were 37% higher during the NEM and 57% higher during the SWM. On average, PP in situ estimates were ∼50% lower than PP VGPM estimates, regardless of monsoon.
The lowest SST in situ (25.0 ± 1.4 • C) and SST sate (25.1 ± 1.2 • C) values were obtained under the prevailing NEM during the cold season. The highest SST in situ (28.9 ± 0.9 • C) and SST sate (28.8 ± 1.0 • C) values were obtained under the prevailing SWM during the warm season ( Table 2). The low SST values recorded during the NEM can be attributed to the combined dynamics of NEM-driven vertical mixing of seawater and reduced solar radiation (Tseng et al., 2005;Zhou et al., 2020). By contrast, the P B opt showed a synchronous but opposite distribution (Figures 2B,C). The (P B opt ) in situ and (P B opt ) sate values computed from SST for the NEM, respectively, ranged from 4.4 to 6.1 and 4.0 to 6.6 mg-C mg-Chl −1 h −1 , and the average value was the same (5.4 ± 0.5 mg-C mg-Chl −1 h −1 ). The (P B opt ) in situ and (P B opt ) sate values for the SWM, respectively, ranged from 4.0 to 4.4 and 4.0 to 6.0 mg-C mg-Chl −1 h −1 , with mean P B opt values of 4.1 ± 0.2 and 4.2 ± 0.2 mg-C mg-Chl −1 h −1 ( Table 2). E 0−in situ and E 0−sate values in the SWM (50 ± 11 and 47 ± 6 Eins m −2 d −1 ) were higher than those in the NEM (33 ± 14 and 36 ± 8 Eins m −2 d −1 ) (see Figures 2D,E and Table 2). From the perspective of variables in the VGPM algorithm, E 0 /(E 0 +4.1) ratios varied little between the in situ measurement and satellite-based observation. The variable of [E 0 /(E 0 +4.1)] thereby seemed not an important parameter to influence the PP level. Z eu was estimated using the expression of -ln(0.01)/K d (Kirk, 1994;Behrenfeld and Falkowski, 1997), which led to a reversal of monsoonal K d and Z eu distributions (see Figures 2F,G). The mean Z eu−in situ was deeper during the SWM than during the NEM (87 ± 10 and 75 ± 14 m, respectively). The Z eu−sate as well as the Z eu−in situ presented similar trends between SWM and NEM (156 ± 17 and 120 ± 30 m, respectively) ( Table 2). Z eu values were shallower in the NEM than in the SWM, due perhaps to reduced light penetration resulting from more abundant biomass in the water column (Chen, 2005;Tseng et al., 2005;Chen et al., 2008). Figure 2H presents long-term temporal variations in Chl s . In situ measurements and satellite-based observations revealed a distinct NEM maximum (0.24 ± 0.16 and 0.20 ± 0.08 mg m −3 , respectively) and a SWM minimum (0.08 ± 0.02 and 0.10 ± 0.02 mg m −3 , respectively). Chl s−in situ and Chl s−sate concentrations were, respectively, 0.07 to 0.58 and 0.09 to 0.71 mg m −3 during the NEM, and were, respectively, 0.04 to 0.11 and 0.05 to 0.43 mg m −3 during the SWM. For the same location, our observations were close to the data reported by Chen (2005) and Tseng et al. (2005). The SWM minimum Chl s values in the TABLE 2 | Summary of euphotic zone PP (PP) and its relevant variables in the current study conducted at the SEATS: SST, P B opt , E 0 , E 0 /(E 0 +4.1), K d , Z eu , Chl s and Z ML .  current study were nearly the same as those reported in the oligotrophic ocean time-series studies at HOT (Hawaii Ocean Time-Series) and BATS (Bermuda Atlantic Time-Series Study) during the summer (0.05 mg m −3 ); however, NEM values at SEATS exceeded the winter maximum at HOT and BATS (0.1 and 0.3 mg m −3 , respectively) (Karl et al., 2003). The high NEM maximum at SEATS may perhaps be explained by an increase in phytoplankton biomass (particularly larger sizes of >3 µm), which was stimulated by the deepening Z ML ( Table 2). The deeper nutrient-rich water was then efficiently transported to the upper surface water under the influence of NEM given that the nutrient-cline depth was shallower in SCS than in other oceans (Chen, 2005;Tseng et al., 2005).
Monthly Variations in PP, SST, E 0 , K d and Chl s The PP VGPM values in this study are in good agreement with PP in situ in terms of amplitude as well as phase. Overall, we observed a maximum PP in situ of 394 ± 190 in January (NEM) and a minimum of 143 ± 70 mg-C m −2 d −1 in August (SWM) ( Figure 3A). The magnitude of PP VGPM was higher than the values obtained using in situ measurements; however, the trend over a span of 12 months was similar, with a pronounced NEM peak 619 ± 113 mg-C m −2 d −1 in January and the lowest SWM value of 230 ± 28 mg-C m −2 d −1 in September. These results are similar to those reported by other researchers for same area using different methods, such as the particulate organic carbon flux re-calculation (NEM: 207, SWM: 149 mg-C m −2 d −1 ) (Chen et al., 1998); on-deck C 14 incubation (NEM: 300-509, SWM: 110-228 mg-C m −2 d −1 ) (Ning et al., 2004;Tseng et al., 2005); on-deck C 13 incubation (NEM: 190-550, SWM: 190-280 mg-C m −2 d −1 ) (Chen, 2005;Chen et al., 2007), and Chl empirical function (NEM: 148-684, SWM: 86-275 mg-C m −2 d −1 ) ; Table 1).
We observed opposing trends between monthly SST and the corresponding P B opt . The highest (P B opt ) in situ values (5.9 mg-C mg-Chl −1 h −1 ) were observed in January and the lowest (4.0 mg-C mg-Chl −1 h −1 ) were observed in August. The (P B opt ) sate values converted from SST sate (5.8 and 4.0 mg-C mg-Chl −1 h −1 in January and June, respectively) were nearly identical to those for the same time slots obtained via in situ measurements (Figures 3B,C). The highest mean monthly E 0−in situ in April (63 ± 0 Eins m −2 d −1 ) was higher than the highest mean monthly E 0−sate (53 ± 3 Eins m −2 d −1 ). Overall, E 0−in situ and E 0−sate levels in December were very similar (∼26 and ∼27 Eins m −2 d −1 , respectively). Based on E 0−in situ , the highest E 0 /(E 0 +4.1) in situ ratio (0.94) occurred in April and the lowest E 0 /(E 0 +4.1) in situ ratio (0.85) occurred in December. E 0 /(E 0 +4.1) sate values of 0.93 in April and 0.86 in December were similar with the E 0 /(E 0 +4.1) in situ values (Figures 3D,E). As for parameterization, it appears that P B opt and E 0 /(E 0 +4.1) values had little influence on VGPM results (PP), regardless of whether the values were obtained from in situ measurements or satellite-based observation.
We did not observe large monsoonal variations in K d−in situ and K d−sate ; however, the highest K d−in situ (0.072 ± 0.009 m −1 ) was obtained in December and the highest K d−sate (0.052 ± 0.009 m −1 ) was obtained in January. These high K d values resulted in a shallower Z eu−in situ (64 ± 8 m) and Z eu−sate (92 ± 19 m), compared to the values converted from the K d−in situ (0.053 ± 0.007 m −1 ) for August (Z eu−in situ = 89 ± 13 m) and the K d−sate (0.029 ± 0.003 m −1 ) for September (Z eu−sate = 162 ± 22 m) (Figures 3F,G). The shallow Z eu observed in December and January can be attributed mainly to an increase in phytoplankton biomass that reduces the light penetration in that region during the NEM (Chen, 2005;Tseng et al., 2005;Chen et al., 2008). Figure 3H presents temporal variations in Chl s−in situ and Chl s−sate . Conspicuously high Chl s−in situ values were observed in December (0.40 ± 0.15 mg m −3 ) and high Chl s−sate values were observed in January (0.29 ± 0.09 mg m −3 ). The high Chl s values observed throughout the NEM were triggered by the monsoonal force, which stimulated phytoplankton photosynthesis, thereby increasing the phytoplankton abundance or enriching the area with Chl s from subsurface waters (Liu et al., 2002;Tseng et al., 2005). The drop in Chl s−in situ and Chl s−sate values to nearly < 0.10 mg m −3 during the SWM has previously been reported in studies focusing on the same region of the SCS (Liu et al., 2002;Chen, 2005;Tseng et al., 2005;Shih et al., 2020a). Similar findings were also recorded in oligotrophic time-series studies, such as HOT and BATS (Karl et al., 2003).

Relationships Among of PP, SST, E 0 , K d and Chl s in in situ Measurements and Satellite-Based Observations
PP in situ was generally lower than PP VGPM ; however, we observed a significantly positive correlation between these two parameters; i.e., slope = 0.70, r 2 = 0.42, p < 0.01 ( Figure 4A). This clearly indicates the feasibility of the VGPM for predictions; however, tuning would be required for the study area. As shown in Figure 4, we identified significantly positive correlations between PP in situ and PP VGPM as well as their respective variables SST, E 0 , K d and Chl s . As indicated by high r 2 and low p values with slopes close to 1, the most significant correlations were found in SST and E 0 (Figures 4C-F) : SST (slope = 0.92, r 2 = 0.92, p < 0.01) and E 0 (slope = 0.62, r 2 = 0.58, p < 0.01). Taken together, these results indicate that SST sate and E 0−sate were the variables most strongly correlated with SST in situ and E 0−in situ , providing the most accurate estimates of PP when using the VGPM. As indicated by the 1:1 lines in Figures 4B,G,H, the parameters with the greatest variability in terms of slope were Chl s (slope = 0.41, r 2 = 0.74, p < 0.01) and K d (slope = 0.59, r 2 = 0.50, p < 0.01). This suggests Chl s−sate and K d−sate could potentially bias PP estimates obtained using the VGPM.
The fact that the regression line between PP VGPM and PP in situ lies above the 1:1 line indicates that PP VGPM estimates exceeded PP in situ . When using the VGPM to estimate PP, the two main variables are P B opt and E 0 /(E 0 +4.1), derived from in situ measurements and satellite-based observations of SST (SST in situ and SST sate ) and E 0 (E 0−in situ and E 0−sate ). Between SST in situ and SST sate as well as between E 0−in situ and E 0−sate , we observed nearly linear relationships (i.e., close  Frontiers in Marine Science | www.frontiersin.org to the 1:1 diagonal). This suggests that P B opt and E 0 /(E 0 +4.1) were not scaling variables governing the magnitude of PP in this study. Influences derived from these two variables [P B opt , E 0 /(E 0 +4.1)] of the both methods on calculated results were less than 1%. By contrast, the correlations between Chl s−in situ and Chl s−sate , and Z eu−in situ and Z eu−sate (converted from K d−in situ and K d−sate , respectively) had a pronounced impact on VGPM PP estimates. If we considered only the difference between Chl s−in situ and Chl s−sate , then PP values estimated using VGPM would be slightly higher (∼ 5%, depended on the given case) than those obtained using in situ measurements. If we considered only the difference between Z eu−in situ and Z eu−sate , then PP values estimated using VGPM would be apparently 1.72-fold higher than those obtained using in situ measurements.

DISCUSSION
Uncertainty in Estimating PP Due to Differences Between Chl s−in situ and Chl s−sate as Well as Z eu−in situ and Z eu−sate : Implications Our analysis of revealed a number of potential uncertainties pertaining to PP estimation using the VGPM algorithm. Most of the discrepancies were due primarily to differences between Z eu−in situ and Z eu−sate as well as between Chl s−in situ and Chl s−sate . Overall, the product of Z eu and Chl s (i.e., the phytoplankton inventory in the euphotic zone), suggests that the base assumption of vertically distributed standing stock biomass in low latitude waters (SCS) may perhaps be erroneous. If so, then it will be necessary to reformulate methods for the prediction of biomass standing stocks when implementing the VGPM algorithm (Ning et al., 2004;Hill and Zimmerman, 2010).
The fundamental structure of the VGPM is based on a relationship between integrated phytoplankton biomass in the euphotic zone and Chl s (Behrenfeld and Falkowski, 1997;Hill and Zimmerman, 2010). Thus, obtaining accurate estimates of PP by comparing PP in situ and PP VGPM results depends on reliable estimates of the integrated phytoplankton biomass in the euphotic zone. However, passive satellites recording the color of the ocean surface are unable to elucidate the situation at arbitrary depths below the surface (Hill and Zimmerman, 2010;Shih et al., 2020b). Contrary to the assumption that Chl decreases exponentially with depth, most observations in the SCS revealed that the subsurface Chl maximum (SCM) was often found at great depths (Liu et al., 2002;Chen, 2005;Shih et al., 2020a,b). This makes it very difficult or even impossible to estimate the integrated biomass in the euphotic zone simply as a product of Z eu and Chl s . Enhancing the reliability of the VGPM requires that we increase the number of PP in situ observations and the corresponding variables in order to improve the correlation between our assumptions pertaining to phytoplankton integrated biomass and actual measurements obtained in the field. This is particularly important in phytoplankton populations, dynamics and assemblages in specific locations under specific conditions (Hill and Zimmerman, 2010;Shih et al., 2015Shih et al., , 2020b. Only by increasing the number of observations and enhancing our analysis of water composition will it be possible to reduce the uncertainty associated with Z eu and Chl s in estimating PP using the VGPM. The mean PP in the euphotic zone (PP/Z eu , mg-C m −3 d −1 ) presented a positive linear relationship between (PP VGPM /Z eu−sate ) and (PP in situ /Z eu−in situ ); i.e., slope: 0.51, r 2 = 0.39, p < 0.01 (Figure 5A). Higher values were observed during the NEM (4.7 ± 2.4 and 3.8 ± 2.2 mg-C m −3 d −1 of PP in situ /Z eu−in situ and PP VGPM /Z eu−sate , respectively) and lower values were observed during the SWM (1.9 ± 0.8 and 1.7 ± 0.4 mg-C m −3 d −1 of PP in situ /Z eu−in situ and PP VGPM /Z eu−sate , respectively). The slope of 0.51 for PP VGPM /Z eu−sate and PP in situ /Z eu−in situ was lower than that of PP VGPM and PP in situ (slope = 0.70), such that most data fell on FIGURE 5 | (A) Relationship between satellite-based observations of (PP VGPM / Z eu−sate ) and in situ measurements of (PP in situ / Z eu−in situ ) ratios; (B) relationship between PP in situ and SST in situ . Red and blue symbols, respectively, indicate observations made during the SW-monsoon and NE-monsoon. The diagonal in panel (A) is 1:1 line. the right side of the 1:1 line. Ratios of PP VGPM /Z eu−sate were nearly 20% lower than those of PP in situ /Z eu−in situ , indicating that Z eu−sate was deeper than Z eu−in situ , particularly during the SWM. This also indicates that the VGPM parameter Z eu−sate indeed substantially affected PP estimates in this study. It has been proposed that satellite-based K d values have to be calibrated against the zenith solar angle during the data processing in accordance with the methods outlined by Lee et al. (2005) and Li et al. (2015). However, users may not to confirm the processing from the downloaded or retrieved satellite-based products of PP and its relevant variables.
The significant linear correlation between (P B opt ) sate and (P B opt ) in situ (slope = 0.97, r 2 = 0.94, p < 0.01; Figure 4D) is simply a reflection of the relationship between the two SST records, resulting from the fact that P B opt is computed as a function of SST (i.e., a seventh-order polynomial, Equation 2-2) (Behrenfeld and Falkowski, 1997). Thus, any potential deviations between SST and P B opt had only a negligible influence on PP estimates obtained using VGPM. Many studies have nevertheless concluded that the accuracy of VGPM-based estimates of PP are poor, when using the 7th order polynomial of SST (Equation 2-2) to calculate the input of P B opt (Mizobata and Saitoh, 2004;Kameda and Ishizaka, 2005;Yamada et al., 2005;Siswanto et al., 2006;Ishizaka et al., 2007;Tang et al., 2008;Tang and Chen, 2016).
Researchers have reported that much of the uncertainty in estimating PP VGPM is related to the computation of P B opt under the effects of phytoplankton physiology, growth conditions, abundance, size, and productivity Kameda and Ishizaka, 2005;Yamada et al., 2005). It has been suggested that P B opt is influenced by SST as well as E 0 and various biological parameters, such as Chl concentration. P B opt represents an optimal daily carbon fixation rate in the water column previously described as a 7th order polynomial of SST, however, when SST exceeds 28.5 • C, P B opt remains fixed at a constant 4 mg-C mg-Chl −1 h −1 (Behrenfeld and Falkowski, 1997;Mizobata and Saitoh, 2004), indicating that P B opt tends to be underestimated when SST exceeds 28.5 • C. In low-latitude regions of the SCS, the constant P B opt mentioned above usually occurred in late spring, summer, and early fall, during which SST in situ and SST sate both exceeded 28.5 • C. This increased the margin between actual PP values and the estimates obtained by inputting P B opt derived using SST in situ or SST sate . The reliability of P B opt estimates obtained using the VGPM seventh-order polynomial SST algorithm is not universally applicable (temporally or spatially). This has prompted oceanographers to tune existing methods or devise new methods for the precise estimation of P B opt , especially for ocean water at low latitudes.
Our results reveal that the satellite-derived primary production (e.g., PP VGPM ) may significantly affect global carbon export flux to deep waters, but what is the overall significance and impact of these PP values on POC fluxes? For example, Dunne et al. (2005)  According to the expression 4-1, the POC fluxes estimated from in situ measurements (PP in situ , SST in situ and Z eu−in situ ) of Hung et al. (2000); Hung and Gong (2007), and Shih et al. (2015) were from 12 to 319 mg-C m −2 d −1 , an average of 20% less than the trap POC fluxes (25-274 mg-C m −2 d −1 ) ( Figure 6). As described above, the inputs of SST in situ and Z eu−in situ to the expression were fixed, the PP in situ was replaced by the PP VGPM , an average difference between estimated and trap POC fluxes (estimated POC flux: 16-781 mg-C m −2 d −1 ) was a factor of 2. It has been suggested that the uncertainty of these POC fluxes is quite large if satellite-based PP is used to estimate carbon sequestrations in oceans. If the discrepancy between in situ measurements and satellite-based observations of PP can be diminished and the reliance on them (e.g., PP VGPM ) can be increased, it is to dedicate the potential importance and goal of the present study.

Impact of Global Warming on in situ Observations
In both the Atlantic and Pacific oceans, increased phytoplankton biomass is observed at low latitudes during the boreal FIGURE 6 | A comparison of the difference between in situ trap measured and estimated POC fluxes. The estimated POC fluxes were computed according to the empirical model expression revealed by Dunne et al. (2005). The dashed line represented the proportional relationship of estimated POC flux and PP. The PP in situ and PP VGPM were also exerted on the expression to estimate POC fluxes (red square and blue triangle, respectively). In situ measurements of PP in situ and trap measured POC fluxes (green circle) were based on Hung et al. (2000); Hung and Gong (2007), and Shih et al. (2015). warm season. Other than equatorial upwelling, there are no other conspicuous seasonal trends in biogeochemical activities (Dandonneau et al., 2004). In the Arabian Sea, enhanced biogeochemical responses are also observed at low latitudes during the SWM (Banse and English, 2000). In this study, long-term and monthly variations in PP in situ and PP VGPM demonstrated influential monsoonal system at the SEATS (Figures 2, 3 and Table 2). The concurrence of low SST, deep Z ML , high Chl s , shallow Z eu , and increased PP suggest that increased phytoplankton biomass or specific phytoplankton communities dominating were triggered during the NEM. It has been reported that the average nitrate+nitrite (N+N) concentration, one of the key nutrients for phytoplankton growth in the SCS, in the MLD in the seasons of NEM (0.1-0.3 µM) is ∼10 times higher than of SWM (∼0.03 µM) (Tseng et al., 2005). Moreover, the inventory of N+N and the depth of nitracline in the NEM (30 ± 19 mmol m −2 and 28-62 m, respectively) have been observed a ∼ 4.5 fold higher and a ∼ 25-50% shallower than those reported in the SWM (7 ± 4 mmol m −2 and 52-82 m, respectively) (Chen, 2005;Shih et al., 2020a). Evidences of abundant nutrient (e.g., N+N) and shallow nitracline depth favoring biological activities, imply that the growth of phytoplankton communities and associated photosynthesis were affected mainly by vertical advection providing inorganic nutrients from deeper waters, under the NEM system prevailing in the SCS (Liu et al., 2002;Bai et al., 2018;Chen et al., 2020;Zhou et al., 2020).
The global decrease in PP is particularly pronounced in high-latitude waters, which lose roughly 2,000 Mt-C y −1 (Mt = 10 12 g), accounting for a 70% decline in carbon fixation via photosynthesis (Gregg et al., 2003). To compare annually reductions in PP in situ (−11 mg-C m −2 d −1 y −1 ) vs. the annual PP in situ (241 mg-C m −2 d −1 ; mean PP in situ of NEM and SWM; Table 1), it indicated that the annual reduction in carbon fixation via photosynthesis was ∼ -5% y −1 . Based on the 200 m isobath boundary of oligotrophic waters in the SCS (2.76 × 10 12 m 2 ; Lin et al., 2003), the e-ratios were 5-16% in the SCS Shih et al., 2019) and the decreased in carbon fixation was roughly 11 Mt-C, thereby accounting for 30-90% of the export production. This suggests a gradual decrease in the efficiency of photosynthetic carbon fixation by phytoplankton. Nevertheless, satellite-based observations do not show the signs of global warming on carbon fixation and sequestration. Figure 5B illustrates the significantly negative relationship between PP in situ and SST in situ at the SEATS site. The slope of −37 mg-C m −2 d −1• C −1 was exceptionally close to the −36 mg-C m −2 d −1• C −1 reported in previous studies Chen et al., 2007Chen et al., , 2008. Scaling factors related to the increase in SST in situ caused by global changes are extremely complex (Sarmiento et al., 1998;Bai et al., 2018). Notwithstanding the complexity of factors governing SST in situ , they can still be used to estimate PP values. The differences between daytime and nighttime SST in situ values were statistically insignificant, during the NEM as well as the SWM (one-tail t test: p = 0.15 and 0.24, respectively (Figure 6). We therefore surmise that SST sampling schedules had no effect on the overall results. Furthermore, we observed a strong statistically significant correlation between SST sate and SST in situ (Figure 4C), indicating the efficacy of SST in estimating PP distributions over a broad horizontal area, regardless of the method used for derivation.
The straightforward relationship between PP in situ and SST in situ is important when seeking to predict PP values and estimate new or export production in euphotic zones. Under the environmental conditions described above, the reduction in carbon fixation due to photosynthesis would be ∼ −15% • C −1 , and the decrease in carbon fixation would be roughly 37 Mt-C, thereby accounting for a 1-to 3-fold quantity of export production. The negative correlation between PP in situ and SST in situ in the current study matched the findings observed in mid-latitude tropical/subtropical regions of the Pacific and Atlantic oceans (Chen, 2000;Tilstone et al., 2009), but differed drastically from those reported in high-latitude regions (Kudryavtseva et al., 2018).
Generally speaking, the sampling resolution of in situ time-series is too low to eliminate temporal uncertainty over all timespans. Fortunately, we can use SST sate to compensate for deficiencies in SST in situ coverage (Bai et al., 2018;Chen et al., 2020). Only one daily SST sate reading can be obtained at any given location; however, it would be perfectly reasonable to substitute that value with one obtained SST in situ . We observed a statistically significant linear relationship between SST sate with SST in situ ( Figure 4C); however, differences between daytime and nighttime SST in situ measurements did not reach the level of significance (Figure 7). This suggests that SST estimates obtained using either method could be used to assess the influence of temperature on biogeochemical phenomena, such as PP.
Asynchronous variations between PP in situ and PP VGPM and their related variables in the VGPM algorithm indicate the following: (1) Satellite-based evaluations depend primarily on the assumption of an exponential decrease in the vertical distribution of PP from the surface to the euphotic zone base. Nonetheless, remote sensing cannot penetrate beyond the surface, and therefore cannot reflect the true vertical distribution of PP at depth Buitenhuis et al., 2013;Shih et al., 2013Shih et al., , 2020b. This assumption is the primary cause for uncertainty between PP values obtained from satellitebased observations (i.e., PP VGPM ) and those derived from in situ measurements (i.e., PP in situ ). (2) When using satellitebased observations to estimate horizontal PP distributions, mathematical extrapolation/interpolation is commonly used to compensate for gaps in data coverage resulting from cloud coverage, heavy rains, rough seas, extreme weather events, natural episodes, suspended particles, and/or chromophoric dissolved organic matter. Thus, this approach cannot reflect "true" or "in situ" biogeochemistry responses to PP in the oceans. (Boyd and Trull, 2007;Shang et al., 2008;Tang et al., 2008;Hung et al., 2009Shih et al., 2020a). For decades, SEATS has been used as a natural laboratory for studies of prolonged environmental changes and reciprocal biogeochemical responses. It is time to tune existing models or develop more reliable models if we are to gain meaningful estimates of biogeochemical phenomena in the oceans. The proposed calibration method aimed at improving PP estimates for the VGPM is expected to enhance our understanding of changes in the SCS.

SUMMARY
Our time-series study (2003 ∼ 2016) at SEATS compared PP estimates based on in situ measurements and those based on the VGPM in the SCS during the NEM and SWM. PP values obtained during the NEM exceeded those obtained during the SWM, which appears to indicate that weather conditions during the cold season are conducive to high PP values. PP in situ values were roughly 50% lower than PP VGPM values, regardless of the season (NEM or SWM). These discrepancies can be attributed to the satellite-based integrated phytoplankton biomass in the euphotic zone. The discrepancies can be derived as the product of Z eu and Chl s , which are two main variables in the VGPM algorithm, especially the impact of difference between in situ and satellite-based Z eu on the magnitude of PP.
The observed overall decrease in PP in situ can be partially explained by an increase in SST in situ . Our results also showed that SST sate could be used to predict horizontal PP distributions over extended time scales, based on our observation of a statistically significant relationship between SST sate with SST in situ . A significantly negative relationship between PP in situ and SST in situ appears to indicate that global changes, such as oceanic warming, could have a negative impact on ocean biogeochemistry in low-latitude regions of the SCS. Nonetheless, further research will be required to assess the influence of global changes on biogeochemical phenomena, particularly in low-latitude waters. The SEATS has been used for decades to assess the sensitivity and resilience of low-latitude oceans to environmental fluctuations. Our analysis of discrepancies between in situ measurements and satellite-based observations could help to guide revisions aimed at enhancing the robustness and reliability of the VGPM in estimating biogeochemical responses. Satellite-based data could be used to expand the spatiotemporal scale of observations and thereby shed light on the actual biogeochemical effects of global environmental changes in low-latitude regions of the SCS.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
Y-YS, F-KS, and C-CH wrote the manuscript with contributions from W-CC. C-YL, J-HT, Y-SW, and C-CL performed the experiments and created the tables and figures. Y-YS, F-KS, C-CH, W-CC, C-YL, J-HT, Y-SW, C-CL and C-YK reviewed and revised the manuscript. All authors listed have made substantial, direct, and intellectual contribution to the work and approved it for publication.