4.2. Carbon isotopic signature in different taxa/groups as niche indicators
Carbon isotopic signature is related to habitat and food sources. Food selectivity also can play a role, particularly as phytoplankton size and shape seem to be important in the seasonality of carbon isotopic signature in general [
63]. The relationship between temperature and carbon isotopic signature is indicative also of the vertical distribution of the zooplankton taxa.
Several studies have supported the concept that carbon isotopic signature is affected by water depth [
46,
50,
61,
70]. In deep lakes, the carbon isotopic signature is strongly influenced by depth, with organisms living in deeper layers characterized by more negative values than those living closer to the surface [
69,
72,
73]. During the production year of a deep lake, SPM and zooplankton
13C-signature increases with water warming and the onset of thermal stratification [
28,
50,
63].
In our study we observed differences in the
13C signature between
Daphnia and adult Cyclopoids, likely due to the differences in food sources and habitat preferences. The changes were quantitative and the correlation between δ
13C
Daphnia and the delta fractionation of cyclopoids (F
cyclopoids) was statistically significant (
Figure 4). Cyclopoids increasingly disengage from the water layer in which
Daphnia feeds and lives when warming and vertical thermal stratification increase. However, the question remains whether food differentiation or habitat changes are more important in the resulting δ
13C signal of the two zooplankton taxa during the warmer months characterised by thermal stratification. However, it seems evident that there is a distinct separation of the cyclopoids from the predatory cladocerans, which closely align with
Daphnia’s isotopic carbon signature, likely inhabiting and moving within more superficial and warmer zones of the pelagic. Another study in a deep lake of the same Italian lake district [
74] drew similar conclusions, with cyclopoids likely occupying a different and deeper habitat (>20 meters depth) than
Daphnia. [
75] conducted an experimental study in an oligotrophic lake in Japan, and found that the population of a cyclopoid copepods (
Cyclops sp.) inhabited the hypolimnion during the summer stratification, while was scattered throughout the whole water column in the winter circulation period. Other studies [
76,
77] found spatial segregation during summer stratification, with daphnids resided above the thermocline whereas copepods at the lower edge of the thermocline.
Another explanation of the increasing variations between
Daphnia and the Cyclopoids could result from their different life history strategies and metabolism. Some authors [
78] have suggested that the δ
13C‰ values of adult Cyclops respond more slowly to changes in the isotopic baseline because they have a lower carbon turnover rate compared with
Daphnia that rapidly respond to changes in the δ
13C‰ values of its phytoplankton food source. Such variation in the isotopic turnover between the herbivorous daphnids and predatory copepods, has been confirmed in experimental studies [
72], and is attributed to distinct life history strategies [
79].
The least negative δ
13C signature of the filter feeding cladoceran
Diaphanosoma brachyurum was likely due to its limited seasonal presence during the warmest months of the year (i.e. usually absent in winter and spring samples) and its tendency to occupy superficial water. This habitat preference of
D. brachyurum has been demonstrated in previous studies in Lake Maggiore [
80], with species biomass increasing during the establishment of thermal stratification in the warmer, upper ten meters of the lake.
Littoral carbon sources have less
13C-depleted signatures [
50,
81,
82,
83,
84]. Seasonality has been demonstrated also in the littoral carbon signature, with the same pattern observed for the pelagic carbon signature, i.e. a transition from more (in winter) to less (in summer) negative δ
13C‰ values. Previous studies on Lake Maggiore [
60,
63,
83] demonstrated that ranges of isotopic signatures measured in the same season were quite distinct for littoral and pelagic habitat, allowing a correct interpretation of data when speculating on migration of zooplankters between different lake zones.
The two particle feeder cladocerans
Daphnia and
Bosmina shared the same δ
13C‰ range and similar values, indicating they share of the same food sources and have similar mobility across lake habitats. The two taxa have, however, distinct feeding modes:
Daphnia is a strict filter feeder, processing food particles through the filtering combs of the third and fourth thoracic limbs [
85]. Whereas the foraging mode of
Bosmina deviates markedly from filter feeding in a way that permits efficiency in the removal of highly edible flagellates, particularly when these items occur at low density [
86]. In our twelve years’ study, on some occasions
Bosmina’s δ
13C-signature was less negative than
Daphnia’s, likely reflecting a mobility through the water column of this small cladoceran and a tendency to feed in more superficial layers than
Daphnia. The relatively low δ
15N-signature of
Bosmina with values approaching zero during the summer in few years, confirms the capability of this zooplankter to exploit nitrogen fixing cyanobacteria [
74,
83,
86,
87]. Cyanobacteria are usually lower in nitrogen isotopic signature than other phytoplankton groups, in particular nitrogen-fixing taxa [
88,
89,
90], found in the phytoplankton community of Lake Maggiore [
63]. The feeding mode of
Bosmina avoids clogging during the growth of filamentous phytoplankton colonies, including those of cyanobacteria.
The range of variation of carbon isotopic signature (Figure 6a) was small for Cyclopoids and Diaptomids, indicating consistency of food sources throughout the studied period and that Cyclopoids and Diaptomids tend to be conservative in their habitat during the year. For Diaptomids, however, a more skewed distribution was observed toward less negative/depleted δ
13C‰ values, suggesting an opportunistic feeding strategy with a discrete mobility across different zones or habitats of the lake. In particular, because less depleted values of δ
13C‰ occurred during the summer periods, the migration of Diaptomids towards more superficial waters was considered [
50].
Within the zooplankton community of Lake Maggiore, a narrower range of variation for nitrogen isotopic signatures (
Figure A3) was found for Cyclopoids; which combined with their small range of carbon isotopic signatures suggests an overall niche specificity of this zooplankton group. Enrichment with respect to the pelagic baseline, however, cannot be calculated, being solidarity with carbon isotopic signature restricted to the unstratified thermal phase.
Figure 6.
Comparison of seasonally decomposed trends of chlorophyll-a (a) and nitrate (N-NO3) (b) with carbon isotopic signatures of selected zooplankton taxa ((c): red line= Bosmina; green line= diaptomids; purple line= Daphnia), during the period 2009-2020.arbon isotopic values vs. cyclopoids’ carbon fractionation (F).
Figure 6.
Comparison of seasonally decomposed trends of chlorophyll-a (a) and nitrate (N-NO3) (b) with carbon isotopic signatures of selected zooplankton taxa ((c): red line= Bosmina; green line= diaptomids; purple line= Daphnia), during the period 2009-2020.arbon isotopic values vs. cyclopoids’ carbon fractionation (F).
The conceptual diagram represents a summary of all the above-described variations in food sources, habitat preferences and seasonal variations of the zooplankton taxa considered in our study in Lake Maggiore. The diagram clearly illustrates how the relative δ
13C‰ values of zooplankton taxa are collocated with respect to each other and to vertical lake depth in the different seasons. In particular, it illustrates the shift of δ
13C‰ of different zooplankton taxa during two very different times of the year: late winter/springtime when water column circulation generally occurs in the lake, and summertime when stratification takes place (
Figure A4).
In large deep lakes at median latitude, during late winter significant water mixing typically occurs, caused by a change in the air temperature and helped by windy conditions, thus increasing homogeneity in water conditions and environmental parameters along the water column (temperature, dissolved oxygen, pH, etc). Phytoplankton is consequently distributed in the water column and in Lake Maggiore the phytoplankton community is mainly represented by diatoms which require a high degree of turbulence to maintain in suspension [
23] and having
13C-depleted values [
46,
91]. At this time of the year, we registered values of δ
13C‰ for zooplankton taxa within in a narrow range, likely as a result of the homogeneity of phytoplankton type and of habitat.
With the advancing of warmer temperatures, water conditions change, establishing thermal stratification and creating a vertical gradient of the main environmental parameters. These conditions create heterogeneity in the lake water column, reflected also in more diverse plankton communities and vertical distribution of groups and biomass. Filamentous green algae, smaller diatoms, cyanobacteria, and other species capable of depth regulation are favored with intensified stratification [
92]. In Lake Maggiore during the summer, cyanobacteria, green algae and diatoms are all relevant groups contributing to the phytoplankton community. A detailed study in Lake Maggiore [
63] revealed a vertical gradient of phytoplankton density and biomass, in three different lake stations, during thermal stratification; recent phytoplankton monitoring with a fluorescent probe (FluoroProbe, BBE Moldaenke, GmbH, Germany) and with high frequency chlorophyll-
a in situ sensors (Turner Design, USA) confirmed different chlorophyll-
a concentrations and phytoplankton groups with water depth [personal comment, data unpublished). In this decadal study, δ
13C‰ values for different zooplankton taxa in summer consistently underwent a clear separation from those of springtime when taxa δ
13C‰ values were similar.
We hypnotize that it is likely because the conditions of thermal stratification during the summer create vertical heterogeneity in both food sources and habitats, allowing partial niche partitions for the different zooplankton taxa.
Although lakes are often considered as homogenous environments/habitats, this concept may be valid for the circulation period but not when vertical thermal stratification is established. Variation in vertical dimensions in temperature, light, dissolved oxygen and nutrients lead to vertical habitat diversification and creates various microhabitats, affecting the vertical distribution of zooplankton taxa with different food and temperature requirements and predation pressure [
76,
93,
94,
95,
96,
97,
98,
99] In our study (referring to our conceptual diagram) the stratification period corresponded to a phase of heterogeneity also for zooplankton taxa vertical distribution, as suggested by the distribution of isotopic signatures during the twelve year study period.
There is growing consensus that lakes, far from being homogeneous environments, offer a rich and dynamic heterogeneity at multiple spatial and temporal scales [
100]. Natural and anthropogenic factors create spatial heterogeneity in lakes, leading to heterogeneity in environmental conditions and resource availability within-lake compartmentalization in food web structure [
101].
4.3. Zooplankton (Daphnia) δ13C and chlorophyll-a: indication of food quality and availability
In our decadal study of Lake Maggiore, chlorophyll-
a concentration was related to zooplankton δ
13C signatures, indicating the importance of phytoplankton standing crop as food availability and utilization for zooplankton. Previous studies on Lake Maggiore reported significant correlations between chlorophyll-
a and sestonic carbon concentration over the long term [
102], suggesting that POM in the lake is mainly of autochthonous origin and likely made up of algal cells, with the seston δ
13C-signature mainly related to phytoplankton. The correlation found in this study between zooplankton δ
13C-signatures and chlorophyll-
a concentration suggests that the phytoplankton biomass is important in determining zooplankton carbon signature, likely representing the main food source for zooplankton in Lake Maggiore. In other freshwater studies [
62] it has been shown that zooplankton tend to utilize mainly algal derived autochthonous carbon during the more productive months. In addition, [
65] found a correlation between the herbivorous zooplankter
Daphnia and chlorophyll-
a in their study of four annual production cycles in a subarctic lake.
Because chlorophyll-
a variability is often a function of primary productivity in freshwaters [
103,
104], the δ
13C seasonality can be interpreted as resulting from changing phytoplankton growth rates [
65]. This finding is in good agreement with other field investigations [
49,
105,
106,
107] and culture experiments [
106,
108,
109] which reveal a strong relationship between growth rate and δ
13C‰ values of aquatic primary producers. Several models predict a rise in δ
13C‰ values with increasing algal growth rate or carbon demand [
31,
108,
109,
110].
Plankton collected during periods of high biomass and primary productivity tends to exhibit the heaviest (i.e., less negative) δ
13C-signatures, due to reduced isotopic fractionation at high cell densities or growth rates, or a switch to utilizing
as CO
2 concentration is depleted at high cell densities or high growth rates [
111,
112,
113]. The relationships are also mediated by the difference in δ
13C-signatures of the phytoplankton groups present in a certain season, since during the summer months Cyanobacteria tend to increase in Lake Maggiore [
63] and this phytoplankton group is shown to possess relatively enriched δ
13C values [
114].
4.4. δ13C and δ15N time series and environmental variables
Decomposition of time series data allowed the detection of seasonality and trends of environmental parameters during the considered period (2010-2020). In particular, the seasonally decomposed trend of chlorophyll-a concentration increased, while the nitrate (N-NO3) trend decreased.
Recent long-term studies on Lake Maggiore [
5] have shown that during the oligotrophication phase of the lake between 1988 and 2005, chlorophyll-
a, and phytoplankton biovolume, significantly decreased down to low concentrations of 2.0-3.0 µg l
-1. However, in the last decade despite nutrient concentrations remaining quite stable and low, chlorophyll-
a concentration slightly increased, with fluctuations observed especially during summer periods. These data were recently confirmed by satellite observations, detecting since 2015 a change in chlorophyll-
a seasonal pattern with less defined spring blooms and a tendency for annual maxima to occur in summer [
115]. In this study, we detected those increases and fluctuations in chlorophyll-
a especially between 2015 and 2020, matching an increasing trend of δ
13C‰ values of zooplankton taxa such as
Bosmina and Diaptomids. Explanation might be found in the new water conditions with increasing algal growth rates or carbon demand, leading to an increase in the δ
13C-signature of phytoplankton as predicted by several models [
31,
108,
109,
110].
In our study, time series decomposition analysis also revealed that δ
15N-signature of the cladocerans
Daphnia and
Bosmina had a decreasing trend during the last decade. A number of studies have demonstrated that the plankton δ
15N-fingerprint is a good indicator of nitrogen sources, nitrogen cycling processes and consumer trophic position in aquatic systems [
116]. Changes in phytoplankton species composition may influence the use of nitrogen sources (N
2 fixation vs. uptake of dissolved combined nitrogen). Under nitrogen limitation, cyanobacterial N
2 fixation is stimulated, resulting in depleted δ
15N‰ values, which may then be transferred to zooplankton through feeding, resulting in a low zooplankton δ
15N-signature [
116]. Although daphnids cannot ingest the largest forms of cyanobacteria, they are capable of feeding on the smaller, decomposing pieces [
118] and on the associated bacteria, so they will still reflect δ
13C-signature of the baseline during cyanobacterial blooms [
78]. The capability of
Bosmina to exploit cyanobacteria has been already mentioned in this paper. Overall, the decreasing trend in δ
15N-signature of both
Daphnia and
Bosmina during the last decade in our results match the conditions of nitrogen limitation and decreasing concentrations in nitrate and total nitrogen in the epilimnion of Lake Maggiore, as revealed by recent long-term studies [
119].
Our twelve-year study of carbon and nitrogen stable isotopes revealed that zooplankton isotope signatures incorporate changes in lake environmental parameters, integrating the trends that have occurred during the last decade, thereby representing candidate sentinels capable of integrating the influence of climate change on the ecosystem of the lake.
Although the nutrient status in Lake Maggiore has reached oligotrophy, climatic change has affected lake thermal dynamics by increasing the water column stability, aggrading resistance of the water column to mixing and diminishing the water mixing depth [
5], all factors which affect the primary production of the lake. The increase of fluctuations and of concentrations in chlorophyll-
a experienced in Lake Maggiore during most recent times, could be the result of the recent changes in the lake which has modified conditions altering lake primary production.
Especially from work on shallow lakes, limnologists are realising that climate change is intensifying the symptoms of eutrophication in freshwaters [
1,
13], and results from direct experiments suggests that warming can increase the symptoms of eutrophication and changes in ecosystems structure. Our results, together with the most recent long-term studies on Lake Maggiore, could be interpreted as signs of eutrophication-like effects on a deep lake: because increased temperatures and water column stability promoted a prolonged period of water stratification and of epilimnion segregation during the warmest months, these new conditions could resemble those in the better studied shallow lakes, with similar final consequences in promoting lake primary production.
Figure 6 compares the seasonally decomposed trend of chlorophyll-a and nitrate with the carbon isotopic signature of selected zooplankton taxa. A clear-cut match is evident between a surge in chlorophyll-a concentration, the decrease in epilimnion nitrates since 2015 and the trend change in carbon isotopic signature of the considered zooplankters.
These conditions suggest that climate change represented the main driving factor of the lake dynamics during the last decade and that its effects have been integrated in the stable isotopes signatures of the zooplankton community, candidates to become sentinels of climatic changes in the lake.
Temperature as a stressor can impact freshwater ecosystems at lower trophic levels (e.g. primary producers) and can cause a cascading effect on other species through the food chain, which in turn affects ecosystem structure and function [
120,
121]. For aquatic ecosystems, the phytoplankton-zooplankton interaction determines the structure of food web, which represents the energy and materials flowing through the population in the system [
122,
123,
124].