Abstract
We developed and validated an empirical model for estimating chlorophyll a concentrations (Chla) in Lake Taihu to generate a long-term Chla and algal bloom area time series from MODIS-Aqua observations for 2003 to 2013. Then, based on the long-term time series data, we quantified the responses of cyanobacterial dynamics to nutrient enrichment and climatic conditions. Chla showed substantial spatial and temporal variability. In addition, the annual mean cyanobacterial surface bloom area exhibited an increasing trend across the entire lake from 2003 to 2013, with the exception of 2006 and 2007. High air temperature and phosphorus levels in the spring can prompt cyanobacterial growth, and low wind speeds and low atmospheric pressure levels favor cyanobacterial surface bloom formation. The sensitivity of cyanobacterial dynamics to climatic conditions was found to vary by region. Our results indicate that temperature is the most important factor controlling Chla inter-annual variability followed by phosphorus and that air pressure is the most important factor controlling cyanobacterial surface bloom formation followed by wind speeds in Lake Taihu.
Similar content being viewed by others
Introduction
Lakes deliver services of enormous global value by acting as indicators and regulators of climate change and by providing crucial resources for humans1,2. However, as human activities and economic development have increased over the last few decades, lake eutrophication has become a serious ecological, environmental and social problem in both industrialized and developing countries3,4,5. Among human-induced changes to lake environments, nutrient over-enrichment (especially of nitrogen and phosphorus) associated with urban, agricultural and industrial development is one of the most notable. This has led to accelerated rates of lake eutrophication6. Lake eutrophication can result in an increase in primary production, which in turn spurs the development of larger and more frequent cyanobacterial blooms7,8. Around the world, numerous eutrophic lakes suffer from severe cyanobacterial blooms, including Lake Taihu9, Lake Chaohu10 and Lake Dianchi11 in China; Lake Erie in the U.S. and Canada5; Lake Winnipeg in Canada12; and Lake Nieuwe Meer in the Netherlands13.
Cyanobacterial blooms can severely stress the ecological structures and functions and aesthetics of lake ecosystems5,14. Large populations of cyanobacteria can increase lake turbidity levels due to high chlorophyll a concentrations (Chla) present at the water surface, thereby decreasing light penetration levels15. In turn, less light is available for submerged aquatic vegetation (SAV), thus suppressing SAV growth and populations16,17. Moreover, severe cyanobacterial blooms can lead to dissolved oxygen depletion at night, in turn decreasing fish populations18. Some cyanobacterial species are known to produce a variety of toxins (toxic peptides and alkaloids)19,20; the ingestion of such cyanotoxins can have serious effects on animal and human health (i.e., liver, digestive and skin diseases; neurological impairment and death)21,22. Therefore, cyanobacterial blooms pose a major threat to the use of lake ecosystems for drinking and irrigation water, fishing and recreational purposes15,20,23.
To manage, control, and treat cyanobacterial blooms, it is crucial to accurately quantify cyanobacterial bloom responses to climate change and nutrients and to identify the mechanisms of these responses. Many recent studies have used field and laboratory data and simulation models to quantitatively and qualitatively assess the role of climate change and nutrients in cyanobacterial surface bloom expansion in inland, coastal and ocean waters4,6,23,24,25,26,27,28,29. However, many of these studies introduce a substantial degree of uncertainty due to the complexities and costs of in situ measurements of cyanobacterial abundance and difficulties associated with generalizing laboratory findings over large spatial and temporal scales30,31. These complex processes interact to produce physical, chemical, and biotic perturbations of lake ecosystem functioning, thus raising the need for new quantitative methods for the detection, quantification, synthesizing and modeling of impacts at broad spatial and temporal resolutions17.
The spatial and temporal coverage provided by satellite remote sensing makes it an attractive tool to use in such studies. Remotely sensed data such as those generated by the Moderate Resolution Imaging Spectroradiometer (MODIS; from the Aqua satellite for 2002-the present) provide opportunities to quantify cyanobacterial biomass and production levels in lake and ocean waters at a fine spatial and temporal resolution. Satellite observations are critical to identifying and removing spatial and temporal aliasing from scarce field measurements when investigating long-term baseline information and trends of cyanobacterial production. To develop a stronger understanding of the effects of human activities and climate change on cyanobacterial blooms, Paerl and Huisman (2009) recommended coupling traditional sampling and remote sensing methods to monitor cyanobacteria-dominated aquatic ecosystems8. Several recent studies have successfully applied long-term remotely sensed data to model the responses of phytoplankton production dynamics to environmental variables in coastal and ocean waters31,32,33,34,35. For example, Vope et al. (2012) used MODIS-derived Chla data for the Mediterranean Sea to establish that phytoplankton populations were responding to physical variables such as sea surface temperatures35, thus identifying mechanisms that link Chla and sea surface temperature. Keith (2014) found a relationship between Chla variations and nutrients in the Neuse and Tar-Pamlico River estuaries, demonstrating that the remote sensing of Chla can support nutrient management32. One widely used proxy for cyanobacterial biomass monitoring is the Chla3,14,36 for inland cyanobacteria-dominated aquatic ecosystems.
Remote sensing techniques have rarely been applied to quantify the responses of cyanobacterial biomass dynamics to climate change and to nutrient variations for inland lake waters. This is partly attributable to the fact that it is challenging to accurately derive Chla (an index of cyanobacterial biomass) values from satellite measurements of optically complex cyanobacteria-dominated inland lakes due to the presence of non-covarying optically active constituents such as colored dissolved matter (CDOM) and organic and inorganic suspended matter37. Several algorithms have been presented that use red to near-infrared bands to estimate Chla values for inland lake waters. These include (1) two-band empirical algorithms (Moses et al. 2009); (2) three- and four-band empirical algorithms36,38,39,40; and (3) maximum peak height (MPH), fluorescence line height (FLH), and maximum chlorophyll index (MCI) algorithms26. Among the satellite sensors with high spatial and temporal resolutions, only the MERIS can generate appropriate data for red to near-infrared band algorithms. MERIS imagery offers a fine spatial (full resolution: 300 m) and spectral resolution and a short revisit time (near daily coverage) and was therefore believed to serve as an ideal satellite image dataset for monitoring Chla in inland lake waters. Unfortunately, the MERIS sensor was discontinued in April of 2012. As an alternative, the MODIS-Aqua sensor was launched in May of 2002; by offering a short revisit time (one day) and relatively high spatial resolution (250 m for the first two bands), it can generate data of a sufficient spatial and temporal resolution for investigating long-term variations in Chla at large scales.
As the third-largest freshwater lake in China, Lake Taihu has a water surface area of 2,338 km2 and a mean water depth of 1.9 m9. 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 (the Yangtze River delta region). The lake supplies water to the approximately 10 million residents of surrounding cities including Wuxi, Suzhou, and Huzhou9. Thus, Lake Taihu’s water quality is vital to local human activities and needs such as drinking, tourism, fishing, and shipping. It also plays a key role in regional ecosystem functioning22. Lake Taihu has become progressively more eutrophic since the 1980 s due to dramatic increases in nutrient loading from urban and agricultural development in its watershed, leading to the frequent formation of cyanobacterial blooms in the spring and summer7,22. Moreover, multi-annual warming trends have worsened these cyanobacterial blooms17. For example, the spring of 2007 cyanobacterial bloom event occurring in this lake caused the contamination of tap water and spurred a drinking water crisis in the city of Wuxi, directly affecting several million people9,22. In our previous study20, we proposed a spectral index derived from the first two bands of MODIS data that can be used to reflect Chla information and to obtain long-term Chla records from MODIS data for Lake Taihu. However, it should be noted that our previous study did not fully develop a Chla estimation model and that only a simple relationship was found between the spectral index and Chla; this is attributable to the fact that the main aim of that study20 was to estimate microcystin concentrations by means of Chla values. Thus, the objectives of the present study are: (1) to complete the development of the MODIS Chla estimation algorithm initially proposed through our previous study20, (2) to generate a long-term Chla record based on MODIS-Aqua observations to help characterize long-term variability trends of cyanobacterial biomass (Chla) and blooms in Lake Taihu, and (3) to determine how cyanobacteria respond to meteorological factors and nutrient variability levels by examining long-term Chla records from MODIS-Aqua data.
Materials and Methods
Description of sampling sites
Field observations of nutrients including total nitrogen (TN) and phosphorus (TP) and of Chla were collected monthly from pre-defined sites across Lake Taihu from 2003 to 2013. Eleven years of continuous monthly sampling at these sampling sites yielded a total of 2,432 water samples. The datasets were collected as part of a long-term monitoring project managed by the Taihu Laboratory for Lake Ecosystem Research (TLLER) of the Nanjing Institute of Geography and Limnology at the Chinese Academy of Sciences. From each sampling site, surface water samples were collected at four water depths (10 cm, 50 cm, 100 cm, and 150 cm), were added to 2 L acid-washed bottles and were then stored on ice while working in the field. All samples were transported to the TLLER on the day of collection. The samples included a broad range of biogeochemical and optical variability levels for Lake Taihu, including water samples with optical properties dominated by strong terrestrial inputs to samples made during strong cyanobacterial bloom events. Further information on the spatial distribution of the sampling sites in Lake Taihu can be found in our previous paper20. To better understand spatial variations in cyanobacterial biomass levels in aquatic environments of Lake Taihu, we divided the lake into six parts: Meiliang Bay, Zhushan Bay, Gonghu Bay, open area, Xukou Bay, and East Lake Taihu (Figure S1).
Measurements of Chla, TN and TP concentrations
We recorded three water quality parameters: Chla, TP, and TN. We used Whatman GF/F fiberglass filters with an average pore size of 0.7 μm to collect algal particles from the water samples. Chla pigments were extracted using 90% ethanol at 80 °C and then were spectrophotometrically analyzed to measure their absorption coefficients at 665 nm and 750 nm. We then calculated Chla concentrations from the absorption coefficients at those two wavelengths41. TP and TN concentrations were measured via combined persulfate digestion42. We fixed phytoplankton samples with Lugol’s iodine solution and sedimented them for 48 h prior to counting them using a microscope.
Image data processing and FAI (Floating algae index) products
The MODIS-Aqua data have been freely available since 2002, and they have a maximum spatial resolution of 250 m (bands 1 and 2) and a very short revisit interval (1 image/day). MODIS-Aqua L-0 data for January of 2003 to of December 2013 (more than 4,000 images) were downloaded from NASA’s Goddard Space Flight Center website (GSFC, http://oceancolor.gsfc.nasa.gov/). As the daily MODIS-Aqua images often included clouds, cloud shadows, or thick aerosols, not all of the images downloaded were used in this study. We selected 1,109 cloud-free images of Lake Taihu for January of 2003 to December of 2013. These images were processed to Level-1 (calibrated spectral radiance) using the SeaDAS software package (version 6.0). Corrections for ozone, water vapor absorption, and molecular (Rayleigh) scattering were then performed following methods proposed through previous studies7, and FAI values were calculated for recognized and delineated floating algae in Lake Taihu using the methods of Hu et al. (2010b) (i.e., FAI > −0.004)7.
Comparisons between MODIS-Aqua data and in situ measurements
To minimize effects of the temporal difference between the field and MODIS-Aqua data, we created a criterion for matching satellite and in situ data to ≤3 hours (the time interval between in situ and corresponding MODIS-Aqua measurements). Our criterion yielded 250 matched pairs of Rayleigh-corrected Rrc and in situ Chla measurements; for these data pairs, we also required that the MODIS images and in situ data were located at the same points in space (the same pixel). The matching samples were distributed across the entire lake and across four seasons, thus representing an overall relationship between MODIS-Aqua and in situ Chla measurements for Lake Taihu. The 250 in situ Chla-Rayleigh-corrected Rrc data pairs were used to develop and validate the Chla estimation model. First, we used half (125) of the 250 in situ Chla-Rrc data pairs for model development and the other half for model validation. We randomly and evenly divided the 250 in situ Chla-Rrc data pairs into two parts ten times. As a result, we produced nine groups of dataset pairs (125–125) to test the stability of the developed model. One group of dataset pairs was created for model development and validation and other nine groups of dataset pairs were used to test stability of the developed model. For this study, we developed an empirical model for the retrieval of Chla from MODIS-Aqua data for Lake Taihu.
Meteorological data
Daily wind speed (m/s), temperature (°C), and atmospheric pressure (hpa) data for 1956 to 2013 were obtained from the Dongshan meteorological station (31°4′, 120°26′E) and were downloaded from the China Meteorological Data Sharing Service System (http://cdc.nmic.cn).
Statistical analysis and accuracy assessment
Statistical Program for Social Sciences (SPSS 17.0) software (version 17.0) was used to perform the statistical analyses. Pearson’s correlation analysis method was used to investigate relationships between the variables. Significance levels are reported as significant (p < 0.05) or not significant (p > 0.05).
We assessed the accuracy of the model developed using relative error (RE), mean absolute percent error (MAPE), and root-mean-square error (RMSE) values between the measured and predicted values using the following equations:



where N is the number of samples, and Ymeasured and Yestimated are the measured and estimated values, respectively.
CARTs (Classification and regression trees) were used to examine relationships between cyanobacterial dynamics and meteorological and nutrient factors to investigate the relative importance of various factors in controlling cyanobacterial dynamics. CARTs can explain the variability in a single dependent variable corresponding to one or more explanatory variables by splitting data recursively based on the most influential independent variable43. CARTs, as efficient tools for extracting key variables and thresholds from a multivariate dataset, have been widely used in a variety of fields such as environmental studies and ecology43. CART analyses were performed in this study using SPSS 17.0 software.
Results
Model development and validation
The spectral index [(EXP(Rrc(645)) − EXP(Rrc(859)))/(EXP(Rrc(645)) + EXP(Rrc(859)))] proposed through our previous study was used to construct the Chla estimation model20. We used linear, logarithmic, exponential, and quadratic functions to model in situ Chla and the spectral index. Among these functions, the linear function offered the highest degree of modeling precision, the highest correlation coefficient (r = −0.85), and the lowest MAPE (24%) and RMSE (12.44 μg/L) values (Fig. 1(a)); the linear function is as follows:

where IndexMODIS = [(EXP(Rrc(645)) − EXP(Rrc(859)))/(EXP(Rrc(645)) + EXP(Rrc(859)))], and Rrc(645) and Rrc(859) are atmospherically Rayleigh-corrected MODIS-Aqua data at 645 nm and 859 nm, respectively.
To assess the performance of the proposed Chla estimation model, we used independent validation data from 125 matched in situ Chla-Rrc data pairs (Chla: 6.56–113.66 μg/L). Without adjusting the spectral index and re-parameterization, the proposed model generally performed well for Chla estimations (Fig. 1(b)). The RE of the model for the validation dataset ranged from 0.4% to 64.5% with a MAPE of 27.1% (RMSE = 15.01 μg/L). The RE values of 40% and 60% of the samples fell below 20% and 30%, respectively. The in situ Chla measurements and Chla values estimated using the proposed model with the normalized spectral index [(EXP(Rrc(645)) − EXP(Rrc(859)))/(EXP(Rrc(645)) + EXP(Rrc(859)))] showed good agreement and a significant linear correlation (p < 0.005; t-test). In addition, the measured and estimated Chla values were evenly distributed along a 1:1 line (Fig. 1(b)). These results suggest that the model developed from the normalized spectral index performs satisfactorily and can be used to derive Chla. We separately used the remaining nine groups of dataset pairs to develop Chla estimation models and to validate the developed model. We found no significant differences in the slopes and intercepts of the nine linear regression equations. The slopes and intercepts varied from −1,450.1 to −1,457.8 and from 53.11 to 72.86, respectively. MAPE values of the nine models for the corresponding validation datasets ranged from 25.5% to 30.6%. This suggests that the developed model is very stable. The model can thus be used to quantify spatial and temporal Chla distributions for Lake Taihu.
Cyanobacterial dynamics
Two indices of Chla and the area of cyanobacterial surface blooms were used to characterize long-term trends in the cyanobacterial dynamics of Lake Taihu. The proposed model (Equation (4)) and the FAI proposed by Hu et al. (2010b)7 were used, respectively, to derive a Chla time series and the cyanobacterial surface bloom areas for this lake from all available MODIS-Aqua data collected from January of 2003 to December of 2013. Several sub-regions of Lake Taihu, including Gonghu Bay, Xukou Bay, and East Lake Taihu, are known to be covered in aquatic plants such as weeds, reeds, and other macrophytes. These aquatic plants can significantly change the reflectance characteristics of a lake’s surface, and thus Chla data and cyanobacterial bloom areas cannot be accurately derived from remote optical satellite data in such areas20. Thus, in the following, an “entire lake” refers to the entire area of Lake Taihu with the exception of these three areas (Gonghu Bay, Xukou Bay, and East Lake Taihu). Chla images of yearly, monthly, and seasonal means for Lake Taihu were developed using the arithmetic means of all MODIS-derived Chla products for between 2003 and 2013.
Overall, the lake water shows relatively high Chla values during all seasons. Chla exhibits strong seasonal variability across the entire lake (Figs 2 and 3). Overall, Chla is substantially higher during the summer (June-August) and autumn (September-November) than in the spring (March-May) and winter (December-February) (p < 0.005; t-test).
MODIS-Aqua derived (2003–2013) monthly mean Chla for Lake Taihu for January to December. The figure was derived from MODIS-Aqua data using ENVI 5.0 (2013) (https://www.harris.com/solution/envi).
Maps of the MODIS-Aqua derived Chla for all four seasons in Lake Taihu generated using the proposed model. The figure was derived from MODIS-Aqua data using ENVI 5.0 (2013) (https://www.harris.com/solution/envi).
Mean seasonal Chla values of the entire lake for spring to winter are 27.86 μg/L (standard deviation (SD) = 4.43 μg/L), 38.07 μg/L (SD = 3.81 μg/L), 37.32 μg/L (SD = 3.68 μg/L), and 27.05 μg/L (SD = 2.43 μg/L), respectively. The highest monthly mean Chla value for the entire lake (42.67 μg/L) was reached in August and the lowest monthly mean Chla value (22.92 μg/L) was reached in February. The highest and the lowest monthly mean Chla values for both Meiliang Bay (with highest and lowest values of 58.42 μg/L and 33.04 μg/L, respectively) and Zhushan Bay (with highest and lowest values of 72.05 μg/L and 36.03 μg/L, respectively) were reached in August and January, respectively. In contrast, for open areas, the highest (September) and lowest (February) monthly mean Chla values (with highest and lowest values of 21.12 μg/L and 40.81 μg/L, respectively) were reached one month later than those of the other two areas.
Substantial interannual variability in Chla from 2003 to 2013 was found (Fig. 4). The lowest annual mean Chla value for the entire lake was 30.81 μg/L in 2009, and the highest Chla value reached was 37.28 μg/L, which occurred in 2007 (Fig. 4).
Annual mean Chla distributions of Lake Taihu for 2003 to 2013. The figure was derived from MODIS-Aqua data using ENVI 5.0 (2013) (https://www.harris.com/solution/envi).
Generally, Chla dynamics in Lake Taihu were characterized by three phases from 2003 to 2013. Chla values increased sharply from 2003 to 2007 at a rate of 2.03 μg/L/year, after which they decreased from 2007 to 2009 (at a rate of 3.2 μg/L/year) and then increased slightly from 2009 to 2013.
Spatially, Chla values of Zhushan Bay were higher than those of the other two sub-regions, and the lowest Chla value was reached in the middle of the open areas (Fig. 5). In Meiliang Bay and Zhushan Bay, Chla values decreased from the inner to the outer parts of the regions (Fig. 5).
Chla spatial distribution for Lake Taihu averaged from all Chla estimates based on MODIS-Aqua data gathered from 2003–2013. The figure was derived from MODIS-Aqua data using ENVI 5.0 (2013) (https://www.harris.com/solution/envi).
For cyanobacterial surface bloom areas, a clear seasonal cycle was observed in Lake Taihu (Fig. 6(a)) whereby cyanobacterial blooms occurred much more often during the summer and autumn and less frequently in the spring and winter. Cyanobacterial surface bloom areas were also significantly higher in the summer and autumn than in the spring and winter (p < 0.005; t-test) (Fig. 6(a)).
We also found substantial interannual variations in the cyanobacterial surface bloom areas (Fig. 6(b)). For the entire lake, we found an obvious difference between 2006 and 2007 and other years. In 2006 and 2007, mean cyanobacterial surface bloom areas were larger than those existing during the other periods. For example, the mean area reached in 2007 was almost twice that reached in 2003. Cyanobacterial surface bloom areas were relatively small between 2003 and 2005, and the areas were larger than those of 2003 to 2005 from 2008 to 2013 but smaller than those of 2006 to 2007. In addition, the annual mean cyanobacterial surface bloom area expanded from 2003 to 2013 (R2 = 0.78, p < 0.005, excluding 2006 and 2007). The annual mean cyanobacterial surface bloom area of Lake Taihu gradually expanded from 115.91 km2 in 2003 to 167.77 km2 in 2013.
Responses to nutrient enrichment and meteorological factors
The Chla and TP time series are somewhat synchronous, showing similar variations except during the period of each year showing low Chla and TP values (Fig. 7(a)). This indicates that Chla responds differently to TP by season. The highest correlation coefficient (r = 0.91; p < 0.005) between Chla and TP occurred in the spring with a simple linear function (Chla = 89.91*TP + 10.14). In contrast, the two parameters were not significantly correlated during the other seasons. However, we found no consistent variation trend and no obvious relationship between Chla and TN (Fig. 7(a))
These results suggest that Chla is only phosphorus-limited in the spring, and they explain variations in Lake Taihu’s spring Chla values. However, the other three seasons showed a weak relationship between Chla and TP, meaning that Chla is not phosphorus-limited during these seasons.
Significantly positive correlations were found between the daily mean MODIS-Aqua derived Chla and air temperature for Lake Taihu (r = 0.43, p < 0.005) (Fig. 8). We note that Pearson’s correlation coefficients between Chla and temperature vary by sub-region: Pearson’s correlation coefficients for Zhushan Bay (r = 0.69, p < 0.005) and Meiliang Bay (r = 0.60, p < 0.005) are higher than those of the open areas (r = 0.39, p < 0.005).
No significant correlation was found between daily mean MODIS-Aqua derived Chla and wind speed. This is likely attributed to other meteorological factors, such as air temperature, that drive the majority of daily variations in cyanobacterial biomass. Previous studies have suggested that wind has spatially variable effects on Lake Taihu44,45, indicating that wind may affect Chla spatial variations independently of other meteorological effects. The spatial relationship between wind effects and Chla for this lake was not investigated, as we did not have access to the spatial distribution of wind fields across Lake Taihu.
Overall, our results suggest that air temperature plays a critical role in the cyanobacterial dynamics of Lake Taihu. More specifically, high air temperatures and low pressure levels prompt cyanobacterial bloom formation. The effects of air temperature and pressure vary depending on the sub-region of Lake Taihu concerned; these variables have a stronger impact on cyanobacterial growth in Zhushan Bay and Meiliang Bay, where wind and wave effects are relatively weak and where nutrient levels are high. However, in open areas with strong wind effects and low nutrient levels, effects of air temperature on cyanobacterial growth are less pronounced. Thus, our results suggest that it is necessary to differentiate between the various sub-regions of Lake Taihu when studying the effects of air temperature and pressure on cyanobacterial growth.
The Pearson’s correlation coefficient between the two parameters was found to be highest in Meiliang Bay (r = 0.87, p < 0.0005) and Zhushan Bay (r = 0.82, p < 0.0005) and lower in open areas (r = 0.62, p < 0.0005). Across the entire lake area, there is a significantly negative correlation between annual cyanobacterial surface bloom areas and atmospheric pressure levels (r = 0.68, p < 0.0005) (Fig. 9), suggesting that cyanobacterial blooms are more likely to occur during periods of low atmospheric pressure in Lake Taihu.
Overall, annual mean cyanobacterial surface bloom areas contracted as wind speeds increased (Fig. 10). However, no significant correlation was found between annual mean cyanobacterial surface bloom areas and wind speed. When we exclude the data for 2006 and 2007, which include anomalously large bloom areas, we find statistically significant negative correlations for Lake Taihu (from 2003 to 2005 and from 2008 to 2013) (r = −0.69, p < 0.0005) and especially for open areas (r = −0.82, p < 0.0005) (Fig. 10). These correlations indicate that low wind speeds favor cyanobacterial surface bloom occurrence. In comparison, correlation coefficients for Meiliang Bay (r = −0.69, p < 0.0005) and Zhushan Bay were found to be lower than those of open areas, indicating that cyanobacterial blooms are less sensitive to wind speed in the two bays than in the open areas.
Our results suggest that low atmospheric pressure and wind speed levels can be hypothesized to facilitate cyanobacterial surface bloom formation. Cyanobacterial blooms in Lake Taihu are generally hindered by high atmospheric pressure and strong winds. However, the responses of cyanobacterial blooms in Lake Taihu to atmospheric pressure and wind speed vary spatially. Cyanobacterial blooms are more sensitive to atmospheric pressure in Meiliang Bay and Zhushan Bay than in open areas; however, wind speed more heavily controls cyanobacterial blooms in open areas than in the two bays.
The effects of nutrient enrichment and meteorological factors on cyanobacterial dynamics
As relationships between cyanobacterial biomass and these environmental factors may be non-linear and may involve complex interactions, we utilized the CART model to further examine effects of the environmental factors on inter-annual Chla and cyanobacterial surface bloom area dynamics. Chla was found to be sensitive to TP and temperature, but cyanobacterial surface bloom areas were found to be sensitive to air pressure and wind speed. Therefore, TP and temperature levels were used as Chla CART model inputs while air pressure and wind speed values were used in the cyanobacterial surface bloom CART model.
A root node (Node-1) that included all of the data (N = 11 years, data = 100%) was used as a starting point (Fig. 11(a)). Temperature was found to be the most important factor affecting inter-annual variability in Chla values in Lake Taihu. The average Chla value for this group was 25.53 μg/L (SD = 1.63 μg/L). The observations were divided (threshold value = 17 °C) further based on high and low temperatures. Under higher temperatures, relatively higher Chla values were observed in Lake Taihu (node = 2, mean = 27.81 μg/L (SD = 1.09 μg/L), percentage = 27.3%). Under lower temperatures, TP values also heavily affected Chla levels in Lake Taihu. A significant increase in TP (≧0.136 mg/L) produced high Chla conditions (node = 5, mean = 25.59 μg/L (SD = 0.51 μg/L), percentage = 36.3%) at this node. In summary, high temperature, high TP and low TP values were found to control inter-annual variations in Chla by approximately 27.3%, 36.4% and 36.4%, respectively, in Lake Taihu from 2003 to 2013. These results indicate that temperature was the most important factor affecting inter-annual Chla variability (followed by TP) in Lake Taihu from 2003 to 2013.
A root node (Node-1) that included all of the data (N = 11 years, data = 100%) was used as a starting point (Fig. 11(b)). Air pressure was found to be the most important factor affecting the inter-annual dynamics of cyanobacterial bloom areas in Lake Taihu. The average cyanobacterial surface bloom area for this group was 154.7 km2 (SD = 27.3 km2). Observations were divided (threshold value = 1,014 hpa) further based on high and low air pressure values. Under lower air pressure conditions, relatively larger cyanobacterial surface bloom areas were observed in Lake Taihu (node = 2, mean = 175.0 km2 (SD = 26.3 km2), percentage = 45.5%). Under higher air pressure conditions, wind speeds also heavily affected the formation of cyanobacterial surface blooms in Lake Taihu. A significant decrease in wind speed (≦2.87 m/s) expanded the cyanobacterial surface bloom area (node = 4, mean = 148.9 km2 (SD = 1.9 km2), percentage = 27.3%) at this node. In summary, low air pressure, low wind speed and high wind speed values controlled inter-annual variations in cyanobacterial surface bloom area values by approximately 45.5%, 27.3% and 27.2%, respectively, in Lake Taihu from 2003 to 2013. These results indicate that air pressure was the most important factor controlling cyanobacterial surface bloom formation (followed by wind speed) in Lake Taihu from 2003 to 2013.
Discussion
Advantages of the 11-year MODIS Cyanobacteria data records
Lake Taihu is managed through an extensive water quality observational network based on a strict nutrient management strategy. However, monthly field sampling can introduce uncertainties into both short- and long-term observations of water quality (e.g., Chla) dynamics due to rapid changes that occur in this lake45. In addition, Lake Taihu is subject to complex interacting dynamics and often experiences changes in aquatic system functioning due to various physical, chemical and biotic drivers9. All of these factors make it necessary to develop a remote sensing approach for monitoring Lake Taihu’s spatial and temporal dynamics in water quality and to quantify its responses to environmental drivers.
Recently, Qi et al. (2014) developed and validated a novel algorithm based on MERIS data for estimating phycocyanin concentrations of cyanobacteria in Lake Taihu46. Using the model and 512 MERIS images, long-term records of phycocyanin data were constructed through their study, which generated valuable information for the study of cyanobacterial biomass dynamics. However, our studies offer more benefits than those of Qi et al. (2014). First, our consideration of many more calibration and validation datasets (250 data pairs) covering the entire area of Lake Taihu and all four seasons ensures the reliability of our Chla estimation model for Lake Taihu. Second, we considered many more high-quality MODIS images (1,109 images) relative to the number of MERIS images collected in Qi et al. (2014) (512 images), and the use of more long-term Chla data should generate more reliable results overall. Second, though MERIS services have be available from April of 2012, MODIS data can be used to continually to monitor Chla in the future.
Based on field observations and laboratory experiments, previous studies have suggested that the frequent cyanobacterial blooms occurring in Lake Taihu have resulted from a combination of environmental factors, including nutrient inputs and climatic conditions9,22,23,46,47,48. However, the results of these studies have largely been based on samples drawn from a limited number of sites and over short time intervals, and so their conclusions are based on data of low spatial and temporal resolutions. It is well known that Lake Taihu’s biogeochemical parameters are characterized by complex spatial and temporal dynamics10,22.
Paerl et al.49 used a series of nutrient addition bioassays to determine that cyanobacterial growth in Lake Taihu is co-limited by nitrogen and phosphorus49; their study also showed that nitrogen and phosphorus limitations change seasonally, with phosphorus limitations generally occurring in the early spring and nitrogen limitations occurring from the summer to autumn (Paerl et al.)49. Similar results on nutrient limitations in Lake Taihu have been drawn by Xu et al.29 using the same approach as that of Paerl et al.49. Clearly, we appear to present results contradictory to those of these previous studies29,50 regarding whether nitrogen limits cyanobacterial growth. The differing time and spatial scales of data used in our study and previous studies45,49 may have caused these differences. In determining whether nitrogen or phosphorus limits cyanobacterial growth, Xu et al.29 performed four in situ nutrient addition experiments in May, July, October, and December of 200829. In the study (Xu et al.)29, each nutrient addition experiment lasted several days (day scale), and three sites in Lake Taihu were used for the four experiments (specific site scale)29. In our study, both Chla and nutrient data were collected for 2003 to 2013 for the entire lake area. Our results are thus based on an examination of the entire lake area on an annual scale, thus revealing the longer-term responses of cyanobacterial dynamics to nutrients at a broader spatial scale.
The following points could explain the seasonality of our findings. Our long-term in situ measurements showed that the ratio values of TN:TP for all seasons from 2003 to 2013 are significantly higher than the Redfield ratio (16:1)51, which has been widely used to identify nitrogen limitations (when TP:TN < 16:1) for Lake Taihu. Second, the spring season is typically considered to be a season of cyanobacterial growth during which cyanobacteria multiply rapidly and require higher levels of phosphorus. In Lake Taihu, TP levels were found to be lower in the spring than during other seasons, when they were high enough not to limit cyanobacterial growth. In addition, ratios of TN:TP in the spring with an average value of 45:1 significantly exceed the Redfield ratio (16:1) (closed to 3 times the Redfield ratio (16:1)), revealing the reasonability of phosphorus limitations in the spring in Lake Taihu. Based on laboratory experiments and field observations, previous studies have suggested that temperature increases can prompt significant cyanobacterial growth4, can cause initial bloom times to occur earlier on23, and can prolong annual bloom periods52. These data, however, due to their inherent limitations, are discrete and sporadic in nature, thus preventing us from conducting a thorough and objective evaluation of climatic driven factors of cyanobacterial dynamics in Lake Taihu. By taking advantage of the MODIS cyanobacterial time series for 2003–2013 generated using a new developed model together with historical records of meteorological data, we offer a comprehensive understanding of the climatic forces driving such cyanobacterial changes. We found that temperature levels can prompt cyanobacterial growth and that cyanobacterial surface bloom patterns are sensitive to wind speeds. The sensitivity of these climatic factors to cyanobacterial dynamics varies regionally. We also made a hypothesis that air pressures may prompt cyanobacterial surface bloom formation in Lake Taihu. Atmospheric pressure has not traditionally been considered a significant factor shaping cyanobacterial surface bloom formation. To date, there are no previous studies that have linked atmospheric pressure to cyanobacterial surface bloom formation. Our results reveal a clear increase in cyanobacterial surface bloom area with a decrease in atmospheric pressure, suggesting that atmospheric pressure can be hypothesized to drive cyanobacterial bloom formation in Lake Taihu. The finding that atmospheric pressure plays a role in cyanobacterial bloom formation contributes a new perspective on the phenomenon of cyanobacterial blooms. Forecasting models for cyanobacterial blooms could benefit from including atmospheric pressure as a predicting parameter. Bubbles resulting from low atmospheric pressure should encourage surface cyanobacterial bloom formation53. Bubbles in water can result from the oversaturation of methane, carbon dioxide, nitrogen, or oxygen. Thus, as atmospheric pressure decreases, gas solubility in water is reduced, in turn leading to gas oversaturation and gas bubble growth.
Implications for future lake management
The effects of meteorological factors on cyanobacterial blooms necessitate a different interpretation of the phenomenon of cyanobacterial blooms. Forecasting models for cyanobacterial blooms could benefit from considering meteorological factors such as air temperature, atmospheric pressure, and wind as prediction parameters. It is noteworthy that from the 1980 s to the present, the ratio of air temperature to atmospheric pressure has significantly increased while wind patterns have declined in the Lake Taihu region (Fig. 12). However, these trends were not apparent from 1956 to 1980 (Fig. 12). Our results suggest that these trends in the ratio of air temperature to atmospheric pressure and wind speed could lead to the generation of more severe cyanobacterial blooms in Lake Taihu. Thus, severe cyanobacterial bloom events in Lake Taihu since the 1980 s could be attributed at least in part to changes in meteorological conditions. As climate change continues to affect the region, we can infer that the severity of cyanobacterial blooms will be enhanced in Lake Taihu.
The long-term Chla data with high spatial and temporal resolutions used in this study are valuable in identifying the synergistic effects of nutrient loading on cyanobacterial growth. Our finding that Lake Taihu is only subject to phosphorus limitations in the spring has implications for controlling cyanobacterial blooms via the reduction of nutrients. This result highlights that severe cyanobacterial blooms may form during the spring seasons of years presenting high TP levels. It can also be inferred that phosphorus reduction during the spring may help control cyanobacterial blooms over short timescales; in contrast, the lack of Chla sensitivity to TN found suggests that efforts to control the magnitude and duration of cyanobacterial blooms should not focus on nitrogen reduction in Lake Taihu. Furthermore, our results suggest that nutrient reduction strategies for controlling cyanobacterial blooms must take into account co-occurring climatic changes that favor cyanobacterial bloom development. More specifically, critical nutrient thresholds above which cyanobacterial blooms are prompted under favorable meteorological conditions (increases in temperature and decreases in atmospheric pressure and wind speed) should be reduced to compensate for more favorable growth conditions50. Moreover, the adjustment of nutrient reduction strategies should be regionally specific, as the sensitivity of cyanobacterial dynamics to these meteorological conditions varies across different sub-regions of Lake Taihu.
The results of this study have significant environmental implications according to the long-term monitoring of cyanobacterial dynamics in Lake Taihu. This is particularly critical in a changing climate, as it is often difficult to elucidate the causal factors of aquatic environmental changes without continuous and long-term water quality monitoring. Our results will help managers and decision-makers account for and modify their strategies for controlling cyanobacterial blooms in response to future climate change and human impacts. We thus recommend that remote sensing approaches be incorporated into future Lake Taihu management systems.
On a broader scale, the approaches and findings of this study may be extended to other lakes in which cyanobacteria dominate such as Lake Chaohu54, Lake Dianchi50, and Lake Erie. Once the Chla estimation model is tuned with local data and once satellite-based Chla are validated, similar long-term Chla series analyses can be conducted with little effort and cost. Climatic warming and anthropogenically enhanced nutrient loading in lakes around the world are also known to be potential drivers of cyanobacterial surface bloom intensity levels50. The integration of satellite-derived data products with other in situ environmental and meteorological data will further aid in elucidating the cause and effect relationships between cyanobacterial dynamics and environmental and meteorological factors.
Conclusions
We used remote sensing techniques to address this issue in Lake Taihu, which is a large, shallow and eutrophic lake. We first developed and validated an empirical model for estimating Chla levels in Lake Taihu based on a spectral index derived from Rrc and in situ data. The proposed model exhibited robust performance for an independent validating dataset, with relative errors (RE) ranging from 0.4% to 64.5% with a mean absolute percent error (MAPE) value of 27.1% (RMSE = 15.01 μg/L) and from 6.56 μg/L to 113.66 μg/L for Chla. Second, long-term Chla and cyanobacterial surface bloom area time series were generated from MODIS-Aqua observations gathered from 2003 to 2013 by means of the proposed model and the floating algae index (FAI). Chla values for Lake Taihu varied considerably in space and time, with higher Chla values found for Zhushan Bay and Meiliang Bay and with lower Chla values found for open areas. Chla values in the summer and autumn were found to be significantly higher than those recorded in the spring and winter (p < 0.005). The annual mean area of cyanobacterial surface bloom events exhibited a clearly increasing trend for the entire lake area from 2003 to 2013 with the exception of 2006 and 2007 (annual mean areas = 3.8791*year-7646.1, p < 0.0005). From 2003 to 2013, the area gradually expanded from 115.91 km2 to 167.77 km2. Third, based on the long-term Chla and cyanobacterial bloom area derived from MODIS-Aqua observations, we quantified the responses of cyanobacterial dynamics to nutrient enrichment and climatic conditions. The results show that high air temperatures and high phosphorus levels in the spring may prompt cyanobacterial growth and that low wind speeds and atmospheric pressure levels may favor cyanobacterial surface bloom formation. The sensitivity of cyanobacterial dynamics to climatic conditions was found to vary regionally, with more sensitivity to air temperature and atmospheric pressure observed in Zhushan Bay and Meiliang Bay and with more sensitivity to wind speed found in open areas.
Additional Information
How to cite this article: Shi, K. et al. Long-term MODIS observations of cyanobacterial dynamics in Lake Taihu: Responses to nutrient enrichment and meteorological factors. Sci. Rep. 7, 40326; doi: 10.1038/srep40326 (2017).
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
References
Matthews, M. W. Eutrophication and cyanobacterial blooms in South African inland waters: 10 years of MERIS observations. Remote Sensing of Environment 155, 161–177 (2014).
M, M.-H. et al. Recent changes in phytoplankton communities associated with rapid regional climate change along the Western Antarctic Peninsula. Science 323, 1470–1473 (2009).
Duan, H., Ma, R. & Hu, C. Evaluation of remote sensing algorithms for cyanobacterial pigment retrievals during spring bloom formation in several lakes of East China. Remote Sensing of Environment 126, 126–135 (2012).
Liu, X., Lu, X. & Chen, Y. The effects of temperature and nutrient ratios on Microcystis blooms in Lake Taihu, China: An 11-year investigation. Harmful Algae 10, 337–343 (2011).
Michalak, A. M. et al. Record-setting algal bloom in Lake Erie caused by agricultural and meteorological trends consistent with expected future conditions. Proceedings of the National Academy of Sciences of the United States of America 110, 6448–6452 (2013).
Paerl, H. W., Gardner, W. S., Mccarthy, M. J., Peierls, B. L. & Wilhelm, S. W. Algal blooms: noteworthy nitrogen. Science 346, 175–175 (2014).
Hu et al. Moderate Resolution Imaging Spectroradiometer (MODIS) observations of cyanobacteria blooms in Taihu Lake, China. Jounral of Geophysical Research (2010).
Paerl, H. W. & Huisman, J. Climate change: a catalyst for global expansion of harmful cyanobacterial blooms. Environmental Microbiology Reports 1, 27–37 (2009).
Qin, B. et al. Cyanobacterial bloom management through integrated monitoring and forecasting in large shallow eutrophic Lake Taihu (China). Journal of Hazardous Materials 287C, 356–363 (2015).
Shi, K., Li, Y., Li, L. & Lu, H. Absorption characteristics of optically complex inland waters: Implications for water optical classification. Journal of Geophysical Research Biogeosciences 118, 860–874 (2013).
Shi, K. et al. Remote estimation of cyanobacteria-dominance ininland waters. Water Research 68, 217–226 (2014).
Schindler, D. W., Hecky, R. E. & Mccullough, G. K. The rapid eutrophication of Lake Winnipeg: Greening under global change. Journal of Great Lakes Research 38, 6–13 (2012).
Jöhnk, K. D. J. et al. Summer heatwaves promote blooms of harmful cyanobacteria. Global Change Biology 14, 495–512 (2008).
Kutser, T., Metsamaa, L., Strömbeck, N. & Vahtmäe, E. Monitoring cyanobacterial blooms by satellite remote sensing. Estuarine Coastal & Shelf Science 67, 303–312 (2006).
Paerl, H. W. et al. Controlling cyanobacterial blooms in hypertrophic Lake Taihu, China: will nitrogen reductions cause replacement of non-N2 fixing by N2 fixing taxa? Plos One 9, e113123–e113123 (2014).
Liu, X., Zhang, Y., Yin, Y., Wang, M. & Qin, B. Wind and submerged aquatic vegetation influence bio-optical properties in large shallow Lake Taihu, China. Journal of Geophysical Research Biogeosciences 118, 713–727 (2013).
Paerl, H. W. & Paul, V. J. Climate change: Links to global expansion of harmful cyanobacteria. Water Research 46, 1349–1363 (2012).
García, R. & Johnstone, R. W. Effects of Lyngbya majuscula (Cyanophycea) blooms on sediment nutrients and meiofaunal assemblages in seagrass beds in Moreton Bay, Australia. Marine & Freshwater Research 57, 155–165 (2006).
Otten, T. G., Xu, H., Qin, B., Zhu, G. & Paerl, H. W. Spatiotemporal patterns and ecophysiology of toxigenic microcystis blooms in Lake Taihu, China: implications for water quality management. Environmental Science & Technology 46, 3480–3488 (2012).
Shi, K. et al. Long-Term Satellite Observations of Microcystin Concentrations in Lake Taihu during Cyanobacterial Bloom Periods. Environmental Science & Technology 49, 6448–6456 (2015).
Otten, T. G. & Paerl, H. W. Phylogenetic Inference of Colony Isolates Comprising Seasonal Microcystis Blooms in Lake Taihu, China. Microbial Ecology 62, 907–918 (2011).
Qin, B., Xu, P., Wu, Q., Luo, L. & Zhang, Y. Environmental issues of Lake Taihu, China. Hydrobiologia 581, 3–14 (2007).
Deng, J. et al. Earlier and warmer springs increase cyanobacterial (Microcystis spp.) blooms in subtropical Lake Taihu, China. Freshwater Biology 59, 1076–1085 (2014).
Elliott, J. A. The seasonal sensitivity of Cyanobacteria and other phytoplankton to changes in flushing rate and water temperature. Global Change Biology 16, 864–876 (2009).
Komatsu, E., Fukushima, T. & Harasawa, H. A modeling approach to forecast the effect of long-term climate change on lake water quality. Ecological Modeling 209, 351–366 (2007).
Matthews, M. W. & Odermatt, D. Improved algorithm for routine monitoring of cyanobacteria and eutrophication in inland and near-coastal waters. Remote Sensing of Environment 156, 374–382 (2015).
O’Neil, J. M., Davis, T. W., Burford, M. A. & Gobler, C. J. The rise of harmful cyanobacteria blooms: The potential roles of eutrophication and climate change. Harmful Algae 14, 313–334 (2013).
Paerl, H. W. & Huisman, J. Blooms like it hot. Science 320, 57–58 (2008).
Xu, H. et al. Controlling harmful cyanobacterial blooms in a hyper-eutrophic lake (Lake Taihu, China): The need for a dual nutrient (N & P) management strategy. Water Research 45, 1973–1983 (2010).
Jutla, A. S., Akanda, A. S. & Islam, S. Satellite remote sensing of space–time plankton variability in the Bay of Bengal: Connections to cholera outbreaks. Remote Sensing of Environment 123, 196–206 (2012).
Li, Y. & He, R. Spatial and temporal variability of SST and ocean color in the Gulf of Maine based on cloud-free SST and chlorophyll reconstructions in 2003–2012. Remote Sensing of Environment 144, 98–108 (2014).
Keith, D. J. Satellite remote sensing of chlorophyll a in support of nutrient management in the Neuse and Tar–Pamlico River (North Carolina) estuaries. Remote Sensing of Environment 153, 61–78 (2014).
Labiosa, R. G., Arrigo, K. R., Genin, A., Monismith, S. G. & Dijken, G. V. The interplay between upwelling and deep convective mixing in determining the seasonal phytoplankton dynamics in the Gulf of Aqaba: Evidence from SeaWiFS and MODIS. Limnology & Oceanography 48, 2355–2368 (2003).
Le, C., Hu, C., English, D., Cannizzaro, J. & Kovach, C. Climate-driven chlorophyll-a changes in a turbid estuary: Observations from satellites and implications for management. Remote Sensing of Environment 130, 11–24 (2013).
Volpe, G., Nardelli, B. B., Cipollini, P., Santoleri, R. & Robinson, I. S. Seasonal to interannual phytoplankton response to physical processes in the Mediterranean Sea from satellite observations. Remote Sensing of Environment 117, 223–235 (2012).
Gons, H. J., Auer, M. T. & Effler, S. W. MERIS satellite chlorophyll mapping of oligotrophic and eutrophic waters in the Laurentian Great Lakes. Annales Chirurgiae Et Gynaecologiae 90 Suppl 215, 29–32 (2008).
Mishra, S. & Mishra, D. R. Normalized difference chlorophyll index: A novel model for remote estimation of chlorophyll-a concentration in turbid productive waters. Remote Sensing of Environment 117, 394–406 (2012).
Dall’Olmo, G. et al. Assessing the potential of SeaWiFS and MODIS for estimating chlorophyll concentration in turbid productive waters using red and near-infrared bands. Remote Sensing of Environment 96, 176–187 (2005).
Gons, H. J. Optical Teledetection of Chlorophyll a in turbid inland waters. Environmental Science & Technology 33, 1127–1132 (1999).
Le, C. et al. A four-band semi-analytical model for estimating chlorophyll a in highly turbid lakes: The case of Taihu Lake, China. Remote Sensing of Environment 113, 1175–1182 (2009).
Jespersen, A. & Christoffersen, K. Measurements of Chl a from phytoplankton using ethanol as extraction solvent. Archiv Fur Hydrobiologie 109 (1987).
Ebina, J., Tsutsui, T. & Shirai, T. Simultaneous determination of total nitrogen and total phosphorus in water using peroxodisulfate oxidation. Water Research 17, 1721–1726 (1983).
Joshi, I. & D’Sa, E. Seasonal variation of colored dissolved organic matter in Barataria Bay, Louisiana, using combined landsat and field data. Remote Sensing 7, 12478–12502 (2015).
Wu, T. et al. Dynamics of cyanobacterial bloom formation during short-term hydrodynamic fluctuation in a large shallow, eutrophic, and wind-exposed Lake Taihu, China. Environmental Science & Pollution Research International 20, 8546–8556 (2013).
Zhang, Y., Shi, K., Liu, X., Zhou, Y. & Qin, B. Lake topography and wind waves determining seasonal-spatial dynamics of total suspended matter in turbid Lake Taihu, China: Assessment using long-term high-Resolution MERIS data. Plos One 9, e98055 (2014).
Ma, J. et al. Environmental factors controlling colony formation in blooms of the cyanobacteria Microcystis spp. in Lake Taihu, China. Harmful Algae 31, 136–142 (2014).
Chen, Y., Fan, C., Teubner, K. & Dokulil, M. Changes of nutrients and phytoplankton chlorophyll- a in a large shallow lake, Taihu, China: an 8-year investigation. Hydrobiologia 506–509, 273–279 (2003).
Chen, Y., Qin, B., Teubner, K. & Dokulil, M. T. Long-term dynamics of phytoplankton assemblages: Microcystis-domination in Lake Taihu, a large shallow lake in China. Journal of Plankton Research 25, 445–453 (449) (2003).
Paerl, H. W. et al. Nutrient limitation dynamics examined on a multi-annual scale in Lake Taihu, China: implications for controlling eutrophication and harmful algal blooms. Journal of Freshwater Ecology 30, 5–24 (2015).
Sun, D., Hu, C., Qiu, Z. & Wang, S. Reconstruction of hyperspectral reflectance for optically complex turbid inland lakes: test of a new scheme and implications for inversion algorithms. Optics Express 23, 718–740 (2015).
Redfield, A. C. The biological control of chemical factors in the environment. Science Progress 11, 150–170 (1960).
Zhang, M., Duan, H. T., Shi, X. L., Yu, Y. & Kong, F. X. Contributions of meteorology to the phenology of cyanobacterial blooms: Implications for future climate change. Water Research 46, 442–452 (2012).
Oliver, R. L. & Walsby, A. E. Direct evidence for the role of light-mediated gas vesicle collapse in the buoyancy regulation of Anabaena flos-aquae (Cyanobacteria). Limnology & Oceanography 29, 879–886 (1984).
Sun, D., Qiu, Z., Li, Y., Shi, K. & Gong, S. Detection of total phosphorus concentrations of turbid inland waters using a remote sensing method. Water Air & Soil Pollution 225, 1122–1132 (2014).
Acknowledgements
This study was supported by grants from the National Natural Science Foundation of China (Nos 41301376, 41325001, 41621002, and 41301022), International Science & Technolog Cooperation Program of China (2014DFG91780), The Provincial Natural Science Foundation of Jiangsu in China (No. BK20141515), Key Research Program of Frontier Sciences, Chinese Academy of Sciences (QYZDB-SSW-DQC016), and Youth Innovation Promotion Association (CAS). The authors would like to thank Hans W. Paerl (University of North Carolina at Chapel Hill) for his suggestions pertaining to the study. The authors thank all members of the Taihu Laboratory for Lake Ecosystem Research (TLLER) for their participation in the field experiments.
Author information
Authors and Affiliations
Contributions
K.S. wrote the main manuscript text. Y.L.Z., X.H.L., K.S. and B.Q.Q. designed and conducted the experiments and collected and analyzed the data. G.G. and G.W.Z. contributed to analyzing the data and writing and editing the manuscript. Y.Q.Z. performed statistical analyses and helped interpret the data. All of the authors reviewed the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Shi, K., Zhang, Y., Zhou, Y. et al. Long-term MODIS observations of cyanobacterial dynamics in Lake Taihu: Responses to nutrient enrichment and meteorological factors. Sci Rep 7, 40326 (2017). https://doi.org/10.1038/srep40326
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/srep40326