Temporal and Spatial Dynamics of Phytoplankton Primary Production in Lake Taihu Derived from MODIS Data

We investigated the long-term variations in primary production in Lake Taihu using Moderate Resolution Imaging Spectroradiometer (MODIS) data, based on the Vertically Generalized Production Model (VGPM). We firstly test the applicability of VGPM in Lake Taihu by comparing the results between the model-derived and the in situ results, and the results showed that a strong significant correlation (R2 = 0.753, p < 0.001, n = 63). Then, VGPM was used to map temporal-spatial distributions of primary production in Lake Taihu. The annual mean daily primary production of Lake Taihu from 2003 to 2013 was 1094.06 ± 720.74 mg·C·m−2·d−1. Long-term primary production maps estimated from the MODIS data demonstrated marked temporal and spatial variations. Spatially, the primary production in bays, especially in Zhushan Bay and Meiliang Bay, was consistently higher than that in the open area of Lake Taihu, which was caused by chlorophyll-a concentrations resulting from high nutrient concentrations. Temporally, the seasonal variation of primary production from 2003 to 2013 was: summer > autumn > spring > winter, with significantly higher primary production found in summer and autumn than in winter (p < 0.005, t-test), primarily caused by seasonal variations in water temperature. On a monthly scale, the primary production exerts a clear character of bimodality, increasing from January to May, decreasing in June or July, and finally reaching its highest value during August or September. Wind is another important factor that could affect the spatial variations of primary production in the large, eutrophic and shallow Lake Taihu.


Introduction
Phytoplankton is at the bottom of the food chain, creating fresh organic matter from inorganic nutrients, carbon dioxide and energy from sunlight, which strongly influence nutrient concentrations and support higher trophic levels such as zooplankton and filter feeder.When either system is employed, phytoplankton creates organic matter from inorganic compounds and carbon dioxide, which is called phytoplankton primary production [1,2].The rate of phytoplankton primary production is a fundamental property of aquatic system and measurements of primary production are critical to our understanding of the carbon cycle [3].Estimation of the phytoplankton primary production is also an important topic in fisheries resource management and global change [4].
One of the major goals of modern biological oceanography is to acquire a better understanding of primary production in various oceanic provinces, with a special emphasis on marine carbon cycling and climate change on regional to global scales [5].As an important part of the global carbon cycle, estimating the primary production in lakes is also the major goals of modern limnology, and important for us to understand the regional ecological environment [4].Actually, lake primary production, which is commonly regarded as the photosynthetic capacity of phytoplankton for the unit volume, is an important parameter to describe ecosystem and environmental characteristics of lake [6].
Traditionally, the primary production of an aquatic environment was obtained from cruise samples.While in situ experiments provided accurate estimates of primary production in small volumes of water, they may not be easily extrapolated to lake-wide estimates [7][8][9].Moreover, these experiments provide an integrated measure of production that is dependent on many variables (e.g., phytoplankton biomass, light, temperature, etc.), thus limiting their predictive value.Early ways to estimate primary production for the great lakes may be biased because of deficiencies in traditional collection and incubation techniques [10][11][12][13][14][15].
Satellite data have been widely used to derive several lake biogeochemical parameters, such as, total suspended matter, chlorophyll a (Chla), and diffuse attenuation coefficient of photosynthetically active radiation (K d (PAR)).Satellite remote sensing has many advantages in estimation of lake biogeochemical parameters.The repeated coverage by remote sensing enables the detection of the temporal and spatial variation, which has proven beneficial for rapidly estimating lake biogeochemical parameters [16].Additionally, satellite remote sensing is more practical and economical than the other monitoring methods and also can be easily integrated into geographic information [17][18][19].However, remote sensing data cannot directly provide information on primary production without support of models [20].Vertically Generalized Productivity Model (VGPM) is the most widely used to estimate primary production due to its minimal input parameters, which has been validated by thousands of in situ measured data over several orders of magnitude.Among these input parameters, Chla and euphotic depth are could be derived from remote sensing data, providing a possibility for estimating primary production using remote sensing data.Many studies have successfully estimated phytoplankton primary production by combining remote sensing data to VGPM for open ocean [21][22][23].However, due to difficulty in retrieval of Chla and euphotic depth in complexly optical waters, there have been few applications of the VGPM to lakes [24][25][26].
As the third-largest freshwater lake in China, Lake Taihu has a water surface area of 2338 km 2 , and a mean depth of 1.9 m [24].This large, shallow lake is located in one of the world's most heavily populated regions, which has experienced rapid economic development in recent years.Lake Taihu serves as a key drinking water resource for the approximately 10 million local residents residing in several large nearby cities, such as Shanghai, Suzhou, Wuxi and Huzhou [27,28]; it also has additional economic functions including tourism, fisheries and shipping [29].Besides, Lake Taihu Basin also contributes 11.6% of China's gross domestic product, despite accounting for only 0.38% of China's total area.With economic development, a lot of wastewater, sewage and polluted water surrounding urban and rural areas were discharged to Lake Taihu, resulting in phytoplankton blooms and then triggering a drinking water crisis in 2007 [30].Phytoplankton are commonly the most important primary producer in lake ecosystems, strongly influencing nutrient concentrations and supporting higher trophic levels such as zooplankton and filter feeders.Phytoplankton primary production is an important indicator of water quality and algal biomass.Therefore, describing the temporal and spatial variations in phytoplankton assemblages and primary production is crucial for understanding the development of eutrophication and ecosystem dynamics in Lake Taihu [31].
Biogeochemical parameters in Lake Taihu are always characterized by highly temporal-spatial dynamics, indicating that estimates of primary production using traditional sampling methods may lead to large uncertainty when we investigate primary production dynamics in the whole lake for a long time period.
In recent years, we have calibrated and validated Chla and euphotic depth estimation model in Lake Taihu, making it possible for us to estimate the phytoplankton primary production based on the VGPM model.Therefore, the aims of ours study were to (1) assess the accuracy and feasibility of VGPM in Lake Taihu using in situ measurement data; (2) estimate the phytoplankton primary production of Lake Taihu from 2003 to 2013 using MODIS-derived data; (3) analyze the temporal and spatial variations in primary production in Lake Taihu and discuss the potential affecting factors.

