SciELO - Scientific Electronic Library Online

 
vol.69 número1Variación corporal paralela en peces de dos ríos costeros del Chocó ecuatorianoÁreas en riesgo de invasión por Procambarus clarkii (Decapoda: Cambaridae) un cangrejo de río introducido en Colombia índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

Compartir


Revista de Biología Tropical

versión On-line ISSN 0034-7744versión impresa ISSN 0034-7744

Rev. biol. trop vol.69 no.1 San José ene./mar. 2021

http://dx.doi.org/10.15517/rbt.v69i1.42731 

Artículo

Effect of climato-environmental parameters on chlorophyll a concentration in the lower Ganga basin, India

Efecto de los parámetros climáticos ambientales sobre la concentración de clorofila a en la cuenca baja del Ganges, India

Soma-Das Sarkar1 

Uttam-Kumar Sarkar1 

Malay Naskar1 

Koushik Roy1  2 

Arun-Kumar Bose1 

Subir-Kumar Nag1 

Gunjan Karnatak1 

Basanta-Kumar Das1 

1. ICAR-Central Inland Fisheries Research Institute, S.N. Banerjee Road, Barrackpore, Kolkata-700120, West Bengal, India; soma.das@icar.gov.in, uksarkar1@gmail.com (*Correspondence), malaynaskar@icar.org.in, koushik.roy.89@gmail.com, arun_bose06@yahoo.co.in; nagsk_67@rediffmail.com, gunjankarnatak87@gmail.com, basantakumard@gmail.com

2. Faculty of Fisheries and Protection of Waters, University of South Bohemia, Ceske Budejovice 37005, Czech Republic.

Abstract

Introduction:

Chlorophyll a concentration proxies the phytoplankton biomass which directly involves in signifying the production functions of aquatic ecosystem. Thus, it is imperative to understand their spatio-temporal kinetics in lotic environment with reference to regional climatic variabilities in the tropical inland waters.

Objective:

In-situ studies were conducted to examine the changes in phytoplankton biomass in lower Ganga basin as influenced by various environmental parameters under regional climatic variability during 2014-2016.

Methods:

Firstly, the most key influential environmental parameters on riverine Chl-a concentration were determined. Then the direct cascading effect of changing climatic variables on key environmental parameters were derived through modeling and quantified probable changes in mean Chl-a concentration in the lower stretch of river.

Results:

Only five environmental parameters namely water temperature, total dissolved solid, salinity, total alkalinity and pH were key factors influencing Chl-a (Multiple R2: 0.638, P < 0.05). Present estimates indicate that if the present rate of regional climatic variability over the last 3 decades (mean air temperature + 0.24 °C, total annual rainfall -196.3 mm) remain consistent over the next three decades (2015-2045), an increase in mean Chl-a by + 170 µgL-1 may likely be expected grossly reaching about 475.94 µg L-1 by the year 2045 or more.

Conclusions:

The present study is first such comprehending a gross hint towards the probable ecosystem response with an alternative model based methodology in data-deficient situations. Subsequently, the output would also be of great benefit for increase water governance and developing strategy protocol for sustainable water management for greater ecosystem services.

Key words: chlorophyll a; climate change; environmental variable; predictive modeling; River Ganga

Resumen

Introducción:

La concentración de clorofila a representa la biomasa de fitoplancton la cual influye directamente en las funciones de producción de los ecosistemas acuáticos. Por lo tanto, es imperativo comprender su cinética espacio-temporal en el ambiente lótico con respecto a las variabilidades climáticas regionales en las aguas continentales tropicales.

Objetivo:

Se realizaron estudios in situ para examinar la influencia de varios parámetros ambientales en la biomasa del fitoplancton en la cuenca baja del Ganges durante 2014-2016.

Métodos:

En primer lugar, se determinaron los parámetros ambientales más influyentes en la concentración de Chl-a fluvial. Luego, el efecto directo en cascada de las variables climáticas sobre los parámetros ambientales clave, mediante el modelado y los cambios en la concentración media de Chl-a en el tramo inferior del río.

Resultados:

Solo cinco parámetros ambientales, entre ellos, temperatura del agua, sólidos disueltos totales, salinidad, alcalinidad total y pH, fueron factores clave que influyeron en la Chl-a (R2 múltiple: 0.638, P < 0.05). Las estimaciones actuales indican que si la tasa actual de variabilidad climática regional durante las últimas 3 décadas (temperatura media del aire + 0.24 °C, precipitación total anual -196.3 mm) permanece constante durante las próximas tres décadas (2015-2045), se presente un aumento en el promedio de la Chl-a en +170 µgL-1 y alcance aproximadamente 475.94 µgL-1 para el 2045 o más.

Conclusiones:

Este estudio presenta una metodología basada en modelos alternativos en situaciones de escasez de datos, la información generada también podría contribuir a mejorar la gobernanza del agua y a desarrollar un protocolo para la gestión sostenible del agua y de esta manera mejorar los servicios ecosistémicos.

Palabras clave: clorofila a; cambio climático; variable ambiental; modelado predictivo; río Ganges

Introduction

Tropical riverine fisheries support millions of people with food, nutrition, livelihood etc. (Vass et al., 2011; Sarkar et al., 2019) and hence contributing indirectly to a large extent to the national economy of the associated countries, especially India. The Ganga river basin has already witnessed anthropogenic pressure due to deforestation, water abstraction, irrational fishing, mining, industrialization, economic development, navigational movements (Das Sarkar et al., 2019a). Changes in the climatic regime such as increased water temperature, uncertain monsoon variation and increased rate of abnormal climatic events in the coastal area have been the recent addition. Synergistically, these affect the ecosystem biota, productivity and ultimately livelihood security (Sharma, Joshi, Naskar, & Das, 2015) of poor fishermen.

