Remote Sensing Observation of Particulate Organic Carbon in the Pearl River Estuary

River estuaries are connectors of terrestrial and marine ecosystems. Riverine particulate organic carbon (POC) is discharged into oceans after a series of biogeochemical reactions in estuaries. Satellite monitoring of POC will improve our understanding of the carbon dynamics of these water bodies. Based on in situ data from four seasonal survey cruises, we developed an algorithm for estimating POC concentrations in the Pearl River Estuary (PRE). Reflectance ratios, Rrs(678)/Rrs(488) and Rrs(748)/Rrs(412), were set as inputs to calculate POC concentration in the PRE. The algorithm was then applied to MODIS/AQUA data to inverse POC concentrations in the PRE from 2002 to 2014. Additionally, sources, impact factors, and seasonal distributions of POC were also investigated. Phytoplankton contributed more to POC in off-shore waters than that in in-shore waters in autumn and spring, but showed the opposite pattern in winter. Under the influence of freshwater from the Pearl River, underwater topography, tides, winds, etc., the seasonal POC concentrations along a specific section, vertical to water depth gradient, decreased in different seasons. These decreases could be described by exponential functions (y = aebx, b < 0). The distribution of POC concentrations in the PRE resulted from complex physical and biogeochemical processes, which can change spatially and seasonally. OPEN ACCESS Remote Sens. 2015, 7 8684