Study Area
Lake Taihu is a crucial water source for the Yangtze River Delta, located between 30 • 56 N-31 • 33 N and 119 • 54 E-120 • 36 E (Figure 1) [30].In this article, Lake Taihu can be divided into six areas based on shoreline geometry, environmental factors, and human activities: Meiliang Bay, Zhushan Bay, Gonghu Bay, Xukou Bay, East Lake Taihu and the open area (Figure 1) [32,33].It should be noted that the East Lake Taihu, Gonghu Bay and Xukou Bay were outside of the scope of the study area, defined as the non-region-of-interest (NROI) for this research, because the retrieved Chla were unusually high in these regions.In fact, these regions were marcophytic-dominated rather than phytoplankton-dominated [34], and the abundant submerged plants resulted in the deviation of retrieved Chla by MODIS data because of similar spectral characteristics with algae.
production of Lake Taihu from 2003 to 2013 using MODIS-derived data; (3) analyze the temporal and spatial variations in primary production in Lake Taihu and discuss the potential affecting factors.

Study Area
Lake Taihu is a crucial water source for the Yangtze River Delta, located between 30°56′N-31°33′N and 119°54′E-120°36′E (Figure 1) [30].In this article, Lake Taihu can be divided into six areas based on shoreline geometry, environmental factors, and human activities: Meiliang Bay, Zhushan Bay, Gonghu Bay, Xukou Bay, East Lake Taihu and the open area (Figure 1) [32,33].It should be noted that the East Lake Taihu, Gonghu Bay and Xukou Bay were outside of the scope of the study area, defined as the non-region-of-interest (NROI) for this research, because the retrieved Chla were unusually high in these regions.In fact, these regions were marcophytic-dominated rather than phytoplankton-dominated [34], and the abundant submerged plants resulted in the deviation of retrieved Chla by MODIS data because of similar spectral characteristics with algae.

Field Data
A total of 63 samples were used to validate the VGPM model and a total of 20 samples were used to analyze the correlation between phytoplankton primary production and nutrients.63 phytoplankton primary production samples were collected between 2005 and 2013 from the two sites (THL04 and THL08).Rates of phytoplankton primary production and respiration were measured using the light-dark bottle incubation method.Within the dark bottle, only respiration can occur; and within the light bottle, both respiration and photosynthesis can occur.From the difference in O2 within the light and dark bottles relative to the initial O2, the rate of primary production can be calculated [35][36][37].The nutrient concentration data in 2013 were collected from 32 sites THL00-THL32 (except THL02) (Figure 1).

Field Data
A total of 63 samples were used to validate the VGPM model and a total of 20 samples were used to analyze the correlation between phytoplankton primary production and nutrients.63 phytoplankton primary production samples were collected between 2005 and 2013 from the two sites (THL04 and THL08).Rates of phytoplankton primary production and respiration were measured using the light-dark bottle incubation method.Within the dark bottle, only respiration can occur; and within the light bottle, both respiration and photosynthesis can occur.From the difference in O 2 within the light and dark bottles relative to the initial O 2 , the rate of primary production can be calculated [35][36][37].

Image Data Description and Processing
The MODIS-Aqua data have been freely available since 2002 and has a maximum spatial resolution of 250 m (band 1 and band 2) and a very short revisit interval (1 image/day).MODIS-Aqua L-0 data from January 2003 to December 2013 (more than 4000 images) were downloaded from NASA's Goddard Space Flight Center (GSFC, http://oceancolor.gsfc.nasa.gov/).Due to the clouds, cloud shadows, or thick aerosols, not all downloaded images were used in this study.We selected 1086 cloud-free images of Lake Taihu from January 2003 to December 2013.These images were processed to Level-1 (calibrated spectral radiance) using the software package SeaDAS (version 6.0).