Primary production in an aquatic system depends on phytoplankton biomass, the food for planktivorous fishes. It is also a parameter hinting overall system productivity. Chlorophyll is the green colouring matter which enables photosynthesis to occur in the presence of light. Out of nine different types of chlorophyll, chlorophyll-a plays a major role in photosynthesis (Ji, 2017) and hence indicates phytoplankton biomass (Gregor & Marsalek, 2004). Although literature surveys indicate several studies on ecology and biodiversity of river Ganges (Manna et al., 2013; Roshith et al., 2018); however, the role of climatic variability and their impact on aquatic ecosystem using different biota are very scanty. Thus, in this rising scenario of global warming it is imperative to understand the impact of climatic variability on primary producers. Microalgae or phytoplankton form a valuable biological parameter in climate models (Ghosal, Rogers, & Wray, 2000). Climate change itself represents a complex combination of stressors, including alterations in temperature (Webb, Hannah, Moore, Brown, & Nobilis, 2008), elevated atmospheric CO2 (Field, Barros, Mach, & Mastrandrea, 2014), and increased frequency and intensity of droughts and extreme flow events (Barnett, Adam, & Lettenmaier, 2005; Milly, Dunne, & Vecchia, 2006). Changes in phytoplankton biomass vis-à-vis primary productivity in aquatic ecosystems resulting from climate change are often attributable to increased water temperatures and/or increased nutrient loading (Wrona et al., 2006). In recent years, the focus has been on predictive and/or trend identification modeling with either isolated or combined influence of an environmental driver or a set of drivers on phytoplankton biomass in rivers. The models developed by Elliot, Thackeray, Huntingford, & Jones (2005), Mooij, Janse, Domis, Hulsmann, & Ibelings (2007), Friberg et al. (2009) are some of the examples where the implications of warming water temperature on riverine phytoplankton biomass have been modeled. Such models formed the basis of the prediction made in Field et al. (2014) regarding the increase of primary productivity and tendency of eutrophication in majority of the global aquatic ecosystems.

Review of literature has revealed that studies on the climate driven changes on primary productivity and/or phytoplankton biomass in rivers have been mostly concentrated in temperate to sub-tropical climate (Klapper, 1991; Friberg et al., 2009; Sipkay, Kiss-Keve1, Vadadi-Fulop, Homoródi, & Hufnagel, 2012; Cloern, Foster, & Kleckner, 2014) and studies on tropically situated Indian rivers is absent. The only study in Indian context was in Bay of Bengal and Arabian sea that too not on the exact hypothesis per se but a related conclusion was drawn from their expert knowledge and extrapolating their results (Chaturvedi, Meghal, & Jasraj, 2013). Authors have hypothesized increased decadal changes of Chl concentration due to high nutrient influx through freshwater river drainage and increase in temperature due to global warming. However, authors have not developed quantitative relationship between Chl concentration and climate-environmental parameters. Although majority have reported an increase in primary production and/or phytoplankton biomass in rivers under future changed climate scenario, but few studies have indicated the opposite i.e. decreased riverine phytoplankton population under changed climatic scenario (Lewandowska & Sommer, 2010; Sommer & Lengfellner, 2008).

Modeling studies often face the problems of lack of required data and access to them which is considered one of the primary hindrances in climate change related modeling research (Porter et al., 2005; Sipkay, Kiss, Vadadi-Fulop, & Hufnagel, 2009). Developing countries like India have dearth of long-term data on riverine chlorophyll concentration or any aquatic system per se; although time series data of climatic variables are well maintained and abundantly available. The present approach explored an innovative and alternative methodology which could generate baseline information on trends and gross quantification of changes in Chl-a concentration in absence of long term data. Therefore, in the absence of long term and continuous data on Chlorophyll-a concentration in River Ganga, it is worthwhile to adopt a model-based framework to identify the key climato-environmental parameters and their isolated individual influences on phytoplankton biomass (Chl-a concentration); thereby, budgeted the Chl-a concentration considering the multiple linear effect of all the key environmental parameters identified. For budgeting, the direct cascading effect of regional climate change (increasing mean air temperature and decreasing rainfall) on influential environmental parameters were first quantified individually and then plugged together to obtain the estimates of changes expected in Chl-a concentration as per the present trend of changing climate within the study area.

In the present study, the photosynthetic pigment ‘Chlorophyll-a’ (henceforth abbreviated as Chl-a) was used as a potential marker that reflects phytoplankton biomass vis-à-vis primary production potential of the aquatic system. With this background, the study was conducted to examine the changes in phytoplankton biomass in lower Ganga basin as influenced by various environmental parameters under regional trend of climatic variability. The objective of this study was to first identify the most influential environmental parameters on riverine chlorophyll-a concentration and to quantify the direct cascading effect of changing climatic variables on key environmental parameters through modeling for deriving the expected changes in Chl-a concentration in the lower Ganga basin.

Materials and methods

Study area: The lower stretch of river Ganga (locally known as river Hooghly) was selected in this study. Monthly water samples were collected in duplicate during October, 2014 to September, 2016 from two different sites viz. Triveni (22°59’411’’ N & 88°24’310’’E) and Godakhali (22°23’629’’ N & 88°07’948’ ’E) as shown in Fig. 1. Both the stations are part of posterior lower stretch of River Ganga. Each station was further sub-divided into two sub-stations (Triveni-1, Triveni-2 and Godakhali-1, Godakhali-2) to obtain representative samples.

Fig. 1 Geographical location of the study area. Blue line indicates the longitudinal course of river Ganga. 

Water sampling, analysis and estimation of Chl-a: Water samples in duplicate were collected from sub-surface region of the water column in 1 Liter HDPE amber bottles. A total of twelve physico-chemical parameters viz. water temperature (°C), dissolved oxygen (DO), free carbon-di-oxide (CO2), pH, electrical conductivity (EC), total hardness, total dissolved solid (TDS), total alkalinity (TA), salinity, Secchi-disc transparency, nitrate and phosphate were analyzed following Standard Methods (APHA, 2012).

Samples for estimation of chlorophyll -a (Chl-a) were collected separately and estimated following APHA (2012). Samples from each replicate were filtered in vacuum with a membrane filters and Millipore membrane filter of pore size 0.45 µm. The residues were wrapped in aluminum foil and transported to laboratory using refrigeration unit in a mobile laboratory van and frozen (-20 °C) immediately upon reaching laboratory. Chl-a pigment was extracted with acetone using a digital laboratory tissue homogenizer (RQT-127AD), Remi make, India. The absorbance of the extracts at wave length 664 nm, 647 nm and 630 nm were recorded in a UV-VIS spectrophotometer. Concentration of Chl-a in water was calculated using following equation: Chl-a = 11.85 (OD 664)-1.54 (OD 647)-0.08 (OD 630), where, Chl-a = concentration of Chl-a in µgL-1, also known as ppb (parts per billion).

Climatic data analysis: Standard “30-years normal” data analysis procedure was used. Daily mean air temperature and precipitation data of 30 years from 1980 to 2015 was obtained from Indian Meteorological Department (IMD). The data corresponding to study location was extracted from 1° x 1° data by using ‘Raster package’ in ‘R’ (R Core Team, 2015). Subsequently monthly means of air temperature and total rainfall was computed from the extracted daily time series data. Finally, annual mean air temperature and total annual rainfall were computed from monthly data. A simple linear trend analysis was carried out to find out the long-term trend. The total rainfall during pre-monsoon, monsoon and post-monsoon) were used to quantify changes in seasonal rainfall pattern (Karnatak et al., 2018).

