Spatial and Temporal Variations of Particulate Organic Carbon Sinking Flux in Global Ocean from 2003 to 2018

: The monitoring of particulate organic carbon (POC) ﬂux at the bottom of the euphotic layer in global ocean using remote sensing satellite data plays an important role in clarifying and evaluating the ocean carbon cycle. Based on the in situ POC ﬂux data, this paper evaluated various estimation models. The global ocean POC ﬂux from 2003 to 2018 was calculated using the optimal model, and its temporal and spatial variation characteristics were analyzed. In general, the annual average of global ocean POC ﬂux is about 8.5–14.3 Gt C yr − 1 for period of 2003–2018. In the spatial dimension, the POC ﬂux in the mid-latitude ocean (30–60 ◦ ) is higher than that in the low-latitude (0–30 ◦ ). The POC ﬂux in Continental Margins with water depth less than 2000 m accounted for 30% of global ocean, which should receive more attention in global carbon cycle research. In the time dimension, the global POC ﬂux decreases year by year generally, but the POC ﬂux abnormally decreases during El Niño and increases during La Niña. In addition, due to global warming, sea ice melting, and bipolar sea area expansion, POC ﬂux in high-latitude oceans (60–90 ◦ ) is increasing year by year.


Introduction
The ocean, an important carbon pool in the world [1,2], has absorbed about one-third of the carbon released by humans [3]. The transport process of carbon from the surface layer to the deep ocean mainly includes two processes, physics and biology, and the biological process is called Biological Carbon Pump (BCP) [4][5][6]. BCP plays an important regulatory role in global atmospheric carbon dioxide content and is also an important indicator of global carbon cycle research [7,8]. Current studies show the strength of BCP decline in response to the global climate change, resulting in the reduction of the ocean carbon storage and the increase of atmospheric CO 2 [9,10]. Therefore, it has great significance to monitor the BCP's efficiency accurately for understanding the global carbon cycle and the disturbance of anthropogenic carbon dioxide emissions [10]. BCP transfers carbon elements from the ocean surface to deeper layers through a series of biological processes. Among them, the vertical migration of particulate organic carbon (POC) plays an important role in this process. Some studies show the sinking flux of POC at the bottom of euphotic layer (POC flux) directly reflects the efficiency of BCP [11].
POC flux is the key to connecting the net primary productivity (NPP) and the seabed sedimentation carbon. Using solar energy and dissolved inorganic nutrients, phytoplankton absorb CO 2 and convert it into POC, most of which is converted back to CO 2 and released to the atmosphere when passes through consumers in the surface ocean, while only a little enters the deep sea water [12]. This process is involved in the deep sea water cycle, which directly leads to the reduction of carbon content in the  [41][42][43]. A total of 285 in situ POC flux samples from 2003 and 2016 were collected at different water depths. After data preprocessing, outlier removal and satellite synchronization matching, 230 POC flux data were used as the verification dataset, accounting for 80.7% of the original data.

POC Flux Estimation Models
The primary production in the oceans that results from allochthonous nutrient inputs to the euphotic zone is called "new production", while the primary production that results from nutrient recycling in the surface waters is called "regenerated production", total primary productivity is the sum of them [23]. Eppley et al. used a large number of measured data to prove that under the premise of stable marine productivity system, the output ratio of POC in the surface layer of the ocean is approximately equal to the ratio of new productivity to total productivity. Therefore, they proposed that the ratio of new productivity to total productivity is the "f-ratio", which links ocean carbon cycle with ocean primary productivity obtained by remote sensing. In fact, the carbon of the upper ocean ecosystem should be stable and balanced from the perspective of long-term and large spatial scales, then the ratio of the POC flux to the total productivity ("e-ratio") is substantially equal to "f-ratio" [26]. This makes the POC flux at the bottom of euphotic layer the product of "e-ratio" and NPP, as a result, many studies focus on the estimation of "e-ratio". Several classical models are listed in Table 1, in which chl tot is the chlorophyll integral value in the euphotic zone and Z eu is the depth of ocean euphotic layer.

