The Relationship between POC Export Efficiency and Primary Production : Opposite on the Shelf and Basin of the Northern South China Sea

Accurate estimation of particulate organic carbon (POC) export efficiency in the euphotic layer is essential to understand the efficiency of the ocean’s biological carbon pump, but field measurements are difficult to conduct and data are sparse. In this study, we investigated the relationship between POC sinking export efficiency and ocean net primary production (NPP) in the euphotic layer of the northern South China Sea (NSCS), with the help of high spatiotemporal coverage satellite-derived NPP. Annual mean POC export efficiency in euphotic zone is 34% for the shelf areas and 24% for the basin of the NSCS in the context of satellite-derived 16-day-composited NPP. Similar to what is generally observed in the global ocean, the POC export efficiency on the shelf areas appears to be strengthened with the increase of NPP. However, in the basin areas, the opposite relationship is observed. That is, the POC export efficiency significantly decreases with the increase of NPP. Seasonal decoupling between NPP and POC export, phytoplankton size structure, grazing by zooplankton, and dissolved organic carbon export might account for the observed negative relationship between the POC export efficiency and NPP in the euphotic layer of basin region. System comparison between shelf and basin would be helpful to promote understanding of the regulation mechanism of POC export in the tropical marginal seas.


Introduction
The ocean is the greatest inventory of actively cycling carbon in the world [1].A key term in the global carbon cycle is the ocean's biological carbon pump, since its magnitude and efficiency have great impact on regulating the concentration of atmospheric carbon dioxide (CO 2 ) [2][3][4].Generally, the biological carbon pump for the open ocean contains the processes of converting inorganic carbon into organic carbon and exporting particulate organic carbon (POC) to the deep ocean, where it can be sequestered on time scales ranging from seasons to centuries [2].On the continental shelf, the lateral transport of POC to the interior of adjacent deep ocean through transcontinental shelf processes, e.g., benthic turbidity current and horizontal particulate plume near the nepheloid layers is also an important carbon fixation mechanism [5].When focused on the euphotic layer, biogenic POC is mainly generated through the photosynthetic process of phytoplankton, and the ratio of biogenic POC export flux to ocean net primary production (NPP) in euphotic zone (namely, the POC export efficiency) is often used to represent the carbon sequestration efficiency of the biological carbon pump [6,7].Understanding the mechanisms that regulate the magnitude and efficiency of POC export from the euphotic layer is crucial for assessing the role of the ocean under global climate change [8], especially in the context of sustained changes in the marine ecosystems [9][10][11][12].
Generally, POC should reach a certain density (through, e.g., zooplankton grazing, formation of aggregates, etc.) before export to the deep ocean, and this process may introduce a time lag between NPP and subsequent POC export [13].Sparse field observation data are insufficient to reveal their relationships and regulation mechanisms.Therefore, several satellite-based algorithms for POC export efficiency of euphotic layer have been established based on satellite-derived NPP, sea surface temperature (SST), sea surface chlorophyll concentration (Chla), etc. [4,14,15].Based on 122 sets of global in situ NPP and corresponding POC export flux data measured at the euphotic depth, Dunne et al. [14] parameterized the POC export efficiency in terms of SST, euphotic layer depth (Z eu ) and NPP (or Chla).Henson et al. [4] defined POC export efficiency with SST, based on globally 306 samples of POC export flux ( 234 Th method) data measured at 100-m depth and 16-day-composited satellite-derived NPP data prior to the 234 Th sampling time.At the global ocean scale, POC export efficiency of euphotic layer was observed to be positively correlated with NPP and negatively correlated with SST [4,14,15], following the pattern that high NPP in the mid-high latitude (SST was low) generated more POC and large size phytoplankton, which produced high POC export flux than that in the low latitude (SST was high).However, studies in the Southern Ocean [16] and the northern Gulf of Mexico [8] showed that POC export efficiency of euphotic layer decreased with the increase of NPP.Using their modeling results, Henson et al. [13] found that both the seawater temperature and fraction of POC export carried by slow sinking particles (less than 3 m•day −1 ) were principal factors in determining POC export efficiency of euphotic layer, which enhanced with increasing NPP in warm regions, while showing an inverse relationship with NPP in relatively cold regions.In addition, Cavan et al. [17] observed that zooplankton grazing could also contribute to the negative correlation between NPP and POC export efficiency of euphotic layer by converting large, fast sinking particles into small, slow sinking particles through sloppy feeding when NPP was high.Therefore, more observation data are needed to understand the regulation mechanisms of POC export efficiency of euphotic layer at a target region.
The South China Sea (SCS) is the largest marginal sea in the world, and one which experiences strong monsoonal influences.Biological activities are affected by monsoon, river plume, upwelling, eddy, typhoon and other meso-scale processes [18][19][20][21].The northern SCS (NSCS), including both shelf (depth between 50 and 200 m) and basin (depth greater than 200 m) areas [22], is an ideal region to study the relationships between climate variability and responses of the biological carbon pump and the mechanisms involved [23].Recently, several studies focused on POC export flux of euphotic layer in the SCS have been carried out [22,[24][25][26][27][28]. Generally, POC export flux of euphotic layer shows significant spatial and temporal differences, with the maximum in winter or autumn, and the minimum in summer in the NSCS [22,[29][30][31].Cai et al. [22] reported that POC stock in euphotic zone and the small phytoplankton (particle size less than 5 µm) could play significant roles in controlling POC export flux on the shelf and basin of the NSCS.However, due to the long incubation periods necessary, there are limited simultaneous observation data of NPP and POC export flux, and the temporal and spatial variation of POC export efficiency of euphotic layer in the NSCS is not yet known.
The objectives of this research included: (1) to reveal the spatiotemporal variation of POC export efficiency in the euphotic layer of NSCS by collecting in situ POC export flux data and match-up satellite-derived NPP which has high spatial and temporal coverage; and (2) to understand the regulation mechanisms of POC export efficiency through comparison analysis of the relationships between POC export efficiency and NPP in the euphotic layer of the shelf and basin areas of the NSCS.
First, we analyzed the applicability of satellite-derived NPP products in the NSCS, and then revealed the spatiotemporal variation of POC export efficiency.Second, we studied the relationship of POC export efficiency and NPP in the euphotic layer of NSCS.The regulation mechanisms of POC export efficiency in the NSCS are discussed in Section 4.

Research Areas
The NSCS includes the Pearl River plume area, shelf and basin areas (Figure 1).Located in the monsoon regime, the NSCS is influenced by northeast monsoon from November to April, and by southwest monsoon from June to September [22,32].The basin-wide circulation in the NSCS displays a large cyclonic gyre, which persists throughout the year [23].The region west of the Luzon Strait is often influenced by the Kuroshio intrusion [33].In summer, the Pearl River plume reaches its maximum and transports large amount of nutrients to the coastal sea and inner shelf, resulting in relatively high NPP in these areas [21,[34][35][36].However, NPP is relatively low in basin areas due to nutrient deficiency caused by stable stratification [37].In winter, the NSCS is controlled by a stronger cyclonic gyre with the strengths of northeast wind and mixed layer depth reaching their peaks.The combined effects of these processes lead to stronger biological activities in winter than in summer in the NSCS.In addition, the biological activities in the NSCS are also affected by typhoons [38], meso-scale eddies [30,39], internal waves [40], etc.Generally, Chla and NPP are high in winter and low in summer in the basin of NSCS, while two peaks of NPP occur in winter and summer at the inner shelf areas [35,41].
In this study, we focused on the shelf (depth 50-200 m) and basin areas (depth > 200 m) of the NSCS.Our research on the shelf was not designed to investigate the terrestrial effects and POC export processes in shallow water.The comparison analyses between shelf and basin areas focused on the POC vertical sinking export, as these two regimes have significant different phytoplankton communities and corresponding environmental and biogeochemical characteristics.The terrestrial effects and lateral export might not be negligible, but have less influences on our research areas based on the following facts.The mixed layer depths of the research regions are generally less than 50 m, while the highest value (50-60 m) occurs in winter [42].By examining the spatial distribution of the climatology monthly average seawater transparency during 2003-2016 (Figure 1), the terrigenous particles were mainly concentrated in coastal waters (depth less than 50 m).According to the research by Cai et al. [22], which reported in situ POC export data of the same cruises as this study, the effect of river diluted water on our research areas was weak based on the measured salinity data.In addition, as reported by Cai et al. [22], the C/N ratio on particles collected from the upper 100 m during the 2009-2011 cruises was about 6.6, which was in excellent agreement with the Redfield ratio of 106:16 for C/N in phytoplankton cells and indicated that biogenic particles were dominant in our study area.Overall, the effects of POC derived from terrestrial and diluted water terrigenous POC were weak in our study regions.On the other hand, Chen [29] observed the lateral transport of POC was weak in the upper 100 m layer in our research areas by comparing the model-derived three-dimensional (considering physical transport item of water based on the Princeton Ocean Model) and one-dimensional (vertical export) POC export fluxes.

In Situ POC Export Flux and Chla Data
Two datasets of in situ POC export flux (EP in_situ ) in the NSCS were gathered from Cai et al. [22] and Chen [29] (Figure 1 and Table 1).The data from Cai et al. [22] were POC export flux measured at 25, 50 and 100 m depth at stations located in the coastal, shelf and basin of the NSCS, respectively.The data from Chen [29] were collected in the basin of the NSCS with POC export flux measured at 100-m depth.Considering the average euphotic layer depth was about 50 and 100 m on the shelf and basin of the NSCS, respectively [22,43], and also to reduce the effects of POC derived from terrigenous and diluted water, only the POC export flux data measured at layers of 50 and 100 m depth were selected for subsequent analysis in this research.We used dataset from Cai et al. [22] (hereafter, the analysis dataset) to investigate the relationships between POC export efficiency and NPP in the euphotic layer of the NSCS.The dataset from Chen [29] (hereafter, the validation dataset) was used to verify the analysis results from the "analysis dataset".In both datasets, POC export flux was measured based on 234 Th/ 238 U disequilibrium method. 234Th is a radiogenic daughter product of the naturally occurring soluble isotope 238 U. 238 U is conservatively distributed in ocean water, whereas 234 Th is easily removed from seawater by particulate adsorption, resulting in a disequilibrium between 238 U and 234 Th.Therefore, the disequilibrium between 238 U and 234 Th could be used to estimate the particles export in the ocean.For observation, 234 Th and 238 U activities at different sampling depths are measured at first, and then the 234 Th export flux at interested depth could be obtained by solving the 234 Th activity balance equation.Based on the measured 234 Th export flux at the interested depth, POC export flux can be calculated using Equation (1), in which POC/ 234 Th is measured using total suspended particles (>1 µm in size).More details about the method can be found in Cai et al. [22].
POC flux = 234 Th flux × POC/ 234 Th, In addition, in situ Chla data were also collected from Cai et al. [22], who measured Chla with high-performance liquid chromatography (HPLC).A volume of 4-16 L of seawater was filtered  We used dataset from Cai et al. [22] (hereafter, the analysis dataset) to investigate the relationships between POC export efficiency and NPP in the euphotic layer of the NSCS.The dataset from Chen [29] (hereafter, the validation dataset) was used to verify the analysis results from the "analysis dataset".In both datasets, POC export flux was measured based on 234 Th/ 238 U disequilibrium method. 234Th is a radiogenic daughter product of the naturally occurring soluble isotope 238 U. 238 U is conservatively distributed in ocean water, whereas 234 Th is easily removed from seawater by particulate adsorption, resulting in a disequilibrium between 238 U and 234 Th.Therefore, the disequilibrium between 238 U and 234 Th could be used to estimate the particles export in the ocean.For observation, 234 Th and 238 U activities at different sampling depths are measured at first, and then the 234 Th export flux at interested depth could be obtained by solving the 234 Th activity balance equation.Based on the measured 234 Th export flux at the interested depth, POC export flux can be calculated using Equation (1), in which POC/ 234 Th is measured using total suspended particles (>1 µm in size).More details about the method can be found in Cai et al. [22].
POC flux = 234 Th flux × POC/ 234 Th, In addition, in situ Chla data were also collected from Cai et al. [22], who measured Chla with high-performance liquid chromatography (HPLC).A volume of 4-16 L of seawater was filtered through a 47-mm (0.7-mm nominal pore size) GF/F filter under a gentle vacuum of less than 150 mm Hg.The filter was folded and stored in liquid nitrogen until the laboratory analysis.More details can be found in Cai et al. [22].

Data of Particle Size Distribution
During the 2009-2011 cruises, inherent optical properties (IOPs) of seawater were measured via a package of bio-optical sensors, including the AC-S Spectrophotometer (the WET Labs) that measured the absorption and attenuation coefficient, the HydroScat-6 Spectral Backscattering Sensor (HS-6) from the HOBI Labs that measured the backscattering coefficient, and the SeaBird CTD (SBE 37-SI MicroCAT) that measured temperature and salinity of seawater.All methods of measurement and calculation of IOPs followed the NASA ocean optics protocols [45].For the AC-S data, pure water calibration values obtained before the sampling period were subtracted, and a temperature-salinity correction was applied by using the Compass2 software (released by the WET Labs specifically for AC-S data processing).Then, a scattering correction was applied to the absorption data by subtracting the mean absorption value at the wavelength between 715 and 735 nm to obtain absorption coefficient (a p ). Particulate scattering coefficient (b p ) was obtained by subtracting the corrected a p from the AC-S attenuation value (c p ).For the HS-6 data, pure water calibration was also conducted before sampling period.Particulate backscattering coefficient (b bp ) was calculated by the corrected HS-6 data (b b ) subtracting the backscattering coefficient of pure seawater (b bw ) from Morel [46].
In practice, inherent optical properties of seawater, such as particulate backscattering ratio (b bp /b p ), can be used to reflect particle size distribution; large particles corresponding to small particulate backscattering ratio, and vice versa [47][48][49].In our study areas, the biogenic particles were dominant [22] and the physical transport of water had a weak effect on the POC export in the upper 100 m layer of the basin region [29], thus we used the profile of particulate backscattering ratio at 442 nm (bbp442/bp442) to study the vertical distribution of POC-related particle size within euphotic zone in the basin areas of NSCS.

Satellite Data
Satellite-derived Chla (mg•m −3 ), photosynthetic available radiation (PAR, Einstein•m −2 •day −1 ), SST ( • C), euphotic layer depth (Z eu , m), and daylength (DL, h) were used to calculate NPP.All satellite data used in the estimation of NPP in this study were daily products with spatial resolution of 1/4 • .
To obtain as many matching data as possible, we used the satellite-derived Chla and PAR data downloaded from the GlobColour Website (http://hermes.acri.fr/index.php?class=archive), which provided a merged product of the Sea-Viewing Wide Field-of-View Sensor (SeaWiFS), Moderate-resolution Imaging Spectroradiometer (MODIS) and Medium-resolution Imaging Spectrometer (MERIS).Chla was produced by using the GSM (Garver-Siegel-Maritorena) algorithm.According to the user guider of GlobColour [50], the uncertainties associated with Chla and PAR products in this merged dataset were about 35% and 7.5%, respectively.Reanalysis SST products were from the Advanced Very High Resolution Radiometer (AVHRR), downloaded from the NOAA Website (http://data.nodc.noaa.gov/pathfinder/Version5.2/).Z eu , defined as the depth where 1% of the surface PAR arrived, was derived using surface Chla [51].DL was calculated using Julian day and latitudes of samplings, and the algorithm can be found at http://orca.science.oregonstate.edu/faq01.php.
where P b opt is the maximum carbon fixation rate within a water column (mg C•mg Chla −1 •h −1 ), and is a function of SST in the VGPM model.In addition, climatology monthly mean seawater transparency distribution in NSCS (from 2003 to 2016) was retrieved by MODIS-derived reflectance, total absorption coefficient and particulate backscattering coefficient at 443 nm based on the semi-analytic algorithm of He et al. [44].The monthly mean SST (11 µ daytime) products from MODIS sensor in the Aqua satellite during the 2009-2011 cruises were also collected as a reference for the spatiotemporal variation of POC export efficiency in Figure 2. Monthly averaged satellite products used in this research were from https://oceancolor.gsfc.nasa.gov/cgi/l3.

Data of POC Export Efficiency in Euphotic Layer
Based on the in situ POC export flux data in the NSCS from Cai et al. [22] and Chen [29], we matched them with satellite-derived NPP data and the corresponding satellite-derived Chla, PAR, SST, and Zeu data.All data are shown in Tables A1 and A2.Geological location (longitude and latitude) of in situ POC sampling stations were matched by the nearest satellite pixel with 1/4° spatial resolution.To make the observation time scale of the satellite-derived NPP and 234 Th method consistent with each other, we used the 16-day-composited satellite data prior to the 234 Th sampling date, as suggested by Henson et al. [4] (see Section 3.1.1for details).In situ POC export efficiency in euphotic layer (EFin_situ) was calculated by the ratio of in situ POC export flux to satellite-derived 16day-composited NPP.POC export efficiency data greater than 1 were excluded.

Models of POC Export Efficiency
In our study, we also compared the applicability of three empirical models established at global ocean scale [4,14,15] and one model developed for the Southern Ocean [16] in the NSCS (Table 2).
Table 2. Typical POC export efficiency models for using satellite data.

Study Areas Formulas Data Sources
Global

Data of POC Export Efficiency in Euphotic Layer
Based on the in situ POC export flux data in the NSCS from Cai et al. [22] and Chen [29], we matched them with satellite-derived NPP data and the corresponding satellite-derived Chla, PAR, SST, and Z eu data.All data are shown in Tables A1 and A2.Geological location (longitude and latitude) of in situ POC sampling stations were matched by the nearest satellite pixel with 1/4 • spatial resolution.To make the observation time scale of the satellite-derived NPP and 234 Th method consistent with each other, we used the 16-day-composited satellite data prior to the 234 Th sampling date, as suggested by Henson et al. [4] (see Section 3.1.1for details).In situ POC export efficiency in euphotic layer (EF in_situ ) was calculated by the ratio of in situ POC export flux to satellite-derived 16-day-composited NPP.POC export efficiency data greater than 1 were excluded.

Models of POC Export Efficiency
In our study, we also compared the applicability of three empirical models established at global ocean scale [4,14,15] and one model developed for the Southern Ocean [16] in the NSCS (Table 2).

Applicability of VGPM Products in the NSCS
The applicability of VGPM algorithm for satellite retrieval of NPP in the NSCS has been demonstrated in previous studies [38,53,54].Based on in situ NPP data from four surveys in the NSCS (including the shelf and basin areas), Hao et al. [53] observed a significant linear relationship between the VGPM-NPP and measured NPP values (r 2 = 0.79, n = 21, p < 0.0001).By collecting in situ NPP data (for the whole SCS in different seasons) from previous researches, Tan and Shi [54] also observed a good positive correlation (R = 0.705, n = 13, p < 0.0001) between VGPM-NPP and in situ NPP data.Xu et al. [55] conducted an inter-comparison of three satellite-derived NPP products estimated by standard VGPM, Eppley-VGPM [56], and CbPM models [57], and found the NPP products estimated by Eppley-VGPM and CbPM cannot reflect seasonal variation characteristics of NPP in the SCS, while the VGPM product can not only characterize the seasonal variation of NPP, but also its average value was closer to previous in situ NPP data than the other two algorithms, making it more suitable as the algorithm for satellite-derived NPP in the basin of SCS.VGPM-NPP was also used in modeling the long-term (1958-2009) variation of POC export flux and structure of the biological carbon pump in the SCS by a three-dimensional physical-biogeochemical coupled model (ROMS-CoSiNE) [31].On the other hand, the results of Huang et al. [58] showed that CbPM-NCP (Carbon-based production model-net community production, equal to CbPM-NPP multiplied by an efficiency factor) were comparable to the oxygen mass balance estimated NCP in the central SCS, whereas the VGPM-NCP (equal to VGPM-NPP multiplied by the same efficiency factor as in calculating CbPM-NCP) were significantly lower.However, Huang et al. [58] also acknowledged that VGPM-NCP showed better performance than CbPM-NCP in representing the seasonal variation of NCP in the SCS.
As one of the main input parameters of the VGPM model, the accuracy of satellite-derived Chla is important for the applicability of VGPM products [52].Figure 2 shows the comparison results between satellite-derived and in situ Chla (Chla in-situ in Table A1) during the 2009-2011 cruises.To eliminate the effect of data missing by the cloud-mask of satellite data, three-day averaged satellite-derived Chla (3d-Chla rs in Table A1) were used in the comparison.Generally, the satellite-derived Chla were significantly correlated with in situ Chla (around the 1:1 line), with the correlation coefficient value of 0.98 and 0.96 for the shelf and basin regions, respectively, although the satellite-derived Chla were slightly lower than the in situ results.
On-deck incubation of NPP is generally carried out around the same time with the 234 Th sampling during 12 or 24 h.However, the residence time of 234 Th is considered to be about 2-20 days [59].This could introduce a mismatch in time scales between 234 Th-derived POC export and the instantaneous NPP measured in situ [4].That means the 234 Th-derived POC export would reflect the NPP situation certain days earlier.According to Ho et al. [60], the residence time of 234 Th in the upper 0-160 m layer was about 0.052 years (20 days) for the basin areas of NSCS.Liu et al. [37] also showed that the peak of POC export flux could lag the peak of NPP by about half a month in the basin of NSCS.Thus, 16-day composited satellite-derived NPP (hereafter, satellite-derived NPP) prior to the in situ measurement were used in subsequent analyses, same as the choice of Henson et al. [4] in their study on estimating POC export efficiency based on satellite-derived SST.

Seasonal Variation of POC Export Efficiency in Euphotic Layer of the NSCS
Based on the results of "analysis dataset", significant seasonal and spatial variability of POC export efficiency in euphotic layer were observed in the NSCS (Figure 3, detailed data are listed in Table A1).In general, the high POC export efficiency was concentrated on the shelf area near the Pearl River plume-dominated region, and decreased toward the interior.High POC export efficiency also existed east of Hainan Island.POC export efficiency was the highest in autumn, followed by the values in spring and summer, and was the lowest in winter.Seasonal mean EPin_situ, EFin_situ and NPP are shown in Table 3.For the shelf areas, the values of NPP were the largest in winter, followed by relatively small values in autumn, and decreased significantly in spring and summer (by nearly two folds).However, both EPin_situ and EFin_situ reached their maxima in autumn.For the basin region, NPP was the largest in winter and the smallest in summer, while EPin_situ was low in both winter and summer.Annual-averaged EFin_situ values were 0.34 and 0.24 for the shelf and basin regions, respectively, consistent with the previous studies in the NSCS [25,29,35].

Correlation between POC Export Efficiency and NPP in Euphotic Layer of the NSCS
Based on the "analysis dataset", EFin_situ in euphotic layer tended to rise with the increase of NPP on the shelf areas (Figure 4a).However, in the basin areas, EFin_situ and NPP showed a significant negative correlation (Figure 4b).While EFin_situ was calculated from NPP, EFin_situ and NPP were not independent variables and their relationships may overlook some of the variability in both NPP and EPin_situ [61].However, by checking the variation of EPin_situ with NPP, we observed that EPin_situ was positively correlated with NPP on the shelf areas (Figure 4c).For the basin areas, EPin_situ tended to Seasonal mean EP in_situ , EF in_situ and NPP are shown in Table 3.For the shelf areas, the values of NPP were the largest in winter, followed by relatively small values in autumn, and decreased significantly in spring and summer (by nearly two folds).However, both EP in_situ and EF in_situ reached their maxima in autumn.For the basin region, NPP was the largest in winter and the smallest in summer, while EP in_situ was low in both winter and summer.Annual-averaged EF in_situ values were 0.34 and 0.24 for the shelf and basin regions, respectively, consistent with the previous studies in the NSCS [25,29,35].

Correlation between POC Export Efficiency and NPP in Euphotic Layer of the NSCS
Based on the "analysis dataset", EF in_situ in euphotic layer tended to rise with the increase of NPP on the shelf areas (Figure 4a).However, in the basin areas, EF in_situ and NPP showed a significant negative correlation (Figure 4b).While EF in_situ was calculated from NPP, EF in_situ and NPP were not independent variables and their relationships may overlook some of the variability in both NPP and EP in_situ [61].However, by checking the variation of EP in_situ with NPP, we observed that EP in_situ was positively correlated with NPP on the shelf areas (Figure 4c).For the basin areas, EP in_situ tended to decrease slightly with the increase of NPP (Figure 4d).This suggests that the relationship between EF in_situ and NPP in the NSCS was caused by the biogeochemical mechanism, not an artifact of plotting EP in_situ /NPP versus NPP [61].decrease slightly with the increase of NPP (Figure 4d).This suggests that the relationship between EFin_situ and NPP in the NSCS was caused by the biogeochemical mechanism, not an artifact of plotting EPin_situ/NPP versus NPP [61].To further check the correlation of POC export efficiency and NPP in euphotic layer of the basin areas, we binned the EFin_situ data based on the NPP rates with the ranges of 0-200, 200-300, 300-400, 400-500, 500-600, 600-700, 700-1000, and 1000-1300 mg C•m −2 •day −1 (Figure 5a).Averaged NPP and EFin_situ values in each bin were used in the following analysis.Wider ranges were used at high NPP values to balance the distribution of data points.The relationship between logarithmic-NPP and POC export efficiency was robust (r = −0.97,p < 0.0001) in the basin areas.To evaluate the applicability of the linear fitting results shown in Figure 5a, POC export efficiency in euphotic layer was re-calculated by in situ POC export flux data from the "validation dataset" and corresponding satellite-derived NPP (see Table A2).The validation results also show an inverse relationship between POC export efficiency and NPP in euphotic layer for the basin areas, with the RMSE (root mean square error) value of 0.15 (Figure 5b).To further check the correlation of POC export efficiency and NPP in euphotic layer of the basin areas, we binned the EF in_situ data based on the NPP rates with the ranges of 0-200, 200-300, 300-400, 400-500, 500-600, 600-700, 700-1000, and 1000-1300 mg C•m −2 •day −1 (Figure 5a).Averaged NPP and EF in_situ values in each bin were used in the following analysis.Wider ranges were used at high NPP values to balance the distribution of data points.The relationship between logarithmic-NPP and POC export efficiency was robust (r = −0.97,p < 0.0001) in the basin areas.To evaluate the applicability of the linear fitting results shown in Figure 5a, POC export efficiency in euphotic layer was re-calculated by in situ POC export flux data from the "validation dataset" and corresponding satellite-derived NPP (see Table A2).The validation results also show an inverse relationship between POC export efficiency and NPP in euphotic layer for the basin areas, with the RMSE (root mean square error) value of 0.15 (Figure 5b).
the linear fitting results shown in Figure 5a, POC export efficiency in euphotic layer was re-calculated by in situ POC export flux data from the "validation dataset" and corresponding satellite-derived NPP (see Table A2).The validation results also show an inverse relationship between POC export efficiency and NPP in euphotic layer for the basin areas, with the RMSE (root mean square error) value of 0.15 (Figure 5b).A1); and (b) results from the "validation dataset" (see Table A2).Open circles are original matching points, and black dots are the averaged data in different NPP bins.Solid line is the regression result in Figure 5a.Error bars show the standard errors of the mean.

Comparison of Different Models of POC Export Efficiency
To contextualize within the field, we compared the calculated EF in_situ (based on the "analysis dataset" and satellite-derived NPP) in the NSCS to the model-derived POC export efficiency results (Equations shown in Table 2, Figure 6).The NPP data used by all models were from the same VGPM-NPP dataset, which could eliminate the systematic deviation of model-derived EF variation trend with VGPM-NPP in Figure 6.dataset" (see Table A2).Open circles are original matching points, and black dots are the averaged data in different NPP bins.Solid line is the regression result in Figure 5a.Error bars show the standard errors of the mean.

Comparison of Different Models of POC Export Efficiency
To contextualize within the field, we compared the calculated EFin_situ (based on the "analysis dataset" and satellite-derived NPP) in the NSCS to the model-derived POC export efficiency results (Equations shown in Table 2, Figure 6).The NPP data used by all models were from the same VGPM-NPP dataset, which could eliminate the systematic deviation of model-derived EF variation trend with VGPM-NPP in Figure 6.  2, with the parameter values given in Table A1.Dashed lines are the 1:1 lines.
For the shelf areas (Figure 6a), the averaged EFin_situ, EFDunne, EFHenson, EFLaws, and EFMaiti were 0.34, 0.11, 0.03, 0.12, and 0.32, respectively, at the selected stations.Generally, EFDunne and EFLaws showed similar trends, but their averaged relative deviations from EFin_situ were high (60% and 53% for EFDunne and EFLaws, respectively).EFMaiti and EFHenson had totally different trends from that of EFin_situ, and they were not suitable to the shelf areas of the NSCS.Located in the tropical-subtropical area with high ocean temperature throughout the year (Figure 3), the solely SST-dependent model developed with the global data in Henson et al. [4] was invalid in the NSCS.As for the equation in Maiti et al. [16], the characteristic on the shelf of the NSCS was also different from that in the Southern Ocean.
For the basin areas (Figure 6b), the averaged EFin_situ, EFDunne, EFHenson, EFLaws, and EFMaiti were 0.25, 0.09, 0.03, 0.11, and 0.35, respectively, at the selected stations.Only EFMaiti had similar trend and value as EFin_situ, and the other three models showed significant disagreement in the basin of NSCS.Comparing with the large range of SST on the global ocean scale, the variation of SST in the NSCS is relatively small (Figure 3); thus, NPP is mainly controlled by the nutrient levels driven by monsoonal winds, the seasonal change of mixed layer depth, etc. [37].According to recent studies, the temperature scaling for autotrophic processes is valid up to 40 °C, but the temperature dependency  2, with the parameter values given in Table A1.Dashed lines are the 1:1 lines.
For the shelf areas (Figure 6a), the averaged EF in_situ , EF Dunne , EF Henson , EF Laws , and EF Maiti were 0.34, 0.11, 0.03, 0.12, and 0.32, respectively, at the selected stations.Generally, EF Dunne and EF Laws showed similar trends, but their averaged relative deviations from EF in_situ were high (60% and 53% for EF Dunne and EF Laws , respectively).EF Maiti and EF Henson had totally different trends from that of EF in_situ , and they were not suitable to the shelf areas of the NSCS.Located in the tropical-subtropical area with high ocean temperature throughout the year (Figure 3), the solely SST-dependent model developed with the global data in Henson et al. [4] was invalid in the NSCS.As for the equation in Maiti et al. [16], the characteristic on the shelf of the NSCS was also different from that in the Southern Ocean.
For the basin areas (Figure 6b), the averaged EF in_situ , EF Dunne , EF Henson , EF Laws , and EF Maiti were 0.25, 0.09, 0.03, 0.11, and 0.35, respectively, at the selected stations.Only EF Maiti had similar trend and value as EF in_situ , and the other three models showed significant disagreement in the basin of NSCS.Comparing with the large range of SST on the global ocean scale, the variation of SST in the NSCS is relatively small (Figure 3); thus, NPP is mainly controlled by the nutrient levels driven by monsoonal winds, the seasonal change of mixed layer depth, etc. [37].According to recent studies, the temperature scaling for autotrophic processes is valid up to 40 • C, but the temperature dependency for zooplankton growth rate (thus, grazing rate) is only reasonable fit up to 24 • C, after which the temperature scaling breaks down and appears to flatten out due to decreased solubility of oxygen at higher temperature [62].Since the SST in shelf and basin regions of NSCS are generally greater than 20 • C throughout the year (Figure 3), the algorithm considering only the influence of temperature may not well characterize the variation of POC export efficiency in euphotic layer of the NSCS.In fact, the relationship between SST and POC export efficiency was also unstable in the dataset of Henson et al. [4] and Laws et al. [15] when SST was higher than 24 • C.

Seasonal Mean POC Export Efficiency and NPP in Euphotic Layer of the Basin Regions
Due to the long incubation periods necessary, we had less synchronous observation dataset of NPP and POC export flux in the NSCS.Therefore, we collected as many published data as possible to check the negative relationship between POC export efficiency and NPP in the basin.We used the in situ seasonal averaged NPP data from previous research [41], which included 30 sets of measurements from more than eight cruises in the NSCS (details can be found in Table A3), for the calculation of seasonal POC export efficiency in euphotic layer of basin regions.Individual field measurements of POC export flux in the basin of the NSCS were also gathered from previous researches (Table 4, only field data of POC export flux at 100-m depth were gathered).Although the inverse relationship between POC export efficiency and NPP in euphotic layer was not statistically significant on seasonal scale (Figure 7), these independent results also indicated that the POC export efficiency was not positively related to NPP in the basin areas of NSCS as suggested in the global models (Table 2).

Seasonal Mean POC Export Efficiency and NPP in Euphotic Layer of the Basin Regions
Due to the long incubation periods necessary, we had less synchronous observation dataset of NPP and POC export flux in the NSCS.Therefore, we collected as many published data as possible to check the negative relationship between POC export efficiency and NPP in the basin.We used the in situ seasonal averaged NPP data from previous research [41], which included 30 sets of measurements from more than eight cruises in the NSCS (details can be found in Table A3), for the calculation of seasonal POC export efficiency in euphotic layer of basin regions.Individual field measurements of POC export flux in the basin of the NSCS were also gathered from previous researches (Table 4, only field data of POC export flux at 100-m depth were gathered).Although the inverse relationship between POC export efficiency and NPP in euphotic layer was not statistically significant on seasonal scale (Figure 7), these independent results also indicated that the POC export efficiency was not positively related to NPP in the basin areas of NSCS as suggested in the global models (Table 2).4) while in situ NPP data were from Chen et al. [41] (Table A3).4) while in situ NPP data were from Chen et al. [41] (Table A3).

Effect of Phytoplankton Community Structure
POC export variability is regulated by a combination of possesses and mechanisms, e.g., ultimate control by particle production and proximate control by the ability to generate sinking materials through phytoplankton size composition and temperature [14,64].Generally, the larger the phytoplankton size is, the more efficient euphotic zone NPP is partitioned into export [7].Larger cells, diatoms for example, have a greater chance of being associated with ballasting materials and the silica cell wall of diatoms is less labile, thus result in greater export efficiency [65].Observation results also show that diatom-dominated phytoplankton communities in productive areas, such as oceans at high latitude and upwelling regions, have high POC export efficiency in euphotic layer; conversely, in lower latitude oligotrophic environments, where diatoms are largely absent, primary production is low and mostly regenerated, resulting in weak POC export efficiency in euphotic layer [65].On the other hand, recent researches also revealed that pico-phytoplankton could dominate POC export through both direct aggregates and indirect grazing export of pico-phytoplankton-derived particles, referred to the studies in an ultra-oligotrophic lake on Vancouver Island, the equatorial Pacific Ocean and the Arabian Sea [66][67][68].In the coastal and inter-shelf areas of the NSCS, it was reported that diatoms dominated the phytoplankton community structure in euphotic zone [22], especially in winter [36].Considering the relatively fast sinking rate accompanied by diatoms and the shallow euphotic depth (about 50 m), micro-phytoplankton may play an important role in regulating the POC export efficiency on the shelf areas of the NSCS.For the basin regions, a weak but intriguing correlation (r ≈ 0.51, p < 0.0001) between POC export flux and the fraction of haptophytes and prasinophytesin (size typically < 5 µm) was reported in Cai et al. [22].Observation of HPLC results during the 2009-2011 cruises also showed that, in the basin areas of NSCS, there was a significant positive correlation (p < 0.001) between POC export efficiency and the fraction of synechococcus and prochlorococcus (pico-phytoplankton) [69].Phytoplankton community structure tended to have different effects on POC export efficiency in euphotic layer of the shelf and basin areas, although the specific regulatory mechanisms are not clear and need further investigation.

Factors in Driving the Inverse Relationship between POC Export Efficiency and NPP in Euphotic Layer of the Basin Regions
According to previous researches [16,17,61,70], the inverse relationship between NPP and POC export efficiency in euphotic layer may be caused by several factors, such as the differences in trophic structure, grazing intensity, recycling efficiency, high bacterial activity and increased DOC export in high NPP areas.Therefore, below we focus some regulation mechanisms that may cause the negative correlation between POC export efficiency and NPP in euphotic layer of the basin regions.

Effect of Zooplankton Grazing
According to the model-simulated results of Liu et al. [37], an average of 70% of the NPP in the euphotic layer of SCS could be removed by zooplankton grazing.Zooplankton grazing would result in POC either passing through the gut and being egested as fecal pellets, being respired as CO 2 or fragmenting into smaller particles through sloppy feeding [17,71].In general, particle size is an important parameter that determines the POC export [72] as small-sized, slow sinking fecal debris produced by small meso-zooplankton may be reutilized either through remineralization or coprophagy within the euphotic layer, while large, fast sinking fecal pellets produced by large copepods significantly contribute to vertical carbon export [73,74].Figure 8 shows the profile distributions of bbp442/bp442 in different seasons in the euphotic layer of the basin regions.Generally, the sea surface bbp442/bp442 values was small in winter and autumn, corresponding to the relatively larger fraction of micro-phytoplankton in surface water in winter and autumn than in the other seasons [36,75].However, for the 100-m layer, the mean particle size was comparably smaller in winter than in the other seasons, with the largest mean value of bbp442/bp442 (0.0247).Thus, it is indicated that fecal pellets from zooplankton grazing appeared to be rapidly mineralized in euphotic zone, or the role of zooplankton grazing was more likely fragmenting large, fast sinking particles into smaller, slower sinking particles through sloppy feeding in winter for the basin areas of NSCS.

Effect of Zooplankton Grazing
According to the model-simulated results of Liu et al. [37], an average of 70% of the NPP in the euphotic layer of SCS could be removed by zooplankton grazing.Zooplankton grazing would result in POC either passing through the gut and being egested as fecal pellets, being respired as CO2 or fragmenting into smaller particles through sloppy feeding [17,71].In general, particle size is an important parameter that determines the POC export [72] as small-sized, slow sinking fecal debris produced by small meso-zooplankton may be reutilized either through remineralization or coprophagy within the euphotic layer, while large, fast sinking fecal pellets produced by large copepods significantly contribute to vertical carbon export [73,74].Figure 8 shows the profile distributions of bbp442/bp442 in different seasons in the euphotic layer of the basin regions.Generally, the sea surface bbp442/bp442 values was small in winter and autumn, corresponding to the relatively larger fraction of micro-phytoplankton in surface water in winter and autumn than in the other seasons [36,75].However, for the 100-m layer, the mean particle size was comparably smaller in winter than in the other seasons, with the largest mean value of bbp442/bp442 (0.0247).Thus, it is indicated that fecal pellets from zooplankton grazing appeared to be rapidly mineralized in euphotic zone, or the role of zooplankton grazing was more likely fragmenting large, fast sinking particles into smaller, slower sinking particles through sloppy feeding in winter for the basin areas of NSCS.

Effect of DOC Export in the Basin
In addition to the export of POC, NPP also results in the export of dissolved organic carbon (DOC) from the upper ocean [76].A large DOC export flux in high NPP season would result in a low POC export efficiency.At HOT (Hawaii Ocean Time-Series) station, DOC production could be about 17-50% of total carbon fixation by phytoplankton [77].Although the model-simulated results of Ma et al. [31] showed that the DOC export could be negligible in the basin of NSCS on annual average scale, with the net advection and diffusion DOC flux was about −0.02 mmol C•m −2 •day −1 (upward), Chen [29] observed that the DOC export could be relatively high in spring and winter.Accordingly, average net DOC diffusion export was about 0.91 mmol C•m −2 •day −1 (downward) in spring of NSCS, which may account for 30% of the POC export flux [29].When in situ total organic carbon (TOC) inventory was used as an approximation of DOC, the seasonal variation of TOC inventory (23.80, 25.20 and 23.60 Tg in spring, autumn and winter, respectively [78]) in the basin of NSCS suggested a relatively strong TOC (hence DOC) export in winter.Given the lack of in situ DOC export data in the basin of NSCS, we attempted to summarize existing results in this section rather than quantitatively assess the impact of DOC export on the inverse relationship between NPP and POC export efficiency in the basin.More field DOC export flux observations are needed to quantify the role of DOC export in driving the inverse relationship between POC export efficiency and NPP in euphotic layer of the basin regions.

Seasonal Decoupling between Primary Production and POC Export
Before organic carbon produced by NPP is converted into sinking flux, it goes through many processes, such as aggregation of individual particles until reaching certain density or size before sinking, or grazing by zooplankton and subsequent excretion as fecal pellets [13].These processes could lead to uncoupling between production and POC export, thus influencing the relationship between POC export efficiency and NPP.Unlike typical oceans in the tropics, the NPP and POC export in the NSCS are affected by monsoons with significant seasonal changes, and bloom occurs in the cold eddy or local upwelling northwest of Luzon in winter [36].Cai et al. [22] reported that the seasonal average Chla stock was almost a factor of two higher in winter than in summer, while the seasonal average POC stock was only about 30% higher in winter than in summer in the basin of NSCS.Phytoplankton blooms might be developing in the NSCS basin during a winter survey and this may lead to an elevated NPP level, but not yet immediately enhanced POC export flux until the following spring, thus resulting in low POC export efficiency in winter in the basin of NSCS.Chen et al. [79] observed that grazed NPP of micro-zooplankton for the basin waters of the NSCS was smaller in winter (73% and 65% in summer and winter, respectively).According to previous research, the peak of POC export flux of euphotic layer could lag the peak of primary production by about half a month in the basin of NSCS [37].In this research, we used the 16-day-composited satellite-derived NPP data prior to the 234 Th sampling time, which can reduce the decoupling effect of NPP and POC export in the euphotic layer to some extent.However, due to the complex processes through the food web, more parameterizations need to be developed.

Conclusions
The carbon sequestration by biological carbon pump plays an important role in the global ocean carbon cycle.However, estimates of the strength and efficiency of the biological carbon pump remain poorly quantified due to the lack of in situ data.In this study, we utilized two datasets of seasonal 234 Th-derived POC export flux at the euphotic depth and satellite-derived 16-day-composited NPP prior to the 234 Th sampling time to examine the correlation between POC export efficiency and NPP in the euphotic layer of NSCS.The results revealed that, on the shelf areas of the NSCS, POC export efficiency in euphotic layer tended to rise with the increase of NPP, and the biological processes associated with large phytoplankton might play an important role in leading to this result.On the other hand, an inverse relationship between POC export efficiency and NPP was observed in euphotic layer of the basin areas.Grazing by zooplankton, DOC export and seasonal decoupling between NPP and POC export in euphotic layer might contribute to the inverse relationship.The global POC export efficiency models are inappropriate in the NSCS.More studies (i.e., POC sink rate, POC inventory, POC composition, etc.) are needed for understanding the evolution of carbon sequestration capacity of the NSCS in the context of climate change.
and Chla is (in situ sea surface Chla, mg•m −3 ) data are from Cai et al. [22].We matched the EP in_situ with the satellite-derived NPP and relative parameters (with the subscript rs).NPP rs (net primary production, mmol C•m −2 •day −1 ) was calculated by using the satellite-derived 16-day-composited SST rs (sea surface temperature, • C), PAR rs (photosynthetically active radiation, Einstein•m −1 •day −1 ), and Chla rs (16d-Chla rs , surface chlorophyll content, mg•m −3 ) prior to 234 Th sampling time.3d-Chla rs is three-day averaged satellite-derived Chla.EF in_situ is the POC export efficiency calculated by EF in_situ = EP in_situ /NPP rs .EF in_situ values greater than 1 were excluded from the analysis.

Sustainability 2018 , 22 Figure 1 .
Figure 1.Spatial distribution of sampling stations.Diamonds and boxes symbols show the POC export flux stations in analysis dataset (2009-2011 cruises) and validation dataset (2004-2005 cruises), respectively.Crosses symbols indicate the IOP stations during the 2009-2011 cruises.Black contours indicate bathymetry at 50 and 200 m, respectively.The background is climatology monthly mean seawater transparency in NSCS from 2003 to 2016 (only results of the corresponding months during 2009-2011 cruises are shown).Seawater transparency was retrieved by MODIS-derived reflectance, total absorption coefficient and particulate backscattering coefficient at 443 nm based on the semianalytic algorithm of He et al. [44].

Figure 1 .
Figure 1.Spatial distribution of sampling stations.Diamonds and boxes symbols show the POC export flux stations in analysis dataset (2009-2011 cruises) and validation dataset (2004-2005 cruises), respectively.Crosses symbols indicate the IOP stations during the 2009-2011 cruises.Black contours indicate bathymetry at 50 and 200 m, respectively.The background is climatology monthly mean seawater transparency in NSCS from 2003 to 2016 (only results of the corresponding months during 2009-2011 cruises are shown).Seawater transparency was retrieved by MODIS-derived reflectance, total absorption coefficient and particulate backscattering coefficient at 443 nm based on the semi-analytic algorithm of He et al. [44].

Sustainability 2018 ,
10, x FOR PEER REVIEW 6 of 22backscattering coefficient at 443 nm based on the semi-analytic algorithm of He et al.[44].The monthly mean SST (11 µ daytime) products from MODIS sensor in the Aqua satellite during the 2009-2011 cruises were also collected as a reference for the spatiotemporal variation of POC export efficiency in Figure2.Monthly averaged satellite products used in this research were from https://oceancolor.gsfc.nasa.gov/cgi/l3.

Figure 2 .
Figure 2. Comparison results between three-day satellite-derived Chla (Chla RS ) and in situ Chla (Chla in-situ ) on the shelf (a) and basin (b) of the NSCS, respectively, during the 2009-2011 cruises.

22 Figure 3 .
Figure 3. Spatiotemporal distributions of POC export efficiency in euphotic layer of the NSCS in four seasons (Table A1).The backgrounds are monthly averaged SST from MODIS/Aqua during the 2009-2011 cruises.Two black lines indicate the 50-and 200-m isobaths which are used to divide the coastal, shelf and basin areas in the NSCS.

Figure 3 .
Figure 3. Spatiotemporal distributions of POC export efficiency in euphotic layer of the NSCS in four seasons (Table A1).The backgrounds are monthly averaged SST from MODIS/Aqua during the 2009-2011 cruises.Two black lines indicate the 50-and 200-m isobaths which are used to divide the coastal, shelf and basin areas in the NSCS.

Figure 4 .
Figure 4. Correlations between POC export efficiency (EFin_situ) and NPP in euphotic layer (top), and between POC export flux (EPin_situ) and NPP (bottom): (a,c) for shelf areas; and (b,d) for basin region.Black lines are the linear fitting results.

Figure 4 .
Figure 4. Correlations between POC export efficiency (EF in_situ ) and NPP in euphotic layer (top), and between POC export flux (EP in_situ ) and NPP (bottom): (a,c) for shelf areas; and (b,d) for basin region.Black lines are the linear fitting results.

Figure 5 .
Figure 5. Relationship between NPP and POC export efficiency (EF) in euphotic layer of the basin regions: (a) results based on the "analysis dataset" (see Table A1); and (b) results from the "validation

Figure 5 .
Figure 5. Relationship between NPP and POC export efficiency (EF) in euphotic layer of the basin regions: (a) results based on the "analysis dataset" (see TableA1); and (b) results from the "validation dataset" (see TableA2).Open circles are original matching points, and black dots are the averaged data in different NPP bins.Solid line is the regression result in Figure5a.Error bars show the standard errors of the mean.

Figure 6 .
Figure 6.Comparison between EFin_situ and model-derived EF on the shelf (a) and basin (b) regions of NSCS.Equations of models are shown in Table2, with the parameter values given in TableA1.Dashed lines are the 1:1 lines.

Figure 6 .
Figure 6.Comparison between EF in_situ and model-derived EF on the shelf (a) and basin (b) regions of NSCS.Equations of models are shown in Table2, with the parameter values given in TableA1.Dashed lines are the 1:1 lines.

Figure 7 .
Figure 7. Seasonal variation of POC export efficiency along with NPP in euphotic layer of the basin areas.Black line is the linear fitting result.In situ POC export flux data were collected from previous studies in the basin areas of NSCS (Table4) while in situ NPP data were from Chen et al.[41] (TableA3).

Figure 7 .
Figure 7. Seasonal variation of POC export efficiency along with NPP in euphotic layer of the basin areas.Black line is the linear fitting result.In situ POC export flux data were collected from previous studies in the basin areas of NSCS (Table4) while in situ NPP data were from Chen et al.[41] (TableA3).

Figure 8 .
Figure 8. Profile distributions of particulate backscattering ratio at 442 nm (bbp442/bp442) in four seasons for the basin areas of NSCS.Black lines are the linear fitting results.See Figure 1 for the locations of optical stations.Different colors indicate profile data for different sampling stations.4.3.2.Effect of DOC Export in the Basin In addition to the export of POC, NPP also results in the export of dissolved organic carbon (DOC) from the upper ocean [76].A large DOC export flux in high NPP season would result in a low POC export efficiency.At HOT (Hawaii Ocean Time-Series) station, DOC production could be about 17-50% of total carbon fixation by phytoplankton [77].Although the model-simulated results of Ma

Figure 8 .
Figure 8. Profile distributions of particulate backscattering ratio at 442 nm (bbp442/bp442) in four seasons for the basin areas of NSCS.Black lines are the linear fitting results.See Figure 1 for the locations of optical stations.Different colors indicate profile data for different sampling stations.

Table 1 .
Details about the in situ data used in this research.

Table 1 .
Details about the in situ data used in this research.

Table 2 .
Typical POC export efficiency models for using satellite data.

Table 4 .
In situ POC export flux data collected from previous studies in the basin areas of NSCS.

Table 4 .
In situ POC export flux data collected from previous studies in the basin areas of NSCS.

Table A2 .
The validation dataset and corresponding satellite-derived products.

Table A2 .
[29].: EP in_situ and 3D_EP in_situ data are from Chen[29].EF rs was retrieved using EP in_situ and satellite-derived 16-day-composited NPP rs according to the algorithm in Figure4b.3D_EPin_situare POC export fluxes considering physical transport item of water based on POM (Princeton Ocean Model).Other parameters are the same as those in TableA1. Note

Table A3 .
Previous report of measured seasonal mean NPP in the basin of NSCS.