Statistical analyses and modelling

Modeling climato-environmental influence on Chl-a concentration: The data was analyzed in multiple steps, employing various statistical tools. In order to identify the key climato-environmental parameters influencing Chl-a concentration in river, a stepwise regression with collinearity diagnostics was applied, designating all environmental parameters as explanatory variables and Chl-a concentration as the response variable. Further the relative importance of each environmental parameter (influencing Chl-a) in the optimized linear model was computed in R software (R Core Team, 2015) by using the computer intensive and most recommended method (Lindeman, Merenda, & Gold, 1980). Finally, looking beyond linearity, non-linear Generalized Additive Model (GAM) (Wood, 2006) was applied to unravel the pattern of individual influence of the key environmental parameters on Chl-a.

Budgeting changes in Chl-a concentration under present rate of climate change: The cascading effect of variations in climatic factors (mean air temperature and rainfall) on the key environmental parameters influencing Chl-a concentration was quantified. To avoid complexity in interpretation, only the direct effect of climatic factors on the key environmental parameters was quantified irrespective of the interactions among other non-key parameters existing in the aquatic ecosystem. This was accomplished through a two-step strategy. Firstly, the relationship between the key environmental parameter and climatic factors (mean air temperature and rainfall) was established by simple regression, taking climatic parameters as independent variable and environmental parameters as dependent variables in order to quantify the nature and degree of change in environmental parameters per unit change of climatic factors, irrespective of spatial differences. Secondly, changes in mean air temperature and rainfall over the sampling stations during last three decades (1980-2015) was assessed using IMD data and the present rate of climate change was assumed to remain true in the next three decades (2015-2045). The present rate of change in mean air temperature and rainfall which may cause direct changes in the key environmental parameters ultimately manifesting into change in Chl-a concentration has been derived empirically by plugging the estimates from regression models developed. Finally, these predictive values of environmental variables were utilized in the optimized regression model for Chl-a to derive predictive changes in Chl-a concentration.

Results

Climate change trend: Analysis of the IMD data on mean air temperature has revealed a unanimous warming trend among the studied sites. During the period of 1980-2015 the mean air temperature has increased by 0.24 °C (Fig. 2). It has also been recorded that the total annual rainfall during the same period has decreased by196.3 mm (Fig. 3).

Fig. 2 Trend of increasing mean air temperature along lower stretch of river Ganga in West Bengal. 

Fig. 3 Decreasing trend of total annual rainfall along lower stretch of river Ganga in West Bengal. 

Modeling climato-environmental influence on Chl-a concentration

Present observations: During the study period, chlorophyll-a (Chl-a) concentration were ranged between 18.92-1 565.45 µgL-1 with an annual mean of 305.94 ± 52.25 µgL-1. The annual distribution pattern of Chl-a encountered is depicted in Digital appendix 1. The inter-site variability (spatial) in Chl-a concentration was statistically insignificant (p = 0.679) as revealed by ANCOVA. Likewise, Chl-a concentration did not vary significantly between seasons pre-monsoon, monsoon and post-monsoon (Pre-M vs. M P = 0.879; Post-M vs. M P = 0.825). The overall descriptive account of physico-chemical and climatic parameters recorded during the study period is presented in Table 1.

Table 1 Descriptive account of habitat parameters recorded in River Ganga during the study period 

Parameters Stretch#
Water temperature (°C) 19.8-34.0 (27.64 ± 0.58)
Dissolved oxygen (ppm) 1.6-10 (4.51 ± 0.28)
Free carbon dioxide (ppm) 0-10 (3.95 ± 0.5)
pH (units) 7.5-8.6 (8.1 ± 0.04)
Total Dissolved Solids (ppm) 113-301 (228.76 ± 6.28)
Conductivity (µmhos) 160.1-432 (320.89 ± 8.87)
Alkalinity (ppm) 40-140 (100.41 ± 2.88)
Salinity (ppm) 80.3-700 (202.86 ± 18.96)
Nitrate (ppm) Trace-1.72 (0.84 ± 0.07)
Phosphate (ppm) Trace-0.56 (0.1 ± 0.07)
Rainfall (mm) 0-882.5 (150.47 ± 33.09)
Chlorophyll-a (µgL-1) 18.92-1565.45 (305.94 ± 52.25)

Data presented in range and mean ± standard error of mean (in parentheses).

The chlorophyll budgeting framework used in the present study therefore took into the account merged dataset of two sampling locations since no spatial variability exists. An attempt was also made to chalk out the seasonal thresholds (cut-offs) of Chl-a concentration (if any) to demarcate signature minima or range of Chl-a concentration characteristic in any season; be it pre-monsoon (Feb-May), monsoon (Jun-Sept) or post-monsoon (Oct-Jan). These thresholds were supposed to be super imposed over thermal and precipitation windows to derive season specific climate optima (temperature and rainfall) within which such thresholds are attained in natural conditions. All this information was assumed to serve better in regional climate change impact assessment over lower Gangetic basin on riverine ecological productivity. Unfortunately, no such season-specific thresholds could be identified and therefore, any further analysis involving seasonal influence on riverine chlorophyll was purposively skipped. Following this, only annual data like annual mean chlorophyll concentration were forecasted through the budgeting framework; not the seasonal changes in Chl-a over next three decades.

Cumulative and overall influence of habitat parameters on Chl-a concentration: The study indicated that among 14 climato-environmental parameters considered (mean air temperature, water temperature, rainfall, secchi disc transparency, pH, dissolved oxygen, free carbon dioxide (CO2), total alkalinity, electrical conductivity, total dissolved solids, phosphate, nitrate and salinity), only five environmental parameters namely water temperature, TDS, salinity, TA and pH were found to be key factors influencing Chl-a concentration in the river as revealed through stepwise regression. The generated optimized model was Chl-a = -5.796 + (0.00521 * Total Alkalinity) + (0.347 * pH) - (0.000389 * Salinity) + (0.00451 * TDS) + (0.0652 * Water temperature), which could explain 63.8 % variability of the chlorophyll-a concentration. Water temperature (P < 0.01), total dissolved solids (P < 0.01), total alkalinity (P < 0.05) and pH (P < 0.05) had significant positive influence on Chl-a concentration while salinity (P < 0.05) had a significant negative influence on the same.