Method Description
After investigating the variability observed in primary production by assembling a dataset of 11,283 14 C-based measurements of daily carbon fixation collected at 1698 oceanographic stations in both open ocean and coastal waters, Behrenfeld and Falkowski discovered a consistent trend in the vertical distribution of primary production and observed that the VGPM well accounted for the trend in normalized primary production estimated from depth integrated primary production [14].VGPM model has been validated by thousands of in situ measured data points on a large scale and in different water areas.Therefore, it has been used widely to accurately estimate primary production.
The equation of VGPM is given below: where PP eu represents primary production from sea surface to the euphotic depth (mg•C•m −2 •d −1 .),P B opt is the maximum optimal rate of carbon fixation of water mass, D irr is the illumination period, Zeu is the euphotic depth, E z represents Photosynthetically Active Radiation(PAR) that is photosynthetic effective radiation at depth Z (mol quanta/m 2 ), E max represents the maximum photosynthetic effective radiation, E opt is the PAR where P B opt is located and β d is the initial slope of curve P-I chlorophyll concentration at depth Z.
After contrasting the results calculated by the VGPM model and ship-measured data, Behrenfeld and Falkowski considered that the model describes 79% of temporal-spatial variety of primary production (n = 10857).After simplification, the model can be presented as: where P B opt is the maximum rate of carbon fixation within a water column (mg C/(mg Chl•h)), E 0 is PAR of the surface photosynthetic effective radiation of the lake, Z eu is the euphotic depth, C opt is Chla concentration where P B opt is located, which can be replaced by remote sensed surface Chla concentration, and D irr is day length in decimal hours.
The model parameters can be acquired through the following methods: (1) Calculation of P B opt is based on the experienced relationship between P B opt and the lake surface temperature (LST, • C) obtained with field investigation data provided by Behrenfeld and Falkowski.
Therefore, P B opt can be commonly considered as a function of the lake surface temperature, which was described by Behrenfeld and Falkowski [12,14] as: where The LST data from Lake Taihu can be calculated from MODIS data.Previous studies have showed that MODIS-derived LST and in situ water temperatures in Lake Taihu were significantly correlated, with a coefficient of determination higher than 0.96 and a root mean square error between 1.2 • C and 1.8 • C [38].Then, the P B opt can be derived from MODIS-derived LST using the Equations ( 3) and ( 4). ( 2) Chla can also be acquired from the MODIS data based on the model which was developed by Shi et al using Rayleigh-corrected reflectance as follows [31]: where Index MODIS = (EXP(R rc 645) − EXP(R rc 859))/(EXP(R rc 645) + EXP(R rc 859)).A significant linear correlation was found between in situ and estimated Chla with the normalized spectral index (R 2 = 0.72, p < 0.005; t-test).The relative error of the model for the validation dataset ranged from 0.4% to 64.5% with a mean absolute percent error of 27.1% (RMSE = 15.01 µg/L).The good agreement showed that the Chla can be used to estimate the primary production of Lake Taihu.
(3) The euphotic depth (Z eu ) in this study was calculated using the relationship between the euphotic depth and the PAR diffuse attenuation coefficient (K d (PAR)) as follows: In this study, the K d (PAR) can be derived using a estimation model developed for Lake Taihu using MODIS-Aqua remote sensing reflectance at the 645 nm band (R rs (645)) (Unpublished data) as follows.
This model has been validated with independent random samples, and indicated a good performance for Taihu K d (PAR) estimation (R 2 = 0.71; p < 0.001, t-test; RMSE = 1.06 m −1 ; MAPE = 25%).It is noted that there will be such a situation arose in which estimation euphotic zone depth exceeded bottom depth and if so, we use the bottom depth to replace the estimation euphotic zone depth [39].
(4) D irr can be calculated based on the time and position of the water column.In this study, northern, southern and central regions of the Lake Taihu were selected to calculate the day length of each day.
(5) E 0 data from Lake Taihu can be calculated according to the following equation from Zhang and Qin [40]: where E 0 is monthly mean total daily PAR and E* is the vapor pressure corrected by air pressure, where P 0 and P are the air pressure at sea level and at the monitoring station, respectively, E is the monthly mean vapor pressure and Q is the monthly mean daily global solar radiation.Q can be calculated according to the percentage of sunshine as follows [41]: where Q is the monthly mean total daily global solar radiation, Q 0 is the monthly mean total daily extra-terrestrial radiation outside the Earth's atmosphere and s 1 is the percentage of actual sunshine expressed relative to the duration of possible sunshine.This model has been validated with independent random samples, and indicated a good performance for Taihu E 0 estimation (R 2 = 0.83, p < 0.001; t-test) With the VGPM model and the Chla data, lake surface photosynthetic available radiation, day length, optimal rate of daily carbon fixation and the euphotic depth, daily primary production in Lake Taihu from 2003 to 2013 has been calculated using the Equation (2).

Statistical Analysis and Accuracy Assessment
Statistical analysis, including calculation of the maximum, minimum, mean, median values and linear and non-linear regressions were performed using the Statistical Program for Social Sciences (SPSS) 17.0 software.The significance level was reported to be significant (p ≤ 0.05) or not significant (p > 0.05).The fitting determination coefficient (R 2 ) was calculated through linear regression using theOriginPro8.5 software.
The accuracy of the phytoplankton primary production estimation was assessed for the two validation datasets by comparing the estimation with the measured primary production values.This comparison was quantified by means of the percent difference (MNB) between the estimated primary production (PP esti ) and analytically measured phytoplankton primary production (PP meas ): Systematic and randomized errors were characterized by the root mean square error (RMSE), and Mean absolute percentage error (MAPE).These metrics were defined as follows:

Validation of VGPM
In theory, the estimated phytoplankton primary production through satellite data and in situ measured phytoplankton primary production should be consistent within a limited period.We set the criterion for matching the satellite data and the in situ observations to ≤12 h (the time interval between the in situ measurements and the corresponding estimated data), according to the hourly variations of the process being measured, to minimize the effects of the temporal difference between the in situ and estimated phytoplankton primary production.Among all collected samples, there were 63 pairs of in situ data and estimated phytoplankton primary production that met the matching criterion.The accuracy of the VGPM was then evaluated through the comparison of the estimated phytoplankton primary production, and in situ phytoplankton primary production is shown in Figure 2. The results show that these values were in good agreement, with a highly significant linear relationship (R 2 = 0.753, RMSE = 384.68mg•C•m −2 •d −1 , MAPE = 40.5%).For all data sets, the minimum and mean of MNB were 1.3% and 41.7%, respectively, between estimated phytoplankton primary production and in situ measurements.The good performance of the VGPM is encouraging and clearly demonstrates that the VGPM could be used for remotely sensed estimations of phytoplankton primary production in the turbid waters of Lake Taihu.

