Estimation of the Particulate Organic Carbon to Chlorophyll- a Ratio Using MODIS-Aqua in the East / Japan Sea, South Korea

: In recent years, the change of marine environment due to climate change and declining primary productivity have been big concerns in the East / Japan Sea, Korea. However, the main causes for the recent changes are still not revealed clearly. The particulate organic carbon (POC) to chlorophyll- a (chl- a ) ratio (POC:chl- a ) could be a useful indicator for ecological and physiological conditions of phytoplankton communities and thus help us to understand the recent reduction of primary productivity in the East / Japan Sea. To derive the POC in the East / Japan Sea from a satellite dataset, the new regional POC algorithm was empirically derived with in-situ measured POC concentrations. A strong positive linear relationship (R 2 = 0.6579) was observed between the estimated and in-situ measured POC concentrations. Our new POC algorithm proved a better performance in the East / Japan Sea compared to the previous one for the global ocean. Based on the new algorithm, long-term POC:chl- a ratios were obtained in the entire East / Japan Sea from 2003 to 2018. The POC:chl- a showed a strong seasonal variability in the East / Japan Sea. The spring and fall blooms of phytoplankton mainly driven by the growth of large diatoms seem to be a major factor for the seasonal variability in the POC:chl- a . Our new regional POC algorithm modiﬁed for the East / Japan Sea could potentially contribute to long-term monitoring for the climate-associated ecosystem changes in the East / Japan Sea. Although the new regional POC algorithm shows a good correspondence with in-situ observed POC concentrations, the algorithm should be further improved with continuous ﬁeld surveys.


Introduction
The East/Japan Sea is a semi-marginal sea located in the northwestern Pacific, and it is bordered by Korea, Japan, and Russia. Many previous studies reported that the East/Japan Sea is a productive region, especially the Ulleung Basin located in the southwestern East/Japan Sea [1][2][3][4][5][6][7]. Recently, not only the changes in environmental conditions but also alterations in biological characteristics were reported in many previous studies [4][5][6]8]. Especially, Joo et al. [5] addressed a significant declining trend of the annual primary production in the Ulleung Basin which is the biological hotspot in the East/Japan Sea. However, although many other studies have observed the changes in the East/Japan Sea, the main causes driving the recent changes remain unclear.
The particulate organic carbon (POC) to chlorophyll-a (chl-a) ratio (POC:chl-a) has been used as an important indicator for ecological and physiological state of phytoplankton. Generally, the concentrations of photosynthetic pigments in phytoplankton cells depend on the environmental factors and the physiological conditions of phytoplankton [9][10][11][12][13][14]. For instance, the POC:chl-a can be increased under high light intensity and low nitrogen supply conditions [1,5,6]. Moreover, the POC:chl-a varies according to the size structure of the phytoplankton community [15][16][17]. Therefore, understanding of spatial and temporal variations in the POC:chl-a can help us to determine the ecological and physiological conditions of a phytoplankton community. However, the existing POC-deriving algorithm using satellite ocean color data [18] is not validated in the East/Japan Sea. Although the algorithm is fully validated in global ocean, it is hard to directly apply to this small regional sea. Thus, the algorithm needs to be calibrated and the accuracy should be evaluated before applying to our study area.
Therefore, this study aims to (1) develop a new regional POC algorithm using an ocean color satellite and (2) investigate spatiotemporal variability of the POC:chl-a in the East/Japan Sea.

Study Area and Sampling
Total chl-a and POC concentrations were measured at 41 stations in the East/Japan Sea (16 stations in 2012, 11 stations in 2013, six stations in 2014, and eight stations in 2015, respectively; Figure 1). Northern and southern regions of the East/Japan Sea were defined as shown in Figure 1 to investigate spatial variations of POC:chl-a. Sampling and analysis of total chl-a and POC concentrations were conducted based on Lee et al. [19]. Water samples for the total chl-a and POC concentrations of phytoplankton were collected from the surface layer using a rosette sampler with Niskin bottles. The collected water samples were immediately filtered on Whatman ® glass microfiber filters (precombusted; Grade GF/F, diameter = 24 mm) using a vacuum pressure lower than 5 in. Hg. The filtered samples were frozen immediately and preserved until analysis at the laboratory. The chl-a concentrations were measured with a precalibrated fluorometer (10-AU, Turner Designs) after extraction in 90% acetone in a freezer for 24 h based on Parsons et al. [20]. The filters for POC concentrations were frozen immediately and preserved for mass spectrometric analysis at the Alaska Stable Isotope Facility of the University of Alaska Fairbanks, USA.