Methods for Model Comparison
The POC flux can be obtained by using these models mentioned above. However, the in situ POC flux data are measured at different water depths, so they cannot be compared directly. Establishing the relationship between POC fluxes at different depth are needed. Researchers have never stopped studying the vertical distribution of POC sedimentation: beginning from the POC sedimentation attenuation model proposed by Martin in 1987 [44], many researchers have studied the vertical distribution characteristics of POC flux [45][46][47][48][49]. However, Martin's model is still deemed as a powerful model to estimate POC flux at different depths when analyzing large-scale POC flux. In this paper we used Martin's model to compare the POC flux calculation results with the in situ POC flux and screen out the optimal estimation method. The formula is as follows: POC(z) refers to the POC flux at the depth of z, z 0 is the reference depth (for example, 150 m at HOT), and -0.858 is the flux attenuation coefficient.
By using the Martin's model, we transformed all estimated POC flux into the equivalent value at the water depth of the in situ data. Based on these equivalent data, we evaluated the accuracy of the above POC flux models. The accuracy of the POC flux models can be assessed by using the following criteria: logarithmic deviation (Bias), logarithmic root mean square error (RMSD), and unbiased root mean square error (uRMSD) proposed by Friedrichs [50]. At the same time, the correlation coefficient (R 2 ) and the average relative error (r.e) were used to calculate the correlation and the average dispersion between the calculated results and the in situ data. The formulas for calculating the evaluation parameters are as follows: where N is the number of sampled data, POC m is the POC flux obtained by the model, and POC f is the in situ POC flux. σ m and σ f are the logarithmic standard deviations of log 10 (POC m ) and log 10 (POC f ).

Evaluation of POC Flux Estimation Models
Common methods for calculating POC flux is to multiply e-ratio and NPP (Equation (7)). The accuracy of the global POC flux is determined by both e-ratio models and different NPP products.
In this paper, we used seven classic e-ratio estimation models listed in Table 1 and three types of NPP products from SAbPM, VGPM, and CbPM to examine the accuracy and stability of different POC flux estimation models, and then chose the best combination for our research. There are three main steps involved: (1) calculating POC flux separately by using different e-ratio models based on remote-sensing derived parameters; (2) converting the calculated POC flux from step 1 to the data with the same water depth of in situ data by using Martin's model; (3) comparing the converted POC flux (step 2) with the in situ data using Equations (2)- (6). Table 2 lists the error evaluation results of different models using different NPP products. From Table 2, we can see Dunne (2005a) model using NPP data from SAbPM and CbPM models have higher accuracy than other e-ratio models using the same two NPP models, while the accuracy of Dunne (2005a) model using SAbPM is better. Although e-ratio models of Baines (1994) and Laws (2011a and 2011b) outperformed the Dunne (2005a) model by using NPP data from VGPM, the accuracy of these e-ratio models are much worse if using NPP Remote Sens. 2019, 11, 2941 5 of 20 data from SAbPM and CbPM models. Therefore, we figured out that e-ratio models of Dunne (2005a) combined with NPP from SAbPM have better stability and accuracy. As the in situ data of high-flux regions are collected from different water depths, some of them are much higher than POC flux in low-flux regions. In order to better reflect the relationship between the in situ POC flux and the estimated results both in high and low flux regions, we converted the in situ data and the inversion data in high-flux regions into the data at the water depth of 150 m based on Martin's model. Then, Figure 1 was plotted to depict the performance of Dunne (2005a) e-ratio model by using different NPP models. It can be seen that the POC flux calculated with SAbPM is closer to the in situ data, with the minimum dispersion and the highest correlation, which is consistent with the error calculation results (R 2 is 0.50 and r.e is 30%) in Table 2.
Based on the analysis above, we chose the Dunne (2005a) e-ratio model and NPP product from SAbPM model to calculate the global monthly average POC flux. Although its average relative error is 30%, which cannot satisfy the detailed study of POC flux, it is good enough to analyze the spatial and temporal characteristics of global ocean POC flux [51].