Temporal-Spatial Distribution of Phytoplankton Primary Production
The daily average data were considered on a monthly, seasonal and annual basis, with the seasons defined as follows: winter, December-February; spring, March-May; summer, June-August; and autumn, September-November [23].

Annual Variation of Phytoplankton Primary Production
The annual means of the MODIS-derived phytoplankton primary production for Lake Taihu are presented in Figure 3, demonstrating inter-annual changes in the phytoplankton primary production from 2003 to 2013.The mean phytoplankton primary production of entire Lake Taihu ranged from 952.84 ± 759.36 to 1259.39 ± 859.71 mg•C•m −2 •d −1 .The mean daily phytoplankton primary production in Lake Taihu (Open area, Meiliang Bay, Zhushan Bay) from 2003 to 2013 is 1094.06 ± 720.74 C/m 2 /d, which is far higher than the average of case I water and some case II water, such as Lake Michigan and Lake Huron [42,43].In detail, the phytoplankton primary production of entire Lake Taihu in 2007 (1259.39  ), revealing that Lake Taihu experienced moderate spatial and inter-annual variations in phytoplankton primary production.It is worth noting that the highest phytoplankton primary production value was recorded in 2007.This was probably due to abnormal weak cold air in farther northern part of the region, which caused high temperature during the winter and spring of 2007, and the abnormal lack of precipitation may have jointly contributed to the monthly variations in 2007 [44].The phytoplankton primary production in Lake Taihu experienced three markedly

Temporal-Spatial Distribution of Phytoplankton Primary Production
The daily average data were considered on a monthly, seasonal and annual basis, with the seasons defined as follows: winter, December-February; spring, March-May; summer, June-August; and autumn, September-November [23].

Annual Variation of Phytoplankton Primary Production
The annual means of the MODIS-derived phytoplankton primary production for Lake Taihu are presented in Figure 3, demonstrating inter-annual changes in the phytoplankton primary production from 2003 to 2013.The mean phytoplankton primary production of entire Lake Taihu ranged from 952.84 ± 759.36 to 1259.39 ± 859.71 mg•C•m −2 •d −1 .The mean daily phytoplankton primary production in Lake Taihu (Open area, Meiliang Bay, Zhushan Bay) from 2003 to 2013 is 1094.06 ± 720.74 C/m 2 /d, which is far higher than the average of case I water and some case II water, such as Lake Michigan and Lake Huron [42,43].In detail, the phytoplankton primary production of entire Lake Taihu in 2007 (1259.39  ), revealing that Lake Taihu experienced moderate spatial and inter-annual variations in phytoplankton primary production.It is worth noting that the highest phytoplankton primary production value was recorded in 2007.This was probably due to abnormal weak cold air in farther northern part of the region, which caused high temperature during the winter and spring of 2007, and the abnormal lack of precipitation may have jointly contributed to the monthly variations in 2007 [44].The phytoplankton primary production in Lake Taihu experienced three markedly

Monthly Variation of Phytoplankton Primary Production
Maps of monthly phytoplankton primary production derived from the VGPM are shown in Figures 6 and 7, demonstrating the inter-monthly changes in the phytoplankton primary production from 2003 to 2013.The highest monthly mean daily phytoplankton primary production in Lake Taihu during the study interval from 2003 to 2013 was measured in August, with an average value of 1675.49± 710.36 mg•C•m −2 •d −1 , nearly 7 times higher than the lowest value of 296.01 ± 130.83 mg•C•m −2 •d −1 in January.Some interesting characteristics demonstrated that there is a clear bimodality in character in the monthly variability in Lake Taihu.Specifically, the monthly distribution of phytoplankton primary production increased from January to May, then decreased in June or July, and finally reached highest value in August or September (Figure 7).

Spatial Variation of Phytoplankton Primary Production
The MODIS-derived phytoplankton primary production data from 2003 to 2013 were averaged to calculate the regional phytoplankton primary production distribution for Lake Taihu.It can be concluded that the decrease of phytoplankton primary production in Lake Taihu from the bay of the lake to the open area is in accordance with the basic distribution principle of the spatial variability (Figures 4 and 6).The area around the estuary such as Meiliang Bay and Zhushan Bay are the regions where the maximum value of phytoplankton primary production was recorded.In the district of the Meiliang Bay and Zhushan Bay, the annual mean daily phytoplankton primary production was 2205.47 ± 1397.98 mg•C•m −2 •d −1 , 1551.39 ± 851.36 mg•C•m −2 •d −1 , respectively, which was nearly 2 times and 1.5 times higher than the lowest value of 1019.25 ± 721.23 was recorded in open area (Table 1).This is in agreement with previous observations of productivity in Lake Taihu [24,26].However, there was no significant spatial variation in the monthly mean daily phytoplankton primary production in winter, especially in January and February, possibly due to the degradation of algae in the whole lake (Figure 6).

