Phenology of Phytoplankton Blooms in a Trophic Lake Observed from

4 days ago - Multiple linear correlation showed that air temperature, wind speed, and ... assessment of phytoplankton phenological shifts and elucidat...
0 downloads 0 Views 1MB Size
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