Introduction
Carbon is found widely in aquatic ecosystems in the form of dissolved organic carbon (DOC), particulate organic carbon (POC), dissolved inorganic carbon (DIC), and particulate inorganic carbon (PIC).In the ocean's euphotic layer, DIC together with nutrients can be used to produce POC through phytoplankton photosynthesis [1].Oceanic POC includes biologically derived debris and autotrophic and heterotrophic microorganisms [1,2].Although POC generally represents only around 4-5 × 10 9 t C of total carbon in marine surface waters, it plays an important role in the global carbon cycle [2].POC is important in regards to its role in sequestering carbon, scavenging associated elements and compounds, and sinking to the bottom as part of "the biological pump" [1][2][3].Moreover, compared with POC stock in the ocean surface, the flux from DIC to POC has been estimated at about 50 × 10 9 t C•yr −1 [2].POC pool variation can reflect substantial changes in the related carbon pools, carbon fluxes, etc.
The concentration of POC in surface waters is affected by biological production, transformation to other carbon pools, downward export, etc. [2,4,5].Temporal and spatial variations in ocean POC concentrations make characterizing it based on data taken from ships or other in situ measurements alone difficult [5].Other technical means are needed.Suspended particles can scatter light back to space, with backscattered light being able to be received by remote sensing sensors.Some researchers have reported on relationships between POC concentration and particle backscattering light [5][6][7], and it has been suggested that POC concentrations in oceanic surface waters can be estimated from ocean color remote sensing data [2,3,[5][6][7][8].However, the spectral variability of the backscattering coefficients of suspended particles increases the difficulty in calculating POC from satellite data.At present, many remote sensing POC algorithms are based on the green band or green-to-blue ratio [3,5,9]; besides, many studies use the decreasing trend of bb with wavelength [5,7,9].In 1999, Stramski et al. [9] firstly published an algorithm for estimating POC concentration in the Southern Ocean from SeaWiFS data based on the correlation between the particle backscattering coefficient, bbp(λ), and POC concentration and the dependence of spectral reflectance, Rrs(λ), on the total backscattering coefficient, bb(λ) [9].Using in situ data from the eastern Atlantic Ocean [5] or field data from cruises across the Atlantic and Pacific [7], further algorithms were developed based on power function fitting between POC concentration and Rrs(490)/Rrs (555) or Rrs(443)/Rrs (555).Correlations between POC concentration and the diffuse attenuation coefficient at 490 nm, K490, or the beam attenuation coefficient due to particles, Cp, were also used to estimate ocean POC concentration from satellite data [8].All these algorithms were designed for retrieval of POC concentrations in open sea surface water or over the euphotic depth [10], where POC is mainly derived from phytoplankton and can even be estimated from Chl-a concentrations [5,11].However, for waters along coastal zones, especially estuaries, POC is derived not only from phytoplankton but also from terrestrial ecosystems [12].The variable particle assemblages (many mineral particles for example) and high concentrations of optically important dissolved organic matter (DOM) in coastal waters mean that open ocean POC algorithms do not perform well [3].Although they account for less than 10% of the world's ocean area, coastal oceans play an important role in the global carbon cycle.On the one hand, approximately 5 × 10 8 t of organic carbon from various terrestrial sources is transported into coastal waters by global rivers every year [13]; on the other hand, non-deltaic shelf deposits account for 45% of organic carbon burial in global oceans, which are thought to be the second largest sink of atmospheric CO2 [14].Son et al. [3] established the Maximum Normalized Difference Carbon Index (MNDCI) to estimate POC concentration in the Gulf of Mexico from remote sensing reflectance.However, the MNDCI does not perform well for highly turbid waters affected by rivers in the East Asian Monsoon region.
The Pearl River Estuary (PRE), which is covered by highly turbid waters, connects the Pearl River and the South China Sea.The POC in the PRE is sourced from both terrestrial and aquatic ecosystems.From November 2013 to August 2014, four seasonal survey cruises were undertaken to determine the sources and dynamics of POC in the PRE.Based on spectrum and POC data from the cruises, we first developed an estimation algorithm for surface POC concentrations in the highly turbid waters of the PRE.We then applied the algorithm to the MODIS/AQUA data to calculate seasonal POC in the PRE; sources, impact factors, and temporal and spatial variations of POC in the PRE were also examined.

Study Area
The PRE is located at 112.93°E-114.29°Eand 21.31°N-22.86°N,ranging from the coastline of the Pearl River Delta to the northern South China Sea (Figure 1).The Pearl River Delta is one of the most prosperous regions in China.Intensive anthropogenic activities, including rapid industrialization and urbanization, have increased organic pollutant discharge into the PRE and eventually to open seas [15].The northern part of Wanshan Archipelago (north of 22° N), also called Lingdingyang Bay, is severely affected by human activities due to a number of large Chinese metropolises are located in this area, including Hong Kong, Shenzhen, and Guangzhou, etc.
The PRE receives freshwater from the Pearl River via eight runoff outlets, which are all located along the western PRE coastline (Figure 1).The eastern side has two deep channels, along which saltwater can intrude into the HuMen outlet at high tide (Figure 1).A counter-clockwise current inside the PRE is present under these conditions [16].The PRE is situated in the subtropical zone, and air temperature above the estuary does not change greatly.However, it is affected by the East Asian Monsoon by, for instance, receiving about 80% of its freshwater from nearby rivers during the wet season [15,17].Ranked according to discharge, the Pearl River is the world's fourteenth largest and China's second largest river [18].Every year, about 9.2 × 10 5 t of dissolved and particulate organic carbon, along with approximately 3.43 × 10 12 m 3 freshwater from the Pearl River, are discharged into the PRE [15,18].Among this, about 4.0 × 10 5 t is dissolved organic carbon (DOC) [15].River water is a main source of chromophoric dissolved organic matter (CDOM) in the PRE, though co-variation between CDOM and DOC is absent [16].

Field Sampling and Laboratory Analyses
Four seasonal cruises were conducted from 18 to 22 November 2013 (autumn cruise), 26 February to 2 March 2014 (winter cruise), 10 to 14 May 2014 (spring cruise), and 24 to 30 August 2014 (summer cruise) (Figure 1).Using Niskin bottles, water samples were collected from 0.5 m below the water surface.The water samples were filtrated immediately through filters to assemble POC, chlorophyll a (Chl-a), and total suspended matter (TSM).Glass fiber filters (Whatman, 0.7 μm, and Φ47 mm) for POC samples had been pre-combusted at 450 °C for 5 h; after filtering, filters with POC were then packed in aluminum foils and stored in a refrigerator with a temperature around −20 °C.Chl-a samples were prepared by filtrating waters through glass fiber filters (Whatman, 0.7 μm, and Φ25 mm) without preprocessing; and filters with Chl-a were packed in aluminum foils and preserved in liquid nitrogen at a temperature around −200 °C.Each TSM sample was assembled by filtrating water through cellulose acetate membrane (Satorius, 0.45 μm, and Φ47 mm), which had been dried at 45 °C for 4 h repeatedly until the difference of weight between the pretreatment and post-treatment was less than 0.1 mg.The filters with TSM were then placed in culture dishes and stored in a refrigerator at −20 °C.The collection of samples followed the protocols of the Joint Global Ocean Flux Study (JGOFS) [19]; more details about Chl-a sampling could also refer to Liu et al. [20].
After returning to the laboratory, the POC samples were measured via the high temperature combustion method (680 °C) using a Total Organic Analyzer (Shimadzu Incorporation) [13].Chl-a samples were measured using a Turner Designs Trilogy TM Fluorometer.Each TSM sample was dried in an oven at 40 °C for 6-8 h repeatedly until the difference of weight between the pretreatment and posttreatment was less than 0.1 mg; the TSM weight was then calculated through the quality difference method [21].Details on Chl-a measurement can be found in Liu et al. [20].

Measurement of Remote Sensing Reflectance
Radiance spectra were measured using a FieldSpec Spectroradiometer (Analytical Spectral Devices Incorporation, ASD, Boulder, CO, USA) with a spectral range of 350 nm to 2390 nm and a resolution of 1 nm.We measured radiance of water (Lws) and radiance of skylight (Lsky) with the probe positioned at an angle of 90°-135° to the plane of incident radiation in order to minimize the effects of sunlight and the ship's shadow.The view angle of the probe for Lws was 30°-45° to the aplomb direction, but 135°-150° for Lsky.After measurements of Lws and Lsky, radiance of the standard gray board (Lp) was measured as similar to Lws.Because of poor sea conditions or occasional factors, spectra were not measured at every designed station.Spectral measurement followed the Ocean Optical Protocols proposed by NASA [21,22].From Lws, Lsky, and Lp, the remote sensing reflectance (Rrs) of water was computed using Equation (1).
where, Lw is the water-leaving radiance; Es(0 + ) is the incident irradiance to the water; r is the reflectance of the air-water interface, r = 0.028 for mixing waters in the PRE; and ρ represents reflectance of the standard gray board, obtained from calibration processes.Figure 2 shows the reflectance spectra ranging from 350 nm to 800 nm, which contain commonly used ocean color bands.Reflectance spectra at stations with high POC concentrations presented absorption peaks around 670 nm and fluorescence peaks around 690 nm.Reflectance spectra of several stations with different POC concentrations measured in the summer cruise are shown in Figure 2.

Figure 2.
Reflectance spectra of water in the PRE.Water-leaving reflectance was calculated using Equation ( 1).Both water-leaving reflectance and POC concentration differed with stations.

Satellite Ocean Color Data
We also obtained daily Level 1B calibrated radiance products (MYD021KM, version 2.16) scanned by the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the AQUA satellite from the Level 1 and Atmosphere Archive and Distribution System (LAADS Web, http://ladsweb.nascom.nasa.gov/).The MODIS/AQUA data have a resolution of 1 km at ocean monitoring bands.All MODIS/AQUA data were observed at about 1:30 pm local time.The near-infrared and shortwave infrared bands (NIR-SWIR) atmospheric correction approach was applied to remove atmospheric effects and calculate water-leaving reflectance from MYD021KM [23].The NIR-SWIR algorithm can derive water-leaving reflectance of open ocean turbid coastal waters at the same time [23].

Algorithm Development
To develop an algorithm for POC concentration estimation, we first calculated equivalent reflectance for eight ocean color bands of MODIS/AQUA data, with central wavelengths at 412 nm, 443 nm, 488 nm, 531 nm, 551 nm, 667 nm, 678 nm, and 748 nm respectively.For a specific band with a central wavelength at λ, equivalent reflectance was computed through Equation (2).
where, requi(λ) is the equivalent reflectance for a band with a central wavelength at λ; f(λ) is the spectral response function, downloaded from OceanColor Website (http://oceancolor.gsfc.nasa.gov/);r(λ) is the reflectance shown in Figure 2; L(λ) is the solar irradiance at mean Earth-Sun distance; λmin is equal to 350 nm, and 800 nm for λmax.
Both particulate and dissolved materials in water influence water-leaving reflectance.An algorithm using the reflectance ratio can decrease the effect of atmospheric correction [5].Many ocean color algorithms have been developed using the reflectance ratio [5,7,21,24].We design our algorithm for POC concentration in the PRE as Equation (3).
In situ data collected in winter and summer cruises, 45 samples in all (red dots in Figure 3a), were used to train the algorithm.To determine the best band combination, correlation coefficient (R 2 ), Mean Square Error (RMSE), Absolute Error (AE), and Relative Error (RE) were employed to evaluate the result of each band combination.Definitions of AE and RE are as per Świrgoń et al. [25] and are defined in Equation ( 4).Based on Equation (3), band tuning and accuracy optimization [26,27] were used to find the optimal bands for λ1, λ2, λ3, and λ4, from the eight ocean color bands of MODIS/AQUA.After testing all possible band combinations, the algorithm achieved a best fitting result with R 2 = 0.865 and AE = 0.0932 mg/L when λ1 = 678 nm, λ2 = 488 nm, λ3 = 748 nm, and λ4 = 412 nm.We also obtained the regression coefficients for the algorithm, with b0 = 0.0078, b1 = 1.3973, and b2 = −1.2397.In other words, surface POC concentrations in the PRE could be estimated through the following equation.
All four bands are commonly used by ocean color researchers: 412 nm, where CDOM has a high absorption coefficient, is often used to inverse CDOM concentration [24]; 488 nm or 748 nm, where suspended particles have a high backscattering coefficient, are often used to inverse TSM concentration [21]; and 678 nm, containing a fluorescence signal from Chl-a, is often used to inverse Chl-a concentration [26,27].

Application to in Situ Data
Parameterized Equation ( 5) was applied to the in situ data collected in autumn and spring cruises, accounting for 26 samples (blue dots in Figure 3a).Figure 3a shows the comparison between in situ POC and modeled POC using the algorithm developed above.For all tested samples, the modeled POC concentrations were greater than the in situ values when in situ POC concentrations were higher than 0.4 mg/L and lower than 0.6 mg/L.Moreover, most modeled POC concentrations were lower than in situ POC concentrations for samples with low in situ POC concentrations (Figure 3a).However, results of applying the developed algorithm to in situ data were acceptable, with a mean AE of 0.103 mg/L and a mean RE of 32.796% for all tested samples.
By matching the in situ POC data with SeaWiFS reflectance, Son et al. [3] developed the MNDCI remote sensing algorithm to estimate water surface POC concentration in the northeastern Mexico Gulf, with almost all POC concentrations used for algorithm development lower than 0.7 mg/L.The MNDCI is shaped by Equation (6).In this paper, 45 samples from the winter and summer PRE cruises were used to compute the coefficients in Equation ( 6).The parameterized equation was then applied to 26 samples from the autumn and spring PRE cruises.Figure 3b shows the comparison between the in situ POC and modeled POC concentrations using the MNDCI.For all tested samples, most modeled POC concentrations were greater than the in situ values for samples with in situ POC concentrations lower than 0.5 mg/L.With a mean AE of 0.118 mg/L and mean RE of 35.81% for all tested samples, the overall estimation accuracy when applying the MNDCI to inverse the PRE POC concentrations was slightly lower than that of the algorithm developed in this paper (Figure 3).It would appear that the algorithm developed by Son et al. [3] was invalid for coastal waters where POC concentrations were greater than 0.8 mg/L.For PRE, all modeled POC concentrations using the MNDCI were lower than 0.7 mg/L when in situ values were higher than 0.8 mg/L (red circles in Figure 3b).MDNDCI where POC represents POC concentration, and requi(412), requi(443), requi(490), and requi(555) are defined in Equation (2).

Application to MODIS/AQUA Data
The new POC algorithm was developed by considering central wavelengths and spectral response functions of ocean color bands of MODIS/AQUA (Section 4).In this section, we applied the new algorithm to estimate POC concentrations in the PRE by inputting reflectance ratios of the MODIS/AQUA data from 2002 to 2014.Table 1 shows the comparisons between in situ POC concentrations and modeled POC concentrations using the newly developed algorithm.Modeled POC in Table 1 denotes the mean value of the 3 × 3 pixels window with central pixel covering the location of in situ data.Due to cloud cover, only nine samples with in situ POC concentrations were found.The nine samples were sampled within four hours prior to or after the passing time of the AQUA satellite.
The mean AE of the modeled POC concentrations from the MODIS/AQUA data was 0.099 mg/L, with a minimum of 0.001 mg/L and maximum of 0.476 mg/L (Table 1), and the mean RE was 19.736%, with a minimum of 0.116% and maximum of 46.825%.When the difference between sampling time and passing time of MODIS/AQUA was large, such as No. 3 in Table 1, AE or RE was large.The RE of the modeled POC concentration was less than 35% for samples collected within three hours prior to or after the passing of the AQUA satellite.On the whole, the developed algorithm could be applied to estimate POC concentration in the PRE based on the MODIS/AQUA data.A good linear relationship was found between the in situ and modeled POC, with an R 2 of 0.855.The Ocean Biology Processing Group (OBPG) recommends using Equation (7) as the standard algorithm to estimate global ocean surface POC concentration [7].For the PRE, however, considerable differences were found in POC using option 2 and option 1 (Equation ( 7)).In addition, POC concentrations obtained using Equation (7) decreased when the POC concentrations in the PRE increased.In other words, the standard OBPG algorithm for POC estimation was not suitable for coastal waters in the PRE.

Particulate Organic Carbon (POC) Distribution from in situ Data
For all 103 samples collected during the four seasonal cruises, the POC concentration ranged from 0.113 mg/L to 1.402 mg/L, with a minimum of 0.113 mg/L in winter and a maximum of 1.402 mg/L in summer (Table 2).In situ data showed there was a wide range of POC concentrations in the PRE (Table 2).All designed stations were covered in the winter and summer cruises (Figure 1).Average POC concentration in the winter cruise was 0.393 mg/L, and in the summer cruise was 0.592 mg/L.Generally, the POC concentration in summer was greater than that from the other three seasonal cruises, especially for stations with salinity lower than 28 PSU (Figure 4).Table 2. Descriptive statistics of in situ POC concentrations (mg/L) in the PRE.Due to rough sea conditions, not all designed stations were covered, especially those located south of 22° N in the autumn and spring cruises (Figure 1).For each cruise, the POC concentrations at different stations varied considerably (Table 2).Figure 4 shows the differences at different stations over the four seasonal cruises.The POC concentrations at stations where salinities were lower than 33 PSU changed more drastically.These stations were closer to the shore and were more affected by freshwater from the Pearl River.However, POC concentrations at stations where salinities were higher than 33 PSU (red circles in Figure 4) were around 0.2 mg/L.However, no clear negative relationship between POC concentration and salinity was observed.

Temporal and Spatial Variations from Satellite-Derived POC
By applying the developed algorithm to the reflectance ratios of MODIS/AQUA, we obtained the distribution of POC concentration in the PRE.We then computed seasonal POC concentration by averaging all daily POC values in four designated months from 2002 to 2014, which included 20 values in January for winter, 16 in April for spring, 24 in July for summer, and 34 in November for autumn.Figure 5 shows the seasonal POC concentrations in the four specific months.Seasonal POC results were used to analyze the seasonality of POC distribution.In Lingdingyang Bay, POC concentrations were higher in the west than those in the east, especially in spring and autumn; however, this distribution pattern was not obvious in summer.In the southern Wanshan Archipelago, POC concentration decreased from the northwest to the southeast, with higher values in near-shore waters; moreover, the area where POC concentrations were larger than 0.3 mg/L was larger in summer than those in the other three seasons.January for winter, April for spring, July for summer, and November for autumn.The white regions are covered by highly turbid waters.Radiance of these regions was beyond the maximum radiance that band 13 or 14 of MODIS/AQUA can record, that is, band saturation occurred.Because of this, we did not obtain valid POC concentrations in the white regions.The pink line shows the vertical to water depth gradient section used to describe spatial variations of POC in the PRE.
To better understand POC distribution in the PRE, we defined a section (pink line in Figure 5) to quantitatively describe temporal and spatial variations of POC concentration in the PRE.This section starts in Lingdingyang Bay (113.61°E, 22.36° N) and ends in southeast of the study area (114.12°E, 21.39° N).The section was perpendicular to the isoline of water depth (Figures 1 and 6).POC concentrations along the section are shown in Figure 6.For all seasons, POC concentrations decreased along the section.POC concentrations decreased more remarkably in northern Wanshan at longitudes smaller than 113.8°E. The decreasing trends observed in all seasons can be described using exponential functions (y = ae bx , b < 0), with R 2 = 0.94 for winter, R 2 = 0.75 for spring, R 2 = 0.96 for summer, and R 2 = 0.96 for autumn.It is worth noting there were POC concentration peaks around the Wanshan Archipelago (~113.8°E) in all seasons except summer.Two regions were also defined in this study to present annual variations in PRE POC concentrations from 2002 to 2014.Region A was located in the Lingdingyang Bay, which was affected obviously by freshwater; and region B was located in the southern Wanshan Archipelago (Figure 1).Annual POC concentrations were calculated by averaging all satellite-derived daily POC values in region A or region B in January, April, July, and November of each year.Annual POC concentrations in region A were much greater than those in region B (Figure 7a).

Sources of POC in the Pearl River Estuary
On a global scale, approximately 200 Tg (Tg = 1 × 10 12 g) of terrestrial POC is poured into marginal oceans by rivers annually [13].POC is one component of TSM, with POC(%TSM) denoting the percentage of POC in TSM [28,29].Previous studies have described negative relationships between POC(%TSM) and TSM concentrations in numerous river systems [30][31][32].We collected POC and TSM data from the world's major rivers and developed a power function (y = 16.191x−0.395 , R 2 = 0.602) to describe the correlation between POC(%TSM) and TSM concentration, which showed that the POC(%TSM) in different rivers differed from each other (Figure 8a).Monsoon rivers, such as the Ganges River, Changjiang River, and Mekong River, were characterized by high POC concentrations but low POC(%TSM).Conversely, non-monsoon rivers, such as the St. Lawrence River, Congo River, and Amurm River, were characterized by low POC concentrations but high POC(%TSM) (Figure 8a).Precipitation and river discharge from monsoon regions occur mainly in the wet season [17,33].High intensity rainfall over monsoon regions in the wet season can lead to severe terrestrial erosion, resulting in an increase in TSM in rivers.When the TSM concentration is high, primary production of phytoplankton can decrease due to light limitation [34].Thus, allochthonous POC becomes the major source of POC, and POC(%TSM) will decrease.This situation shows the opposite pattern in nonmonsoon rivers.Of course, basin characteristics might also influence POC(%TSM), such as what observed in the Nile River and Yellow River (not included).Based on our four seasonal cruises, Figure 8b shows the relationship between POC(%TSM) and TSM concentration in the PRE, as also described through the power function (y = 21.630x−0.637 , R 2 = 0.756).The PRE receives freshwater from the Pearl River.However, the POC(%TSM) in the PRE differed from that in the Pearl River.POC(%TSM) in the PRE ranged from 1.262% to 50.377%, with an average of 11.111%, while annual POC(%TSM) in the Pearl River was only 2% (Figure 8).This considerable difference might be due to additional POC sourced from phytoplankton in the more clear waters of PRE [12].Estuaries are the mixing zone between freshwater and seawater, and thus sources of POC in estuaries are complicated.Terrestrial POC comes from terrestrial ecosystems and is produced primarily through photosynthesis of vegetation, while aquatic POC comes from aquatic ecosystems and is produced primarily through photosynthesis of phytoplankton [12,36,37].Many methods have been used to explore sources of POC in estuaries, including the elemental ratio of carbon to nitrogen (C/N), stable isotope composition (δ 13 C), and stable isotope signature of particulate nitrogen (δ 15 N), etc. [38].Phytoplankton contributed significantly to POC in the PRE.The positive correlation between Chl-a and POC concentration was fitted through exponential function, with R 2 = 0.541 and p < 0.001 (Figure 9a).In accordance with Ni et al. [15], we also found a positive correlation between TSM and POC concentration in the PRE (Figure 9b).Around the outlets of HongQiMen and HEngMen, Liu et al. [12] reported that waters had high TSM and low Chl-a concentration.Therefore, POC concentrations in these waters appear to be majorly affected by POC transported by the Pearl River from land.Using δ 13 C, Liu et al. [12] also investigated the mixed terrestrial and aquatic sources of POC in the PRE and found that in summer, POC in the north and middle region was primarily sourced from terrestrial ecosystems (64%), but POC in the southern region was mainly from aquatic ecosystems (68%) [12].
Although the POC/Chl-a ratio might vary over a broad range in different ocean areas, POC/Chl-a can be used to explore sources of POC [39].The influence of Chl-a on POC shows a constant value when data are constrained to regions with less terrestrial input [5].Stramska et al. [5] estimated the POC concentration from Chl-a in the North Polar Atlantic based on a linear equation.Similar to open seas, the POC concentrations in off-shore areas of the PRE, where waters were deeper than 18 m, were little affected by freshwater from the Pearl River, especially in seasons with small discharge [16].There was a strong linear relationship between Chl-a and POC concentration in these waters (Figure 10).According to Eppley et al. [39], the linear regression slope of POC on Chl-a by weight provides information on the carbon/Chl-a ratio of phytoplankton, while the intercept provides the non-algal POC.POC concentrations in off-shore areas of the PRE can also be estimated from Chl-a.As seen in Figure 10, phytoplankton with 1 μg of Chl-a contained about 0.129 mg of POC in off-shore waters.Based on this, POC = 0.231 mg/L at Chl-a = 1 μg/L in off-shore waters.For regional PRE, we first assumed that phytoplankton communities were similar.We then estimated the phytoplankton contribution to POC in the PRE from the Chl-a data collected from the seasonal cruises (Table 3).In autumn and spring, phytoplankton contributed more to POC in off-shore waters than that in in-shore waters; however, phytoplankton contributed more to POC in in-shore waters than that in off-shore waters in winter.This phenomenon was likely the result of smaller riverine POC flux, together with smaller discharge, from the Pearl River in winter [15,17].We assumed that phytoplankton communities were similar in regional PRE so as to calculate the contribution of phytoplankton to POC semi-quantitatively.However, spatial variability of the phytoplankton community structure might exert an effect on POC/Chl-a variations.Based on only one cruise in April 2007, He et al. [40] also reported that POC and Chl-a had a significant positive correlation in the freshwater region of the PRE (including part of the Pearl River), with the POC/Chl-a regression slope being only 0.0234 by weight.Therefore, the contribution of phytoplankton to POC in in-shore waters might have been overestimated.

Impact Factors on POC Distribution in the PRE
After entering marginal oceans, riverine POC will experience a series of biogeochemical reactions.Part of riverine POC will be decomposed by bacteria, zooplanktons, viruses, etc., [35,41].In addition, different sized POC will be partitioned in estuaries and further partitioned when transported along-shore or across-shelf [14].Large POC particles are prone to settle down to the seafloor; however, tiny POC particles might be transported to the continental shelf, or even open seas, by currents [32].Except for riverine POC, salinity-induced flocculation of dissolved organic carbon (DOC) to POC might increase POC concentration in coastal oceans [41].Moreover, phytoplankton production in some parts of the river plume will be generally high because of rich available nutrients, buoyant freshwater, and adequate light penetration [35].However, the location of maximum phytoplankton production along the salinity gradient might change with seasons and estuary conditions [35].Of course, benthic resuspension, which occurs primarily in winter in the PRE, can also increase POC concentrations in coastal waters [1,12].Thus, the distribution of POC in coastal oceans is strongly influenced by all these major processes.
Seawater can intrude into the PRE from the east, while most freshwater from the Pearl River is discharged into the western PRE.Water masses mixing processes can influence the spatial distribution of POC concentration in the PRE.Seawater with lower POC concentrations can prevent freshwater from spreading to the east (Figure 5).Like in many other estuaries, freshwater, usually with higher POC concentrations, will be diluted by seawater [1,3,42].Therefore, POC concentrations in the west were generally greater than those in the east (Figure 5).In the PRE, POC concentration decreased with distance from the Pearl River outlets (Figures 5 and 6).In situ data in the different seasons also showed that POC concentrations decreased with increasing salinity (Figure 4).Based on more than 10 years of data, the exponential function (y = ae bx , b < 0) also described the decrease in seasonal POC with increasing distance from the coastline (Figure 6).Along the specific section (Figure 5), seasonal POC in different seasons were approximate to each other in the east of 113.95°E. The POC concentration in region B was less affected by freshwater than that in region A, and did not change significantly from 2002 to 2014 (Figure 7).
Influenced by the East Asia Monsoon, discharge from the Pearl River differs greatly with season.About 80% of freshwater is discharged into the PRE during the wet season (April-September), and only 20% occurs during the dry season [17].While freshwater brings terrestrial POC into the PRE [12,15], the additional transported nutrients can also trigger algal growth [35].Under the influences of freshwater, seasonal POC concentrations from 2002 to 2014 along the section (Figure 5) in summer, with higher discharge, were generally greater than those in other seasons (Figures 5 and 6).However, high turbidity in near-shore waters can also prevent phytoplankton from growing due to light limitation [16,43].This might be the reason why POC concentrations in near-shore waters were lower in summer than those in autumn and winter (Figure 6).Furthermore, discharge in summer might also differ greatly with days.For example, discharge from the Pearl River on 6 July, measured at Gaoyao, Shijiao, and Boluo hydrometric stations (http://xxfb.hydroinfo.gov.cn), was greater than that on 31 July [17].Similarly, POC concentrations in the northern PRE on 6 July 2011, were apparently greater than those on 31 July 2011 (Figure 11).Temporal distribution of POC differed greatly with discharge changing in an identical season, with higher POC concentration observed following greater discharge.During the wet season, the plume of the Pearl River flows eastward and offshore under the impact of southwest wind [17,44].Therefore, POC concentrations south of Hong Kong in summer were greater than those in other seasons (Figure 5).However, the plume is advected westward due to northeast wind during the dry season, so the eastward expansion of high POC concentrations is much smaller in winter and spring [17,44].The high POC concentrations around the Wanshan Archipelago might be due to high phytoplankton production in low velocity water or POC discharged from surrounding islands.Other factors that influence phytoplankton growth might also have influenced POC distribution in the PRE, including coastal upwelling [45] and air temperature [1,2,6].

Conclusions
Under the influence of the East Asian Monsoon, most freshwater from the Pearl River occurs during the wet season.The POC (%TSM) in a monsoon river differs greatly to that in non-monsoon river.In this paper, an algorithm was developed to estimate surface POC concentrations in the PRE using remote sensing reflectance ratios.The algorithm was applied to MODIS/AQUA data to monitor POC concentration in the PRE dynamically.Importantly, it could describe high POC concentrations (>0.8 mg/L) in near-shore waters.
In the PRE, phytoplankton contribution to POC differed in different waters and changed with the seasons.Phytoplankton contributed more to POC in in-shore waters than that in off-shore waters in winter, but showed the inverse pattern in autumn and spring.Water masses mixing influenced the distribution of POC in the PRE.The temporal distribution of POC in the PRE was negatively affected by discharge from the Pearl River.Even with more POC from the Pearl River, POC concentration was lower at high flow on account of light limitation to phytoplankton growth.Moreover, POC concentrations in off-shore waters did not change significantly in the past decade; but POC concentrations in in-shore waters fluctuated with discharge from the Pearl River.Considering spatial characteristics, decreasing seasonal POC along a specific vertical to water depth gradient section could be fitted by exponential function (y = ae bx , b < 0) in all seasons.
Rich nutrients brought by river might result in phytoplankton flourish and increase POC concentration in estuary.Terrestrial and aquatic POCs in estuary experience a series of complicated physical and biogeochemical reactions before being transported far away.Reflectance ratios of the developed algorithm use common-used bands of ocean color sensors, such as geostationary ocean color imager (GOCI), medium-resolution imaging spectrometer (MERIS), visible infrared imaging radiometer suite (VIIRS), etc.Based on these global covering satellite data, we can monitor POC dynamic changes in highly turbid waters of worldwide estuaries using the developed algorithm.Furthermore, POC concentration calculated from satellite data and water flux estimated through numerical model simulation can be combined to compute the POC flux from coastal ocean to open sea.All in all, the developed algorithm can be applied to explore sources, distributions, and dynamic changes of POC in estuaries using ocean color satellite data.Based on this POC information, we can explore POC budgets in global estuaries and improve management policies.

Figure 3 .
Figure 3. Comparisons between in situ POC and modeled POC concentrations using the algorithm developed in this study (a) and by Son et al. [3] (b). AE and RE are defined in Equation (4).

Figure 4 .
Figure 4. Changes in POC concentrations with salinity (Practical Salinity Units, PSU) in different cruises.Salinity was recorded by a multi-parameter water quality monitor (SBE 19Plus, Sea-Bird Electronics, Inc.: Bellevue, WA, USA).

Figure 5 .
Figure 5. Seasonal POC concentrations (mg/L) derived from the MODIS/AQUA data from 2002 to 2014 using the proposed algorithm.Each season is represented by a specific month:January for winter, April for spring, July for summer, and November for autumn.The white regions are covered by highly turbid waters.Radiance of these regions was beyond the maximum radiance that band 13 or 14 of MODIS/AQUA can record, that is, band saturation occurred.Because of this, we did not obtain valid POC concentrations in the white regions.The pink line shows the vertical to water depth gradient section used to describe spatial variations of POC in the PRE.

Figure 6 .
Figure 6.Seasonal POC concentrations (mg/L) along the section defined in Figure 5.The Wanshan Archipelago is located around 113.8° E. Each season is presented by a specific month: January for winter, April for spring, July for summer, and November for autumn.

Figure 7 .
Figure 7. Annual satellite-derived POC in the PRE in different years: (a) Annual POC concentrations in region A and B, which are defined in Figure 1; (b) relationship between annual POC in the PRE and discharge of the Pearl River.Pearl River discharge was calculated from the daily discharge of three control hydrometric stations, Boluo, Shijiao, and Gaoyao.

Figure 9 .
Figure 9. Relationship between POC concentration in the PRE and ocean color parameters: (a) Chl-a; (b) TSM.Chl-a data are not available for the summer cruise.

Figure 10 .
Figure10.Correlation between POC concentrations and Chl-a concentrations for stations located in off-shore areas, where waters were deeper than 18 m (Figure1).Summer data are not included.

Figure 11 .
Figure 11.Satellite-derived POC concentrations in summer with different discharge from the Pearl River on (a) 6 July 2011; (b) 31 July 2011.Valid POC concentrations in the white regions were not obtained due to cloud cover, contamination by sun-glint, band saturation of MODIS, etc.

Table 1 .
Comparison results between in situ and modeled POC concentrations using the algorithm developed in this study.AE and RE are defined in Equation (4).Sampling time is the local time when MODIS/AQUA was passing. No.

Table 3 .
Contribution of phytoplankton to POC in different seasons and water types.Chl-a data was not available for the summer cruise.Off-shore waters contained stations with waters deeper than 18 m.