Monthly Variation of Phytoplankton Primary Production
Maps of monthly phytoplankton primary production derived from the VGPM are shown in Figures 6 and 7, demonstrating the inter-monthly changes in the phytoplankton primary production from 2003 to 2013.The highest monthly mean daily phytoplankton primary production in Lake Taihu during the study interval from 2003 to 2013 was measured in August, with an average value of 1675.49± 710.36 mg•C•m −2 •d −1 , nearly 7 times higher than the lowest value of 296.01 ± 130.83 mg•C•m −2 •d −1 in January.Some interesting characteristics demonstrated that there is a clear bimodality in character in the monthly variability in Lake Taihu.Specifically, the monthly distribution of phytoplankton primary production increased from January to May, then decreased in June or July, and finally reached highest value in August or September (Figure 7).

Spatial Variation of Phytoplankton Primary Production
The MODIS-derived phytoplankton primary production data from 2003 to 2013 were averaged to calculate the regional phytoplankton primary production distribution for Lake Taihu.It can be concluded that the decrease of phytoplankton primary production in Lake Taihu from the bay of the lake to the open area is in accordance with the basic distribution principle of the spatial variability (Figures 4 and 6).The area around the estuary such as Meiliang Bay and Zhushan Bay are the regions where the maximum value of phytoplankton primary production was recorded.In the district of the Meiliang Bay and Zhushan Bay, the annual mean daily phytoplankton primary production was 2205.47 ± 1397.98 mg•C•m −2 •d −1 , 1551.39 ± 851.36 mg•C•m −2 •d −1 , respectively, which was nearly 2 times and 1.5 times higher than the lowest value of 1019.25 ± 721.23 mg•C•m −2 •d −1 that was recorded in open area (Table 1).This is in agreement with previous observations of productivity in Lake Taihu [24,26].However, there was no significant spatial variation in the monthly mean daily phytoplankton primary production in winter, especially in January and February, possibly due to the degradation of algae in the whole lake (Figure 6).

Correlation between the Primary Production with Chla and Nutrient Concentrations
Previous studies have suggested that phytoplankton biomass is the most important factor affecting the temporal variations of phytoplankton primary production [45][46][47][48][49][50][51][52].The linear relationship between Chla and phytoplankton primary production demonstrates that the Chla is strongly positively correlated with phytoplankton primary production in Lake Taihu (Figure 8).Deng reported that Chla was low in winter and early spring, and increased in late spring to a

Correlation between the Primary Production with Chla and Nutrient Concentrations
Previous studies have suggested that phytoplankton biomass is the most important factor affecting the temporal variations of phytoplankton primary production [45][46][47][48][49][50][51][52].The linear relationship between Chla and phytoplankton primary production demonstrates that the Chla is strongly positively correlated with phytoplankton primary production in Lake Taihu (Figure 8).Deng reported that Chla was low in winter and early spring, and increased in late spring to a maximum in late summer and early autumn; the algae bloom began in May or June, and lasted until October.This finding is consistent with the recorded temporal variations in phytoplankton primary production in Lake Taihu [53].The spatial variation in phytoplankton primary production was also consistent with the variations in Chla; the highest productivity was observed near Zhushan Bay, where Chla was also the highest among all the study areas (Figure 9).maximum in late summer and early autumn; the algae bloom began in May or June, and lasted until October.This finding is consistent with the recorded temporal variations in phytoplankton primary production in Lake Taihu [53].The spatial variation in phytoplankton primary production was also consistent with the variations in Chla; the highest productivity was observed near Zhushan Bay, where Chla was also the highest among all the study areas (Figure 9).The phenomenon in which both higher phytoplankton primary production and higher Chla appeared in the bay is related to nutrient distribution.The nutrients mainly affect phytoplankton biomass to control the level of the phytoplankton primary production of phytoplankton.In fact, previous studies have showed that total nutrients, especially the total phosphorus, play a key role in determining phytoplankton in Lake Taihu [25,54,55].Higher total nitrogen and total phosphorus concentrations at ZhuShan Bay and MeilLiang Bay resulted in a marked increase in Chla.As a highly eutrophic lake, high concentrations of nutrients in Lake Taihu provide good conditions for phytoplankton to grow.The high phytoplankton biomass keeps phytoplankton primary production at a very high level, and the spatial distribution of nutrients can affect the spatial variation in phytoplankton primary production of phytoplankton through affecting the Chla.Some studies have shown that the spatial distribution of nutrient concentrations in Lake Taihu was consistent with the variations in Chla [56,57].Based on the above analysis, it can be found that the influence of maximum in late summer and early autumn; the algae bloom began in May or June, and lasted until October.This finding is consistent with the recorded temporal variations in phytoplankton primary production in Lake Taihu [53].The spatial variation in phytoplankton primary production was also consistent with the variations in Chla; the highest productivity was observed near Zhushan Bay, where Chla was also the highest among all the study areas (Figure 9).The phenomenon in which both higher phytoplankton primary production and higher Chla appeared in the bay is related to nutrient distribution.The nutrients mainly affect phytoplankton biomass to control the level of the phytoplankton primary production of phytoplankton.In fact, previous studies have showed that total nutrients, especially the total phosphorus, play a key role in determining phytoplankton in Lake Taihu [25,54,55].Higher total nitrogen and total phosphorus concentrations at ZhuShan Bay and MeilLiang Bay resulted in a marked increase in Chla.As a highly eutrophic lake, high concentrations of nutrients in Lake Taihu provide good conditions for phytoplankton to grow.The high phytoplankton biomass keeps phytoplankton primary production at a very high level, and the spatial distribution of nutrients can affect the spatial variation in phytoplankton primary production of phytoplankton through affecting the Chla.Some studies have shown that the spatial distribution of nutrient concentrations in Lake Taihu was consistent with the variations in Chla [56,57].Based on the above analysis, it can be found that the influence of The phenomenon in which both higher phytoplankton primary production and higher Chla appeared in the bay is related to nutrient distribution.The nutrients mainly affect phytoplankton biomass to control the level of the phytoplankton primary production of phytoplankton.In fact, previous studies have showed that total nutrients, especially the total phosphorus, play a key role in determining phytoplankton in Lake Taihu [25,54,55].Higher total nitrogen and total phosphorus concentrations at ZhuShan Bay and MeilLiang Bay resulted in a marked increase in Chla.As a highly eutrophic lake, high concentrations of nutrients in Lake Taihu provide good conditions for phytoplankton to grow.The high phytoplankton biomass keeps phytoplankton primary production at a very high level, and the spatial distribution of nutrients can affect the spatial variation in phytoplankton primary production of phytoplankton through affecting the Chla.Some studies have shown that the spatial distribution of nutrient concentrations in Lake Taihu was consistent with the variations in Chla [56,57].Based on the above analysis, it can be found that the influence of phytoplankton primary production of phytoplankton in Lake Taihu by nutrients is highly significant, and the concentration and spatial distribution of nutrients are major factors responsible for the level and spatial variations in phytoplankton primary production in Lake Taihu.Additionally, the temporal variations of phytoplankton primary production in Lake Taihu appear to show a correlation with the concentrations of nutrients, especially the total phosphorus concentration (Figure 10).Therefore, one of the reasons for the temporal-spatial variations of phytoplankton primary production in Lake Taihu is the temporal-spatial variation in Chla, originated by different nutrient concentrations.
Remote Sens. 2017, 9, 195 13 of 18 phytoplankton primary production of phytoplankton in Lake Taihu by nutrients is highly significant, and the concentration and spatial distribution of nutrients are major factors responsible for the level and spatial variations in phytoplankton primary production in Lake Taihu.Additionally, the temporal variations of phytoplankton primary production in Lake Taihu appear to show a correlation with the concentrations of nutrients, especially the total phosphorus concentration (Figure 10).Therefore, one of the reasons for the temporal-spatial variations of phytoplankton primary production in Lake Taihu is the temporal-spatial variation in Chla, originated by different nutrient concentrations.
Figure10.Time series of total nitrogen and total phosphorusin 2013, and the correlation between phytoplankton primary production and total phosphorus.