Inversion Results of POC Flux
There are two main steps involved in calculating the POC flux in our research: 1) Calculating the global oceanic POC output ratio (namely e-ratio) at the bottom of the ocean euphotic layer using the Dunne (2005a) e-ratio model and the monthly average Chl, SST, and Z product spanning from 2003 to 2018 (download from https://oceandata.sci.gsfc.nasa.gov/); 2) Calculating the global monthly POC flux by multiply the e-ratio value and NPP from the SAbPM model (http://geodoi.ac.cn/WebCn/doi.aspx?Id=1280).

Inversion Results of POC Flux
There are two main steps involved in calculating the POC flux in our research: (1) Calculating the global oceanic POC output ratio (namely e-ratio) at the bottom of the ocean euphotic layer using the Dunne (2005a) e-ratio model and the monthly average Chl, SST, and Z eu product spanning from 2003 to 2018 (download from https://oceandata.sci.gsfc.nasa.gov/); (2) Calculating the global monthly POC flux by multiply the e-ratio value and NPP from the SAbPM model (http://geodoi.ac.cn/WebCn/doi.aspx?Id=1280). Figure 2 is the spatial distribution of the global average e-ratio at the bottom of euphotic layer from 2003 to 2018. From it we can see that the e-ratio is less than 0.1 in the low-latitude ocean, while the e-ratio is higher than 0.3 in the high-latitude ocean. The e-ratio in coastal ocean is generally larger than that in deep ocean. 1) Calculating the global oceanic POC output ratio (namely e-ratio) at the bottom of the ocean euphotic layer using the Dunne (2005a) e-ratio model and the monthly average Chl, SST, and Z product spanning from 2003 to 2018 (download from https://oceandata.sci.gsfc.nasa.gov/); 2) Calculating the global monthly POC flux by multiply the e-ratio value and NPP from the SAbPM model (http://geodoi.ac.cn/WebCn/doi.aspx?Id=1280).    From it we can see that the e-ratio is less than 0.1 in the low-latitude ocean, while the e-ratio is higher than 0.3 in the high-latitude ocean. The e-ratio in coastal ocean is generally larger than that in deep ocean. Figure 3 depicts the spatial distribution of the annual average POC flux from 2003 to 2018. It can be seen that global annual POC flux has different distribution characteristics at different latitudes: annual average POC flux in the low-latitude ocean near the equator is below 30 g C per square meter per year, while POC flux of the ocean above 45° in both hemispheres is as high as 70 g C per square meter per year, or even more than 100 in continental costal ocean.
According to our estimates, the global average annual POC flux is approximately 11 Gt C per year, which is close to the estimated results of most models except from Henson and Lutz (Table 3) [10,26,27,52,53]. Considering the relative error of the model (30%, in Table 2), the range of the annual average is between 8.5 and 14.3 Gt C per year. Table 3. Annual average of global total POC flux in several models. According to our estimates, the global average annual POC flux is approximately 11 Gt C per year, which is close to the estimated results of most models except from Henson and Lutz (Table 3) [10,26,27,52,53]. Considering the relative error of the model (30%, in Table 2), the range of the annual average is between 8.5 and 14.3 Gt C per year. Table 3. Annual average of global total POC flux in several models.

Estimated Sources Global Total POC Flux (Gt C yr −1 ) Instructions
The result of this paper 8. The entry of POC into deep seawater directly leads to a reduction in the carbon content of surface seawater, which further regulates the CO 2 flux and promotes the exchange of CO 2 at the air-sea interface [12]. Thence, we also conducted research on the spatial distribution of global sedimentation of organic carbon at ocean bottom (F POC ). First, we obtained global ocean depth data (https://www.bodc.ac.uk/) with a spatial resolution of 1 rad from the British Ocean Data Center (BODC) and resampled it to 9 km. Then based on the global POC flux data, the global ocean depth data and Martin's model, F POC from 2003 to 2018 were acquired (shown in Figure 4). From Figure 4, we can see the deposition flux of organic carbon in offshore seabed is obviously higher than that in deep ocean.  Figure 4). From Figure 4, we can see the deposition flux of organic carbon in offshore seabed is obviously higher than that in deep ocean.

Spatial and Temporal Variations of POC Flux
In order to explore the spatial and temporal distribution of total NPP, POC flux, and F , we calculated the integration of each monthly parameter respectively over each 1° latitude band from 2003 to 2018. The spatial and temporal variations distribution map of three parameters with latitude as y-axis and time as x-axis are depicted in Figure 5.

Spatial and Temporal Variations of POC Flux
In order to explore the spatial and temporal distribution of total NPP, POC flux, and F POC , we calculated the integration of each monthly parameter respectively over each 1 • latitude band from 2003 to 2018. The spatial and temporal variations distribution map of three parameters with latitude as y-axis and time as x-axis are depicted in Figure 5

Spatial and Temporal Variations of POC Flux
In order to explore the spatial and temporal distribution of total NPP, POC flux, and F , we calculated the integration of each monthly parameter respectively over each 1° latitude band from 2003 to 2018. The spatial and temporal variations distribution map of three parameters with latitude as y-axis and time as x-axis are depicted in Figure 5.

Spatial Variation
From Figure 5 we can see the POC flux and NPP show very similar spatial distribution characteristics. Both NPP and POC flux have peaks in the northern and southern hemispheres, but the peak in the southern hemisphere is significantly higher than that in the northern hemisphere. The high peak of POC flux in the region of 40-50°S can be attributed to high value of NPP formed by the westerly drift, high value of e-ratio, and also the large proportion of ocean area in the southern hemisphere and so on. In addition, the POC flux and NPP in both hemispheres show obviously seasonal variation due to the sun's direct point moving back and forth between the north and south tropics. For example, in the winter, POC flux and NPP are at their low peaks in the northern

Spatial Variation
From Figure 5 we can see the POC flux and NPP show very similar spatial distribution characteristics. Both NPP and POC flux have peaks in the northern and southern hemispheres, but the peak in the southern hemisphere is significantly higher than that in the northern hemisphere. The high peak of POC flux in the region of 40-50 • S can be attributed to high value of NPP formed by the westerly drift, high value of e-ratio, and also the large proportion of ocean area in the southern hemisphere and so on. In addition, the POC flux and NPP in both hemispheres show obviously seasonal variation due to the sun's direct point moving back and forth between the north and south tropics. For example, in the winter, POC flux and NPP are at their low peaks in the northern hemisphere while POC flux and NPP are at their high peaks in the southern hemisphere, which is bathed in summer sunlight and heat. From Figure 5 we also found that although the high temperature and high humidity of the equatorial sea lead to a relatively higher NPP value, the total POC flux is significantly lower due to the low e-ratio in the low-latitude ocean (30 • S-30 • N), which accounts for only 30%-50% of POC flux in the mid-latitude ocean.
In order to further explore the spatial distribution characteristics of these parameters, we divided the world into several regions from two perspectives: six regions according to latitude (low-latitude (0-30 • ), mid-latitude (30-60 • ), and high-latitude (60-90 • ) in different hemispheres, respectively) and two regions according to depth of ocean (Continental Margins and Deep Ocean). In this paper, we divided water depth <= 2000 m as the Continental Margins (total area is about 4.8 × 10 7 km 2 ) and the water depth >2000 m as the Deep Ocean (total area is about 31 × 10 7 km 2 ) according to the classification criteria of Liu and Muller-Karger [54][55][56]. Table 4 presents the statistics about the proportion of NPP, POC, and F POC in different regions. From Table 4, we found that the percentage of NPP at low-latitude accounts for about 56% due to the sufficient illumination and vast sea area, while POC flux in the same region only accounts for 28.7% owing to the low e-ratio. The total NPP (25.4%) and POC flux (37.0%) in the mid-latitude ocean in the southern hemisphere are higher than those in the northern hemisphere, which can be attributed to the lager sea area in the southern hemisphere. In the high-latitude ocean, the percentage of NPP and POC flux are quite low, accounting for only 4% and 10%, due to the limitation of illumination conditions and smaller ocean area. Though the NPP and POC flux in Deep Ocean account for 81.9% and 70.5% respectively, the proportion of NPP and POC flux in the Continental Margins are also important enough (18.1% and 29.5%).
When POC settles to the bottom of the ocean, the change of F POC between different latitudes is not obvious except for the middle and high latitudes in the northern hemisphere ( Figure 5). As shown in Table 4, the percentages of F POC in the northern hemisphere in generally higher than that in the southern hemisphere and this pattern is particularly evident in high latitudes. The percentages of F POC in the high-latitude ocean of the northern and southern hemisphere are 14.6% and 0.4%, respectively. For the Continental Margins, the proportion of F POC is very high (over 90%) due to the shallower water depth and higher POC flux, while in the Deep Ocean, deeper sea depth and longer sedimentation period greatly reduce the proportion of F POC (less than 10%).

Temporal Variation
From Figure 5 we can observe not only the spatial but also the temporal distribution characteristics of NPP, POC flux, and F POC . NPP and POC flux show a downward trend overall, especially after 2008. For example, from 2003 to 2018 the NPP in the equatorial sea area is reduced from 0.1 Pg C deg −1 mon −1 to 0.07 or even lower, and the peak's color of NPP in the southern hemisphere also becomes shallower from deep red which represents 0.12 Pg C deg −1 mon −1 . Unlike NPP, the POC flux shows a significant downward trend in the mid-high latitude ocean of both two hemispheres. From 2003 to 2018, the peaks of POC flux in the southern hemisphere decline from 0.03 Pg C deg −1 mon −1 to about 0.025, while peaks in the northern hemisphere decrease from 0.02 to about 0.015. Due to the very low e-ratio in the equatorial ocean, the difference of POC flux in the equatorial region is not obvious. As for the F POC , the slowness of the sedimentation process and the long period of the deep-sea change (hundred years) make the inter-annual variation of F POC very weak [57]. Therefore, we cannot observe the obvious change of F POC . Figure 6 is the trend graph for global NPP and POC flux from 2003 to 2018. In this figure, the monthly NPP and POC flux shows an overall decreasing trend. The slope of the NPP trend line is -0.0028, while the POC flux's is -0.00033. Therefore, we determined that the POC flux is decreasing slower than NPP which can be attributed to the low e-ratio (less than 0.7). Although the general trend of POC flux is declining year by year, there are some abnormal years. For example, the black line in Figure 6 has many abnormal bulges, which is much higher than the trend line. In particular, the monthly total POC flux increased abnormally during 2010-2011 and 2016-2017.

POC Flux In Continental Margins
Compared with the open oceans, the biological productivity and POC flux of the Continental Margins is sufficiently higher due to the higher rates of nutrient supply through upwelling, riverine inputs, and terrestrial runoff [20]. Some studies show that Continental Margins contribute to 15%-20% of global primary production [58]. According to our calculations, Continental Margins, with water depth less than 2000 m account for only 13% of global ocean surface area, but they contribute to approximately 18% of global NPP and 30% of global POC flux. Moreover, F of Continental Margins accounts for more than 90% of the global ocean. These results are evidence to show the importance of Continental Margins in oceanic carbon cycling.
More and more researchers are using remote sensing data to study POC flux in Continental Margins. They found that using satellite monitoring with high spatiotemporal resolution and carefully developed algorithms can greatly enhance the spatiotemporal observation capabilities of

POC Flux in Continental Margins
Compared with the open oceans, the biological productivity and POC flux of the Continental Margins is sufficiently higher due to the higher rates of nutrient supply through upwelling, riverine inputs, and terrestrial runoff [20]. Some studies show that Continental Margins contribute to 15%-20% of global primary production [58]. According to our calculations, Continental Margins, with water depth less than 2000 m account for only 13% of global ocean surface area, but they contribute to approximately 18% of global NPP and 30% of global POC flux. Moreover, F POC of Continental Margins accounts for more than 90% of the global ocean. These results are evidence to show the importance of Continental Margins in oceanic carbon cycling.
More and more researchers are using remote sensing data to study POC flux in Continental Margins. They found that using satellite monitoring with high spatiotemporal resolution and carefully developed algorithms can greatly enhance the spatiotemporal observation capabilities of POC flux on continental shelf and estuaries [59,60]. Despite the tremendous efforts made by researchers in the development of remote sensing algorithms, estimating the POC flux in coastal areas from satellite measurements still remains a challenge [60]. Among them, it is mainly limited by the low accuracy of remote sensing products and the complexity of carbon cycle in Continental Margins [61,62]. Except for remote sensing data, many researchers have also studied the POC flux of the Continental Margins by other means [42,[63][64][65]. These researchers studied the POC flux in the Continental Margins using in situ data and simulation analysis methods, and their results also have explained the importance of the change of POC flux in Continental Margins from different aspects.
Therefore, we should strengthen the study of NPP and POC flux in Continental Margins. Remote sensing data with higher spatial-temporal resolution can provide more accurate observations of NPP and POC flux. Biochemical models focus on the observations and targeted research of plankton community structure, POC composition and sinking behavior, and particulate aggregation ballasting effects [53], which can provide precise theory and process analysis. Through the effective combination of these technologies, the monitoring of POC flux in Continental Margins will be more accurate and reliable.

Anomalies During the Annual Decreasing of POC Flux
In Section 3.3.2, we noticed some abnormal phenomena in the downward trend of POC flux over years. In order to investigate these outliers, we depicted the relationship between monthly POC flux and month average POC flux between 2003 and 2018 ( Figure 7). In Figure 7, the downward trend is very obvious. Before 2008, most of monthly POC flux is above the average, but after 2008, most POC flux is below the average. We have found several special periods in this downtrend, marked with ellipses in the POC histogram, where the red ellipse indicates an abnormal increase of POC flux, and the blue ellipse indicates an abrupt decrease. Based on the existing research and related data, we suspected this unusual phenomenon is related to the El Niño-La Niña phenomenon. In Section 3.3.2, we noticed some abnormal phenomena in the downward trend of POC flux over years. In order to investigate these outliers, we depicted the relationship between monthly POC flux and month average POC flux between 2003 and 2018 ( Figure 7). In Figure 7, the downward trend is very obvious. Before 2008, most of monthly POC flux is above the average, but after 2008, most POC flux is below the average. We have found several special periods in this downtrend, marked with ellipses in the POC histogram, where the red ellipse indicates an abnormal increase of POC flux, and the blue ellipse indicates an abrupt decrease. Based on the existing research and related data, we suspected this unusual phenomenon is related to the El Niño-La Niña phenomenon.  Figure 7, which proves that the change of POC flux is indeed affected by the climatic event. This also indirectly confirms our conjecture.
The El Niño event may affect NPP by altering the nutrient supply in the equatorial Pacific. During El Niño, the weakening of the equatorial trade wind reduces the equatorial upwelling in the central and eastern equatorial Pacific [66]. The weakened upwelling reduces the nutrient supply of phytoplankton and consequently decreases primary production below the climatological mean [67].  Figure 7, which proves that the change of POC flux is indeed affected by the climatic event. This also indirectly confirms our conjecture.
The El Niño event may affect NPP by altering the nutrient supply in the equatorial Pacific. During El Niño, the weakening of the equatorial trade wind reduces the equatorial upwelling in the central and eastern equatorial Pacific [66]. The weakened upwelling reduces the nutrient supply of phytoplankton and consequently decreases primary production below the climatological mean [67]. As a result, the El Niño event led to a gradual decline in POC flux in the central and eastern equatorial Pacific. Iriart and González suggested that the high surface temperature anomaly during the El Niño episode of 1997-1998 led to lower phytoplankton biomass and primary production, and reduced export of POC [68,69]. In addition, according to studies by Behrenfeld, the El Niño phenomenon will enhance the stratification of water bodies in low-latitude ocean, making it difficult for nutrients to be replenished to surface seawater [70]. These will reduce carbon dioxide in the ocean and lead to a decrease in primary productivity of phytoplankton, and the decrease of NPP directly affects the decline of POC flux. Fagan also demonstrated a decrease of POC concentration during the El Niño period from 2014 to 2015 [71]. The El Niño and La Niña events generally alternate. In contrast to the El Niño event, the La Niña event leads to an increase in particulate organic matter and a decrease in sea surface temperature [72]. Especially after the end of the strong El Niño period, phytoplankton will erupt and the POC flux will also increase.
El Niño and La Niña affect not only tropical Pacific ecosystems, but changes in global marine ecosystems by influencing high and low pressure systems, winds, and precipitation. For example, La found the mean NPP decreased from 789 to 493 mg C m −2 d −1 between 2010 and 2013 in the Amundsen Sea polynya, Antarctica. However, at the same time, the NPP and POC flux showed an upward trend from 2010 to 2011 (La Niña) [73]. From observations of the Ocean Station in the northeastern subarctic Pacific between 1998 and 2000, M.S. Lipsen found that POC flux decreased during 1998 (El Niño) and increased during 1999 (La Niña) [74]. The results of these studies in different oceans are consistent with ours, and also prove the impact of El Niño on global POC flux.
Although our results confirmed that POC flux was closely related to the El Niño-La Niña phenomenon, it is not certain that this relationship applies to all ocean regions around the world. More detailed research is needed to demonstrate the attribution analysis of POC flux and NPP changes in regional ocean and whether they have a clear relationship with the El Niño-La Niña phenomenon. In future research, we will determine the effects of ocean currents, nutrients, and other parameters on regional ocean POC flux changes.

Increasing Trend of POC Flux in High-Latitude Ocean
According to our estimation results, the global POC flux is indeed decreasing year by year. The analysis of long time POC flux verifies Stramska's conjecture about the weakening of BCP [37], which is also consistent with the prediction by most models [9, 36,75]. In addition, both HOT and BATS are time series observation stations that provide in situ POC flux from 2003 to 2016. Their variation trends are good indicators of the global POC flux. Therefore, this article attempted to draw a chronological scatter plot using in situ data from these two sites. As shown in Figure 8, it can be seen from the trend lines and the equations that the POC fluxes of both sites are decreasing, which further proves that the trend analysis of the global POC flux over time based on remote sensing data observations is correct.
increase of grazing pressure and other loss processes caused by climate warming is an important additional factor leading to the reduction of NPP and POC flux in the future [77]. However, after running four marine ecosystem model under the RCP8.5 scenario and analyzing the results, Laufkötter showed that e-ratio changes and NPP changes are equally important for changes in POC flux [53]. It reminds us to pay attention to the impact of e-ratio when using the remote sensing data to analyze the change of POC flux, which is also the focus of subsequent research.  At present, the cause for POC flux decline has not come to an agreement. A lot of work has been done on analyzing NPP as one of the main drivers for POC flux in models [9, 36,76]. For example, Marinov suggested that the reduction of nutrient supply caused by stratification in low-latitude ocean is the main driver of NPP and POC flux changes. Additionally, Laufkötter showed that the increase of grazing pressure and other loss processes caused by climate warming is an important additional factor leading to the reduction of NPP and POC flux in the future [77]. However, after running four marine ecosystem model under the RCP8.5 scenario and analyzing the results, Laufkötter showed that e-ratio changes and NPP changes are equally important for changes in POC flux [53]. It reminds us to pay attention to the impact of e-ratio when using the remote sensing data to analyze the change of POC flux, which is also the focus of subsequent research.
It can be seen from our research that the global POC flux is decreasing over time obviously, but the variation characteristic of different regions is unclear. Therefore, we divided the global ocean into three parts according to different latitudes, and performed statistical analysis on the change of POC flux in different regions. We plotted the statistical results in time series as shown in Figure 9. The POC flux of the low-and mid-latitude ocean shows a declining trend consistent with the global POC flux. It can be seen from the slope of the trend line that the POC flux in low-latitude ocean has the fastest decline rate, which is related to the wide ocean area, trade wind, ocean current, and other factors in the equatorial region. From Table 4, total POC flux of the low-latitude ocean only accounts for 28.7% of global POC flux, lower POC output and higher temperature make low-latitude ocean an important "source area" of global ocean CO 2 [78]. The mid-latitude ocean is more stable and less affected by other environmental variables, so its reduction rate is lower than that in low-latitude regions. However, the total POC flux of mid-latitude ocean is higher than that of other regions, because the mid-latitude ocean has more continental shelves, and the e-ratio and NPP is becoming higher when closer to the continental shelves. The e-ratio up to 63.6% makes the POC flux as high as 61% in the mid-latitude ocean of both hemispheres, far higher than the proportion of NPP (40%), which makes the mid-latitude area an important "sink area" of global ocean CO 2 [79]. It can be seen from our research that the global POC flux is decreasing over time obviously, but the variation characteristic of different regions is unclear. Therefore, we divided the global ocean into three parts according to different latitudes, and performed statistical analysis on the change of POC flux in different regions. We plotted the statistical results in time series as shown in Figure 9. The POC flux of the low-and mid-latitude ocean shows a declining trend consistent with the global POC flux. It can be seen from the slope of the trend line that the POC flux in low-latitude ocean has the fastest decline rate, which is related to the wide ocean area, trade wind, ocean current, and other factors in the equatorial region. From Table 4, total POC flux of the low-latitude ocean only accounts for 28.7% of global POC flux, lower POC output and higher temperature make low-latitude ocean an important "source area" of global ocean CO [78]. The mid-latitude ocean is more stable and less affected by other environmental variables, so its reduction rate is lower than that in low-latitude regions. However, the total POC flux of mid-latitude ocean is higher than that of other regions, because the mid-latitude ocean has more continental shelves, and the e-ratio and NPP is becoming higher when closer to the continental shelves. The e-ratio up to 63.6% makes the POC flux as high as 61% in the mid-latitude ocean of both hemispheres, far higher than the proportion of NPP (40%), which makes the mid-latitude area an important "sink area" of global ocean CO [79].
Different to the mid-and low-latitude ocean, the POC flux of high-latitude ocean shows an increasing trend over time. Considering to the proportion of POC flux in high-latitude ocean of both hemispheres (8.7% and 1.5% in Table 4), the Arctic ocean is the main driver to this trend. Since NPP is an important factor affecting POC flux, we first studied the changes of NPP in the Arctic Ocean. The possible trend of NPP in Arctic Ocean over time has received a lot of attention in recent years. Arrigo and Van Dijken first reported a statistically significant increase in annual NPP of the entire Arctic Ocean, which had increased by 20% from 1998 to 2009 [80]. Arrigo extended the NPP time series to the year 2012 and demonstrated an increase of 30% [81]. Li pointed out that the Arctic's NPP increased from 2003 to 2016, and reached a maximum of 525.74 mg m d in 2016 [83]. Moreover, they found a strong correlation between the increasing trend of open water area and the growing NPP in Arctic Ocean [80,82].
With the decline of global NPP, the NPP in the Arctic Ocean has shown an upward trend, which has also led to an increase in POC flux. We suspect that this is related to global warming. Global warming and polar ice melting are generally recognized, researchers noted that the extent of Arctic sea ice has rapidly declined in recent years [83,84]. The IPCC Special Report on the Ocean and Cryosphere in a Changing Climate (SROCC: https://www.ipcc.ch/srocc/home/) indicated that the global ocean has warmed unabated since 1970 and the rate of ocean warming has more than doubled since 1993. Over the last decades, global warming has led to widespread shrinking of the cryosphere, Different to the mid-and low-latitude ocean, the POC flux of high-latitude ocean shows an increasing trend over time. Considering to the proportion of POC flux in high-latitude ocean of both hemispheres (8.7% and 1.5% in Table 4), the Arctic ocean is the main driver to this trend. Since NPP is an important factor affecting POC flux, we first studied the changes of NPP in the Arctic Ocean. The possible trend of NPP in Arctic Ocean over time has received a lot of attention in recent years. Arrigo and Van Dijken first reported a statistically significant increase in annual NPP of the entire Arctic Ocean, which had increased by 20% from 1998 to 2009 [80]. Arrigo extended the NPP time series to the year 2012 and demonstrated an increase of 30% [81]. Li pointed out that the Arctic's NPP increased from 2003 to 2016, and reached a maximum of 525.74 mg m −2 d −1 in 2016 [82]. Moreover, they found a strong correlation between the increasing trend of open water area and the growing NPP in Arctic Ocean [80,82].
With the decline of global NPP, the NPP in the Arctic Ocean has shown an upward trend, which has also led to an increase in POC flux. We suspect that this is related to global warming. Global warming and polar ice melting are generally recognized, researchers noted that the extent of Arctic sea ice has rapidly declined in recent years [83,84]. The IPCC Special Report on the Ocean and Cryosphere in a Changing Climate (SROCC: https://www.ipcc.ch/srocc/home/) indicated that the global ocean has warmed unabated since 1970 and the rate of ocean warming has more than doubled since 1993. Over the last decades, global warming has led to widespread shrinking of the cryosphere, and the rates of ice loss from the Greenland and Antarctic ice sheets are increasing (SROCC). NASA (https://www.nasa.gov/) also announced that Arctic sea ice reached its annual minimum on 18 September 2019, the analysis of satellite data by NSIDC and NASA shows that the Arctic sea ice sheet shrank to 4.15 million km 2 , surpassing 2007 and 2016 as the second lowest in satellite observation history. From 1984 to 2019, Arctic sea ice shrank rapidly at a visible rate, which has not slowed down at present. As Arctic sea ice melt, the ocean area in the high-latitude has increased year by year, leading to an increase of NPP and POC flux in the high-latitude ocean. Therefore, it is not surprising that there has been an increase in POC flux of high-latitude ocean.

Conclusions
The ocean plays an important role in regulating the carbon dioxide content in the atmosphere, and the transfer of carbon from the surface water to the deep ocean directly affects the CO 2 flux between the ocean and the atmosphere. Therefore, the monitoring of the global ocean POC flux is of great significance for the study of global carbon cycle. Based on the in situ POC flux data, this paper compared the existing classical models of estimating POC flux using remote sensing satellite data. Then, the optimal e-ratio model and the most appropriate NPP products were selected to calculate the global ocean POC flux from 2003 to 2018. The total POC flux is between 8.5 to 14.3 Gt C per year, which is consistent with these model results from other researchers. Based on the analysis of temporal and spatial variations of POC flux, three conclusions are drawn in this paper: Firstly, in terms of spatial distribution, the Continental Margins with water depth less than 2000 m account for less than 1/7 of the global ocean area, but the POC flux accounts for about 30% of the global POC flux. Therefore, we should combine remote sensing observations, biochemical models, and some effective methods to strengthen the monitoring and evaluation of POC flux in Continental Margins.
Secondly, POC flux decreases from 2003 to 2018, which is consistent with NPP, but there are also some abnormal increases and decreases. This paper found that the occurrence of this mutation is very consistent with the appearance of the El Niño-La Niña phenomenon. Some scholars have analyzed the effect of a single El Niño-La Niña phenomenon on POC flux and confirmed this relationship. Our results further prove that the El Niño-La Niña phenomenon has an important impact on global POC flux over a long period of time.
Thirdly, the melting of polar sea ice caused by global warming enlarges the area of ocean in the high-latitudes, which makes POC flux in high-latitude ocean increase from 2003 to 2018.
In summary, the reduction of POC flux in the ocean will inevitably lead to an increase of CO 2 in the atmosphere, thereby aggravating the greenhouse effect. This paper suggests that measures to deal with greenhouse gases can include two aspects: one is to reduce anthropogenic greenhouse gas emissions, and the other is to protect the marine and terrestrial ecosystems and restore their ability to fix carbon, nitrogen, and other elements. Remote sensing technology can provide effective technical support for these measures, and we will conduct more in-depth research in subsequent studies.
Author Contributions: Z.T. designed the research, F.X. and Z.T. carried out the experiments and prepared the manuscript. All authors contributed to the scientific content, the interpretation of the results, and manuscript revisions.