Human-Induced Sharp Salinity Changes in the World’s Largest Hypersaline Lagoon Bay Sivash (Crimea) and Their Effects on the Ecosystem

: Lakes and lagoons play an important role worldwide, and salinity ﬂuctuations signiﬁcantly affect their ecosystems. Bay Sivash, the world’s largest hypersaline water body, underwent a sharp change in salinity, induced by the closing of the North Crimean Canal. To monitor a shift in the ecosystem, a study was carried out from 2014 to 2020 at 15 sites of the lagoon. Since the closure of the canal, the average salinity increased from 22 g L − 1 (2013) to 94 g L − 1 (2020). Suspended solids and dissolved organic matter also increased. When salinity increased above 50 g L − 1 , the number of taxa signiﬁcantly decreased; this was a negative linear relation. The increase in salinity signiﬁcantly changed the structure of zooplankton and benthos. The most dramatic changes occurred with the salinity increase from 25 to 70 g L − 1 . Chironomidae larvae numbers began to increase greatly in the ecosystem of the bay, and since 2014, they have rapidly increased their contribution to the abundance of benthos and plankton. The concentration of benthic–planktonic species increased in plankton, in particular, in Harpacticoida and Chironomidae. At salinity above 80–90 g L − 1 , nauplii and adult brine shrimp appeared to become abundant in plankton and benthos. The transit of the ecosystem to a new alternative state occurred.