Correlation between the Primary Production and Lake Water Temperature
Although influenced by nutrient concentrations, the temporal variations of phytoplankton primary production were also likely affected by water temperature.Chla and phytoplankton primary production peaks appeared during different months (Figure 8b), suggesting that the water temperature was likely to be responsible for this difference.The P B opt is potentially the most important variable in the estimation of phytoplankton primary production, and one of the principal approaches for estimating P B opt is to define predictive relationships between P B opt and lake surface temperature by a high-order polynomial.The suitable water temperature for phytoplankton growth was approximately 20 °C [34]; a temperature that is too high or too low will suppress the phytoplankton photosynthesis rate, therefore maximal P B opt values occur from18-23 °C (Figure 11).In January, the water temperature was the lowest, which corresponded to a very low value for P B opt , greatly inhibiting phytoplankton primary production.As temperatures increased, the phytoplankton primary production dropped briefly due to the value of inhibition of P B opt result from the high water temperature of 23.6 °C in June and 26.6 °C in July.In August and September, although the temperature remained high, the phytoplankton primary production began to gradually rise to its maximum due to the explosive growth of the phytoplankton biomass.In November and December, as with the decrease of phytoplankton biomass, the photoperiod and P B opt , phytoplankton primary production began to drop quickly.Based on the above analysis, the water temperature in Lake Taihu seem to influence the temporal distribution by affecting the maximum optimal rate of carbon fixation of water mass.Additionally, the mutual movement between phytoplankton biomass and P B opt was the intrinsic c motivation of the temporal variation mechanism of phytoplankton primary production in Lake Taihu.

Correlation between the Primary Production and Lake Water Temperature
Although influenced by nutrient concentrations, the temporal variations of phytoplankton primary production were also likely affected by water temperature.Chla and phytoplankton primary production peaks appeared during different months (Figure 8b), suggesting that the water temperature was likely to be responsible for this difference.The P B opt is potentially the most important variable in the estimation of phytoplankton primary production, and one of the principal approaches for estimating P B opt is to define predictive relationships between P B opt and lake surface temperature by a high-order polynomial.The suitable water temperature for phytoplankton growth was approximately 20 • C [34]; a temperature that is too high or too low will suppress the phytoplankton photosynthesis rate, therefore maximal P B opt values occur from18-23 • C (Figure 11).In January, the water temperature was the lowest, which corresponded to a very low value for P B opt , greatly inhibiting phytoplankton primary production.As temperatures increased, the phytoplankton primary production dropped briefly due to the value of inhibition of P B opt result from the high water temperature of 23.6 • C in June and 26.6 • C in July.In August and September, although the temperature remained high, the phytoplankton primary production began to gradually rise to its maximum due to the explosive growth of the phytoplankton biomass.In November and December, as with the decrease of phytoplankton biomass, the photoperiod and P B opt , phytoplankton primary production began to drop quickly.Based on the above analysis, the water temperature in Lake Taihu seem to influence the temporal distribution by affecting the maximum optimal rate of carbon fixation of water mass.Additionally, the mutual movement between phytoplankton biomass and P B opt was the intrinsic c motivation of the temporal variation mechanism of phytoplankton primary production in Lake Taihu.