Satellite Datasets
We obtained the chl-a and remote sensing reflectance (Rrs) data from the MODIS (Moderate Resolution Imaging Spectroradiometer) onboard the satellite Aqua provided by the OBPG (Ocean Biology Processing Group at NASA Goddard Space Flight Center; https://oceandata.sci.gsfc.nasa.gov/ MODIS-Aqua/). We used the Level-3 daily composite datasets covering the East/Japan Sea from July 2002 to December 2018 at 4-km of spatial resolution [21,22].
The POC derivation algorithm for the East/Japan Sea was empirically derived with our in-situ observation data. Based on the Stramski et al. [18], a power-law relationship between a blue-to-green band ratio of Rrs and POC were used to estimate POC concentrations. The equation for the algorithm is expressed below: where a and b are constants. Constants were determined empirically by regression analysis with our in-situ dataset. The input wavelength for the green band can be replaced with available band between 547 and 565 nm. In this study, Rrs(443) and Rrs(547) which are available bands of MODIS were used as input wavelengths for blue-to-green band ratio for the POC algorithm.
POC:chl-a ratios were calculated by dividing our estimated POC concentrations by remotely sensed chl-a concentrations. Monthly composited data for POC and POC:chl-a were obtained by averaging daily data for each month.

POC Algorithm Derivation
In-situ measured POC concentrations ranged from 84.07 to 713.69 mg m −3 , and the average was 262.93 ± 205.82 mg m −3 . The blue-to-green band Rrs ratio were extracted from MODIS-Aqua monthly composite datasets. To determine two constants in the POC algorithm, a and b, the curve fitting using nonlinear regression between Rrs ratio and in-situ POC concentration was conducted ( Figure 2a).  From this result, the constants a and b were determined as 295.7 and -1.028, respectively. Based on these two constants, the POC algorithm was derived as follows: The determination coefficient (R 2 ) and Spearman's correlation coefficient for the relationship From this result, the constants a and b were determined as 295.7 and −1.028, respectively. Based on these two constants, the POC algorithm was derived as follows: The determination coefficient (R 2 ) and Spearman's correlation coefficient for the relationship between Rrs ratio and in-situ POC concentrations were 0.8017 and −0.861, respectively (Figure 2a). The POC concentrations derived from the regional model also showed a strong linear relationship with in-situ measured POC concentrations (R 2 = 0.6579), and most of the satellite-derived POC were plotted within 95% prediction bounds (Figure 2b).
The new POC algorithm showed lower RMSE and bias in comparison to the existing algorithm [18] in the East/Japan Sea (Figure 3). The RMSE and bias of the new algorithm were 115.37 and −17.43, respectively, and those of Stramski et al. [18] were 161.47 and −97.29, respectively. From this result, the constants a and b were determined as 295.7 and -1.028, respectively. Based on these two constants, the POC algorithm was derived as follows: The determination coefficient (R 2 ) and Spearman's correlation coefficient for the relationship between Rrs ratio and in-situ POC concentrations were 0.8017 and -0.861, respectively (Figure 2a). The POC concentrations derived from the regional model also showed a strong linear relationship with in-situ measured POC concentrations (R 2 = 0.6579), and most of the satellite-derived POC were plotted within 95% prediction bounds (Figure 2b).
The new POC algorithm showed lower RMSE and bias in comparison to the existing algorithm [18] in the East/Japan Sea (Figure 3). The RMSE and bias of the new algorithm were 115.37 and -17.43, respectively, and those of Stramski et al. [18] were 161.47 and -97.29, respectively. The POC concentrations were derived for the East/Japan Sea from our new regional algorithm ( Figure 4). The climatological monthly mean POC (January 2003-December 2018) showed a seasonal variation of the POC concentration in the East/Japan Sea (Figure 4). Generally, the POC concentrations were relatively higher during spring season, and the lowest POC concentrations were observed during summer season.
The POC concentrations were derived for the East/Japan Sea from our new regional algorithm (Figure 4). The climatological monthly mean POC (January 2003-December 2018) showed a seasonal variation of the POC concentration in the East/Japan Sea (Figure 4). Generally, the POC concentrations were relatively higher during spring season, and the lowest POC concentrations were observed during summer season.

POC:chl-a
The climatological monthly distribution of the POC:chl-a (January 2003-December 2018) in the East/Japan Sea showed strong seasonal variations ( Figure 5). In contrast to the seasonal pattern of POC, relatively lower POC:chl-a ratios were observed during spring and autumn compared to those during winter and summer.

POC:chl-a
The climatological monthly distribution of the POC:chl-a (January 2003-December 2018) in the East/Japan Sea showed strong seasonal variations ( Figure 5). In contrast to the seasonal pattern of POC, relatively lower POC:chl-a ratios were observed during spring and autumn compared to those during winter and summer.   The ranges of the mean POC:chl-a in the northern and southern East/Japan Sea were 169.6-528.4 and 172.4-549.3, respectively ( Figure 6). Domains for the two regions are shown in Figure 1. The average of POC:chl-a in the northern and southern regions of the East/Japan Sea were 377.6 ± 80.4 and 388.1 ± 69.7, respectively. No statistically significant difference in POC:chl-a was observed between the northern and southern regions (t-test, p > 0.05).
. The POC:chl-a showed a strong seasonal variation (Figure 7). Climatological monthly mean POC:chl-a showed the lowest values during April (270.3 ± 74.7 and 261.9 ± 60.9 for the northern and southern East/Japan Sea, respectively), and the highest values were observed during August (461.0 ± 18.1) and July (451.5 ± 23.6) for the northern and southern East/Japan Sea, respectively (Figure 7). The The POC:chl-a showed a strong seasonal variation (Figure 7). Climatological monthly mean POC:chl-a showed the lowest values during April (270.3 ± 74.7 and 261.9 ± 60.9 for the northern and southern East/Japan Sea, respectively), and the highest values were observed during August (461.0 ± 18.1) and July (451.5 ± 23.6) for the northern and southern East/Japan Sea, respectively (Figure 7). The POC:chl-a during spring and autumn were significantly lower than summer and winter (t-test, p < 0.01).
Remote Sens. 2020, 12, x FOR PEER REVIEW 7 of 10 POC:chl-a during spring and autumn were significantly lower than summer and winter (t-test, p < 0.01).

New regional POC algorithm
Since the POC algorithm reported by Stramski et al. [18] was developed for the eastern South Pacific and eastern Atlantic Oceans, we validated and modified the algorithm to derive a suitable POC model for the East/Japan Sea. The previously reported POC algorithm by Stramski et al. [18] tends to underestimate the POC concentrations in the East/Japan Sea (Figure 3). However, our modified regional POC algorithm showed significantly improved accuracy in the East/Japan Sea, and the POC concentrations estimated from the regional algorithm derived in this study showed a strong

New Regional POC Algorithm
Since the POC algorithm reported by Stramski et al. [18] was developed for the eastern South Pacific and eastern Atlantic Oceans, we validated and modified the algorithm to derive a suitable POC Remote Sens. 2020, 12, 840 7 of 10 model for the East/Japan Sea. The previously reported POC algorithm by Stramski et al. [18] tends to underestimate the POC concentrations in the East/Japan Sea (Figure 3). However, our modified regional POC algorithm showed significantly improved accuracy in the East/Japan Sea, and the POC concentrations estimated from the regional algorithm derived in this study showed a strong linear relationship with in-situ measured POC concentrations and its linear regression line located near the 1:1 line (Figure 3).
However, only 41 data points were used to derive the POC model in this study, and in fact, the number of data points could not be considered sufficient. Nevertheless, the modified POC model showed a good correspondence with field measured POC concentrations. If continuous field observations and the calibration and validation of the algorithm are conducted, the POC algorithm with a better performance would be derived.

Spatial and Temporal Variability of the POC:chl-a
The seasonal variation of the POC:chl-a in the East/Japan Sea would be closely related with the physiological conditions of phytoplankton communities. Generally, phytoplankton tend to accumulate more carbon in their cells under the high light intensity and nutrient depleted conditions [9,13,14,23,24]. Phytoplankton increase their chl-a contents to maximize light absorption under the low light condition [11,25,26]. Moreover, many previous studies reported that the carbon to chl-a ratio of phytoplankton decreases with increased growth rate [14,[27][28][29][30]. The size structure of the phytoplankton community also can affect the carbon to chl-a ratio [15][16][17]. Small cell-sized phytoplankton such as flagellates usually show higher carbon to chl-a ratio than large cell-sized phytoplankton such as diatoms [15][16][17]. However, the daily carbon uptake rate by phytoplankton tends to be lowered when the productivity contribution of picoplankton to the total primary production is high [31,32]. It suggests that the investigation of POC:chl-a in the East/Japan Sea can provide some potential clues on the recent changes of primary productivity in the East/Japan Sea.
In general, there are two phytoplankton blooms per year in the East/Japan Sea; spring bloom and fall bloom [4,5,[33][34][35]. Signals of the blooms were also observed in the climatological monthly mean distribution of the POC (Figure 3). The strongest bloom occurs during spring season while the weaker blooms appears in fall. Both spring and fall blooms are mainly caused by the massive growth of diatoms [3,4,33,35,36]. During spring and fall blooms, the size distribution of phytoplankton cells would shift from smaller to larger size due to the rapid growth of diatoms. Consequently, lower POC:chl-a during spring and fall might be caused by the bloom mainly driven by the growth of diatoms. However, these suggestions are only our hypothesis based on several previous studies. To understand the temporal variation of the POC:chl-a in the East/Japan Sea, further research with field observations is needed.
On the other hand, the timing of spring bloom showed a difference in the northern and southern regions of the East/Japan Sea (Figures 3 and 4). POC concentrations were highest in April in the southern region and highest in May in the northern region. Additionally, the distribution of POC:chl-a during April and May appeared to be the opposite of POC concentrations. Other previous studies have also observed a similar spatial distribution of chl-a in the East/Japan Sea with satellite datasets [37,38]. Those spatiotemporal distributions of POC concentrations and POC:chl-a suggest that the spring bloom occurs earlier in the southern regions of the East/Japan Sea.

Summary and Conclusions
In this study, the regional POC algorithm for the East/Japan Sea was derived empirically using in-situ measured POC and MODIS-Aqua satellite datasets. In-situ measured POC concentrations at the 41 stations in the East/Japan Sea were used to calibrate and validate our new regional POC algorithm. The power-law relationship between POC concentration and blue-to-green band of remote sensing reflectance, Rrs(443)/Rrs(547), was used to derive the algorithm based on Stramski et al. [18] (Figure 2a).