In terms of the relative importance of individual environmental parameters (Fig. 4) in influencing Chl-a concentration, water temperature (33.13 % of R2) was the most important while place effect was least important (3.4 % of R2).The relative importance of environmental parameters influencing Chl-a concentration of the river can be summarized as: water temperature > alkalinity > TDS > pH > salinity > place-effect; depicted in Fig 4. The ‘place effect’ diagnostics in the model reveal that Chl-a concentration in the river had no significant spatial difference between two sites. More specifically, there is no significant spatial difference between the two sites of River Ganga in terms of environmental parameters-Chl-a relationship.

Fig. 4 Relative importance of key environmental parameters for explaining variability in chlorophyll a concentration. 

Individual influence pattern of habitat parameters on Chl-a concentration: Individual influence(s) of identified key environmental parameters on Chl-a concentration have been discussed in descending order of their importance.

Temperature: GAM models showed that Chl-a concentration in rivers can be expected to initially increase with rise in water temperature (Fig. 5) but may reach a plateau beyond 30 °C (Deviance explained: 12.2 %). It is emphasized that the present model might be useful mostly for pattern identification and may not be optimum for prediction purposes due to short period of study. The prediction belt in the model shows Chl-a concentration may either rise or fall when water warms up beyond 34 °C which requires further study to draw an inference. A gross statement can be made on the fact that although global warming manifested into warmer water temperatures round the year may initially result in increased phytoplankton population but in the long run this effect might not prove beneficial. The study hints at possibility of increased phytoplankton population in the coming years due to enhanced thermal regime mostly during winter. Water temperature was also found to be the most important driver for influencing Chl-a concentration and has been given special emphasis during model interpretation and discussion in climate change perspective; assumed as primary driver.

Fig. 5 Influence of water temperature on chlorophyll a concentration. 

Total Alkalinity: A progressive increase in Chl-a concentration with increasing alkalinity of river water has been predicted by GAM model (Deviance explained: 39.6 %) (Fig. 6). Furthermore, not much change in phytoplankton population can be expected initially when the total alkalinity values remain below 100 ppm while as the total alkalinity in river increases beyond 100 ppm a steep response in phytoplankton biomass is noticeable due to a progressive increase in Chl-a concentration. Alkalinity values > 100 ppm exists in the lower stretch of River Ganga especially throughout the winter followed by early (February to March) and midsummer (March to April). An increasing trend in mean air temperature and decreasing rainfall in the area of study indicates a likely increase in average alkalinity values of river water in the coming years manifested into increased phytoplankton biomass.

Fig. 6 Influence of total alkalinity on chlorophyll a concentration. 

Total dissolved solids (TDS): No particular trend could be observed for change in Chl-a concentration with increase in TDS owing to the multiple peaks in the GAM model but prediction belt in the model shows extreme TDS > 300 ppm and < 175 ppm may be considered as the critical limit beyond which steep fall in Chl-a concentration can be expected (Deviance explained: 72.9 %) (Fig. 7). Eliminating the extremes of TDS, a moderate value of 210-230 ppm TDS may be considered optimum for healthy Chl-a concentration. This implies that high TDS in river during intense summer and diluted TDS due to high rainfall are both detrimental for phytoplankton population.

Fig. 7 Influence of TDS on chlorophyll a concentration. 

pH: Similarly, for pH-Chl-a relationship indicated that water pH might have a linear influence as the non linear influence seems to be weak (Deviance explained: 5.29 %) (Fig. 8). From the model it can be observed that response of phytoplankton biomass with changes in pH of river water alone is very subtle and the model does not permit any conclusion in this regard. The purpose of the model generated is only for pattern identification.

Fig. 8 Influence of water pH on chlorophyll a concentration. 

Salinity: The influence of salinity on Chl-a concentration showed that salinity may have a linear influence owing to the near straight curve symmetry and broad prediction belt of the model at higher salinities implying low confidence (Deviance explained: 3.43 %) (Fig. 9). It is highlighted that the model might be useful mostly for pattern identification and may not be most favorable for prediction purposes. From the model it is evident that response of phytoplankton population in river with changes in salinity alone is very subtle. It was further observed that the phytoplankton population may likely remain stable as the salinity in the river increases up to around 600 ppm under increased temperature and reduced rainfall scenario. Beyond 600 ppm salinity a probable reduction in Chl-a concentration vis-à-vis phytoplankton biomass might occur, which can happen especially during summer.

Fig. 9 Influence of salinity on chlorophyll a concentration. 

Budgeting changes in Chl-a concentration under present rate of climatic variability

Budgeting results: Change in mean annual temperature over last three decades was considered in the budgeting framework but not the minimum or maximum air temperature data (Winslow, Read, Hansen, Rose, & Robertson, 2017). Likewise, changes in total annual rainfall over last three decades were considered but not the seasonal, minima or maxima. Accordingly, the forecasting was done on mean annual Chl-a basis (elaborated above).

A trend of 0.24 °C rise in mean air temperature in West Bengal may reflect in increased mean water temperature by + 0.20 °C (R2 = 0.843, Fig. 10) in the coming years. This alone may lead to an increase in mean Chl-a concentration by + 10µgL-1in the estimates of multiple linear model (Multiple R2 = 0.684) as deduced from the linear models. Likewise, a 0.20 °C increase in mean water temperature (R2 = 0.377) and 196.3 mm decrease in rainfall (R2 = 0.577) may result in a net gain of average TDS by + 26.21 ppm which may alone lead to increased mean Chl-a concentration by +110 µgL-1 (Multiple R2 = 0.684) in the coming years. There is a possibility that, 196.3 mm decrease in rainfall over the area may probably bring an increase in average salinity by + 68.51 ppm (R2 = 0.261, low confidence) and average total alkalinity by + 9.03 ppm (R2 = 0.292, low confidence) in the coming years that may individually manifest changes in mean Chl-a concentration by -20 µgL-1 (Multiple R2 = 0.684) and +50 µgL-1 (Multiple R2 = 0.684), respectively. Lastly, a 0.20 °C increase in mean water temperature (R2 = 0.421) and 196.3 mm decrease in rainfall (R2 = 0.335) may be reflected through a net gain in average pH by + 0.09 units that may alone increase mean Chl-a concentration by + 30 µgL-1 (Multiple R2 = 0.684) in the coming years.

Fig. 10 Air-water temperature regression in the river. 