Introduction
Among ongoing global changes, salinity modification in different water bodies is one of the most influential, affecting the structure, functioning, and dynamics of aquatic ecosystems [1][2][3][4]. This modification may have various causes (natural and anthropogenic) as well as different directions, either increasing or decreasing [5][6][7][8]. The same global change can lead to an increase in salinity in some water bodies and a decrease in others. For example, air temperatures increase, leading to an acceleration of melting of mountain glaciers, which causes a salinity decrease in lakes that are fed from these glaciers [6]; however, in other lakes, it may cause a salinity increase by increasing evaporation [2]. Anthropogenic influences also may strongly influence the salinity of water bodies, with a wide variety of ecosystems responding to them [9][10][11][12][13].
Each water body has individual abiotic and biological characteristics and differing external natural and anthropogenic factors, which result in each water body being unique. Understanding the effect of changes in salinity on an aquatic ecosystem requires a combination of knowledge of the uniqueness of the particular ecosystem, overlaid with known general patterns. To this end, data from salinity changes in a variety of water bodies should be accumulated to obtain general patterns. While there are some accumulated data on the response to fast or slow salinity changes on various aquatic systems, there is not sufficient collected data to use as a tool to understand, anticipate, and mitigate those Before the construction of the North Crimean Canal (1963)(1964)(1965)(1966)(1967)(1968)(1969)(1970)(1971)(1972)(1973)(1974)(1975), the semi-closed shallow lagoon was hypersaline (mean salinity of about 140 g L −1 reaching >200 g L −1 in the southern part) and highly productive. Due to scarcity of fresh water in Crimea (the largest peninsula in the Black Sea), the North Crimean Canal was constructed (total length of 465 km, 375 km in Crimea) to supply water from the River Dnieper into Crimea. This started the development of irrigated agriculture on 4000 km 2 of land, and drainage water from that area began discharging mostly into Bay Sivash. In 1985, 521 million m 3 and 109 million m 3 of drainage waters were discharged from the territories of Crimea and the Kherson region, respectively [12]. As a result, after the start of the operation of the North Crimean Canal, a decrease in salinity began in the Sivash, which continued in various areas until 1997; the average salinity decreased from 141.00 g L −1 (data from 1955) to 22.60 g L −1 in 1989, and to a minimum value of 17.00 g L −1 in 1997-2009 [15]. The Sivash ecosystem gradually transformed into a new ecosystem of brackish waters [12,[16][17][18]. The transformation of the ecosystem was accompanied by the appearance of large populations of commercial fish in the bay, that came from the Sea of Azov, due to which, fishing began to intensively develop in the bay. Due to political reasons, the Ukrainian government decided to stop the supply of Dnieper water via the North Crimean Canal, and did so in April 2014. The discharge of freshwater to the lagoon almost stopped, and salinity started to increase gradually, resulting in a new ecosystem shift. To monitor this shift, a comprehensive study was carried out from 2014 until 2020. Partial results which were obtained before June 2018, and were analyzed and published previously [4,12,[18][19][20].
The proposed goals of this work are the following: (1) to present new unpublished data obtained from 2016 (June) to 2020 (December), summarize, and analyze all these data; (2) find the general patterns in the different occurred changes (means and ranges) and the correlation between them; (3) prove or disprove the assumption that salinity itself can be a good predictor of ecosystem shift character.

Materials and Methods
Bay Sivash (from 45 • 17 to 46 • 12 N, from 34 • 30 to 35 • 28 E) is separated from the Sea of Azov by the sandy narrow Arabat Spit (from 270 m to 8 km in width, length 112 km) and occupies an area of around 2560 km 2 ( Figure 1). To the north, the narrow Henichesk Strait connects the lagoon with the sea; it has been described previously in several papers [4,12,18,20]. It is a shallow water body with a depth of 0.1 to 2.5 m. Few small rivers enter the lagoon, and no river drains the water from it. The water balance in the lagoon before 2014 was mainly determined by the water from the North Crimean canal, precipitation, exchange with the Sea of Azov, and evaporation. Flowing freshwater from the North Crimean canal was 10 times more than that from all the small rivers. The volume of water from the small rivers was approximately 15 times lower than average annual precipitation. Before the canal construction, floating mats of the filamentous green alga Cladophora siwaschensis C. Meyer occupied huge areas in the bay (Figure 2c). Currently, floating mats have formed again, occupying a large area and reaching high biomass (up to 2.5-3.0 kg m −2 in wet weight). In these mats, there is an abundance of animals of different taxa. Hypoxic and even anoxic conditions often exist below these floating mats. After construction of the canal, changes in Bay Sivash were accompanied by shifts on its shores, where large areas were occupied by belts of the common reed, which, in some places, reached 100 meters wide. With the salinity increase, the common reed began to die, and its belt width decreased significantly ( Figure 2d). Data used in this paper were collected during eight expeditions to Sivash in the period from June 2016 to December 2020; sampling sites and their coordinates are given in Figure 1 and Table 1. The methods used were the same as in previous works [4,12,18,20]. The soft-bottom samples were collected by benthic tubes (diameter of 5 cm) in duplicate at depths of 0.2-0.6 m. Quantitative zooplankton samples were taken by filtration of 100-150 L of water through a plankton net with a mesh size of 110 µm. All samples were fixed by 4% formaldehyde. Samples were processed in a Bogorov chamber with binocular MBS-9, Olympus SZ61TR. The concentration of total suspended solids (TSS), as well as dissolved organic matter (DOM), was measured with the hydrobiophysical multiparametric complex 'Condor' (made by Akvastandart-Yug, Sevastopol, Russia, http://ecodevice.com.ru/ecodevice-catalogue/ multiturbidimeter-kondor, accessed date 1 December 2021). Salinity and temperature were also determined by manual refractometer Kellong WZ212 and an electronic thermometer PHH-830, respectively.
All data were subject to standard statistical treatments: mean values, the regression equation parameters, coefficients of variation (CV), correlation (R), and determination (R 2 ). These were calculated in the standard programs MS Excel 2007 and STATISTICA 6.0. The confidence level of correlation coefficients (p) was also evaluated [21].

Abiotic Parameters
Average values and indices of their variability of the studied parameters in different periods of observation are given in Table 2.
The average salinity and the range of its spatial differences constantly increased during the study. At the northernmost point close to the strait connecting the lagoon with the sea, salinity did not change significantly, for example, in December 2020, it ranged from 30 to 41 g L −1 . At the southernmost point, salinity varied significantly from 82 to 112 g L −1 during the same period. The direction and speed of wind could modify salinity distribution, given the general tendency for salinity to increase from north to south. The TSS concentration showed a large spatial heterogeneity, which was even higher during the period of strong autumn-winter winds. In the summer months, a significant correlation (p < 0.01) was noted between TSS and salinity, for example, in May 2018, the relationship had the following form (R = 0.913, p = 0.0005): where: TSS-total suspended solids, mg L −1 ; S-salinity, g L −1 . A similar significant dependence was noted in other summer periods, for example in June 2019 ( Figure 3). In the autumn-winter period, this dependence was not so pronounced, and the main influence on the distribution of suspended matter was exerted by winds. As an example, there was a northeast wind during sampling in November 2018 influencing TSS spatial distribution in the bay. TSS and DOM concentrations were significantly (p = 0.001) lower near the windward shore than near the leeward shore. TSS concentrations were on average 22.2 mg L −1 (CV = 1.085) near the windward shore and only 5.7 mg L −1 (CV = 0.660) opposite, and for DOM were 13.95 mg L −1 (CV = 1.011) and 4.66 mg L −1 (CV = 0.324), respectively. The heterogeneity of the spatial distribution of DOM was less pronounced than the TSS; however, significant positive dependence on salinity was also noted, which was seen, for example, in June 2019 ( Figure 3). There was a close positive correlation in the spatial distribution between DOM and TSS in all observation periods. For example, in June 2019, it looked like the following (R = 0.922, p = 0.0005): where DOM and TSS are measured in mg L −1 . The data showed the existence of a general trend of increasing DOM concentration from spring to autumn, and a slight increase in DOM from 2018 to 2020.

Zooplankton
Generalized data on changes in the zooplankton structure are given in Table 3. During the period under consideration, there were significant changes in the taxonomic structure of zooplankton. In June 2016, Copepoda (Acartia (Acartiura) margalefi Alcaraz, 1976) still dominated in numbers, accounting for 86% of the abundance of zooplankton.
In June and November 2018, copepods accounted for only 7% and 2% of the total zooplankton abundance, respectively. In later samples, their contribution to the total number did not even reach 1%. In the period from 2018 to 2020, Artemia sp. (Anostraca) or Harpacticoida dominated. At the same time, during this period, the macro-taxonomic diversity of plankton decreased. At a salinity of up to 90 g L −1 , representatives of 7 classes of animals were present in plankton; from 70 to 112 g L −1 , 6 classes were represented, and at a salinity of 138 g L −1 , only 3 classes. In the period from June 2016 to September 2020, the total abundance of zooplankton increased, which can be approximated as follows (R = 0.862, p = 0.005): where: N t -where the average total number of zooplankton at the time of sampling, ind. m −3 ; t-the number of months since the beginning of the study (June 2016 was taken as first month).  Between September and December 2020, the total number of zooplankton decreased by 4 times, which was probably due to winter conditions. If the mid-winter point (December 2020) is excluded, then a significant relationship is revealed between the total abundance of zooplankton and salinity (R = 0.893, p = 0.005), as follows: where: N t -where the average total number of zooplankton at the time of sampling, ind. m −3 ; S-salinity, g L −1 .
In each of the study periods, in the spatial aspect, the total abundance of zooplankton, and its constituent groups did not unambiguously depend on salinity, TSS, and DOM. However, in all periods of the study, the maximum total abundance was found in the salinity range from 80 to 110 g L −1 . This was because the number of active stages of Artemia sharply increased at salinity above 75-80 g L −1 . For example, in July 2020, in the salinity range from 41 to 86 g L −1 , the average number of active stages of Artemia was 238 ind. m −3 (SD = 384, CV = 1.612), and in the salinity range from 90 to 128 g L −1 was 75,271 ind. m −3 (SD = 82,637, CV = 1.098). The differences are significant (p = 0.001). Winds, especially in the autumn-winter period, had a significant impact on the distribution of zooplankton, this effect was not the same for its components. In November 2018, as an example, average Artemia cyst abundance was significantly higher near the leeward shore (by 73 times) than near the windward shore (p = 0.0001). Artemia nauplii on average were also significantly (p = 0.05) more abundant (by 2.8 times) near the leeward shore than near the opposite shore. Another trend was found for adult Artemia, when more females and males were observed near the windward shore (average 718 and 1700 ind. m −3 , respectively) than near the leeward shore (553 and 493 ind. m −3 , respectively). Cyst abundance did not correlate with salinity near the windward shore, and there was a positive significant correlation near the leeward shore (R = 0.982; p = 0.001), as follows: where: N c -the number of cysts, S-salinity, g L −1 .

Zoobenthos
The variability in the composition of macrotaxa in zoobenthos is presented in Table 4. Macro-taxonomic diversity in the benthos from May to November 2018 increased from 11 to 16 taxa, and then gradually decreased to 10 (December 2020).
There is a significant correlation between the total number of animal and protist macrotaxa and salinity in the period from November 2018 to December 2020 (R = 0.982, p = 0.005): where: K b -the total number of animal and protist macrotaxa; S-salinity, g L −1 . The total abundance of zoobenthos increased in the period from May 2018 to June 2019 and then began to gradually decrease. The temporal dependence of the average abundance on salinity was bell-shaped rather than linear. At the same time, the maximum abundance was observed in the salinity range from 70 to 80 g L −1 . Table 4 also shows the maximum salinities at which different taxa were encountered. The total abundance of zoobenthos is negatively statistically significantly related to salinity and TSS, for example in June 2019 ( Figure 4). The abundance of individual groups of animals did not correlate with salinity. Changes in different groups did not correlate with each other. Of macrobenthos organisms at salinity above 40-60 g L −1 , only larvae of Chironomidae (Baeotendipes noctivagus (Kieffer, 1911)) were usually present, the share of which, in the total abundance, in some samples exceeded 80%, averaging in the lagoon from 7.5 to 29% for different observation periods. A linear relationship between the number of larvae and salinity was not observed, but the maximum proportion of larvae in the total number of zoobenthos usually was within the salinity range from 80 to 120 g L −1 . Taking into account the fact that chironomid pupae were massively present in plankton only in May, this was the time when the adults emerge. At the beginning of summer, the number of chironomid larvae in the benthos is significantly higher than in November-December. Occasionally, a high abundance of other macrobenthos organisms Gastropoda (Hydrobia acuta (Draparnaud, 1805)), Amphipoda (Gammarus aequicauda (Martynov, 1931)), and Polychaeta (Protodrilius sp., Polydora cornuta Bosc, 1802,) were observed in single sites, and Oligochaeta were found only once in high numbers (Table 4). During the days with a strong northeasterly wind (November 2018), zoobenthos average abundance near the leeward and windward coasts were 106,345 ind. m −2 (CV = 0.956) and 167,979 ind. m −2 (CV = 0.524), respectively. Abundance was higher near the windward shore at 1.7 times. A similar trend was observed also separately for Chironomidae and Harpacticoida. In Gastropoda, in November 2018, the average abundance was 1010 ind. m −2 (CV = 1.126) and 2706 ind. m −2 (CV = 1.545) near windward and leeward shores, respectively.

Zooplankton vs. Zoobenthos
As the comparison of Tables 3 and 4 showed, the most common species of different taxa cannot be classified as purely planktonic or benthic species, since they made a significant contribution to both zooplankton and zoobenthos. The trend where the total number of zooplankton and zoobenthos organisms decreased with increasing salinity was not statistically significant. The ratio of the total abundance of zooplankton and zoobenthos (abundance of zoobenthos/abundance of zooplankton) significantly changed in the bay with an increase in the proportion of zooplankton in the total abundance; the changes that took place can be described (R = 0.903, p = 0.01) as follows: where: N b /N z -an abundance of zoobenthos/an abundance of zooplankton; t-the number of months since the beginning of the study (June 2016 was taken as the first month). Harpacticoida, as a rule, were simultaneously present in plankton and benthos, but their planktonic share in the total abundance significantly increased with increasing salinity, as, for example, in June 2019 ( Figure 5). This dependence can be described (R = 0.780, p = 0.025) as follows: Y = 0.0002 e 0.109 S , (8) where: Y-the proportion of planktonic Harpacticoida in their total plankton and benthos abundance; S-the salinity, g L −1 . Other factors can significantly modify this dependence; therefore, it did not manifest itself as clearly in all periods as in July 2019. A similar trend was noted for the larvae of Chironomidae, but it is not significant. The observations indicate that the bottom oxygen deficiency significantly affected the distribution of larvae between plankton and benthos. Strong winds had a significant impact on the distribution of the total pool of animals in the bay between plankton and the bottom. For example, in November 2018, near the leeward coast, on average, 28% (CV = 1.399) of all recorded animals and protists were found in plankton, and only 7% near the windward coast (CV = 0.883).
The analyses of Tables 3 and 4 showed ( Figure 6) that there is a significant negative dependence of the number of classes of animals in plankton and/or benthos on the maximum salinity at which they were found in the lagoon (R = 0.969, p = 0.001), as follows: where: K s -the number of animal classes in Bay Sivash, S-the salinity, g L −1 .

Discussion
The use of new and published data [12,[18][19][20] showed that since the closure of the North Crimean Canal, salinity, TSS, and DOM increased in Bay Sivash ( Figure 7). As was shown before, the content of TSS positively correlates with salinity in different water bodies of Crimea [12,22] and different regions [23][24][25][26][27]. There are different reasons for the observed correlation of increases in salinity, TSS, and DOM, and the main one in Bay Sivash may be a positive correlation between liquid density with salinity that resulted in a lower sedimentation rate of single-sized particles [12]. Additionally, the disappearance of reed beds may be responsible for increasing the mineral particle flow into the bay [12]. Dramatically decreased reed beds of P. australis do not prevent erosion of small clay cliffs on the western coast of the lagoon now. Currently, muddy water stains near those cliffs have been observed regularly [4]. Accelerated cliff erosion may also contribute to the TSS increase in the Bay Sivash. The increase in DOM concentration can probably be explained by the TSS increase and/or the addition to primary production due to the development of floating mats of filamentous green algae Cladophora [4,12,19]. The intensive development of Cladophora in hypersaline water bodies can increase their primary production by more than 10 times [3,28]. In November 2018, the abundance of benthic microalgae related positively to an increase in salinity and TSS [19]. Another reason for an increase in DOM is the increased release of exopolymers by algae into the water. Under increased salinity, a larger amount of primary production is excreted into the environment as exo-metabolites [29][30][31][32]. However, changes in the bay environment caused by an increase in salinity are not limited to changes in the concentration of TSS and DOM. A cascade of interrelated changes in the physical and chemical parameters of the environment was initiated. For example, the heat capacity of a solution decreases linearly with increasing salinity [12]. With a decrease in heat capacity, the temperature regime changes, in particular, less heat is required to heat a unit volume of water per 1 • C. Consequently, the same amount of solar radiation heating the water in the bay will increase more saline water to a higher temperature. The calculation, according to the known formula [12], showed it can heat water in the lagoon up by 2-5 • C more than in the period before the closure of the North Crimean channel. A wider range of daily temperature fluctuations now occur. This effect of a salinity increase on the water temperature was observed previously in the Dead Sea and hypersaline lakes of Crimea and India [12,33]. The solubility of oxygen in water exponentially decreases with a salinity increase [12,34]. A temperature increase acts the same way. Currently, in different seasons and different parts of the bay, according to calculations, the oxygen concentration should be 20-50% lower than before the start of the salinity increase. An increase in salinity, as well as changes in other physical and chemical parameters, shift the geochemical background of the lagoon ecosystem, including the solubility of various salts, as well as direction and character of some reactions [35][36][37][38]. All this may often increase the concentration of toxic substances adding more extremality to the environment.
With an increase in salinity in the lagoon, the massive development of floating mats of green filamentous algae of Cladophora began, which was also shown for other Crimean hypersaline water bodies with a significant increase in their total primary productivity [3,28]. This often led to the formation of hypoxic and anoxic conditions at the bottom, which often led to a decrease (to zero) in the diversity and number of benthic animals [4,12]. At the same time, Cladophora filaments are intensively overgrown with various microalgae, seen also in Bay Sivash [19,32], significantly improving the food supply of animals and leading to an increase in their total number. Normally, benthic animals move to live in the Cladophora floating mats [4,12].
The decrease in the number of taxa of different levels with an increase in salinity above 50 g L −1 is a well-known phenomenon, and the dependence at the global level has been quantitatively characterized [39]. Using the data of that article, the authors calculated the following equation (R = 0.987, p = 0.001): where: K g -the number of animal classes at the global level; S-the salinity, g L −1 . Calculated Equations (9) and (10) showed that, despite a much larger number of classes at all salinities at the global level, the general trend is the same in both cases. In the case of Bay Sivash, the number of classes decreased with an increase in salinity from 50 to 150 g L −1 by 3.6 times, and at the global level, with the same rise in salinity, the number of animal classes decreases 1.5 times.
As new data and the previously published articles cited above show, with an increase in the salinity of the environment in Bay Sivash significant changes in the structure of zooplankton and benthos occurred, similar to those observed in similar situations in different water bodies of Crimea and other regions [3,28,40,41]. Summarizing all the available information on these changes, the following generalized trends can be distinguished: (1) In Bay Sivash, the most dramatic changes occurred with an increase in salinity in the range from 25 to about 70 g L −1 , i.e., in the period from 2014 to 2017. The same can be said about other similar well-known cases; (2) Changes in different taxa occurred unevenly, although a general tendency towards a decrease in taxonomic diversity was present in all of them; (3) With an increase in salinity, the proportion of calanoid copepods in the total abundance of zooplankton decreased-they practically disappeared from plankton at salinity above 80-100 g L −1 ; (4) At salinity above about 50-60 g L −1 , Chironomidae larvae began to massively develop in the ecosystem of the bay, and since 2014, they have rapidly increased their contribution to the abundance of both benthos and plankton. Their larvae are actively consumed by both fish and invertebrates, for example, amphipods [42]. Therefore, it can be confidently assumed that at lower salinities, when there are many predators, they are quickly consumed; (5) The concentration of benthic-planktonic species increases in plankton, in particular, Harpacticoida and Chironomidae; (6) At salinity above 80-90 g L −1 , nauplii and adult brine shrimp appeared and became abundant in plankton and benthos. Below this salinity, fish, amphipods, and ostracods are usually present in water bodies, which eat up all brine shrimp at the naupliar stage [43]. At higher salinity, as in our case, brine shrimp begin to dominate in plankton. It should be noted that the appearance and dominance of Artemia in water bodies leads to a significant restructuring of intra-ecosystem relationships [40,44]. In particular, the concentration of phytoplankton sharply decreases, and the transparency of the water increases significantly. This may be one of the factors responsible for the high concentrations of benthic microalgae in Bay Sivash [19]. It is interesting to note that when salinity in reservoirs begins to decrease, brine shrimp persist in plankton approximately down to a salinity of 50-55 g L −1 , and when salinity begins to increase in reservoirs, adult brine shrimp appear at salinity above 80 g L −1 [40,45]. This covers the hysteresis effect, which is common for system oscillations between alternative states [45,46]; (7) The variability of different taxa were practically not correlated with each other. There was also no clear relationship between changes in plankton and benthos; (8) The correlation between changes in abundance of different taxa and variations in environmental factors, if they were observed in some cases, was very weak. In the salinity range from 40-50 g L −1 to about 120-140 g L −1 , a significant role in the dynamics of zooplankton was probably played not by salinity itself, but by a complex of all interacting factors. In this interaction, an important role is played by unaccounted factors, internal and external, such as the strength and direction of the wind, the precipitation regime, and the drift of water masses from the Sea of Azov. This is typical for similar lagoons and lakes [47][48][49][50]. An equally important role is played by the factor of chance, which, for example, may be associated with the feeding of birds that have arrived episodically in the lagoon or extreme climatic events. The important role of the factor of chance in the dynamics of the biotic structure of ecosystems is now widely recognized [51][52][53].
The general conclusion following from all this is that salinity by itself is not a reliable predictor of possible ecosystem changes in the range from about 25 to 120-140 g L −1 . Salinity becomes a hard environmental filter for the existence of certain animal species at its value above 120-140 g L −1 [54]. Consequently, the ecosystem outcomes of anthropogenically induced changes of salinity in saline water bodies are poorly predictable and should therefore be avoided. The transit of ecosystems to new alternative states is always uncertain [55], and new data confirm this conclusion. That is why it is difficult to adapt the management of the external environment to such shifts in ecosystems [8].
A comparison can be made between the rates and ranges of salinity changes in two periods when it decreased after the canal was built and when it increased after the closure of the canal. When the canal was built, salinity decreased from 140 to 17-25 g L −1 , i.e., approximately by 7 times in 40 years [15], and when the canal was closed, it increased from 22 to 94 g L −1 , i.e., about by 4 times over 6 years. The calculation showed that the relative rate of change (g L −1 year −1 ) during the period of salinity decline was about 6 times less than during the period of increase. Drainage water from the canal led to a significant decrease in salinity also in the hypersaline lake Kyzyl-Yar, where the period of salinity decline in the lake lasted 30 years [45]. At the same time, the mineralization of water in the bottom deposits declined significantly slower and did not finish after 45 years. This fact showed that the speed of salinity changes in saline water bodies is determined not only by water budget but also by features of the bottom sediments as one of the aquatic ecosystem components.
Transformations occurred two times-not only in the lagoon, but also in surrounding territory and the avifauna. During the first shift, reed vegetation had occupied large areas, which stabilized the higher water level in the lagoon with reducing its seasonal fluctuations, and the bird composition changed [12]. Ciconiiformes birds, which were absent in the hypersaline stage, became common in the Bay Sivash area [56]. In 2014-2018, after the canal closing, the total diversity and abundance of visually observe birds decreased. In 2020-2021, after increasing Artemia and chironomid larvae abundance, bird presence increased, and new species were present. The most interesting case is Phoenicopterus roseus Pallas, 1811, which started to nest in Bay Sivash in 2017 [57]. Adult and young greater flamingo Ph. roseus were observed in Bay Sivash on 29 September 2021 ( Figure 8) (photo by G. Prokopov). As is known, flamingos prefer to use such hypersaline water bodies with a similar set of invertebrates in various regions, in particular, in South America [58][59][60]. The introduction of flamingos into Bay Sivash can lead to a new round of ecosystem changes as was previously noted in other water bodies [61]. To conserve the biodiversity of the unique hypersaline lagoon and birds, further monitoring of changes in its ecosystem is required. Author Contributions: Conceptualization, N.S.; methodology, N.S. and E.A.; field investigation and sampling, N.S., E.A. and A.L.; sample processing, species identification, E.A., E.K. and T.R.; data formal analysis, N.S. and E.A.; writing-original draft preparation, N.S.; writing-review and editing, N.S., E.A., E.K., T.R. and A.L.; supervision, E.A.; project administration, E.A.; funding acquisition, E.A. All authors have read and agreed to the published version of the manuscript.
Funding: Sample processing, species identification, data analysis, and manuscript writing was supported by the Russian Science Foundation, grant number 18-16-00001; the part concerning the long-term study of the lagoon environment was conducted in the framework by the state assignment of A.O. Kovalevsky Institute of Biology of the Southern Seas of RAS, number 121041500203-3; TSS and DOM measurements were conducted in the framework by the state assignment of Marine Hydrophysical Institute of RAS, number 0555-2021-0003.
Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: All data used in this study are available upon request from the corresponding author.