The Effects of Wind-Driven Sediment Resuspension on the Spatial Distribution of Phytoplankton Primary Production
Wind-induced sediment resuspension can affect phytoplankton primary production by influencing euphotic depth and nutrient availability.To investigate the effect of wind-driven sediment resuspension on the spatial distribution of phytoplankton primary production in Lake Taihu, we performed a correlated analysis between the MODIS-derived phytoplankton primary production and Kd(PAR).The results (Figure 12) showed the existence of a negative linear correlation between the phytoplankton primary production and the Kd(PAR) in the open area, indicating that the process of primary production in the open area was easily disturbed by wind-derived waves.This result is in agreement with our finding that lower phytoplankton primary production was found in the open area compared to the other regions of Lake Taihu.The dominant wind directions in Lake Taihu were found to be NNW and ESE, resulting in longer wind fetches for northern and western parts of Lake Taihu.Longer wind fetches yield stronger wind forces and thus more sediment resuspension in the northern and western regions of Lake Taihu.In the littoral zones and bays, such as Zhushan Bay and Meiliang Bay, the wind fetches are shorter than in the open area.Therefore, wind-driven sediment resuspension, expressed using disturbance index, could explain the spatial variations in the phytoplankton primary production in Lake Taihu.

The Effects of Wind-Driven Sediment Resuspension on the Spatial Distribution of Phytoplankton Primary Production
Wind-induced sediment resuspension can affect phytoplankton primary production by influencing euphotic depth and nutrient availability.To investigate the effect of wind-driven sediment resuspension on the spatial distribution of phytoplankton primary production in Lake Taihu, we performed a correlated analysis between the MODIS-derived phytoplankton primary production and K d (PAR).The results (Figure 12) showed the existence of a negative linear correlation between the phytoplankton primary production and the K d (PAR) in the open area, indicating that the process of primary production in the open area was easily disturbed by wind-derived waves.This result is in agreement with our finding that lower phytoplankton primary production was found in the open area compared to the other regions of Lake Taihu.The dominant wind directions in Lake Taihu were found to be NNW and ESE, resulting in longer wind fetches for northern and western parts of Lake Taihu.Longer wind fetches yield stronger wind forces and thus more sediment resuspension in the northern and western regions of Lake Taihu.In the littoral zones and bays, such as Zhushan Bay and Meiliang Bay, the wind fetches are shorter than in the open area.Therefore, wind-driven sediment resuspension, expressed using disturbance index, could explain the spatial variations in the phytoplankton primary production in Lake Taihu.

The Effects of Wind-Driven Sediment Resuspension on the Spatial Distribution of Phytoplankton Primary Production
Wind-induced sediment resuspension can affect phytoplankton primary production by influencing euphotic depth and nutrient availability.To investigate the effect of wind-driven sediment resuspension on the spatial distribution of phytoplankton primary production in Lake Taihu, we performed a correlated analysis between the MODIS-derived phytoplankton primary production and Kd(PAR).The results (Figure 12) showed the existence of a negative linear correlation between the phytoplankton primary production and the Kd(PAR) in the open area, indicating that the process of primary production in the open area was easily disturbed by wind-derived waves.This result is in agreement with our finding that lower phytoplankton primary production was found in the open area compared to the other regions of Lake Taihu.The dominant wind directions in Lake Taihu were found to be NNW and ESE, resulting in longer wind fetches for northern and western parts of Lake Taihu.Longer wind fetches yield stronger wind forces and thus more sediment resuspension in the northern and western regions of Lake Taihu.In the littoral zones and bays, such as Zhushan Bay and Meiliang Bay, the wind fetches are shorter than in the open area.Therefore, wind-driven sediment resuspension, expressed using disturbance index, could explain the spatial variations in the phytoplankton primary production in Lake Taihu.

Conclusions
The efficient and accurate of the distribution and variation of phytoplankton primary production is important for ecosystem and water resource management of lake.In this study, we addressed one technical challenge and understanding of the long-term primary production distribution variations in Lake Taihu, from 2003 to 2013.The technical challenge is that the VGPM model is a good model for phytoplankton primary production estimation in Lake Taihu (R 2 = 0.753, RMSE = 384.68mg•C•m −2 •d −1 ).It holds considerable promise and could be useful for primary production estimation in large lakes.
The annual mean daily phytoplankton primary production was 1094.06±720.74mg•C•m −2 •d −1 in Lake Taihu from 2003 to 2013.The spatial pattern in daily mean phytoplankton primary production was similar to that of Chla, decreasing from bay to open area.The phytoplankton primary production exhibited typical seasonal variability over the entire Lake Taihu, and there is a clear bimodality characteristic of phytoplankton primary production in the monthly variability of Lake Taihu.
Temporal-spatial variations in phytoplankton primary production appeared to be mainly driven by Chla variations, and the spatial distribution of nutrients affected the spatial variation of phytoplankton primary production of through affecting Chla.High nutrient concentrations determined the high Chla and then resulted in the high phytoplankton primary production of Lake Taihu.Other contributing factors include water temperature, surface PAR, photoperiod and euphotic depth responsible for the spatial variation of the phytoplankton primary production.In addition, as a large, shallow lake, the wind-driven sediment resuspension plays an important role in spatial distribution of phytoplankton primary production, especially in the open area.

