Spatiotemporal Variation in Particulate Organic Carbon Based on Long-Term MODIS Observations in Taihu Lake , China

In situ measured values of particulate organic carbon (POC) in Taihu Lake and remote sensing reflectance observed by three satellite courses from 2014 to 2015 were used to develop an near infrared-red (NIR-Red) empirical algorithm of POC for the Moderate Resolution Imaging Spectroradiometer (MODIS-Aqua) satellite image. The performance of the POC algorithm is highly consistent with the in situ measured POC, with root mean square error percentage (RMSPs) of 38.9% and 31.5% for two independent validations, respectively. The MODIS-derived POC also shows an acceptable result, with RMSPs of 53.6% and 61.0% for two periods of match-up data. POC from 2005 to 2007 is much higher than it is from 2002 to 2004 and 2008 to 2013, due to a large area of algal bloom. Riverine flux is an important source of POC in Taihu Lake, especially in the lake’s bank and bays. The influence of a terrigenous source of POC can reach the center lake during periods of heavy precipitation. Sediment resuspension is also a source of POC in the lake due to the area’s high dynamic ratio (25.4) and wind speed. The source of POC in an inland shallow lake is particularly complex, and additional research on POC is needed to more clearly reveal its variation in inland water.


Introduction
Inland freshwater has an important role in the regional and global carbon cycle despite its small fraction of the earth's surface relative to the land and ocean [1][2][3][4][5][6][7].In the global carbon cycle, inland freshwater is not only a pivotal point for carbon transport between terrestrial and marine (about 0.9 Pg C y −1 ) locations, but it also serves to bury terrestrially derived carbon (about 0.98 Pg C y −1 ).This buried carbon will be emitted to the atmosphere as CO 2 (0.75 Pg C y −1 ) during the mineralization process of carbon and finally stored in the sediment (0.23 Pg C y −1 ) [4,8].Particulate organic carbon (POC) in inland waters is an important carbon pool that has a clear effect on dissolved organic and inorganic carbon due to the degradation and mineralization [2,9,10].Little change in the POC of inland water will indicate that there are substantial changes in related carbon fluxes, even though POC represents only a small portion of the total carbon reservoir on a global scale [11][12][13][14].A detailed understanding of the distribution and long-term patterns of POC in inland water is therefore important to understanding the carbon cycle in inland water and its effect on the global carbon cycle [15][16][17][18].
Large POC-retrieval algorithms for application in remote sensing have been developed because they have the advantage of providing synchronous and timely observations over a large area [12,14,[19][20][21][22][23][24][25][26].Some of these POC algorithms, which are based on the relationships between POC and optical properties (backscattering coefficient, beam attenuation and diffuse attenuation coefficient), performed very well in open ocean waters [19,20,22,[27][28][29][30].An empirical algorithm based on the blue and green bands was also used to estimate POC and showed high accuracy in oceans [12,25,31].However, these POC algorithms fail in inland waters due to the much more complex optical properties [32].This significantly limits the gathering of large-scale and long-term properties of POC distribution over time in inland water bodies.Several new POC retrieval algorithms focusing on inland turbid and eutrophic water were developed and applied to remote sensing satellite data [32][33][34].Most of these POC retrieval algorithms for inland water are based on the near infrared and red bands due to the low impact of tripton and colored dissolved organic matter [35].Nevertheless, long-term observation of POC in inland water is still scarce due to the limitations of the validated atmospheric correction method and POC algorithms for specific satellite images (such as Moderate-resolution Imaging Spectroradiometer (MODIS)) [33].It is thus difficult to evaluate the historical levels and changing trends of POC in inland water.The validated POC algorithm and atmospheric correction method for satellite images, which have a long-term observation dataset, are useful for assessing of the carbon status of inland waters from remote sensing [32][33][34].The amount of POC data for inland water is still insufficient for the development of validated POC algorithms, although several studies have described the large-scale distributions and short-term variations of POC in the inland water [32,33,36].Thus, much is still not well understood about POC in inland water, and field measurements and POC-retrieval algorithms are needed to capture its changing trends.
The storage of POC in Taihu Lake is highly variable because of the strong influence of terrestrial carbon input, sedimentary releases from sediment resuspension and local biological products.This study develops an NIR-Red band ratio empirical algorithm of POC, based on the in situ POC measurement data and remote sensing reflectance (R rs ) in Taihu Lake.This empirical algorithm was then applied to satellite data from MODIS to obtain the long-term distribution of POC.The primary objective is to reveal long-term spatial and temporal variation in the POC in Taihu Lake during the period 2002-2014, and to evaluate the factors influencing this variation.

Study Area and Sampling
Taihu Lake is a eutrophic lake with the third largest lake area (2338 km 2 ) in China; the lake is located downstream of the Yangtze River.The lake faces crucial issues of eutrophication and the consequent algal bloom [37,38].The concentration of Dissolved Organic Carbon (DOC) in the lake can reach 7.21 mg/L with mean value of 4.51 ± 0.89 mg/L (mean value ± standard deviation) and the concentration of POC can reach 18.07 mg/L with mean value of 5.45 ± 6.28 mg/L from in situ measurements [38].The optical-property and water-constituent measurements were made during five cruises over a period from 2014 to 2015 (namely, 5 August 2014, 24 October 2014, 21 May 2015, 10 June 2015 and 15 July 2015).The coverage of sampling stations was evenly distributed across the entire lake (Figure 1).The total number of sampling stations with POC data was 197 (this study dataset), and 143 samples were randomly selected to calibrate the POC-retrieval algorithm; the remaining 54 samples (validation data) were used to validate the POC algorithm.The match-up points for the validation of atmospheric correction were 129 during 2006 to 2014 from Huang [39].

Optical Property Measurements
Remote sensing reflectance (Rrs) was obtained by measuring the radiance of water, sky and gray panels.The detailed steps and calculation method follows methods described by Huang [39] and Lee [40].The absorption coefficient of pigment phytoplankton (aph), no-pigment matter (aNAP) and chromophoric dissolved organic matter (aCDOM) were measured using a quantitative membrane filter technique with spectrophotometer [41].

Water Constituent Measurement
Water samples were collected from surface water (30 cm) and stored in Niskin bottles in ice.The water was filtered through pre-combusted (450 °C for 4 h) Whatman glass fiber filters (GF/F) under a low vacuum.A filter with suspended particles was dried at 55 °C for 8 hours, and then treated with 0.01 mol/L HCl to remove inorganic carbon.After acidification treatment, the filter was dried again and weighed.POC concentration (mg/L) was determined by the combustion of the filter in a carbon elemental analyzer (TOC-L-SSM-5000A, Shimadzu Corporation) [42].
The concentration of chlorophyll-a (Chl-a) was measured by a spectrophotometric method with hot alcohol (90%) [41].The concentration of phycocyanin (PC) was determined by spectrofluorimetry [43].The concentrations of total suspended particulate matter (SPM), organic suspended matter (OSM), and inorganic suspended matter (ISM) were detected following the gravimetric method [41].

MODIS Data Processing
Medium-resolution satellite data (with 250 m spatial resolution) from MODIS Aqua Level 0 images from June 2002 to 2014, with no or minimal cloud cover, were chosen and downloaded from the US NASA Goddard Space Flight Center (GSFC, http://oceancolor.gsfc.nasa.gov).SeaDAS 6.4 was used to produce the radiometric and geometric correction of acquired MODIS Aqua Level 0 data.The corresponding MODIS cloud mask products (MOD35) were used to perform the cloud masking.
Satellite angle information (solar zenith angle, solar azimuth angle, satellite zenith angle, and satellite azimuth angle), ozone, and water vapor content from SeaDAS 6.4 were used as input parameters in the atmospheric correction.The atmospheric correction method was based on the hypothesis that water reflectance at near-infrared band is zero [44] and may generate large errors in inland turbid water due to the high reflectance at near-infrared wavelength [45,46].An improved land target-based atmospheric correction method [47,48] was selected to derive the reflectance from

Optical Property Measurements
Remote sensing reflectance (R rs ) was obtained by measuring the radiance of water, sky and gray panels.The detailed steps and calculation method follows methods described by Huang [39] and Lee [40].The absorption coefficient of pigment phytoplankton (a ph ), no-pigment matter (a NAP ) and chromophoric dissolved organic matter (a CDOM ) were measured using a quantitative membrane filter technique with spectrophotometer [41].

Water Constituent Measurement
Water samples were collected from surface water (30 cm) and stored in Niskin bottles in ice.The water was filtered through pre-combusted (450 • C for 4 h) Whatman glass fiber filters (GF/F) under a low vacuum.A filter with suspended particles was dried at 55 • C for 8 hours, and then treated with 0.01 mol/L HCl to remove inorganic carbon.After acidification treatment, the filter was dried again and weighed.POC concentration (mg/L) was determined by the combustion of the filter in a carbon elemental analyzer (TOC-L-SSM-5000A, Shimadzu Corporation) [42].
The concentration of chlorophyll-a (Chl-a) was measured by a spectrophotometric method with hot alcohol (90%) [41].The concentration of phycocyanin (PC) was determined by spectrofluorimetry [43].The concentrations of total suspended particulate matter (SPM), organic suspended matter (OSM), and inorganic suspended matter (ISM) were detected following the gravimetric method [41].

MODIS Data Processing
Medium-resolution satellite data (with 250 m spatial resolution) from MODIS Aqua Level 0 images from June 2002 to 2014, with no or minimal cloud cover, were chosen and downloaded from the US NASA Goddard Space Flight Center (GSFC, http://oceancolor.gsfc.nasa.gov).SeaDAS 6.4 was used to produce the radiometric and geometric correction of acquired MODIS Aqua Level 0 data.The corresponding MODIS cloud mask products (MOD35) were used to perform the cloud masking.
Satellite angle information (solar zenith angle, solar azimuth angle, satellite zenith angle, and satellite azimuth angle), ozone, and water vapor content from SeaDAS 6.4 were used as input parameters in the atmospheric correction.The atmospheric correction method was based on the hypothesis that water reflectance at near-infrared band is zero [44] and may generate large errors in inland turbid water due to the high reflectance at near-infrared wavelength [45,46].An improved land target-based atmospheric correction method [47,48] was selected to derive the reflectance from the MODIS-Aqua data over Taihu Lake.This method involves three simple steps: first, the selection of dense dark vegetation (DDV) pixels from lakeshore and island areas; second, the retrieval of the aerosol optical thickness at 550 nm (AOT 550 ) for the selected DDV pixels; and third, the retrieval of the remote sensing reflectance over water.The DDV pixels were selected from 5-km buffer along the lakeshore and island of Taihu Lake.The mixed pixels were eliminated via a checking process with 3 × 3 uniformity windows.The DDV pixel was reserved only when all nine pixels (3 × 3 windows) were preliminary DDV pixels.DDV pixels with NDVI were chosen to derive AOT 550 .AOT 550 for each pixel in Taihu Lake was from spatial interpolation of retrieved AOT 550 for DDV pixels.For detailed process steps can refer to Liu [48].

Statistical Analysis and Accuracy
Data were analyzed to calculate the spatial and temporal average value of POC for different lake segments and times (monthly), and the temporal average value of POC (seasonal and annual) for the whole lake.We compiled an Interactive Data Language (IDL) program to complete the statistical analysis and map the box figures.The accuracy of the retrieval algorithm was evaluated by root mean square error percentage (RMSP), which expresses the normalized difference between measured and retrieved value to measured value.It can be calculated using Equation (1).
where C POC-r and C POC-m are the retrieved and measured concentrations of POC, respectively.

Calibration and Validation of NIR-Red Empirical Algorithm of POC
Our POC-retrieval algorithm relies on the empirical relationship between surface POC concentration and remote sensing reflectance.The reflectance at a specific wavelength was determined using an optimized iterative algorithm between POC and band ratio of remote sensing reflectance [39,49].The optimal spectral bands (λ 2 and λ 1 ) were determined by the iteration of the MODIS band's setting.Finally, the NIR-Red band ratio (R rs (859)/R rs (645)) empirical algorithm for POC was confirmed.The calibration result shows that the exponential growth function fits the co-variation between POC and R rs (859)/R rs (645) well (Figure 2A).The NIR-Red band ratio in the empirical POC algorithm is very similar to the current algorithm of chlorophyll-a concentration in inland water (or turbid case II water) [39,[49][50][51][52].This also occurs in the ocean, where the blue-green band ratio in the POC algorithm is very similar to the current algorithm of ocean chlorophyll-a concentration [53].This could be caused by any number of POC-containing particles, which include both pigmented phytoplankton and organic detritus, and have higher absorption in the blue spectral region.However, in inland water, the POC-containing particles, not just particles of phytoplankton, but also total suspended particulate matter, have specific fluorescence and scattering properties in the red and near-infrared spectral region [52,54,55].All calibration data in the POC retrieval algorithm were in the 95% confidence line of the linear regression, and the best regression equation in the 95% prediction line can be expressed as Equation (2): N = 143, R 2 = 0.73, p < 0.0001 The validation result for validation data indicates that the POC-retrieval algorithm can measure POC successfully with relatively high accuracy.The RMSP of validation data is 38.9% (Figure 2B).However, this model will underestimate extremely high POC (>30 mg/L), and the underestimated error will increase with high POC.
Remote Sens. 2017, 9, 624 5 of 19 To further verify the performance of the POC algorithm to MODIS satellite data, the match-up points were used to validate the atmospheric correction (match-up points for the Rrs) and POC-retrieval results (match-up points for the Rrs and POC).The match-up MODIS satellite image was selected with the condition that the time interval for matching MODIS and in situ Rrs was less than three hours (similar to Cédric et al. [56]).The match-up pixel is consistent with the coordinates of the sampling site.The performance of this improved land target-based atmospheric correction method is acceptable for turbid water or eutrophic water, aside from some blue and green bands (Figure 3A).The atmospheric correction results indicate that the accuracy of Rrs (645) retrieved by MODIS is much better than that of Rrs (859).The RMSP of Rrs (645) and Rrs (859) are 22.6% and 40.2%, respectively.The RMSP of band ratio of Rrs (645) and Rrs (859) is 29.6%.The accuracy of MODIS-derived Rrs (859)/Rrs (645) is better than that of Rrs (859), but worse than that of Rrs (645) (Figure 3B).The validation of match-up points for the POC-retrieval algorithm shows that the MODIS-derived POC is highly consistent with the in situ measured data.The RMSP between MODIS-derived POC and POC measured on 5 August 2014 and 24 October 2014 is 53.6% and 61.0%, respectively (Figure 3C).The validation of the atmospheric correction and POC-retrieval algorithm shows the applicability of the algorithm to MODIS satellite data.To further verify the performance of the POC algorithm to MODIS satellite data, the match-up points were used to validate the atmospheric correction (match-up points for the R rs ) and POC-retrieval results (match-up points for the R rs and POC).The match-up MODIS satellite image was selected with the condition that the time interval for matching MODIS and in situ R rs was less than three hours (similar to Cédric et al. [56]).The match-up pixel is consistent with the coordinates of the sampling site.The performance of this improved land target-based atmospheric correction method is acceptable for turbid water or eutrophic water, aside from some blue and green bands (Figure 3A).The atmospheric correction results indicate that the accuracy of R rs (645) retrieved by MODIS is much better than that of R rs (859).The RMSP of R rs (645) and R rs (859) are 22.6% and 40.2%, respectively.The RMSP of band ratio of R rs (645) and R rs (859) is 29.6%.The accuracy of MODIS-derived R rs (859)/R rs (645) is better than that of R rs (859), but worse than that of R rs (645) (Figure 3B).The validation of match-up points for the POC-retrieval algorithm shows that the MODIS-derived POC is highly consistent with the in situ measured data.The RMSP between MODIS-derived POC and POC measured on 5 August 2014 and 24 October 2014 is 53.6% and 61.0%, respectively (Figure 3C).The validation of the atmospheric correction and POC-retrieval algorithm shows the applicability of the algorithm to MODIS satellite data.To further verify the performance of the POC algorithm to MODIS satellite data, the match-up points were used to validate the atmospheric correction (match-up points for the Rrs) and POC-retrieval results (match-up points for the Rrs and POC).The match-up MODIS satellite image was selected with the condition that the time interval for matching MODIS and in situ Rrs was less than three hours (similar to Cédric et al. [56]).The match-up pixel is consistent with the coordinates of the sampling site.The performance of this improved land target-based atmospheric correction method is acceptable for turbid water or eutrophic water, aside from some blue and green bands (Figure 3A).The atmospheric correction results indicate that the accuracy of Rrs (645) retrieved by MODIS is much better than that of Rrs (859).The RMSP of Rrs (645) and Rrs (859) are 22.6% and 40.2%, respectively.The RMSP of band ratio of Rrs (645) and Rrs (859) is 29.6%.The accuracy of MODIS-derived Rrs (859)/Rrs (645) is better than that of Rrs (859), but worse than that of Rrs (645) (Figure 3B).The validation of match-up points for the POC-retrieval algorithm shows that the MODIS-derived POC is highly consistent with the in situ measured data.The RMSP between MODIS-derived POC and POC measured on 5 August 2014 and 24 October 2014 is 53.6% and 61.0%, respectively (Figure 3C).The validation of the atmospheric correction and POC-retrieval algorithm shows the applicability of the algorithm to MODIS satellite data.

Annual Scale
MODIS R rs data were used to extract the long-term temporal and spatial distribution of POC in Taihu Lake via a POC-retrieval algorithm (Equation ( 2)).The annual mean value of POC for the entire lake shows that its spatial variation is highly significant.We should note that the POC retrieval result is invalidated in the East Lake and East Bay due to the NIR signal from submerged vegetation [57] and in surrounding land border pixels due to the heavy land border effect.The distribution of POC along the lake's bank is high, and the lake's bays, such as Meilaing, Gonghu and Zhushan Bay, show high levels of POC (East Bay should be excluded due to the invalidated POC retrieval result).The POC distribution in the center of the lake is relatively low (Figure 4).This may be explained by the fact that there is a clear influence of POC from terrain on the lake and the phytoplankton is flourishing in the bay and at the shore.The annual variation in the POC distribution shows that the distribution of POC from 2005 to 2007 was relatively more homogeneous than it was from 2002 to 2004 and from 2008 to 2013 (reflected in the difference between the central lake and coastal water).This may be caused by the long period and large area of algal bloom in the lake from 2005 to 2007 [38].The diffusion and migration of algae from bay and coastal waters obviously increases the POC in the center of the lake.

Annual Scale
MODIS Rrs data were used to extract the long-term temporal and spatial distribution of POC in Taihu Lake via a POC-retrieval algorithm (Equation ( 2)).The annual mean value of POC for the entire lake shows that its spatial variation is highly significant.We should note that the POC retrieval result is invalidated in the East Lake and East Bay due to the NIR signal from submerged vegetation [57] and in surrounding land border pixels due to the heavy land border effect.The distribution of POC along the lake's bank is high, and the lake's bays, such as Meilaing, Gonghu and Zhushan Bay, show high levels of POC (East Bay should be excluded due to the invalidated POC retrieval result).The POC distribution in the center of the lake is relatively low (Figure 4).This may be explained by the fact that there is a clear influence of POC from terrain on the lake and the phytoplankton is flourishing in the bay and at the shore.The annual variation in the POC distribution shows that the distribution of POC from 2005 to 2007 was relatively more homogeneous than it was from 2002 to 2004 and from 2008 to 2013 (reflected in the difference between the central lake and coastal water).This may be caused by the long period and large area of algal bloom in the lake from 2005 to 2007 [38].The diffusion and migration of algae from bay and coastal waters obviously increases the POC in the center of the lake.

Monthly Scale
The monthly mean POC value for thirteen years was calculated from MODIS-derived daily POC data using the proposed algorithm.The monthly variation in POC shows that the highest POC

Monthly Scale
The monthly mean POC value for thirteen years was calculated from MODIS-derived daily POC data using the proposed algorithm.The monthly variation in POC shows that the highest POC value appeared in the summer, and the relatively low value occurred in the winter (Figure 5).The high POC level in January is mainly distributed in the southwestern portion of the lake; the level decreased in February, but the relatively low POC in the north of the lake (Meilaing Bay, Gonghu Lake and the northwestern portion of the lake) in January increased significantly.The POC in the north rose dramatically from Match to April, due to the high level of algal growth.This trend of increasing POC in the northern portion of the lake was sustained until August.The POC in the coastal waters increased as well.After August, however, the POC was diffused from north to south in the lake.The POC in the center and the south part of the lake increased significantly.In general, the POC in the lake's north increased from January to August and then decreased from August to December.The variation trend in the lake's south almost maintains smoothly high concentration of POC.
Remote Sens. 2017, 9, 624 7 of 19 value appeared in the summer, and the relatively low value occurred in the winter (Figure 5).The high POC level in January is mainly distributed in the southwestern portion of the lake; the level decreased in February, but the relatively low POC in the north of the lake (Meilaing Bay, Gonghu Lake and the northwestern portion of the lake) in January increased significantly.The POC in the north rose dramatically from Match to April, due to the high level of algal growth.This trend of increasing POC in the northern portion of the lake was sustained until August.The POC in the coastal waters increased as well.After August, however, the POC was diffused from north to south in the lake.The POC in the center and the south part of the lake increased significantly.In general, the POC in the lake's north increased from January to August and then decreased from August to December.The variation trend in the lake's south almost maintains smoothly high concentration of POC.

Monthly POC Variation for Each Segment of Taihu Lake
Monthly variation in POC generally increases from January to August, and then decreases from September to December (Figure 6).This trend is probably caused for the most part by two conditions.First, the growth of algae is enhanced with rising temperatures (endogenetic input).Second, the terrestrial input is significantly increased with increasing precipitation (terrestrial input).Both temperature and precipitation have a strong seasonality in Taihu Lake, in which the variation trend is to first increase then decline.However, the asynchronism between temperature and precipitation produced a POC variation trend that is not perfectly consistent with the monthly variation in temperature or precipitation.It is obvious that there is a temporal variation cycle of the high POC concentration (75% red line in Figure 6).Sediment resuspension may be another important factor affecting the inconsistency of POC variation, temperature and precipitation [58].

Monthly POC Variation for Each Segment of Taihu Lake
Monthly variation in POC generally increases from January to August, and then decreases from September to December (Figure 6).This trend is probably caused for the most part by two conditions.First, the growth of algae is enhanced with rising temperatures (endogenetic input).Second, the terrestrial input is significantly increased with increasing precipitation (terrestrial input).Both temperature and precipitation have a strong seasonality in Taihu Lake, in which the variation trend is to first increase then decline.However, the asynchronism between temperature and precipitation produced a POC variation trend that is not perfectly consistent with the monthly variation in temperature or precipitation.It is obvious that there is a temporal variation cycle of the high POC concentration (75% red line in Figure 6).Sediment resuspension may be another important factor affecting the inconsistency of POC variation, temperature and precipitation [58].The different influences of endogenetic and terrestrial input to the POC distribution significantly affected the temporal POC variation in each lake segment.The eastern bay and eastern portion of the lake are extremely limpid with large areas of submerged aquatic plants, and no algal blooms occurred in these areas [57].The discharge of pollutants in these lake segments is very limited.The POC-retrieval results in the East Bay might be affected by the submerged plants.However, the significant seasonal variation in POC in the east bay allows for a clear view of the growth of aquatic plants (Figure 7 East Bay).This means that the atmospheric-corrected MODIS satellite data have high sensitivity to plant variation (or chlorophyll-a).The POC variation in the eastern portion of the lake is primarily caused by the aquatic plants, and terrestrial input of POC may be important as well.
The monthly POC variation in the lake's center is not clear (Figure 7 central lake), although several high POC values appeared in summer (August).This indicates that the POC from algae is not the dominant component of POC in the central lake area.This POC might result primarily from sediment resuspension [58].The other monthly POC variations (Meilaing Bay, Gonghu Bay and the lake's northwest) are highly significant (Figure 7).It is clear that the dominant component of POC in these areas derives primarily from algae, with the exception of 2013.This is primarily because of the large patches of algal bloom that occur frequently there [38,59], although riverine flux in this area is also significant.The monthly POC variation in the southwestern portion of the lake is not significant.However, different factors influence POC variation in here than in the central lake area.Large areas of algal blooms in the southwest are rare because of the hydrodynamic conditions common in open areas, although it is true that small areas of algal bloom happen frequently there [38].In addition, this area is in the Tiaoxi Basin, and large quantities of terrestrial POC input were discharged into this lake segment [60].The different influences of endogenetic and terrestrial input to the POC distribution significantly affected the temporal POC variation in each lake segment.The eastern bay and eastern portion of the lake are extremely limpid with large areas of submerged aquatic plants, and no algal blooms occurred in these areas [57].The discharge of pollutants in these lake segments is very limited.The POC-retrieval results in the East Bay might be affected by the submerged plants.However, the significant seasonal variation in POC in the east bay allows for a clear view of the growth of aquatic plants (Figure 7 East Bay).This means that the atmospheric-corrected MODIS satellite data have high sensitivity to plant variation (or chlorophyll-a).The POC variation in the eastern portion of the lake is primarily caused by the aquatic plants, and terrestrial input of POC may be important as well.
The monthly POC variation in the lake's center is not clear (Figure 7 central lake), although several high POC values appeared in summer (August).This indicates that the POC from algae is not the dominant component of POC in the central lake area.This POC might result primarily from sediment resuspension [58].The other monthly POC variations (Meilaing Bay, Gonghu Bay and the lake's northwest) are highly significant (Figure 7).It is clear that the dominant component of POC in these areas derives primarily from algae, with the exception of 2013.This is primarily because of the large patches of algal bloom that occur frequently there [38,59], although riverine flux in this area is also significant.The monthly POC variation in the southwestern portion of the lake is not significant.However, different factors influence POC variation in here than in the central lake area.Large areas of algal blooms in the southwest are rare because of the hydrodynamic conditions common in open areas, although it is true that small areas of algal bloom happen frequently there [38].In addition, this area is in the Tiaoxi Basin, and large quantities of terrestrial POC input were discharged into this lake segment [60].

Comparisons with Previous Model
The remote estimating algorithms of POC in previous studies are based on the co-variation of POC and optical properties (particulate backscattering and attenuation coefficients) and Rrs [30,32,61].The working POC algorithms in ocean water, e.g., Rrs(443)/Rrs(555) and (Rrs(555) − Rrs(443))/(Rrs(555) + Rrs(443)), have been successfully applied to estimate POC due to the dominant part of phytoplankton [12,25,31,46].These algorithms showed poor performance, with high estimation error (Mod1 and Mod2 in Table 1), but the mean value of relative error (RE) absolute value (M|RE|) performed well (Mod1) (Figure 8).This indicates that the blue-green band ratio [Rrs(443)/Rrs(555)] can work well in some cases with high Chl-a and low ISM and aCDOM.The POC will be significantly overestimated when the water has relative high ISM and low aCDOM, but POC will be underestimated when the water is dominated by OSM.The Mod2 and Mod3 provided poor performance in Taihu Lake (Table 1).

Comparisons with Previous Model
The remote estimating algorithms of POC in previous studies are based on the co-variation of POC and optical properties (particulate backscattering and attenuation coefficients) and R rs [30,32,61].The working POC algorithms in ocean water, e.g., R rs (443)/R rs (555) and (R rs (555) − R rs (443))/(R rs (555) + R rs (443)), have been successfully applied to estimate POC due to the dominant part of phytoplankton [12,25,31,46].These algorithms showed poor performance, with high estimation error (Mod1 and Mod2 in Table 1), but the mean value of relative error (RE) absolute value (M|RE|) performed well (Mod1) (Figure 8).This indicates that the blue-green band ratio [R rs (443)/R rs (555)] can work well in some cases with high Chl-a and low ISM and a CDOM .The POC will be significantly overestimated when the water has relative high ISM and low a CDOM , but POC will be underestimated when the water is dominated by OSM.The Mod2 and Mod3 provided poor performance in Taihu Lake (Table 1).    1.

Suitability and Uncertainty of the POC Model
The water constitutions of each sampling point were classified into four grades according to the RE of our empirical model: >50%, 25%-50%, −25%-25% and <−25% (Table 2).The influence of CDOM on the overestimation of POC is significant; the biggest RE in this dataset is 174% with 6.45 m −1 of aCDOM(440) and 1.36 mg/L of POC (SPM, 68.13; ISM, 55.87 mg/L).However, this is not the only influencing factor on POC retrieval accuracy.POC is underestimated with increases in Chl-a, OSM and POC (mean values of Chl-a, OSM and POC are listed in Table 2).The water constituents for the accurately estimated POC, with RE ranging from −25% to 25%, cover a wide range of Chl-a (7.5-343.7 μg/L), SPM (10.9-121.1 mg/L), ISM (3.9-111.1 mg/L) and aCDOM(440) (0.4-8.2 m −1 ).However, this does not mean the POC retrieval model can be used when the water is covered by algae matting (such as when Chl-a > 431.5 μg/L with 40.55 mg/L of POC in our dataset).The SPM  1.

Suitability and Uncertainty of the POC Model
The water constitutions of each sampling point were classified into four grades according to the RE of our empirical model: >50%, 25%-50%, −25%-25% and <−25% (Table 2).The influence of CDOM on the overestimation of POC is significant; the biggest RE in this dataset is 174% with 6.45 m −1 of a CDOM (440) and 1.36 mg/L of POC (SPM, 68.13; ISM, 55.87 mg/L).However, this is not the only influencing factor on POC retrieval accuracy.POC is underestimated with increases in Chl-a, OSM and POC (mean values of Chl-a, OSM and POC are listed in Table 2).The water constituents for the accurately estimated POC, with RE ranging from −25% to 25%, cover a wide range of Chl-a (7.5-343.7 µg/L), SPM (10.9-121.1 mg/L), ISM (3.9-111.1 mg/L) and a CDOM (440) The optical activity matter, Chl-a and SPM, which can be estimated analytically by remote sensing reflectance, are the main sources of POC in inland water.Field-measured data presented a high correlation between POC and Chl-a (POC = 59.59 × Chl-a + 1.0147, R 2 = 0.66) (Figure 9A), but a low correlation between POC and SPM (POC = 0.0855 × SPM + 1.38, R 2 = 0.18) (Figure 9B).However, the relationship between POC and SPM can be separated into two parts, i.e., low POC/SPM (0.0471 ± 0.0236, Figure 9B, in black ellipse) and high POC/SPM (0.2433 ± 0.1817, Figure 9B, in gray ellipse).The relationships between POC and SPM are POC = 0.0235 × SPM + 0.8214 (R 2 = 0.42, low POC/SPM) and POC = 0.1354 × SPM + 3.2159 (R 2 = 0.44, high POC/SPM).Obviously, the POC could be over-or underestimated when the dominant matter of SPM shifts between ISM and OSM even POC is high correlated to SPM.The validated relationship between POC and Chl-a or SPM may be a basic condition to accurately estimate POC [32,33,61].The dominant effect of Chl-a or SPM on POC varied with algal growth, sediment resuspension, river transport, etc.The relationship between POC and Chl-a (or SPM) presented in this dataset is relatively complex [32,33].
The optical activity matter, Chl-a and SPM, which can be estimated analytically by remote sensing reflectance, are the main sources of POC in inland water.Field-measured data presented a high correlation between POC and Chl-a (POC = 59.59 × Chl-a + 1.0147, R 2 = 0.66) (Figure 9A), but a low correlation between POC and SPM (POC = 0.0855 × SPM + 1.38, R 2 = 0.18) (Figure 9B).However, the relationship between POC and SPM can be separated into two parts, i.e., low POC/SPM (0.0471 ± 0.0236, Figure 9B, in black ellipse) and high POC/SPM (0.2433 ± 0.1817, Figure 9B, in gray ellipse).The relationships between POC and SPM are POC = 0.0235 × SPM + 0.8214 (R 2 = 0.42, low POC/SPM) and POC = 0.1354 × SPM + 3.2159 (R 2 = 0.44, high POC/SPM).Obviously, the POC could be over-or underestimated when the dominant matter of SPM shifts between ISM and OSM even POC is high correlated to SPM.The validated relationship between POC and Chl-a or SPM may be a basic condition to accurately estimate POC [32,33,61].The dominant effect of Chl-a or SPM on POC varied with algal growth, sediment resuspension, river transport, etc.The relationship between POC and Chl-a (or SPM) presented in this dataset is relatively complex [32,33].The independent datasets from Duan et al. [33] and Jiang et al. [32], for which POC is either correlated to SPM [33] or Chl-a [32], were used to evaluate the applicability of using this empirical model to estimate POC with different dominant water constituents (Figure 10).The validation result from the Chl-a dominated POC [32] suggests that this empirical model estimate POC well, except for at two points, i.e., one was overestimated 169.6% and another was underestimated −52.9%The independent datasets from Duan et al. [33] and Jiang et al. [32], for which POC is either correlated to SPM [33] or Chl-a [32], were used to evaluate the applicability of using this empirical model to estimate POC with different dominant water constituents (Figure 10).The validation result from the Chl-a dominated POC [32] suggests that this empirical model estimate POC well, except for at two points, i.e., one was overestimated 169.6% and another was underestimated −52.9% (measured 0.738 mg/L, estimated 1.99; measured 1.530 mg/L, estimated 0.721).The validation result from the SPM dominated POC [33] shows that most of the samples were underestimated (mean value: −26.96%), the biggest RE was −63.36%.This is inconsistent with the result in Figure 8, where POC is overestimated when POC is low.This is primary because OSM is the main component of SPM for in Duan and coworkers' [33] data, but the main component is ISM for our data.Thus, this empirical model will overestimate POC in the ISM-dominated water and underestimate POC in OSM-dominated water.This may conform to the result in Figure 9B.(measured 0.738 mg/L, estimated 1.99; measured 1.530 mg/L, estimated 0.721).The validation result from the SPM dominated POC [33] shows that most of the samples were underestimated (mean value: −26.96%), the biggest RE was −63.36%.This is inconsistent with the result in Figure 8, where POC is overestimated when POC is low.This is primary because OSM is the main component of SPM for in Duan and coworkers' [33] data, but the main component is ISM for our data.Thus, this empirical model will overestimate POC in the ISM-dominated water and underestimate POC in OSM-dominated water.This may conform to the result in Figure 9B.The effect of atmospheric correction on the POC algorithm also should be considered in addition to the dominant water constituents of POC.Chl-a will be overestimated when water reflectance is very low in the NIR wavelength [47].Similarly, low POC will be overestimated due to overestimated reflectance, which should be relatively low at NIR wavelengths.The traditional dark pixel atmospheric correction model makes it hard to get accurate reflectance at NIR wavelengths in Taihu Lake due to the highly turbid and eutrophic water.Getting accurate reflectance from satellites is still difficult although we used an improved land target-based atmospheric correction method (Figure 3A).However, some noise caused by the atmosphere can be weakened by the band ratio [35].This is consistent with the result in Figure 3B and indicates that the NIR-red band-ratio model in this paper can reduce some errors systematically caused by atmosphere [35].Consequently, a perfect atmospheric correction method is not crucial to our POC model.

Probable Source of POC
Taihu Lake is a shallow inland eutrophic lake with high productivity.The POC in the lake has three likely sources: phytoplankton, sediment resuspension and terrestrial discharge.The Chl-a, which is representative of the living phytoplankton biomass, is mostly present in phytoplankton (such microcystis aeruginosa, scenedesmus obliquus and diatom) in the lake [62].The ratio of POC and Chl-a (POC/Chl-a) is a good index for evaluating the influence of terrestrial POC input to total POC in the lake [63][64][65].The POC/Chl-a in Taihu Lake ranges from 8.63 to 718.65, with a mean value of 96.47 ± 92.58 (Figure 11A).The POC/Chl-a is in the range of fresh phytoplankton ranges, between 33 and 100 [66]; this is much higher than that of sub-Antarctic water (50-60) [67], Great Wall Cove and Ardley Cove (74.62 ± 27.97) [65] and the East China Sea (64) [68].However, it is much lower than that of dystrophic water (250-2500) [69], the Pacific Ocean (178 ± 30) [70], the Atlantic Ocean (253) and other ocean subsystems (Table 2 in Fabiano et al. [71]; Arrigo et al. [72]; Table1 in Hadjimitsis et al.[73]).The effect of atmospheric correction on the POC algorithm also should be considered in addition to the dominant water constituents of POC.Chl-a will be overestimated when water reflectance is very low in the NIR wavelength [47].Similarly, low POC will be overestimated due to overestimated reflectance, which should be relatively low at NIR wavelengths.The traditional dark pixel atmospheric correction model makes it hard to get accurate reflectance at NIR wavelengths in Taihu Lake due to the highly turbid and eutrophic water.Getting accurate reflectance from satellites is still difficult although we used an improved land target-based atmospheric correction method (Figure 3A).However, some noise caused by the atmosphere can be weakened by the band ratio [35].This is consistent with the result in Figure 3B and indicates that the NIR-red band-ratio model in this paper can reduce some errors systematically caused by atmosphere [35].Consequently, a perfect atmospheric correction method is not crucial to our POC model.

in Hadjimitsis et al. [73]).
POC/SPM is an indicator of the particle organic carbon fraction of the SPM [64,74].The POC/SPM in Taihu Lake ranges from 0.02 to 0.72, with a mean value of 0.11 ± 0.13.This is similar to the inland water of the Chena River in Alaska, where the mean POC/SPM value is 0.158 ± 0.074 [74].The value for Taihu Lake is much higher than it is for the Tana River, a highly turbid river in Kenya (0.042 above reservoir and 0.046 below dam, due to the negligible contribution of phytoplankton), [75,76] and for sediment in Taihu Lake itself [58].Unsurprisingly, it is much less than that of clear ocean water (Pacific entrance, 0.19; Cabo Pilar area, 0.28) [71] and extremely eutrophic reservoirs (Masinga Reservoir, 0.326 ± 0.036) [75].
Low POC/Chl-a and high POC/SPM indicate a strong influence from phytoplankton on POC and little effect from suspended particles.The ratio of POC/Chl-a and POC/SPM in Taihu Lake indicates that both phytoplankton and suspended particles are the factors that most influence POC distribution.The correlation analysis of POC and water quality (Pearson's correlation coefficient, r) also demonstrates that phytoplankton (high r of POC, PC and Chl-a) is an important source of POC.Suspended particle (SPM and OSM) is also a primary source, deriving mostly from terrestrial discharge and sediment resuspension [77][78][79].This is because the SPM level is one thousand times higher than Chl-a (1886 ± 1443 times), although the r of POC and phytoplankton is slightly greater than the r of POC and suspended particle (Figure 11B).
Remote Sens. 2017, 9, 624 13 of 19 POC/SPM is an indicator of the particle organic carbon fraction of the SPM [64,74].The POC/SPM in Taihu Lake ranges from 0.02 to 0.72, with a mean value of 0.11 ± 0.13.This is similar to the inland water of the Chena River in Alaska, where the mean POC/SPM value is 0.158 ± 0.074 [74].The value for Taihu Lake is much higher than it is for the Tana River, a highly turbid river in Kenya (0.042 above reservoir and 0.046 below dam, due to the negligible contribution of phytoplankton), [75,76] and for sediment in Taihu Lake itself [58].Unsurprisingly, it is much less than that of clear ocean water (Pacific entrance, 0.19; Cabo Pilar area, 0.28) [71] and extremely eutrophic reservoirs (Masinga Reservoir, 0.326 ± 0.036) [75].
Low POC/Chl-a and high POC/SPM indicate a strong influence from phytoplankton on POC and little effect from suspended particles.The ratio of POC/Chl-a and POC/SPM in Taihu Lake indicates that both phytoplankton and suspended particles are the factors that most influence POC distribution.The correlation analysis of POC and water quality (Pearson's correlation coefficient, r) also demonstrates that phytoplankton (high r of POC, PC and Chl-a) is an important source of POC.Suspended particle (SPM and OSM) is also a primary source, deriving mostly from terrestrial discharge and sediment resuspension [77][78][79].This is because the SPM level is one thousand times higher than Chl-a (1886 ± 1443 times), although the r of POC and phytoplankton is slightly greater than the r of POC and suspended particle (Figure 11B).

Factors Influencing the Source of POC
The source of POC in Taihu Lake is complex because of its high productivity, riverine flux and sediment resuspension.The spatio-temporal distributions of these influencing factors for each lake segment aggravate the inconsistent findings of the POC source.Temperature, wind speed and precipitation greatly affect algal growth, sediment resuspension and riverine flux.These three parameters were therefore used to represent the influence of algae, sediment resuspension and riverine flux on POC.The high intensity of precipitation indicating high riverine flux is consistent with the high POC for all lake segments (blue line in Figure 12) [34].The only low-intensity precipitation that influenced POC was in Meiliang Bay (red line in Figure 12).High temperature is not perfectly consistent with high POC for all lake segments.It most often influences just Meilaing Bay and Gonghu Bay (green line in Figure 12), although it sometimes influences all lake segments, because of the large algal bloom area that occurred in 2007 (green-dash line in Figure 12) [38].The pattern of influence for wind speed indicates that the contribution of sediment resuspension to POC is extremely complex.

Factors Influencing the Source of POC
The source of POC in Taihu Lake is complex because of its high productivity, riverine flux and sediment resuspension.The spatio-temporal distributions of these influencing factors for each lake segment aggravate the inconsistent findings of the POC source.Temperature, wind speed and precipitation greatly affect algal growth, sediment resuspension and riverine flux.These three parameters were therefore used to represent the influence of algae, sediment resuspension and riverine flux on POC.The high intensity of precipitation indicating high riverine flux is consistent with the high POC for all lake segments (blue line in Figure 12) [34].The only low-intensity precipitation that influenced POC was in Meiliang Bay (red line in Figure 12).High temperature is not perfectly consistent with high POC for all lake segments.It most often influences just Meilaing Bay and Gonghu Bay (green line in Figure 12), although it sometimes influences all lake segments, because of the large algal bloom area that occurred in 2007 (green-dash line in Figure 12) [38].The pattern of influence for wind speed indicates that the contribution of sediment resuspension to POC is extremely complex.
The correlation analysis of POC for each lake segment and temperature, precipitation and wind speed shows that the change in POC in Taihu Lake is quite similar to the change in the center and northwest of the lake (Pearson's correlation coefficients are 0.90 and 0.91, respectively).POC in the northwest and southwest of the lake is similar to POC in the center.In Gonghu Lake, POC is similar to Meiliang Bay (Table 3).The high correlation coefficients for temperature in Meiliang Bay and Gonghu Bay indicate that the POC in these two bays is highly affected by algal growth.The correlation coefficients of wind speed and POC for each lake segment demonstrate that the sediment resuspension has a significant effect on POC in the center and southwest of the lake.The correlation coefficients for precipitation and POC in each lake segment show that riverine flux has a strong influence on POC in Meiliang Bay and Gonghu Bay, as well as in the northwest and southwest of the lake.The riverine flux was found to slightly affect POC in the center of the lake.
Remote Sens. 2017, 9, 624 14 of 19 The correlation analysis of POC for each lake segment and temperature, precipitation and wind speed shows that the change in POC in Taihu Lake is quite similar to the change in the center and northwest of the lake (Pearson's correlation coefficients are 0.90 and 0.91, respectively).POC in the northwest and southwest of the lake is similar to POC in the center.In Gonghu Lake, POC is similar to Meiliang Bay (Table 3).The high correlation coefficients for temperature in Meiliang Bay and Gonghu Bay indicate that the POC in these two bays is highly affected by algal growth.The correlation coefficients of wind speed and POC for each lake segment demonstrate that the sediment resuspension has a significant effect on POC in the center and southwest of the lake.The correlation coefficients for precipitation and POC in each lake segment show that riverine flux has a strong influence on POC in Meiliang Bay and Gonghu Bay, as well as in the northwest and southwest of the lake.The riverine flux was found to slightly affect POC in the center of the lake.Table 3. Correlation coefficients of particulate organic carbon (POC) to T, W, P and autocorrelation coefficients of POC for each lake segment (T is temperature (°C), P is precipitation (mm), and W is wind speed (m/s)); the number of samples in the statistical analysis is 197.

Conclusions
Based on the in situ measurements of particulate organic carbon (POC) and remote-sensing reflectance, an NIR-Red empirical POC algorithm was developed to observe the long-term variation of POC from MODIS-Aqua.The algorithm was found to provide acceptable validation results, with 38.9% RMSP for the validation data and 31.5% RMSP for the independent validation data.The application of the empirical algorithm to the atmospheric corrected MODIS images showed that the algorithm is highly consistent with the in situ measured POC of match-up points.
The temporal and spatial variation of POC and the analysis of factors that influence POC distribution led to some interesting findings.High POC levels were primarily distributed in the lake near the shore and bays of Meiliang, Zhushan and Gonghu, because of the influence of both riverine flux and algal growth.Sediment resuspension is another important source of POC that significantly affected its distribution during the research period when there was high wind speed.Heavy precipitation will increase the influence of terrigenous POC on the POC distribution in the bays, and this influence can even reach the center of the lake.Algal bloom (overgrowth) will dramatically increase POC.

Figure 1 .
Figure 1.Study area: (A) map of the land-use cover in the Taihu basin in 2000; and (B) a picture showing the sampling points at Taihu Lake for the cruises of 5 August 2014, 24 October 2014, and 21 May 2015' the sample sites of 10 June 2015 and 15 July 2015 are same as 21 May 2015.The segments of the lake include Meiliang Bay (MLB); Gonghu Bay (GHB); East Bay (EB); and the central (CL), northwestern (NWL), southwestern (SWL) and eastern (EL) parts of the lake.

Figure 1 .
Figure 1.Study area: (A) map of the land-use cover in the Taihu basin in 2000; and (B) a picture showing the sampling points at Taihu Lake for the cruises of 5 August 2014, 24 October 2014, and 21 May 2015' the sample sites of 10 June 2015 and 15 July 2015 are same as 21 May 2015.The segments of the lake include Meiliang Bay (MLB); Gonghu Bay (GHB); East Bay (EB); and the central (CL), northwestern (NWL), southwestern (SWL) and eastern (EL) parts of the lake.

Figure 2 .
Figure 2. Calibration and validation of near infrared-red (NIR-Red) band ratio algorithm to the particulate organic carbon (POC) in Taihu Lake: (A) calibration results; and (B) validation results from data entirely independent of the calibration data.

Figure 3 .
Figure 3. (A) Atmospheric correction result for the comparison of in situ measured remote sensing reflectance (Rrs)(dot dashed line) and Moderate Resolution Imaging Spectroradiometer (MODIS) -derived Rrs (dot solid line), the dotted line with the same color is the pair of measured and corrected Rrs, and the high reflectance at 859 nm is induced by extremely high amounts of phytoplankton; (B) the atmospheric correction result for the comparison of in situ measured Rrs and MODIS-derived Rrs at specific bands (645 nm and 859 nm); and (C) validation results for the MODIS-derived particulate organic carbon (POC) of match-up points.

Figure 2 .
Figure 2. Calibration and validation of near infrared-red (NIR-Red) band ratio algorithm to the particulate organic carbon (POC) in Taihu Lake: (A) calibration results; and (B) validation results from data entirely independent of the calibration data.

Figure 2 .
Figure 2. Calibration and validation of near infrared-red (NIR-Red) band ratio algorithm to the particulate organic carbon (POC) in Taihu Lake: (A) calibration results; and (B) validation results from data entirely independent of the calibration data.

Figure 3 .
Figure 3. (A) Atmospheric correction result for the comparison of in situ measured remote sensing reflectance (Rrs)(dot dashed line) and Moderate Resolution Imaging Spectroradiometer (MODIS) -derived Rrs (dot solid line), the dotted line with the same color is the pair of measured and corrected Rrs, and the high reflectance at 859 nm is induced by extremely high amounts of phytoplankton; (B) the atmospheric correction result for the comparison of in situ measured Rrs and MODIS-derived Rrs at specific bands (645 nm and 859 nm); and (C) validation results for the MODIS-derived particulate organic carbon (POC) of match-up points.

Figure 3 .
Figure 3. (A) Atmospheric correction result for the comparison of in situ measured remote sensing reflectance (R rs ) (dot dashed line) and Moderate Resolution Imaging Spectroradiometer (MODIS) -derived R rs (dot solid line), the dotted line with the same color is the pair of measured and corrected R rs , and the high reflectance at 859 nm is induced by extremely high amounts of phytoplankton; (B) the atmospheric correction result for the comparison of in situ measured R rs and MODIS-derived R rs at specific bands (645 nm and 859 nm); and (C) validation results for the MODIS-derived particulate organic carbon (POC) of match-up points.

Figure 4 .
Figure 4. Annual variations in particulate organic carbon (POC) in Taihu Lake.Annual mean value of POC is the mean value for all images each year, from 2002 to 2014.The POC retrieval result in East Lake and East Bay is invalidated and should not be treated as a real distribution of POC (masked by stripe in first map).

Figure 4 .
Figure 4. Annual variations in particulate organic carbon (POC) in Taihu Lake.Annual mean value of POC is the mean value for all images each year, from 2002 to 2014.The POC retrieval result in East Lake and East Bay is invalidated and should not be treated as a real distribution of POC (masked by stripe in first map).

Figure 5 .
Figure 5. Annual variations in particulate organic carbon (POC) in Taihu Lake.Annual mean value of POC is the mean value for all images each year, from 2002 to 2014.The POC retrieval result in East Lake and East Bay is invalidated and should not be treated as a real distribution of POC (masked by stripe in the first map).

Figure 5 .
Figure 5. Annual variations in particulate organic carbon (POC) in Taihu Lake.Annual mean value of POC is the mean value for all images each year, from 2002 to 2014.The POC retrieval result in East Lake and East Bay is invalidated and should not be treated as a real distribution of POC (masked by stripe in the first map).

Figure 6 .
Figure 6.Monthly particulate organic carbon (POC) variations from June 2002 to April 2014 for the mean value of Taihu Lake.The POC retrieval result in East Lake and East Bay is invalidated and did not participate in the analysis of monthly POC variations.The boundary of the box is the 25% and 75% position of POC, the line in the box is the median of POC, and the boundary of the line beside the box is the minimum and maximum value of POC.

Figure 6 .
Figure 6.Monthly particulate organic carbon (POC) variations from June 2002 to April 2014 for the mean value of Taihu Lake.The POC retrieval result in East Lake and East Bay is invalidated and did not participate in the analysis of monthly POC variations.The boundary of the box is the 25% and 75% position of POC, the line in the box is the median of POC, and the boundary of the line beside the box is the minimum and maximum value of POC.

Figure 7 .
Figure 7. Monthly particulate organic carbon (POC) variations from June 2002 to April 2014 for the mean value of Taihu Lake.The POC retrieval result in East Lake and East Bay is invalidated and did not participate in the analysis of monthly POC variations.The boundary of the box is the 25% and 75% position of POC, the line in the box is the median of POC, and the boundary of the line beside the box is the minimum and maximum value of POC.

Figure 7 .
Figure 7. Monthly particulate organic carbon (POC) variations from June 2002 to April 2014 for the mean value of Taihu Lake.The POC retrieval result in East Lake and East Bay is invalidated and did not participate in the analysis of monthly POC variations.The boundary of the box is the 25% and 75% position of POC, the line in the box is the median of POC, and the boundary of the line beside the box is the minimum and maximum value of POC.

Figure 8 .
Figure 8.Comparison of measured and retrieved particulate organic carbon (POC) for previous POC estimation algorithms.The information of Mod1, Mod2, Mod3, Mod4 and Mod5 are listed in Table1.

Figure 8 .
Figure 8.Comparison of measured and retrieved particulate organic carbon (POC) for previous POC estimation algorithms.The information of Mod1, Mod2, Mod3, Mod4 and Mod5 are listed in Table1.

Figure 9 .
Figure 9. Relationship between particulate organic carbon (POC) and optical activity matter: (A) Chlorophyll-a (Chl-a); and (B) Suspended particulate matter (SPM) for different RE ranges.

Figure 9 .
Figure 9. Relationship between particulate organic carbon (POC) and optical activity matter: (A) Chlorophyll-a (Chl-a); and (B) Suspended particulate matter (SPM) for different RE ranges.

Figure 10 .
Figure 10.Evaluation of the empirical model applicability for different dominant water constituents.The data include 21 samples (red circle) from Duan et al. [33] and 14 samples (blue box) from Jiang Duan et al.[32]

Figure 10 .
Figure 10.Evaluation of the empirical model applicability for different dominant water constituents.The data include 21 samples (red circle) from Duan et al. [33] and 14 samples (blue box) from Jiang Duan et al. [32].

Figure 12 .
Figure 12.Time-series of T, W, and P to particulate organic carbon (POC) from June 2002 to April 2014 for each lake segment (T is temperature (°C), P is precipitation (mm), and W is wind speed (m/s)).

Figure 12 .
Figure 12.Time-series of T, W, and P to particulate organic carbon (POC) from June 2002 to April 2014 for each lake segment (T is temperature ( • C), P is precipitation (mm), and W is wind speed (m/s)).

Table 1 .
Summary of error statistics for particulate organic carbon (POC) algorithms.All the parameters in the algorithms were re-calibrated by POC and remote sensing reflectance (R rs ) data collected from Taihu Lake.
[33] ± 0.29 R 2 is the linear determined coefficient between measured and estimated POC, RE is the relative error, calculated by RE = (POC r − POC m )/POC m , M|RE| is the mean value of the absolute value of RE, RMSP is root mean square error percentage and STDEV is standard deviation.For Mod4*, R rs (1240) in Duan and coworkers'[33]algorithm is out of our FieldSpec spectroradiometer (ASD, 350-1050 nm range) measurement range, thus R rs (900) was used instead of R rs (1240) in this study.

Table 1 .
Summary of error statistics for particulate organic carbon (POC) algorithms.All the parameters in the algorithms were re-calibrated by POC and remote sensing reflectance (Rrs) data collected from Taihu Lake.
[33]is the linear determined coefficient between measured and estimated POC, RE is the relative error, calculated by RE = (POCr − POCm)/POCm, M|RE| is the mean value of the absolute value of RE, RMSP is root mean square error percentage and STDEV is standard deviation.For Mod4*, Rrs(1240) in Duan and coworkers'[33]algorithm is out of our FieldSpec spectroradiometer (ASD, 350-1050 nm range) measurement range, thus Rrs(900) was used instead of Rrs(1240) in this study.

Table 2 .
Water constituents for different accuracy ranges retrieved by particulate organic carbon (POC) estimation algorithm.The value in the bracket is the mean value of each water constituent in the range above it.

Table 2 .
Remote Sens. 2017, 9, 624 11 of 19 and ISM ranges for different grades of RE with insignificant distinction suggest that this empirical model can be used in highly turbid water (ISM can reach 89.27 mg/L with median Chl-a of 66.96 μg/L and median POC of 5.35 mg/L).Water constituents for different accuracy ranges retrieved by particulate organic carbon (POC) estimation algorithm.The value in the bracket is the mean value of each water constituent in the range above it.

Table 3 .
Correlation coefficients of particulate organic carbon (POC) to T, W, P and autocorrelation coefficients of POC for each lake segment (T is temperature ( • C), P is precipitation (mm), and W is wind speed (m/s)); the number of samples in the statistical analysis is 197.