Inference of Budgeting: Keeping in view the rise in mean air temperature (+ 0.24 °C) and decrease in rainfall (196.3 mm) during the last 3 decades (1980-2015) over the lower stretch of river Ganga passing through West Bengal, India (Fig. 1, Fig. 2), we assume that this present rate of climate change might remain true in the next three decades (2015-2045) . Based on this assumption and quantified cascading effect of such climatic change on the key environmental parameters influencing Chl-a concentration, our model summarizes that a gradual increase in the existing mean Chl-a concentration (305.94 ± 52.25 µg L-1) by approximately + 170 µgL-1 (ppb) (Multiple R2 = 0.6839) in the lower Ganga basin may probably occur over the next three decades (2015-2045).

Discussion

Importance of Chl-a mapping and environmental influence: Chl-a indicates health of water in aquatic ecosystem (Boyer, Christopher, Peter, & David, 2009). It also provides trophic status based on primary productivity of the ecosystem (Das Sarkar et al., 2020) and measures trophy index (Bbalali, Hoseini, Ghorbani, & Kordi, 2013). Ecohydrological interaction studies in natural river system showed that physical variables such as, water discharge, suspended solids and turbidity regulate the development of phytoplankton biomass critically in comparison with nutrients (Salmaso & Braioni, 2008).

In the present study, only five environmental parameters namely water temperature, TDS, salinity, TA and pH were found to be key factors influencing Chl-a concentration in the river. The nature of their individual influence was found to be in conformity with published literatures. Rising temperature has aided to increase growth and productivity of algae which may be boosted by the higher enzymatic activity (Nweze & Ude, 2013). The changing temperature and precipitation pattern over lower Gangetic stretch are in the line with Rathore, Attri, & Jaswal (2013), Paul & Birthal (2015), Naskar, Roy, Karnatak, Nandy, & Roy (2017). High nutrient influx through freshwater river drainage in lower Ganga has supported rising trend of chlorophyll concentration. Similar positive impact of elevated water temperature for short term on phytoplankton biomass productivity and photosynthetic activity is also suggested tropical rice fields (Roger & Kulasooriya, 1980) and in subtropical waters (Willis, Chuang, Orr, Beardall, & Burford, 2019). A positive correlation was found between alkalinity and TDS with Chl-a concentration (P < 0.01) in the present study. The study further denotes a positive correlation (P < 0.05) between Chl-a and pH which may be due to higher phytoplankton photosynthesis and biological productivity. Simultaneously, time series models also showed influence of nutrients and pH on phytoplankton biomass for shorter time span (Jeong, Kim, Whigham, & Joo, 2003). Salinity (P < 0.05) ingression shows a significant negative influence on the Chl-a concentration in the lower Ganga which may be due to inability in generating turgor pressure to produce intracellular osmolarity needed with increasing salinity (Mitra, Zaman, & Raha, 2014). However, present study of identifying the relative influence of environmental parameters on Chl-a concentration shows water temperature (33.13 % of R2) as the most significant variable followed by TA (29.5 % of R2) and TDS (23.1 % of R2). Rainfall plays a vital role in mixing of environmental variables and also responsible for physical dispersal of phytoplankton in wetlands of tropical semiarid region (Brasil et al., 2020), hence resulting into considerable reduction in chlorophyll concentration. At the same time, rainfall did not contribute as influential climatological parameter in the present study. Taken together, it could be interesting to note that this study has also sketched out the relative importance of climato-environmental parameters influencing Chl-a concentration in inland open waters which is unavailable in most of the studies consulted. Therefore, the present information on relative importance may also be regarded as one of the recent additions towards better understanding of phytoplankton response to environmental dynamics in riverine ecosystems.

Relatively high dependence on water temperature and relevance for climate research: The prime importance of water temperature besides other environmental parameters on dictating phytoplankton production or phytoplankton biomass in lotic systems such as rivers is well documented (Schabhuttl et al., 2013; Wrona et al., 2006). Climate change is altering the physical, chemical and biological characteristics of freshwater habitats (Hartmann et al., 2013) including inland fish and fisheries (Carlson & Lederman, 2016). The importance of temperature over other water quality parameters in influencing Chl-a concentration observed during the present study has both biological and ecological basis. According to Cloern et al. (2014), phytoplankton growth rate is regulated by water temperature, nutrient concentrations and forms, and the amount and quality of photosynthetically available radiation (PAR).

From the ecological perspective, temperature regulates ecosystem functioning directly by influencing primary production (Dallas, 2008). Studies have shown that increase temperature is accompanied by a shift in the dominance of phytoplankton groups from diatoms (< 20 °C) to green algae (15 to 30 °C) to blue-green algae (> 30°C) (DeNicola, 1996). According to DeNicola (1996), phytoplankton biomass increases with temperature from approximately 0 to 30 °C and decreases at 30 to 34 °C; while species diversity increases from 0 to 25 °C and decreases at temperatures > 30 °C. This was further confirmed in a lab-based experiment by Schabhuttl et al. (2013) where it was observed that temperature is among the major determinants to influence phytoplankton growth rates positively but higher temperatures resulted in higher fractions of blue green algae over diatoms and green algae.

In the present study, although a gradual increase in Chl-a concentration was observed between 20-30 °C but a plateau was reached between 30-34 °C instead of directly decreasing. However, the prediction belt in our temperature-chlorophyll predictive model also hinted the probability of either decrease or increase in Chl-a concentration beyond 30 °C which might be attributed to the role of other key environmental parameters that significantly influence phytoplankton biomass as well as in the presence of increasing temperature, as concluded in Sipkay et al. (2012). Several modeling methodologies have also forecasted increase of phytoplankton biomass in rivers under increasing temperatures (Klapper, 1991; Elliot et al., 2005; Komatsu, Fukushima, & Shiraishi, 2006; Mooij et al., 2007; Sipkay et al., 2012). Similarly, predictive models have also forecasted the broad-spectrum increase in picophytoplankton with rise in temperature in subtropical gyres (Agusti, Lubián, Moreno-Ostos, Estrada, & Duarte 2019). However, few studies also report otherwise i.e. phytoplankton biomass in rivers to decrease with rising temperature (Sommer & Lengfellner, 2008; Lewandowska & Sommer, 2010).