Figure 1 .
Figure 1.Maps showing location of Lake Taihu and the sampling sites.The in situ primary production were measured at two sites, THL04 (31.44583°N, 120.18883°E) and THL08 (31.24816°N, 120.10062°E), from 2005 to 2013.The in situ nutrient concentration in 2013 including total nitrogen and total phosphorus were collected from 32 sites THL00-THL32 (except THL02).

Figure 1 .
Figure 1.Maps showing location of Lake Taihu and the sampling sites.The in situ primary production were measured at two sites, THL04 (31.44583 • N, 120.18883 • E) and THL08 (31.24816 • N, 120.10062 • E), from 2005 to 2013.The in situ nutrient concentration in 2013 including total nitrogen and total phosphorus were collected from 32 sites THL00-THL32 (except THL02).

Figure 2 .
Figure 2. Relationship between measured and estimated primary production by VGPM.

Figure 2 .
Figure 2. Relationship between measured and estimated primary production by VGPM.

Figure 3 .
Figure 3. Annual variations of phytoplankton primary production in Lake Taihu from 2003 to 2013.

Figure 4 .
Figure 4. Maps of the MODIS-derived mean daily phytoplankton primary production for the four seasons in Lake Taihu generated using the VGPM.

Figure 3 .
Figure 3. Annual variations of phytoplankton primary production in Lake Taihu from 2003 to 2013.

Figure 3 .
Figure 3. Annual variations of phytoplankton primary production in Lake Taihu from 2003 to 2013.

Figure 4 .
Figure 4. Maps of the MODIS-derived mean daily phytoplankton primary production for the four seasons in Lake Taihu generated using the VGPM.

Figure 4 .
Figure 4. Maps of the MODIS-derived mean daily phytoplankton primary production for the four seasons in Lake Taihu generated using the VGPM.

Figure 5 .
Figure 5. Seasonal variation of phytoplankton primary production in each region based on VGPM estimation from 2003 to 2013 in (a) Zhanshan Bay; (b) Meiliang Bay; (c) the open area and (d) the entire lake.

Figure 5 .
Figure 5. Seasonal variation of phytoplankton primary production in each region based on VGPM estimation from 2003 to 2013 in (a) Zhanshan Bay; (b) Meiliang Bay; (c) the open area and (d) the entire lake.

Figure 6 .
Figure 6.Monthly mean daily phytoplankton primary production of Lake Taihu from January to December.

Figure 6 .
Figure 6.Monthly mean daily phytoplankton primary production of Lake Taihu from January to December.

Figure7.
Figure7.Time series of monthly phytoplankton primary production in Lake Taihu from 2003 to 2013.

Figure 7 .
Figure 7. Time series of monthly phytoplankton primary production in Lake Taihu from 2003 to 2013.

Figure 8 .
Figure 8.(a) Relationship between Chla and MODIS derived primary productivity of Lake Taihu from 2003-2013 (b) Monthly mean MODIS derived primary productivity and Chla of Lake Taihu from 2003-2013.

Figure 9 .
Figure 9. (a) Chla distribution in Lake Taihu averaged from all Chla estimates from the MODIS data gathered from 2003-2013 and (b) spatial PP distribution in Lake Taihu based on all PP estimates from the MODIS data gathered from 2003-2013.

Figure 8 .
Figure 8.(a) Relationship between Chla and MODIS derived primary productivity of Lake Taihu from 2003-2013 (b) Monthly mean MODIS derived primary productivity and Chla of Lake Taihu from 2003-2013.

Figure 8 .
Figure 8.(a) Relationship between Chla and MODIS derived primary productivity of Lake Taihu from 2003-2013 (b) Monthly mean MODIS derived primary productivity and Chla of Lake Taihu from 2003-2013.

Figure 9 .
Figure 9. (a) Chla distribution in Lake Taihu averaged from all Chla estimates from the MODIS data gathered from 2003-2013 and (b) spatial PP distribution in Lake Taihu based on all PP estimates from the MODIS data gathered from 2003-2013.

Figure 9 .
Figure 9. (a) Chla distribution in Lake Taihu averaged from all Chla estimates from the MODIS data gathered from 2003-2013 and (b) spatial PP distribution in Lake Taihu based on all PP estimates from the MODIS data gathered from 2003-2013.

Figure 10 .
Figure 10.Time series of total nitrogen and total phosphorusin 2013, and the correlation between phytoplankton primary production and total phosphorus.

Figure 11 .
Figure 11.Correlations between water temperature and maximum carbon fixation rate and phytoplankton primary production.

Figure12.
Figure12.Correlations between PAR diffuse attenuation coefficient and phytoplankton primary production for the regions of (a) the entire Lake Taihu; (b) the open area; (c) Meiliang Bay and (d) Zhushan Bay.

Figure 11 .
Figure 11.Correlations between water temperature and maximum carbon fixation rate and phytoplankton primary production.

Figure 11 .
Figure 11.Correlations between water temperature and maximum carbon fixation rate and phytoplankton primary production.

Figure12.
Figure12.Correlations between PAR diffuse attenuation coefficient and phytoplankton primary production for the regions of (a) the entire Lake Taihu; (b) the open area; (c) Meiliang Bay and (d) Zhushan Bay.

Figure 12 .
Figure 12.Correlations between PAR diffuse attenuation coefficient and phytoplankton primary production for the regions of (a) the entire Lake Taihu; (b) the open area; (c) Meiliang Bay and (d) Zhushan Bay.