Subscriber access provided by UNIV OF NEW ENGLAND ARMIDALE
Characterization of Natural and Affected Environments
Phenology of Phytoplankton Blooms in a Trophic Lake Observed from Long-Term MODIS Data Kun Shi, Yunlin Zhang, Yibo Zhang, Na Li, Boqiang Qin, Guangwei Zhu, and Yongqiang Zhou Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b06887 • Publication Date (Web): 18 Feb 2019 Downloaded from http://pubs.acs.org on February 18, 2019
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 19
Environmental Science & Technology
1
Phenology of Phytoplankton Blooms in a Trophic Lake Observed from Long-Term
2
MODIS Data
3
Kun Shi 1, 2*, 3, Yunlin Zhang1, 2, Yibo Zhang1, 2, Na Li1, 2, Boqiang Qin1, 2, Guangwei Zhu1, 2,
4
Yongqiang Zhou1, 2
5 6
1
7
Environment, Nanjing Institute of Geography and Limnology, Chinese Academy of Sciences,
8
Nanjing 210008, China
Taihu Laboratory for Lake Ecosystem Research, State Key Laboratory of Lake Science and
9 10
2
3
University of Chinese Academy of Sciences, Beijing 100049, China
CAS Center for Excellence in Tibetan Plateau Earth Sciences, Beijing 100101, China
11 12 13 14 15
* Corresponding author. Email:
[email protected] 16
Phone: (+86) -25-86882174, Fax: (+86) 25-57714759
17
University of Chinese Academy of Sciences, Beijing 100049, China
18
1
ACS Paragon Plus Environment
Environmental Science & Technology
19
ABSTRACT
20
Phytoplankton phenology critically affects elements biogeochemical cycles, ecosystem
21
structure, and productivity. However, our understanding about the phenological process and
22
driving mechanism is still very limited due to the shortage of long-term observation data. We
23
used all available daily MODIS-Aqua data from 2003 to 2017 to determine bloom start dates
24
(BSDs) in a typical trophic lake (Lake Taihu) and investigate how phytoplankton BSDs respond
25
to climate change and environmental factors. The results indicate that BSDs have advanced
26
29.9 days for the entire Lake Taihu from 2003 to 2017. Spatially, an earlier phytoplankton
27
bloom was recorded in the northern bays and the littoral regions than in the center open water.
28
Air temperature, wind speed, and N/P ratio (N: total nitrogen; P: total phosphorus) were three
29
important factors affecting phytoplankton phenology. Multiple linear correlation showed that
30
air temperature, wind speed, and N/P ratio in spring could explain 59.9% variability of BSDs
31
for Lake Taihu. This study provides a quantitative assessment of phytoplankton phenological
32
shifts and elucidates the interrelationship between phenology parameters and environmental
33
factors, thus improving our understanding on the potential impact of climate change and
34
eutrophication on lake ecosystems. The starting earlier and lasting longer of phytoplankton are
35
consistent with the expected effects of climate warming on aquatic ecosystem in recent decades,
36
which will bring new challenges for algal bloom management in eutrophic Lake Taihu.
37
Key words: Phytoplankton phenology, Remote sensing, Climate warming, MODIS, Lake
38
Taihu
39
2
ACS Paragon Plus Environment
Page 2 of 19
Page 3 of 19
Environmental Science & Technology
40
TOC
41
3
ACS Paragon Plus Environment
Environmental Science & Technology
42
INTRODUCTION
43
Global change has been proven to influence the temporal-spatial variation and magnitude
44
of several climatic conditions, such as air temperature, rainfall, and wind. Meanwhile, it exerts
45
a series of extensive and profound effects on the global ecosystem 1 - 3. Phenological processes
46
and changes are the critical response mechanisms of ecosystem to climate change. Therefore,
47
extensive studies have been conducted on the phenology in terrestrial and aquatic systems 4 - 6.
48
However, changes in the temperature, precipitation, and wind speed regimes are neither
49
spatially nor temporally uniform, implying spatially heterogeneous responses of ecosystems to
50
climatic change 3. However, a comprehensive and profound understanding of the phenological
51
process at individual, species, population, community and ecosystem levels is still very limited,
52
and more data, evidences and studies are needed to support ecosystem management.
53
Traditionally, a time series for monitoring the phenology of terrestrial and aquatic
54
ecosystems is built by regularly recording the timing and collect in situ samples in certain sites
55
5 - 8.
56
over intervals of a week to two weeks during the stratified period and monthly during the
57
unstratified period since 1962, which was used to address the impacts of climate change on the
58
phenology of lake processes 6. Remotely sensed data provide strong evidence and promote
59
significant progresses for monitoring phenological processes of terrestrial and aquatic systems
60
from local, regional to global scales since the mid-1980s because of their synoptic coverage,
61
repeated temporal sampling and low cost that satellite observations afford 9 - 11. Several studies
62
have demonstrated different remote sensing water color products deriving from SeaWiFS,
63
MODIS, and MERIS images can be used to examine the timing, magnitude, and duration of
64
phytoplankton blooms 12 - 14.
For example, the limnological data have been obtained at a central site of Lake Washington
65
Aquatic ecosystems and especially phytoplankton are highly sensitive to environmental
66
and climate changes 15 - 17. Recently, several studies have demonstrated that eutrophication and
67
climate warming are anticipated to prompt harmful algal blooms in various aquatic ecosystems
68
over the world 18, 19. Accurately quantifying phytoplankton dynamics (phenology) is essential
69
for improving our understanding in the food web, fishery resources and carbon cycle with
70
different trophic levels in aquatic ecosystems 5. Previous studies indicated that the timing of
71
phytoplankton blooms advanced much faster than that of plants on land under when facing 4
ACS Paragon Plus Environment
Page 4 of 19
Page 5 of 19
Environmental Science & Technology
72
climate warming 20, 21. For example, the annual phytoplankton bloom maximum has advanced
73
by up to 50 days from 1997 to 2009 as a consequence of changes in seasonal ice cover which
74
may have significant consequences for the Arctic food web structure and carbon cycling 9.
75
As one of the most biologically productive regions in the world, Lake Taihu is
76
characterized by frequent algal blooms 22 - 24. Climate change is likely to prompt the occurrence
77
of harmful algal blooms and ecological disasters 25 - 27. For example, a combination of excessive
78
nutrient loading and extreme climatic conditions induced a severe cyanobacterial bloom of
79
Lake Taihu in 2007, resulting that approximately 2, 000, 000 people did not have drinking water
80
for more than seven days 24. Characterization of the phytoplankton phenology in Lake Taihu
81
has important implications for the understanding the roles of climate change and eutrophication
82
in controlling the phytoplankton phenology and instigating ecological disasters.
83
Here we combined a long-term series of remotely sensed chlorophyll a (Chla) with climate,
84
hydrology and water chemistry data in Lake Taihu to address how BSDs have shifted from
85
2003 to 2017 and determine which climatic, hydrological or water chemistry factors mainly
86
influence the initial time of phytoplankton blooms. The results will have important implications
87
in the monitoring, early warning and management of algal bloom to mitigate the disaster effects
88
in eutrophic Lake Taihu.
89
MATERIALS AND METHODS
90
Study Area
91
As a large and shallow water body located in the flood plain in Yangtze River Delta, the
92
water quality of Lake Taihu has been seriously affected by eutrophication and phytoplankton
93
blooms
94
industry, agriculture, tourism activities. Taihu Lake is a unique and complex lake ecosystem
95
that provides an indispensable habitat for a variety of species, as well as diverse natural and
96
economic resources to surrounding cities 23. Over the last four decades, however, Lake Taihu
97
has a changed from a dominant position in macrophyte to a phytoplankton state, and is affected
98
by anthropogenic pressure 28. Excessive nitrogen and phosphorus nutrients promote the growth
99
of phytoplankton, leading to the proliferation of surface cyanobacterial blooms 26.
100 101
23, 24.
It receives a large amount of pollutants from various sources, mainly including
MODIS-Derived Daily Chlorophyll a Concentration (Chla) Data MODIS-Aqua data have a maximum spatial resolution of 250 meters (bands 1 and 2) and 5
ACS Paragon Plus Environment
Environmental Science & Technology
Page 6 of 19
102
a very short revisit interval (1 image/day), which have been available free of charge since 2002.
103
We downloaded more than 5,000 MODIS-Aqua L-0 images acquired from January 2003 to
104
December
105
http://oceancolor.gsfc.nasa.gov/). MODIS-Aqua images often suffer the impacts from clouds,
106
cloud shadows, or thick aerosols, meaning that not all the acquired images could be useful for
107
this study. After visual examination of these images, we selected 1,401 high-quality images to
108
process to Level-1 (calibrated spectral radiance) data with the SeaDAS software package
109
(version 7.4). The atmospherically Rayleigh-corrected MODIS-Aqua data (Rrc) were then
110
derived following the methods detailed by previous study 26.
111
2017
from
NASA's
Goddard
Space
Flight
Center
website
(GSFC,
Using 250 matched data pairs of in situ Chla – MODIS Rrc, Shi et al. (2017) developed
112
and validated a statistically emperical model for deriving Chla for Lake Taihu
113
1545.3*[(EXP(Rrc(645)) − EXP(Rrc(859)))/(EXP(Rrc(645)) + EXP(Rrc(859)))]+69.346. For the
114
validation dataset, the relative errors of the model were from 0.4% to 64.5% with a mean-
115
absolute-percent-error of 27.1% (root-mean-square error = 15.01 μg/L), suggesting that the
116
model performed well in Lake Taihu. Here, we applied this model to the selected MODIS-Aqua
117
images (atmospherically Rayleigh-corrected MODIS-Aqua products) to derive daily Chla data
118
for Lake Taihu from January 2003 to December 2017 26.
119
Derivation of Phenological Metrics
26:
Chla=-
120
Phenological metrics are defined to represent the critical characteristics of the
121
phytoplankton annual cycle. They were selected to focus on the phytoplankton biomass and on
122
the time of the main growing period. The extraction of BSDs generally has two steps:
123
smoothing daily Chla time series and defining BSDs. The process of data smoothing is critical
124
for accurately deriving BSDs from MODIS-derived Chla data because the smoothing can
125
remove noise and abnormal data. Smoothing of the MODIS-derived daily Chla data was
126
performed using the TIMESAT software, which was developed originally for retrieving and
127
mapping phenology features by use of terrestrial vegetation indices’ maps. This software can
128
also be applied to Chla time series of lakes which exhibit the regularly cycling characteristics
129
as terrestrial vegetation indices do 12, 29. The TIMESAT software integrated three methods for
130
smoothing: Double Logistic, asymmetric Gaussian and Savitzky – Golay filtering 29. Palmer et
131
al. (2015) 12 suggested that asymmetric Gaussian filtering may be better for inland waters and 6
ACS Paragon Plus Environment
Page 7 of 19
132
Environmental Science & Technology
thus we selected this method for processing our MODIS-derived Chla data in this study.
133
The BSDs are usually defined in two ways. One way is to use the date at which the input
134
mapped value (Chla) rises above a pre-defined percentage of the detected peak; the other uses
135
the date at which the mapped value rises above a pre-defined Chla value 12. The two approaches
136
to determining BSDs, agree upon each other in determining the large-scale pattern of
137
phytoplankton phenology 12. Previous studies for assessing BSDs in pelagic oceanic and inland
138
define the BSDs as the dates when Chla rising above background median concentrations + 5%
139
12, 30, 31.
140
Meiliang Bay, Zhushan Bay, Open area and entire Lake Taihu for the full 2003 - 2017 time
141
series + 5% were found to range between 20.03 and 33.43 μg/L, between 22.05 and 37.11 μg/L,
142
between 13.42 and 23.96 μg/L, and between 14.53 and 25.12 μg/L, respectively. It should be
143
noted that the median values varied significantly across different regions of Lake Taihu. This
144
indicates that the TIMESAT software was not appropriate for deriving BSDs in Lake Taihu as
145
the values defining bloom start timing cannot be assigned on a per-pixel basis. We calculated
146
yearly median + 5% values for a pixel position from 2003 to 2017 for Lake Taihu and assigned
147
these median + 5% values as threshold of BSDs. Subsequently, we could derive BSDs values
148
and map the spatial distribution of BSDs.
149
Chemical Data
Here, we adapted the similar approach to defining BSDs. Median Chla values of
150
In situ pH, dissolved oxygen, Secchi disc depth (SDD) etc. were measured during the
151
sampling cruise and other chemical parameters including total nitrogen (TN), total phosphorus
152
(TP), Chla were measured by collecting samples for analysis in Taihu Laboratory for Lake
153
Ecosystem Research (TLLER).These monthly or seasonal in situ data were collected from 2003
154
to 2017 via many ship surveys at pre-defined sampling stations distributed in Meiliang Bay,
155
Zhushan Bay and the open water. We used a standard 30-cm diameter Secchi disk to measure
156
SDD. TN and TP concentrations were determined using alkaline potassium persulfate digestion,
157
followed by absorbance analysis, using a Shimadzu UV–2550PC spectrophotometer.
158
Meteorological Data
159
To investigate the ecological response of phytoplankton phenology to meteorological
160
conditions, daily meteorological data including air temperature, wind speed; sunshine duration
161
and precipitation were downloaded from the China Meteorological Data Sharing Service 7
ACS Paragon Plus Environment
Environmental Science & Technology
162
System (http://cdc.nmic.cn). These meteorological data were measured at Dongshan station
163
(31.06o N, 120.43o E). The water temperature at 0.5 m depth in the open water along TLLER
164
was measured three times/day at 8:00, 14:00 and 20:00 since 1992. In order to exhibit the long-
165
term trend and maintain the consistency of observation of air and water temperature, the data
166
from 2003 to 2017 were used in this study.
167
Statistical Approach
168
We used the Statistical Program for Social Sciences (SPSS) 20.0 software to investigate
169
the relations between the variables. Multiple linear regression was performed to assess the joint
170
effects of key meteorological, chemical factors on BSDs. As the temporal trends of
171
environmental factors and BSDs were hypothesized to be non-linear and the generalized
172
additive model (GAM) provided a flexible and effective technique for modeling a nonlinear
173
time-series, the trend for air temperature, water temperature, sunshine duration, precipitation,
174
wind speed, SDD, TN, TP, N/P ratio (N: total nitrogen; P: total phosphorus), pH, and BSDs
175
was fitted using the GAM procedures (“mgcv” function) in R software (Version 3.0.0.).
176
RESULTS
177
Temporal and Spatial Patterns of BSDs
178
The Chla concentration derived from the daily MODIS images between 2003 and 2017
179
demonstrated apparent seasonal cycles, with the yearly peaks observed in July–August (Fig. 1).
180
In addition, a marked Chla concentration peak is recorded in 2017 but not for Meiliang Bay.
181
Although this peak was not present for Meiliang Bay, the average Chla concentration of this
182
year was higher than that of other years. The spatial distribution of phytoplankton BSDs (day
183
of the year) from 2003 to 2017 in Lake Taihu was shown in Fig. 2. There are significant
184
temporal and spatial differences in phytoplankton BSDs (Fig. 1). Overall, the northern bays
185
such as Zhushan Bay and Meiliang Bay and the coastal regions showed earlier phytoplankton
186
bloom than the center open water (Figs. 2 and 3). The mean BSDs from 2003 to 2017 are 106.6,
187
111.5 and 121.7 in Meiliang Bay, Zhushan Bay, and Open area, respectively. Generally, the
188
mean start day of the dominant phytoplankton blooms varied with latitude. Obviously, latitude
189
is not the controlling factor of the fine spatial differences of phytoplankton BSDs in Lake Taihu.
190
Higher nutrients and water temperature in northern bays and coastal regions may be the
191
potential controlling factors. 8
ACS Paragon Plus Environment
Page 8 of 19
Page 9 of 19
Environmental Science & Technology
192
193
194
195 196
Figure 1 Long time series of Chla derived from MODIS-Aqua image data in Lake Taihu from
197
2003 to 2017 and the corresponding BSDs
9
ACS Paragon Plus Environment
Environmental Science & Technology
198
199
200
201
202
10
ACS Paragon Plus Environment
Page 10 of 19
Page 11 of 19
Environmental Science & Technology
203
Figure 2 Maps of phytoplankton BSDs based on MODIS 250 m spatial resolution daily Chla
204
data in Lake Taihu from 2003 to 2017
205 206
Long Term Trend of BSDs
207
Fig. 3 shows the long term trend of BSDs and GAM fitting line for Meiliang Bay, Zhushan
208
Bay, Open area and entire Lake Taihu. Both linear and GAM fitting show that BSDs have
209
significantly advanced for all the lake regions (p ≤ 0.05) especially for the past 15 years. BSDs
210
have advanced 48.8, 40.0, 29.4, and 32.0 days in Meiliang Bay, Zhushan Bay, Open area and
211
entire Lake Taihu based on the linear fitting. Therefore, BSDs advance more in northern lake
212
bays with high nutrients than in Open area.
213 214 215
Fig. 3 Long term trend of BSDs in Lake Taihu from 2003 to 2017 Affecting Factors
216
Our efforts detect the relationships between abiotic factors and phytoplankton bloom
217
timing and size to elucidate the potential affecting factors. The impacts of climatic and water
218
chemical factors on phytoplankton phenology were examined using the linear fitting approach.
219
BSDs were found to correlate with air temperature, water temperature, wind speed, sunshine
220
duration, precipitation, Secchi disc depth, total nitrogen, total phosphorus, N/P ratio and pH, 11
ACS Paragon Plus Environment
Environmental Science & Technology
Page 12 of 19
221
but exhibited different spatial and seasonal patterns. Overall, high and significantly positive
222
correlations are found between BSDs and spring mean wind speed, and between BSDs and
223
spring mean N/P ratio, but negative correlations are expectedly found between BSDs and spring
224
mean air temperature (also water temperature), and between BSDs and spring mean pH (Fig.
225
4). Therefore, these climatic variables in spring are the key predictors of BSDs in Lake Taihu.
226 227
Fig. 4 Linear correlations between BSDs and key meteorological, chemical factors. (a): spring
228
mean air temperature; (b): spring mean wind speed; (c): spring mean N/P ratio; (d): spring
229
mean pH. The “Entire”, “Meiliang”, “Zhushan”, and “Open” in the legend stand for “entire
230
Lake Taihu”, “Meiliang Bay”, “Zhushan Bay”, “Open area”, respectively.
231
DISCUSSION
232
Advantages of MODIS Data Derived Phytoplankton Phenology
233
Although 15-year remote sensing data are not as long those by many traditional site
234
specific observations of phytoplankton phenology which generally exceed 40 years
235
However, there are at least two marked advantages for MODIS data derived phytoplankton
236
phenology. Our MODIS remote sensing data enable to observe day-by-day fine evolution of
237
phytoplankton phenology monthly or seasonal in situ observation. At present, sampling at a
238
few pre-defined stations (e.g. 32 stations carried out by TLLER) at monthly or, seasonal 12
ACS Paragon Plus Environment
6 - 8.
Page 13 of 19
Environmental Science & Technology
239
frequency is generally insufficient for phytoplankton phenology analysis in Lake Taihu because
240
the temporal and spatial resolution of these data is too low to derive BSDs. The increase of
241
observation frequency can reduce the uncertainty of derived phytoplankton phenology
242
parameters. In addition, reliable phytoplankton phenology assessments are restricted by the
243
spatial coverage of in situ observation which may miss patchy and transient phytoplankton
244
blooms. By use of remote sensing data, we found that phytoplankton phenology markedly
245
advanced in Meiliang Bay and Zhushan Bay as compared to the Open area, which was not
246
observed in previous phytoplankton phenology studies of Lake Taihu
247
phytoplankton phenology features should be very promising method of phenology study and
248
be an important continued direction of water remote sensing.
249
Roles of Climate Change and Eutrophication in Controlling Phytoplankton Phenology
22, 25, 27.
This
250
A study using Chla concentration time series during 1998–2015 over a global grid showed
251
the timing and size of phytoplankton blooms have changed on both regional and global scales
252
32.
253
small to large lakes 6, 7, 12, 33. The BSDs' advance of 2.1 days per year resulting from linear fitting
254
to the 15-year daily remote sensing observations for the entire Lake Taihu further confirmed
255
this fact. Many previous studies attributed BSDs advances to climate warming because limited
256
nutrient data make it difficult to assess the impact of nutrient concentration on BSDs. In fact,
257
in addition to climate warming inducing an increased water temperature and thermal
258
stratification, BSDs also depend upon several variables, such as light and nutrient availability,
259
wind-driven vertical mixing, grazing and predating , which influence population growth and
260
loss rates of phytoplankton 7, 8, 33, 34. Indeed, our statistical analysis showed that spring mean air
261
temperature, spring mean wind speed and spring mean N/P ratio were three key meteorological
262
and chemical factors which are responsible for the significant advance of BSDs. This is in
263
agreement with the observations by many previous studies on significantly increased air
264
temperature but wind decreased speed 28, 35. Meanwhile, N/P ratio significantly decreased from
265
2003 to 2017 and more and more reach to 16:1. The value is the most favorable for
266
cyanobacteria blooms 36. The pH value of Lake Taihu was found to have high linear correlation
267
to BSDs but pH should not be the cause of phytoplankton blooms because many studies have
268
shown that phytoplankton blooms resulted in high pH
Meanwhile, BSDs have been shown to advance in lake ecosystems from shallow to deep and
37, 38.
Therefore, we conclude that the
13
ACS Paragon Plus Environment
Environmental Science & Technology
Page 14 of 19
269
BSDs' advance is not only determined by climate warming but also nutrient availability. The
270
same observation was also reported in other studies. For example, a 58-year site specific
271
observation study for four lakes of English Lake District catchment indicated that the effect of
272
climate warming was not consistent on the three phytoplankton taxa with some showing
273
advanced BSDs but others delayed BSDs. On the contrast, soluble reactive phosphorus
274
concentration, had a more consistent effect upon the phenology of all taxa in nearly all lakes 7.
275
It is worth mentioning that we have considered various possible drivers here, but the influence
276
of parasitism and predation cannot be addressed owing to the lack of more long term data
277
records.. Therefore, future research and additional statistical matrices will help separate climate
278
vs. non-climate drivers for phytoplankton phenology.
279
Implications for Algal Bloom Management
280
The objectives of algal bloom management are to protect drinking water and fisheries,
281
minimize economic and ecosystem losses, and protect public health. Therefore, the most critical
282
step of algal bloom management is to develop an effective monitoring method and a robust
283
forecast model. Normally, remote sensing, in situ monitoring, and algal bloom forecast
284
ensuring drinking water supply were performed from May to October every year
285
study results clearly demonstrated that the BSDs of Lake Taihu have significantly advanced
286
and been driven by climate and non-climate factors. For example, a massive algal bloom event
287
began two months earlier than previously documented for Microcystis blooms in Taihu in 2007
288
due to an unusually warm spring, which caused a serious drinking water crisis 24. Therefore,
289
advanced BSDs means more material, manpower and finical resources are needed to manage
290
this lake especially in warm spring. In addition, advanced BSDs resulting from climate
291
warming could offset the already time consuming and costly management strategies that are
292
designed to control algal blooms and achieve lake restoration goal. Therefore, more strictly
293
nutrient control and reduction strategies are needed to mitigate the pollution by sewage, waste
294
discharges, fertilizer use, agricultural non-point and industrial point sources. In short, climate
295
warming poses enormous challenges for lake managers, requiring that traditional lake
296
restoration techniques and nutrient control have to be practiced and enforced strictly.
297 14
ACS Paragon Plus Environment
35, 39.
for Our
Page 15 of 19
298
Environmental Science & Technology
ACKNOWLEDGMENTS
299
This study was jointly funded by the National Natural Science Foundation of China (grants
300
41771472 and 41621002), the Key Research Program of Frontier Sciences of the Chinese
301
Academy of Sciences (QYZDB-SSW-DQC016), Jiangsu 333 Talents Program (BRA2018092),
302
the Key Program of the Chinese Academy of Sciences (ZDRW-ZS-2017-3-4), “Strategic
303
Priority Research Program” of the Chinese Academy of Sciences (XDA19070301), and Youth
304
Innovation Promotion Association (CAS) (2017365). We are grateful to all staff of TLLER for
305
assisting with sample collection, experiment measurement, and data analysis of long-term water
306
quality observations.
307 308
AUTHORS INFORMATION
309
Corresponding Author
310
Phone: (+86) -25-86882174; Email:
[email protected] 311
Notes
312
The authors declare no competing financial interest.
313 314
15
ACS Paragon Plus Environment
Environmental Science & Technology
315 316 317 318 319
REFERENCES (1) Pörtner, H. O.; Farrell, A. P. Physiology and climate change. Science 2008, 322 (5902), 690-692. (2) Parmesan, C. Ecological and evolutionary responses to recent climate change. Annu. Rev. Ecol. Evol. S. 2006, 37 (1), 637-669.
320
(3) Gian-Reto, W.; Eric, P.; Peter, C.; Annette, M.; Camille, P.; Beebee, T. J. C.;
321
Jean-Marc, F.; Ove, H. G.; Franz, B. Ecological responses to recent climate change. Nature
322
2002, 416 (6879), 389-395.
323
(4) Cleland, E. E.; Isabelle, C.; Annette, M.; Mooney, H. A.; Schwartz, M. D.
324
Shifting plant phenology in response to global change. Trends Ecol. Evol. 2007, 22 (7),
325
357-365.
326 327 328 329
(5) Martin, E.; Richardson, A. J. Impact of climate change on marine pelagic phenology and trophic mismatch. Nature 2004, 430 (7002), 881-884. (6) Winder, M.; Schindler, D. E. Climatic effects on the phenology of lake processes. Global Change Biol. 2010, 10 (11), 1844-1856.
330
(7) Feuchtmayr, H.; Thackeray, S. J.; Jones, I. D.; Ville, M. D.; Fletcher, J.; James,
331
B.; Kelly, J. Spring phytoplankton phenology – are patterns and drivers of change
332
consistent among lakes in the same climatological region? Freshwater Biol. 2012, 57 (2),
333
331-344.
334
(8) Thackeray, S.; Jones, I.; Maberly, S. Long-term change in the phenology of
335
spring phytoplankton: species-specific responses to nutrient enrichment and climatic
336
change. J. Ecol. 2010, 96 (3), 523-535.
337 338
(9) Kahru, M.; Brotas, V.; Manzano-Sarabia, M.; Mitchell, B. G. Are phytoplankton blooms occurring earlier in the Arctic? Global Change Biol. 2011, 17 (4), 1733-1739.
339
(10) White, M. A.; Beurs, K. M. D.; Didan, K.; Inouye, D. W.; Richardson, A. D.;
340
Jensen, O. P.; O'Keefe, J.; Zhang, G.; Nemani, R. R.; Leeuwen, W. J. D. V.
341
Intercomparison, interpretation, and assessment of spring phenology in North America
342
estimated from remote sensing for 1982-2006. Global Change Biol. 2009, 15 (8), 2335-
343
2359.
344
(11) Zhang, X.; Friedl, M. A.; Schaaf, C. B.; Strahler, A. H.; Hodges, J. C. F.; Gao, 16
ACS Paragon Plus Environment
Page 16 of 19
Page 17 of 19
Environmental Science & Technology
345
F.; Reed, B. C.; Huete, A. Monitoring vegetation phenology using MODIS. Remote Sens.
346
Environ. 2003, 84 (3), 471-475.
347
(12) Palmer, S. C. J.; Odermatt, D.; Hunter, P. D.; Brockmann, C.; Présing, M.;
348
Balzter, H.; Tóth, V. R. Satellite remote sensing of phytoplankton phenology in Lake
349
Balaton using 10 years of MERIS observations. Remote Sens. Environ. 2015, 158, 441-
350
452.
351
(13) Sasaoka, K.; Chiba, S.; Saino, T. Climatic forcing and phytoplankton phenology
352
over the subarctic North Pacific from 1998 to 2006, as observed from ocean color data.
353
Geophys. Res. Lett. 2011, 38 (15), 165-176.
354
(14) Fernando, G. T.; Ricardo, A. Seasonality of North Atlantic phytoplankton from
355
space: impact of environmental forcing on a changing phenology (1998-2012). Global
356
Change Biol. 2014, 20 (3), 698-712.
357 358
(15) Ove, H. G.; Bruno, J. F. The impact of climate change on the world's marine ecosystems. Science 2010, 328 (5985), 1523-1528.
359
(16) Tranvik, L. J.; Downing, J. A.; Cotner, J. B.; Loiselle, S. A.; Striegl, R. G.;
360
Ballatore, T. J.; Dillon, P.; Finlay, K.; Fortino, K.; Knoll, L. B. Lakes and reservoirs as
361
regulators of carbon cycling and climate. Limnol. Oceanogr. 2009, 54 (6), 2298-2314.
362
(17) Williamson, C. E,; Saros, J. E,; Vincent, W. F,; Smol, J. P. Lakes and reservoirs
363
as sentinels, integrators, and regulators of climate change. Limnol. Oceanogr. 2009, 54 (6),
364
2273-2282.
365 366
(18) Huisman, J.; Codd, G. A.; Paerl, H. W.; Ibelings, B. W.; Verspagen, J.; Visser, P. M. Cyanobacterial blooms. Nat. Rev. Microbiol. 2018, 16, 471-483.
367
(19) Paerl, H. W.; Huisman, J. Blooms like it hot. Science 2008, 320 (5872), 57-58.
368
(20) Wolkovich, E. M.; Cook, B. I.; Allen, J. M.; Crimmins, T. M.; Betancourt, J. L.;
369
Travers, S. E.; Pau, S.; Regetz, J.; Davies, T. J.; Kraft, N. J. B. Warming experiments
370
underpredict plant phenological responses to climate change. Nature 2012, 485 (7399),
371
494-497.
372
(21) Poloczanska, E. S.; Brown, C. J.; Sydeman, W. J.; Kiessling, W.; Schoeman, D.
373
S.; Moore, P. J.; Brander, K.; Bruno, J. F.; Buckley, L. B.; Burrows, M. T.; Duarte, C. M.;
374
Halpern, B. S.; Holding, J.; Kappel, C. V.; O’Connor, M. I.; Pandolfi, J. M.; Parmesan, C.; 17
ACS Paragon Plus Environment
Environmental Science & Technology
375
Schwing, F.; Thompson, S. A.; Richardson, A. J. Global imprint of climate change on
376
marine life. Nat. Clim. Change 2013, 3, 919-925.
377
(22) Duan, H.; Ma, R.; Xu, X.; Kong, F.; Zhang, S.; Kong, W.; Hao, J.; Shang, L.
378
Two-decade reconstruction of algal blooms in China's Lake Taihu. Environ. Sci. Technol.
379
2012, 43 (10), 3522-3528.
380 381
(23) Qin, B.; Xu, P.; Wu, Q.; Luo, L.; Zhang, Y. Environmental issues of Lake Taihu, China. Hydrobiologia 2007, 581 (1), 3-14.
382
(24) Qin, B.; Zhu, G.; Gao, G.; Zhang, Y.; Wei, L.; Paerl, H. W.; Carmichael, W. W.
383
A drinking water crisis in Lake Taihu, China: Linkage to climatic variability and lake
384
management. Environ. Manage. 2010, 45 (1), 105-112.
385
(25) Deng, J.; Qin, B.; Paerl, H. W.; Zhang, Y.; Ma, J.; Chen, Y. Earlier and warmer
386
springs increase cyanobacterial (Microcystis spp.) blooms in subtropical Lake Taihu,
387
China. Freshwater Biol. 2014, 59 (5), 1076-1085.
388
(26) Shi, K.; Zhang, Y.; Zhou, Y.; Liu, X.; Zhu, G.; Qin, B.; Gao, G. Long-term
389
MODIS observations of cyanobacterial dynamics in Lake Taihu: Responses to nutrient
390
enrichment and meteorological factors. Sci. Rep. 2017, 7, 40326; DOI: 10.1038/srep40326.
391
(27) Zhang, M.; Duan, H.; Shi, X.; Yu, Y.; Kong, F. Contributions of meteorology to
392
the phenology of cyanobacterial blooms: Implications for future climate change. Water
393
Res. 2012, 46 (2), 442-52.
394
(28) Zhang, Y.; Boqiang, Q.; Guangwei, Z.; Kun, S.; Zhou, Y. Profound changes in
395
the physical environment of Lake Taihu from 25 years of long-term observations:
396
implications for algal bloom outbreaks and aquatic macrophytes loss. Water Resour. Res.
397
2018, 54 (7), 4319-4331.
398 399
(29) Jönsson, P.; Eklundh, L. TIMESAT—a program for analyzing time-series of satellite sensor data . Comput. Geosci.-UK. 2004, 30 (8), 833-845.
400
(30) Cole, H.; Henson, S.; Martin, A.; Yool, A. Mind the gap: The impact of missing
401
data on the calculation of phytoplankton phenology metrics. J. Geophys. Res. –Oceans
402
2012, 117 (C8); DOI: 10.1029/2012JC008249.
403 404
(31) Racault, M. F.; Quéré, C. L.; Buitenhuis, E.; Sathyendranath, S.; Platt, T. Phytoplankton phenology in the global ocean. Ecol. Indic. 2012, 14 (1), 152-163. 18
ACS Paragon Plus Environment
Page 18 of 19
Page 19 of 19
Environmental Science & Technology
405
(32) Friedland, K. D.; Mouw, C. B.; Asch, R. G.; Ferreira, A. S. A.; Brady, D. C.
406
Phenology and time series trends of the dominant seasonal phytoplankton bloom across
407
global scales. Global Ecol. Biologeogr. 2018, 27 (5), 551-569.
408
(33) Meis, S.; Thackeray, S.; Jones, I. Effects of recent climate change on
409
phytoplankton phenology in a temperate lake. Freshwater Biol. 2010, 54 (9), 1888-1898.
410
(34) Henson, S. A.; Cole, H. S.; Hopkins, J.; Martin, A. P.; Yool, A. Detection of
411
climate change-driven trends in phytoplankton phenology. Global Change Biol. 2018, 24
412
(1), e101-e111.
413
(35) Zhang, Y.; Shi, K.; Liu, J.; Deng, J.; Qin, B.; Zhu, G.; Zhou, Y. Meteorological
414
and hydrological conditions driving the formation and disappearance of black blooms, an
415
ecological disaster phenomena of eutrophication and algal blooms. Sci. Total. Environ.
416
2016, 569-570, 1517-1529.
417
(36) Havens, K. E.; James, R. T.; East, T. L.; Smith, V. H. N:P ratios, light limitation,
418
and cyanobacterial dominance in a subtropical lake impacted by non-point source nutrient
419
pollution. Environ. Pollut. 2003, 122 (3), 379-390.
420
(37) Sandrini, G.; Tann, R. P.; Schuurmans, J. M.; van Beusekom, S. A.; Matthijs, H.
421
C.; Huisman, J. Diel variation in gene expression of the CO2-concentrating mechanism
422
during a harmful cyanobacterial bloom. Front. Microbiol. 2016, 7, 551; DOI:
423
10.3389/fmicb.2016.00551.
424
(38) Xie, L. Q.; Xie, P., .; Tang, H. J. Enhancement of dissolved phosphorus release
425
from sediment to lake water by Microcystis blooms--an enclosure experiment in a hyper-
426
eutrophic, subtropical Chinese lake. Environ. Pollut. 2003, 122 (3), 391-399.
427
(39) Qin, B.; Li, W.; Zhu, G.; Zhang, Y.; Wu, T.; Gao, G. Cyanobacterial bloom
428
management through integrated monitoring and forecasting in large shallow eutrophic
429
Lake Taihu (China). J. Hazard. Mater. 2015, 287 (2), 356-363.
430
19
ACS Paragon Plus Environment