Implications of regional climatic variability on riverine Chl-a: collating the present observations: The changing climate has implications on thermal regime and nutrient availability in aquatic ecosystems through changes in temperature and precipitation (flow) pattern, that is manifested into altered phytoplankton biomass or primary productivity (Webb & Nobilis, 2007; Schabhuttl et al., 2013; Das Sarkar et al., 2019b). According to Whitehead et al. (2012), in many rivers the phytoplankton blooms in summers are a feature of river ecology and the frequency and intensity of these may increase under future climatic scenario. On the other hand, study suggests less pronounced for tropical phytoplankton communities already acclimatized to warmer temperatures round the year (Schabhuttl et al., 2013). Situations of higher water temperatures and lower flow have been attributed to enhance phytoplankton growth in rivers (Verdonschot et al., 2010; Whitehead et al., 2012).

The present study has also hinted evidences of increased phytoplankton biomass viz. Chl-a concentration under increasing temperature and decreasing rainfall over the region. Additionally, estimates indicate that if the present rate of climatic variability over the last 3 decades (mean air temperature + 0.24 °C and total annual rainfall -196.3 mm) remain consistent over the next three decades (2015-2045) as well, an increase in mean Chl-a concentration by +170 µgL-1 may likely be expected. The resulting mean Chl-a concentration may grossly arrive at around 475.94 µgL-1 over the next three decades. Further repetitive time-series investigations are required to validate this estimate. It remains to be a researchable issue whether the dynamics of this change will follow a direct or alternating waveform pattern. The present study also could not conclude whether the degree of quantified change of phytoplankton biomass in river system by 2045 or more could be tagged as progressively eutrophic or lie within the mesotrophic range; difficult to probe without flow data. Nonetheless, our observations with short-term data are in line with several others who have forecasted increased phytoplankton biomass and/or primary production in rivers due to climate change using long-term data (Wrona et al., 2006; Friberg et al., 2009; Verdonschot et al., 2010; Sipkay et al., 2012; Whitehead et al., 2012). This increase in phytoplankton biomass projected for lower Ganga basin under future climatic scenario can be partly explained through the conceptual mechanism discussed in Whitehead et al. (2012) and inferences drawn by Sipkay et al. (2012). The limitation of DO in lotic systems is not a serious issue due to its flowing nature, the observed trend of increased Chl-a concentration under future climatic scenario may only be beneficial from fisheries point of view (increased abundance and biomass of fish food organisms) if, (a) the phytoplankton succession does not get significantly dominated by blue green algae (tropical waters at low risk as explained (Schabhuttl et al., 2013) and, (b) water stress does not become evident resulting into eutrophication.

The present study pioneers the ecological research in Indian rivers which conceptualized a model based framework and identified key climato-hydrochemical parameters for water managers dealing with complex water ecosystems. The study signifies the model based assessment of the impact of climatogical variables on phytoplankton biomass which stands as an emerging topic for environmental biologist and ecologists. Developing countries like India have dearth of long-term data on riverine chlorophyll concentration or any aquatic system per se; although time series data of climatic variables are well maintained and abundantly available. The present methodology explored an alternative approach to generate baseline information on trends and gross quantification of changes in aquatic Chl-a concentration under future climatic scenario -in absence of long term data. In the present study, the photosynthetic pigment Chl-a was used as a potential marker that reflects phytoplankton biomass vis-à-vis primary production potential of the aquatic system. Looking at the multitude of environmental influence on dictating riverine chlorophyll concentration, it is quite difficult to pinpoint any climate change related hypotheses in open waters. The present study may not be a key estimate of the future climate change scenario but certainly a gross hint of the probable ecosystem response and would fit for developing managerial strategies for water managers which augment the sustainable management practices for stakeholders.

Ethical statement: authors declare that they all agree with this publication and made significant contributions; that there is no conflict of interest of any kind; and that we followed all pertinent ethical and legal procedures and requirements. All financial sources are fully and clearly stated in the acknowledgements section. A signed document has been filed in the journal archives.

Acknowledgments

Authors are thankful to the Indian Council of Agricultural Research, New Delhi for financial help under project National Innovations on Climate Resilient Agriculture (NICRA). We duly acknowledge the support and cooperation of M. Prabhakar, PI NICRA, ICAR-CRIDA, Hyderabad for all technical facilities. We acknowledge the efforts rendered by Sanjeev Kumar Sahu, ICAR-CIFRI, Barrackpore for schematic representation of the study area using GIS tools.

References

Agusti, S., Lubián, L.M., Moreno-Ostos, E., Estrada, M., & Duarte, C.M. (2019). Projected changes in photosynthetic picoplankton in a warmer subtropical ocean. Frontiers in Marine Science, 5, 506. [ Links ]

APHA. (2012). Standard methods for the examination of water and wastewater (22nd Ed.). Washington, USA: American Public Health Association. [ Links ]

Barnett, T.P., Adam, J.C., & Lettenmaier, D.P. (2005). Potential impacts of a warming climate on water availability in snow-dominated regions. Nature, 438, 303-309. DOI: 10.1038/nature04141 [ Links ]

Bbalali, S., Hoseini, S.A., Ghorbani, R., & Kordi, H. (2013). Relationships between Nutrients and Chlorophyll a Concentration in the International Alma Gol Wetland, Iran. Journal of Aquaculture Research and Development, 4, 173. DOI: 10.4172/2155-9546.1000173 [ Links ]

Boyer, J.N., Christopher, R.K., Peter, B.O., & David, T.R. (2009). Phytoplankton bloom status: Chlorophyll a biomass as an indicator of water quality condition in the southern estuaries of Florida, USA. Ecological Indicators, 9, 56-67. DOI: 10.1016/j.ecolind.2008.11.013 [ Links ]

Brasil, J., Santos, J.B., Sousa, W., Menezes, R.F., Huszar, V.L., & Attayde, J.L. (2020). Rainfall leads to habitat homogenization and facilitates plankton dispersal in tropical semiarid lakes. Aquatic Ecology, 54(1), 225-241. [ Links ]

Carlson, A.K., & Lederman, N.J. (2016). Climate Change and Fisheries Education. Fisheries, 41, 411-412. DOI: 10.1080/03632415.2016.1182510 [ Links ]

Chaturvedi, N., Meghal, S.A., & Jasraj, Y. (2013). Is there impact of climate change on biological productivity in the Indian Ocean? Indian Journal of Geo-Marine Sciences, 42, 50-57. [ Links ]

Cloern, J.E., Foster, S.Q., & Kleckner, A.E. (2014). Phytoplankton primary production in the world’s estuarine-coastal ecosystems. Biogeosciences, 11, 2477-2501. DOI: 10.5194/bg-11-2477-2014 [ Links ]

Dallas, H. (2008). Water temperature and riverine ecosystems: An overview of knowledge and approaches for assessing biotic responses, with special reference to South Africa. Water SA, 34, 393-404. [ Links ]

Das Sarkar, S., Naskar, M., Gogoi, P., Raman, R.K., Manna, R.K., Samanta, S. … Das, B.K. (2019a). Impact assessment of barge trafficking on phytoplankton abundance and Chl a concentration, in River Ganga, India. Plos One, 14(9), p.e0221451. [ Links ]

Das Sarkar, S., Sahoo, A.K., Gogoi, P., Raman, R.K., Munivenkatappa, M.H., Kumari, K., … Das, B.K. (2019b). Phytoplankton biomass in relation to flow dynamics: the case of a tropical river Mahanadi, India. Tropical Ecology, 60(4), 485-494. DOI: 10.1007/s42965-019-00048-7 [ Links ]

Das Sarkar, S., Sarkar, U.K., Lianthuamluaia, L., Ghosh, B.D., Roy, K., Mishal, P., & Das, B.K. (2020). Pattern of the state of eutrophication in the floodplain wetlands of eastern India in context of climate change: a comparative evaluation of 27 wetlands. Environmental Monitoring and Assessment, 192(3), 1-12. [ Links ]

DeNicola, D. (1996). Periphyton responses to temperature at different ecological levels. In R.J. Stevenson, M.L. Bothwell, & R.L. Lowe (Eds.), Algal Ecology in Freshwater Benthic Ecosystems (pp. 57-68). New York, USA: Academic Press. [ Links ]

Elliot, J.A., Thackeray, S.J., Huntingford, C., & Jones, R.G. (2005). Combining a regional climate model with a phytoplankton community model to predict future changes in phytoplankton in lakes. Freshwater Biology, 50, 1404-1411. [ Links ]

Field, C.B., Barros, V.R., Mach, K.J., & Mastrandrea, M.D. (2014). Technical summary. In C.B. Field, V.R. Barros, D.J. Dokken, K.J. Mach, & M.D. Mastrandrea (Eds.), Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part A: Global and Sectoral Aspects (pp. 35-39). Cambridge, United Kingdom: Cambridge University Press. [ Links ]

Friberg, N., Christensen, J.B., Olafsson, J.S., Gislason, G.M., Larsen, S.E., & Lauridsen, T.L. (2009). Relationship between structure and function in streams contrasting in temperature: possible impacts of climate change on running water ecosystems. Freshwater Biology, 54, 2051-2206. DOI: 10.1111/j.1365-2427.2009.02234.x [ Links ]

Ghosal, S., Rogers, M., & Wray, A. (2000). The turbulent life of phytoplankton. Center for Turbulence Research, Proceedings of the Summer Program. Retrieved from https://www.researchgate.net/profile/Sandip_Ghosal/publication/23831690_The_Turbulent_Life_of_Phytoplankton/links/00b4951fa3e956d1b8000000/The-Turbulent-Life-of-Phytoplankton.pdfLinks ]

Gregor, J., & Marsalek, B. (2004). Freshwater phytoplankton quantification by chlorophyll a: a comparative study of in vitro, in vivo and in situ methods. Water Research, 38, 517-522. DOI: 10.1016/j.watres.2003.10.033 [ Links ]

Hartmann, D.L., Klein-Tank, A.M.G., Rusticucci, M., Alexander, L.V., Brönnimann, S., & Charabi, Y. (2013). Observations: atmosphere and surface. In T.F. Stocker, D. Qin, G.K. Plattner, M. Tignor, S.K. Allen, & J. Boschung (Eds.), Climate change 2013: the physical science basis (pp. 159-254). Cambridge, United Kingdom: Cambridge University Press. [ Links ]

Jeong, K.S., Kim, D.K., Whigham, P., & Joo, G.J. (2003). Modelling Microcystis aeruginosa bloom dynamics in the Nakdong River by means of evolutionary computation and statistical approach. Ecological Modelling, 161, 67-78. [ Links ]

Ji, Z.G. (2017). Water quality and eutrophication. In Z.G. Ji (Ed.), Hydrodynamics and Water Quality: Modeling Rivers, Lakes, and Estuaries (2nd Ed., pp. 161-270). New Jersey, USA: John Wiley and Sons Inc. [ Links ]

Karnatak, G., Sarkar, U.K., Naskar, M., Roy, K., Gupta, S., Nandy, S.K., … Verma, V.K. (2018). Understanding role of climatic and environmental parameters in gonadal maturation and spawning periodicity of spotted snakehead, Channa punctata (Bloch, 1793) in a tropical floodplain wetland, India. Environmental Biology of Fishes, 101, 595-607. DOI: 10.1007/s10641-018-0722-6 [ Links ]

Klapper, H. (1991). Control of eutrophication in Inland waters. Sussex, United Kingdom: Ellis Horwood Ltd. West. [ Links ]

Komatsu, E., Fukushima, T., & Shiraishi, H. (2006). Modeling of P-dynamics and algal growth in a stratified reservoir mechanisms of the P-cycle and interactions between water and sediment. Ecological Modelling, 197, 331-349. [ Links ]

Lewandowska, A., & Sommer, U. (2010). Climate change and the spring bloom: a mesocosm study on the influence of light and temperature on phytoplankton and mesozooplankton. Marine Ecology Progress Series, 405, 101-111. [ Links ]

Lindeman, R.H., Merenda, P.F., & Gold, R.Z. (1980). Introduction to Bivariate and Multivariate Analysis. Illinois, USA: Glenview. [ Links ]

Manna, R.K., Satpathy, B.B., Roshith, C.M., Naskar, M., Bhaumik, U., & Sharma, A.P. (2013). Spatio-temporal changes of hydro-chemical parameters in the estuarine part of the river Ganges under altered hydrological regime and its impact on biotic communities. Aquatic Ecosystem Health & Management, 16(4), 433-444. [ Links ]

Milly, P.C.D., Dunne, K.A., & Vecchia, A.V. (2006). Global pattern of trends in stream flow and water availability in a changing climate. Nature, 438, 347-350. DOI: 10.1038/nature04312 [ Links ]

Mitra, A., Zaman, S., & Raha, A.K. (2014). Phytoplankton cell volume and diversity in Indian Sundarbans. Indian Journal of Geo-Marine Sciences, 43, 1914-1921. [ Links ]

Mooij, W.M., Janse, J.H., Domis, L.N., Hulsmann, S., & Ibelings, B.W. (2007). Predicting the effect of climate change on temperate shallow lakes with the ecosystem model PC Lake. Hydrobiologia, 584, 443-454. [ Links ]

Naskar, M., Roy, K., Karnatak, G., Nandy, S.K., & Roy, A. (2017). Quantifying climate change induced threats to wetland fisheries: a stakeholder-driven approach. Environment Development and Sustainability, 20, 2811-2830. DOI: 10.1007/s10668-017-0018-6 [ Links ]

Nweze, N.O., & Ude, B.O. (2013). Algae and physico-chemical characteristics of Adani rice field, Enugu State, Nigeria. IOSR Journal of Pharmacy and Biological Sciences, 8(5), 12-18. [ Links ]

Paul, R.K., & Birthal, P.S. (2015). Investigating rainfall trend over India using the wavelet technique. Journal of Water and Climate Change, 7, 353-364. DOI: 10.2166/wcc.2015.079 [ Links ]

Porter, J., Arzberger, P., Braun, H.W., Bryant, P., Gage, S., Hansen, T., … Williams, T. (2005). Wireless sensor networks for ecology. Bioscience, 55, 561-572. DOI: 10.1641/0006-3568(2005)055[0561:WSNFE]2.0.CO;2 [ Links ]

R Core Team. (2015). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Retrieved from https://www.R-project.org/Links ]

Rathore, L.S., Attri, S.D., & Jaswal, A.K. (2013). State level climate change trends in India. New Delhi, India: Indian Meteorological Department (IMD). [ Links ]

Roger, P.A., & Kulasooriya, S.A. (1980). Blue-Green Algae and Rice. Laguna, Philippines: International Rice Research Institute. [ Links ]

Roshith, C.M., Meena, D.K., Manna, R.K., Sahoo, A.K., Swain, H.S., Raman, R.K., Sengupta, A., & Das, B.K. (2018). Phytoplankton community structure of the Gangetic (Hooghly-Matla) estuary: Status and ecological implications in relation to eco-climatic variability. Flora, 240, 133-143. [ Links ]

Salmaso, N., & Braioni, M.G. (2008). Factors controlling the seasonal development and distribution of the phytoplankton community in the lowland course of a large river in Northern Italy (River Adige). Aquatic Ecology, 42, 533-545. DOI: 10.1007/s10452-007-9135-x. [ Links ]

Sarkar, U.K., Naskar, M., Srivastava, P.K., Roy, K., Sarkar, S.D., Gupta, S., … Karnatak, G. (2019). Climato-environmental influence on breeding phenology of native catfishes in River Ganga and modeling species response to climatic variability for their conservation. International Journal of Biometeorology, 63(8), 991-1004. [ Links ]

Schabhutt, S., Hingsamer, P., Weigelhofer, G., Hein, T., Weigert, A., & Striebel, M. (2013). Temperature and species richness effects in phytoplankton communities. Oecologia, 171, 527-536. DOI: 10.1007/s00442-012-2419-4 [ Links ]

Sharma, A.P., Joshi, K.D., Naskar, M., & Das, M.K. (2015). Inland Fisheries and Climate Change: Vulnerability and Adaptation options. Kolkata, India: Central Inland Fisheries Research Institute. [ Links ]

Sipkay, C.S., Kiss, K.T., Vadadi-Fulop, C.S., & Hufnagel, L. (2009). Trends in research on the possible effects of climate change concerning aquatic ecosystems with special emphasis on the modeling approach. Applied Ecology and Environmental Research, 7, 171-198. [ Links ]

Sipkay, C., Kiss, K.T., Vadadi-Fulop, C.S., Homoródi, R., & Hufnagel, L. (2012). Simulation modeling of phytoplankton dynamics in a large eutrophic river, Hungary-Danubian Phytoplankton Growth Model (DPGM). Biologia, 67, 323-337. DOI: 10.2478/s11756-012-0004-2 [ Links ]

Sommer, U., & Lengfellner, K. (2008). Climate change and the timing, magnitude, and composition of the phytoplankton spring bloom. Global Change Biology, 14, 1199-1208. [ Links ]

Vass, K.K., Das, M.K., Tyagi, R.K., Katiha, P.K., Samanta, S., Shrivastava, N.P., … Gopal, B. (2011). Strategies for Sustainable Fisheries in the Indian Part of the Ganga Brahmaputra River Basins. International Journal of Ecology and Environmental Sciences, 37(4), 157-218. [ Links ]

Verdonschot, P.F.M., Hering, D., Murphy, J., Jahnig, S.C., Rose, N.L., Graf, W., …Sandin, L. (2010). Climate Change and the Hydrology and Morphology of Freshwater Ecosystem. In M. Kernan, R.W. Battarbee, & B. Moss (Eds.), Climate Change Impacts on Freshwater Ecosystems (pp. 1079-1081). New Jersey, USA: Blackwell Publishing. [ Links ]

Webb, B.W., Hannah, D.M., Moore, R.D., Brown, L.E., & Nobilis, F. (2008). Recent advances in stream and river temperature studies. Hydrological Processes, 22, 902-918. DOI: 10.1002/hyp.6994 [ Links ]

Webb, B.W., & Nobilis, F. (2007). Long-term changes in river temperature and the influence of climatic and hydrological factors. Hydrological Sciences Journal, 52, 74-85. [ Links ]

Whitehead, P.G., Battarbee, R.W., Crossman, J., Elliott, J.A., Wilby, R., Monteith, D.T., & Kernan, M. (2012). River and lake water quality: future trends. NERC/Centre for Ecology & Hydrology. Retrieved from http://nora.nerc.ac.uk/17812/Links ]

Willis, A., Chuang, A.W., Orr, P.T., Beardall, J., & Burford, M.A. (2019). Subtropical freshwater phytoplankton show a greater response to increased temperature than to increased pCO2. Harmful Algae, 90, 101705. [ Links ]

Winslow, L.A., Read, J.S., Hansen, G.J.A., Rose, K.C., & Robertson, D.M. (2017). Seasonality of change: Summer warming rates do not fully represent effects of climate change on lake temperatures. Limnological Oceanography, 62, 2168-2178. DOI: 10.1002/lno.10557 [ Links ]

Wood, S.N. (2006). Generalized Additive Models: An Introduction with R. Boca Ratón, USA: Chapman & Hall/CRC. [ Links ]

Wrona, F.J., Prowse, T.D., Reist, J.D., Hobbie, J.E., Levesque, L.M.J., & Vincent, W.F. (2006). Climate Change Effects on Aquatic Biota, Ecosystem Structure and Function. Ambio, 35, 359-369. [ Links ]

Received: July 03, 2020; Revised: October 15, 2020; Accepted: October 19, 2020

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License