Assessment of River-Sea Interaction in the Danube Nearshore Area (Ukraine) by Bioindicators and Statistical Mapping

: There is a lack of understanding of the main drivers that form the picture of biological communities of transitional waters in deltaic ecosystems under the inﬂuence of terrestrial sources. Analysis of hydrochemical parameters in relation to phytoplankton communities in the Ukrainian part of the Danube coastal zone (in August 2018) is the focus of current work. The goal was to identify patterns in the distribution of environmental parameters (salinity and nutrients) in the area of the shipping channel through the Bystry arm, as well as to assess the state of water quality. The ecological bioindicators approach using modern statistical methods, and ecological mapping shows sufﬁcient achievements in interpreting the results. The indicators of salinity (mesohalobes) had better describe the character of the transportation of fresh riverine waters than salinity gradient. The composition of 35 indicator phytoplankton species corresponds to 3 and 4 water quality classes in the coastal zone. High N:P ratios showed an imbalance in the ecosystem as an indicator of production and destruction processes. Statistical maps of the indicator species distribution revealed the river current’s inﬂuence on the nearshore water mass. Ecological maps of surface and bottom variables show various environmental impacts resulting from dredging in the shipping channel and excavated soil dumping. Canonical correspondence analysis (CCA) and statistical maps revealed two pools of factors with oppositely directed effects on phytoplankton: salinity, on one hand, and nutrients, on the other. Miozoa and Chlorophyta have an opposite interaction with salinity and oxygen and can be ecosystem change indicators in further analysis for the studied area. The biomass of phytoplankton was presented for 63.2% by diatoms (523.1 mg m − 3 ) and for 33.6% by Miozoa (303.2 mg m − 3 ). The NT:PT from 34:1 to 47:1 and from 27:1 to 72:1 in the surface and bottom layers, respectively, which indicate the imbalance of the ecosystem with respect to nitrogen.


Introduction
The Danube River plays a great role for the Black Sea by forming the hydrochemical regime in the northwestern shelf, which is the biggest shallow area in the sea with the high stocks of biological production. The Danube, together with the Dnieper and Dniester rivers, make up to 70% of the input in the northwestern part of the Black Sea [1][2][3]. Understanding the processes of interaction of one of the largest rivers of Europe with the sea is of great importance.
The Danube nearshore area represents a transitional river-sea contact zone. It is formed due to dynamic interaction, mixing, and transformation of the river's water masses and the sea. The Ukrainian part of the nearshore area covers the region from the edge of the Kiliya Arm (Kiliya Delta of the Danube River) to the seawaters' border with a salinity of 17‰. The nearshore outer border is conditional because its location depends on the runoff rate of the Danube River and the wind regime [1,2]. The vertical distribution of waters at the nearshore area is characterized by sharp salinity changes and temperature, especially in spring and early autumn. It is noted a wide variability and fractal in the distribution of hydrochemical and biological parameters. They are associated with the seasonal variability of river runoff in the nearshore, which is the primary source of nutrients that ensure the seasonal development of production processes [2]: the special hydrologicalhydrochemical and hydrobiological regimes in the transition waters of Danube nearshore results in increased biological productivity. The Danube water flowing to the seaside occupies its upper layer up to 5 m. As a rule, these are seaside areas located directly near the mouths of large delta's arms. With a distance from the mouth, the river water gradually transforms into seawater, passing through three successive zones: the principal (salinity up to 10‰), frontal (10-12‰), and final (12-17‰) [3].
There is only one navigable arm in the Ukrainian part of the Kiliya Danube Delta-the Bystryi Arm. To maintain the Danube-Black Sea channel's working depths through the Bystryi Arm, periodically dredging operations are carried out. These operations' impact is overlapped with the river-sea interactions that form the main background in these transition waters. The salinity of the nearshore water masses at a depth of more than 5 m makes 16-18‰. The dumping area is located at a distance of 8 km from the delta edge, in the final zone. Since 2004, the Institute of marine biology of the NAS of Ukraine organizes environmental monitoring to determine the impact of these manipulations on the ecotone marine-river ecosystem [4,5]. However, there is still a lack of a general description of the ecosystem functioning and understanding of how river-sea interactions influence biological communities.
The studies of 2004-2011 describe phytoplankton's biodiversity in the Ukrainian part of the Danube nearshore area with 409 taxa [6,7]. The most remarkable species richness is typical for diatoms (155 species) and green (122) algae. They are followed by dinophytes (70) and cyanobacteria (34). A small number of taxa represents golden, euglenoid, and cryptophyte algae. Compared to the Romanian part of the Danube Delta, 957 species of microalgae are indicated, excluding the section of green ones, which include another 526 species [8].
The monitoring of water bodies, especially assessing the anthropogenic pressure, implies following the European directives' requirements (Water Framework Directive and Marine Strategy Framework Directive) [9,10], using different approaches to assessing the quality of the aquatic environment [11,12]. One of the most effective methods for water quality assessment is bioindication based on species composition shifts as a response of biota to the environmental changes and influence from land base sources [13][14][15].
In the Black Sea and many other marine water areas, species indicators' elaboration suitable for environmental monitoring is in the initial stage [16][17][18]. The level of nutrients and chlorophyll concentrations is widely used as an indicator of the eutrophication assessment [19]. In the northwestern part of the Black Sea and western part of the Azov Sea, the bioindication method was started for benthic species [20][21][22][23]. Colleagues from Ukraine, working within the framework of the development of a national marine strategy, are actively introducing a morphological and functional approach to assess the state of the environment using indices of the specific surface area of macrophytes and microalgae [24][25][26]. Integrated indices and comparative statistics, using various structural indicators of marine and estuarine ecosystems, is widespread [27,28].
Modern statistical programs and mapping have already been tested for many freshwater objects in different climatic zones [13][14][15][29][30][31]. However, applying these approaches to assessing the state of marine ecosystems has not yet received widespread use due to the limited number of indicator species [20,21,32]. The combined approach of studying the species composition of algae and statistics in the coastal area has good examples of its application [28]. We believe it rational to apply this complex of methods for the Danube nearshore area-a dynamically changing ecosystem, using phytoplankton's structural indicators and the content of nutrients. In addition, we think it is essential to contribute to the development of the database on marine species indicators that could be the base for assessing the ecological class of the water areas.
We hypothesize that it is possible to reveal the Danube Delta waters' influence on the coastal ecosystems of the Black Sea using bioindication of phytoplankton communities and statistical methods. The work aimed to assess the hydrochemical state and phytoplankton community of the Ukrainian part of the Danube nearshore area in the summer employing statistical analysis and ecological mapping. This work will be the base for the long-term analysis to identify the patterns in the distribution of environmental parameters and phytoplankton in the studied region.

Materials and Methods
The materials used in the study were obtained during environmental monitoring at the Ukrainian part of the nearshore area of the Danube Delta in August 2018 ( Figure 1). Water samples for the determination of hydrochemical parameters and phytoplankton were simultaneously taken in the surface and bottom layers using a Niskin bathometer. Determination of hydrochemical parameters: salinity, dissolved oxygen and% of its saturation, and pH value were carried out onboard the vessel. Dissolved inorganic nitrogen (DIN = NH 4 + + NO 2 − + NO 3 − ), dissolved organic nitrogen (DON), total nitrogen (NT), dissolved inorganic phosphorus (DIP = PO 4 3− ), dissolved organic phosphorus (DOP), total phosphorus (PT), and silicic acid (SiO 3 2− ), total suspended matter (TSS) in water samples were determined in a stationary laboratory of the Institute of Marine Biology of the NAS of Ukraine by standard methods accepted in international practice [33].

Hydrological and Hydrochemical Parameters
The distribution of hydrochemical parameters at the Danube nearshore area in August 2018 was characterized by significant spatial variability, which is associated with both the two-layer structure of the waters and the development of production-destruction processes ( Table 2 and Table A1; Figures 2-6).  Phytoplankton samples were fixed in 4% neutral formaldehyde and concentrated by sedimentation method to a volume of 40-80 mL. Sample processing was carried out using a Konus Biorex microscope at a magnification of ×160-640. Cells were counted in a 0.05 mL counting chamber. To calculate the biomass, we used the "true" cell volumes calculated on the basis of measurements of the linear dimensions of the cells [34]. For species identification and nomenclature, we used the following sources [35][36][37][38][39].
The index of saprobity was chosen as it reflects the organic pollution and includes not only diatoms but other algae that are rather represented in studied river-mouth communities [40,41]. Indices saprobity (S) were calculated based on identified species for each community and quantitative investigations of phytoplankton as: where S-saprobity index of the algal community; s i -species-specific saprobity index; a i -species abundance. The ecosystem state index WESI was calculated for each community based on the classification ranks of the nitrate-nitrogen content and the index of saprobity of the community [13]. The classes of water quality were given according to [31] ( Table 1). The sites that were included in the analysis are distinguished according to the tasks of ecological monitoring: dredging operations on the offshore approach channel to the mouth of the Bystryi Arm with depths of 0-10 m (sites DR0, 1, 2, 3); dumping area of dredging soil with depths of 23-25 m (sites DUM 1, 2, 3, 4); background areas to the north and south of the mouth of the Bystryi Arm and dumping area (sites B1, 2, 3, 4) with depths of 17-23 m. Table 1. Determination of water quality (EU color codes) based on saprobity index (S) according to [31]. Environmental mapping was carried out in the Statistica 12.0 program according to the parameter values and geographic coordinates of each site. Canonical correspondence analysis (CCA) was done with the CANOCO Program 4.0 [42], calculation of similarity was doing as the network analysis in JASP on the botnet package in R Statistica package of [43].

Hydrological and Hydrochemical Parameters
The distribution of hydrochemical parameters at the Danube nearshore area in August 2018 was characterized by significant spatial variability, which is associated with both the two-layer structure of the waters and the development of production-destruction processes (Tables 2 and A1; Figures 2-6). ency (Secchi) is demonstrated in Figure 6a,b showing the transition of particles in the surface layers from the river to sea. The range of variability of dissolved oxygen in the surface layer was 93-139% of saturation, pH-8.39-9.12 (Figure 6c,d). In the near-bottom horizon at station B2, the water was not enough saturated with oxygen (41%). The ratio NT: PT varied from 34:1 to 47:1 and from 27:1 to 72:1 in the surface and bottom layers, respectively, in different parts of the nearshore area in August 2018 (Table  2, Figure 7e,f).           The distribution of total suspended substance (TSS) and gradient of water transparency (Secchi) is demonstrated in Figure 6a
Between the diatoms, Cyclotella caspia Grunow 1878, C. meneghiniana Kützing, Pseudonitzschia delicatissima (Cleve) Heiden, Pseudosolenia calcar-avis (Schultze) B.G.Sundström were met frequently. Between the Miozoa Prorocentrum cordatum (Ostenfeld, 1901) Dodge, Prorocentrum micans Ehrenberg. We also noted the diatoms Proboscia alata rarely-just on three sites, in the dredging and dumping area and were presented with two species (Table 3). Out of the 66 species found, 35 were indicative (Table 3), with almost half of the species having an indicative role only in relation to salinity and habitat. Our bioindication results show that most species are planktonic (54.3%) and plankto-benthic (45.7%). More than half of the species (54.5%) were marine, and 27.3% oligohalobes-indifferents. The highest values of the abundance and biomass of phytoplankton in surface layers were registered in the zone of dumping with maximum on-site DUM1 (Table 4, Figure 2a,b). Biomass and abundance of phytoplankton in other regions were distributed homogeneously. In the bottom layers, the situation is different, showing the maximum growth of phytoplankton abundance in the river mouth area (sites DR0, DR1, DR2), whereas the maximum biomass was noted in site DUM4. This was due to the presence of dinoflagellates (Miozoa) which results in high biomass because of big cell volume.
The abundance of phytoplankton in the surface layers of the dredging area was twice higher than on the bottom one, whereas, in other areas, these discrepancies were much more considerable ( Table 4). The diatoms dominated (523.1 mg m −3 ), their contribution to the total phytoplankton biomass made 63.2%. Biomass of Miozoa made 303.2 mg m −3 or 33.6%, Chlorophyta-1.7 mg m −3 (0.2%), Cyanobacteria-8.2 mg m −3 (1.0%). Other groups of microalgae met rarely; their biomass did not exceed 1.5%.
In Figure 8a, the indicators of salinity-mesohalobes (mh) that survive in brackish waters demonstrate the character of the distribution of riverine freshwaters in the nearshore area. The presence of mesohaline (mh) species demonstrates dividing the nearshore area into three water masses to the south from the mouth, the north one, and the south one. Species-indicators of class 3 of water quality (Figure 8b) and nitrogen-autotrophic taxa (ate) (Figure 8c) represent the spreading of riverine waters within the marine water. Maximums of plankto-benthic species (Figure 8d) are also associated with the character of water flow from the river. Waters from the dumping area with a bigger amount of plankto-benthic species are transported with intensive water flow to the background sites.

Indices
Analyzing the distribution of the index of saprobity, we see the input of the organic nitrogen and phosphorus from the area of dredging (  (Table 4), which corresponds in average to the class 3 of water quality (beta-mesosaprobity), but in the mouth of the river it increases up to class 4. The data on the water quality class are given in Table 5. The index of the negative impact on the aquatic ecosystem (WESI) turned out to be much higher than 1, which indicates the ability of the ecosystem of the Danube nearshore area to counteract factors negative for photosynthesis. However, in Figure 7c,d, it can be seen that the lowest index values are associated with dredging at both surface and bottom layers.   (Table 4), which corresponds in average to the class 3 of water quality (beta-mesosaprobity), but in the mouth of the river it increases up to class 4. The data on the water quality class are given in Table 5.  The index of the negative impact on the aquatic ecosystem (WESI) turned out to be much higher than 1, which indicates the ability of the ecosystem of the Danube nearshore area to counteract factors negative for photosynthesis. However, in Figure 7c,d, it can be seen that the lowest index values are associated with dredging at both surface and bottom layers.

Indices
Moreover, the maximum of quantitative characteristics of the phytoplankton on the surface horizon coincides with the level of saprobity index and WESI index. The phytoplankton distribution on the bottom layer in its turn coincides with the saprobity index that represents the two spots of microalgae growth: in the river mouth area and in the zone where the soil from the dredging area was dumped.

Diagrams and Clusters
The relation of the quantitative development of phytoplankton with the main environmental parameters on the surface and in the bottom layer is shown in Figure 9a,b. For both horizons, we mention two vectors of the distribution patterns of the factors: the effect of salinity on phytoplankton was opposite to the nutrients. Index of saprobity (SLA) and WESI demonstrate the influence of environmental parameters (organic and toxic contamination) on the phytoplankton community, so they are included in the CCA. ronmental parameters on the surface and in the bottom layer is shown in Figure 9a,b. For both horizons, we mention two vectors of the distribution patterns of the factors: the effect of salinity on phytoplankton was opposite to the nutrients. Index of saprobity (SLA) and WESI demonstrate the influence of environmental parameters (organic and toxic contamination) on the phytoplankton community, so they are included in the CCA.
CCA analysis showed that the most indicative groups of surface phytoplankton were dinoflagellates (Miozoa), which density increased at the surface in low oxygen concentrations and higher ammonium (N-NH4). The opposite groups of phytoplankton were cyanobacteria and Euglenophyta that are managed by low nutrients and higher oxygen, salinity and SLA. On the surface layer, two more vectors were presented: ammonium (N-NH4) and the number of species (No Spec) opposite to transparency (Secchi) (Figure 9a). The growth of diatoms in the surface horizon matches with the green algae and regulated by nutrients contrary to dinoflagellates. Figure 9b demonstrates that in the bottom layer, the environmental parameters influencing the production of phytoplankton biomass were divided into two clusters, 1) salinity, transparency (Secchi), WESI and 2) nutrients and oxygen. Green and diatom algae, as well as cyanobacteria, react most actively to factor group 2. At the same time, dinoflagellates and Cryptophyta developed in conditions of increasing salinity, transparency (Secchi), WESI and decreasing nutrients and oxygen. The suspended substance (TSS) parameter was united with nutrients both on the surface and bottom layer.  CCA analysis showed that the most indicative groups of surface phytoplankton were dinoflagellates (Miozoa), which density increased at the surface in low oxygen concentrations and higher ammonium (N-NH 4 ). The opposite groups of phytoplankton were cyanobacteria and Euglenophyta that are managed by low nutrients and higher oxygen, salinity and SLA. On the surface layer, two more vectors were presented: ammonium (N-NH 4 ) and the number of species (No Spec) opposite to transparency (Secchi) (Figure 9a). The growth of diatoms in the surface horizon matches with the green algae and regulated by nutrients contrary to dinoflagellates. Figure 9b demonstrates that in the bottom layer, the environmental parameters influencing the production of phytoplankton biomass were divided into two clusters, 1) salinity, transparency (Secchi), WESI and 2) nutrients and oxygen. Green and diatom algae, as well as cyanobacteria, react most actively to factor group 2. At the same time, dinoflagellates and Cryptophyta developed in conditions of increasing salinity, transparency (Secchi), Land 2021, 10, 310 13 of 18 WESI and decreasing nutrients and oxygen. The suspended substance (TSS) parameter was united with nutrients both on the surface and bottom layer. Figure 10 shows the similarity of studied parameters on monitoring sites (p < 0.05). First, we see that in surface layer background sites have similarities both to dredging and to dumping sites. This explains their different geographical location. The data are divided into 4 separate clusters (Figure 10a); two of them are closer to one another-these are dredging areas and two background sites (B3, B4). They can be united in one cluster because of light links between them and due to joint hydrodynamic and hydrological conditions under the influence of the river discharge (Figure 1). Clusters 3 and 4 in Figure 10a show the similarity of the dumping area and two other background sites (B1, B2) that are located in much more deep areas of the nearshore area of the Danube. The strong positive links in the surface horizon were mentioned for DUM4 and B2; DR1 and DR0.
Land 2021, 10, x FOR PEER REVIEW 13 of 18 Figure 10 shows the similarity of studied parameters on monitoring sites (p < 0.05). First, we see that in surface layer background sites have similarities both to dredging and to dumping sites. This explains their different geographical location. The data are divided into 4 separate clusters ( Figure 10a); two of them are closer to one another-these are dredging areas and two background sites (B3, B4). They can be united in one cluster because of light links between them and due to joint hydrodynamic and hydrological conditions under the influence of the river discharge ( Figure 1). Clusters 3 and 4 in Figure 10a show the similarity of the dumping area and two other background sites (B1, B2) that are located in much more deep areas of the nearshore area of the Danube. The strong positive links in the surface horizon were mentioned for DUM4 and B2; DR1 and DR0.
However, on the bottom horizon, we note a little bit different situation: B1 and DR3 sites were united with the dumping sites in one cluster (Figure 10b). This distribution is explained by the depth differences between DR2 and DR3, which is critical for phytoplankton in the bottom layer. The strong positive links on the bottom horizon were noted for the sites DR3 and DUM4; B3 and B4; DR1 and DR0.  Table 2. Sets of stations in the legend. The line thickness between stations reflects the correlation value; green is positive, red is negative. Clusters marked by dashed lines and numbered.
Thus, we can see significant differences between the stations in their depth and type of anthropogenic impact. In general, the comparative statistical analysis showed differences in the strength and direction of the environmental impacts for the surface and the bottom. On the surface, it is possible to trace the flow from the impact of dredging on the background stations, whereas at the bottom, the background stations are well separated from the dredging and dumping sites, representing more stable conditions in these areas of the nearshore.

Discussion
The bioindication approach and ecological mapping help to reveal the main vectors of the interaction of environmental parameters and phytoplankton community in a very dynamic region of the Danube nearshore area. Application of these methods, based on nutrient concentrations and species bioindicators, was successfully tested for many water bodies and revealed the spatial dynamics and seasonal fluctuations of pollution in fresh and brackish waters [13,14,31,44]. For the analysis in the present study, we used 35 indicator species, which is 53% out of all found species in the study area of river-sea interaction. This testifies the necessity to find more indicator species in transition waters and the relevance of further research. It is extremely important to conduct a more detailed study of the species diversity in the Danube nearshore area in order to identify the species that form the signal about certain changes.  Table 2. Sets of stations in the legend. The line thickness between stations reflects the correlation value; green is positive, red is negative. Clusters marked by dashed lines and numbered.
However, on the bottom horizon, we note a little bit different situation: B1 and DR3 sites were united with the dumping sites in one cluster (Figure 10b). This distribution is explained by the depth differences between DR2 and DR3, which is critical for phytoplankton in the bottom layer. The strong positive links on the bottom horizon were noted for the sites DR3 and DUM4; B3 and B4; DR1 and DR0.
Thus, we can see significant differences between the stations in their depth and type of anthropogenic impact. In general, the comparative statistical analysis showed differences in the strength and direction of the environmental impacts for the surface and the bottom. On the surface, it is possible to trace the flow from the impact of dredging on the background stations, whereas at the bottom, the background stations are well separated from the dredging and dumping sites, representing more stable conditions in these areas of the nearshore.

Discussion
The bioindication approach and ecological mapping help to reveal the main vectors of the interaction of environmental parameters and phytoplankton community in a very dynamic region of the Danube nearshore area. Application of these methods, based on nutrient concentrations and species bioindicators, was successfully tested for many water bodies and revealed the spatial dynamics and seasonal fluctuations of pollution in fresh and brackish waters [13,14,31,44]. For the analysis in the present study, we used 35 indicator species, which is 53% out of all found species in the study area of river-sea interaction. This testifies the necessity to find more indicator species in transition waters and the relevance of further research. It is extremely important to conduct a more detailed study of the species diversity in the Danube nearshore area in order to identify the species that form the signal about certain changes.
The level of trophicity, which is determined by the concentration of nutrients and their transformation from one form to another, is the determining factor in the development of the phytoplankton assemblage in the studied region ( Figure 9). The river-sea interaction results in the stratification of water masses: a surface layer with a salinity of 0.5-12‰ and a bottom layer formed by waters of marine genesis with a salinity of 16-18‰. This is an important aspect in such transition water areas located in the zone of river runoff influence. The formation of pelagic communities is determined by the complex hydrodynamic structure in the nearshore [2]. The main impact is the formation of a hydrofront-a zone of maximum horizontal gradients of water salinity. Its basic features are stable temperature stratification of water masses in summer and disturbance of diffuse exchange between the surface and bottom layers of the nearshore with the formation of bottom hypoxia. In addition, a great influence in the deltas is made by the rate of circulation of biogenic compounds of nitrogen and phosphorus with the participation of algal communities during production and destruction processes [45][46][47].
A peculiar feature of the hydrochemical regime of the Danube nearshore area is the imbalance of the ecosystem in terms of the content of the main biogenic elements-nitrogen and phosphorus, which is in line with previous research [2]. The balanced ecosystems are characterized with the ratio NT: PT = 7.2:1 (the standard ratio for an organic matter of plankton) [48]. The ratio of the mean values of NT and PT in the surface (44:1) and bottom (62:1) layers from the studied region also indicate the imbalance of the ecosystem with respect to nitrogen, which is associated with an excessive supply of mineral nitrogen compounds (nitrates) with river runoff (Figure 6).
Concentrations of phosphor compounds (DIP, DOP) and pH level testify to the higher trophic state in the region. It is known that DIP concentrations above 40 µg P L −1 are typical for hypereutrophic waters [31]. The pH values are not usually considered to the trophic level, but in the cited classification system [31], it is given that high pH usually reflected the high trophicity of a waterbody. In the bottom water masses with higher salinity, there is no intensive accumulation of DOP and DON, which is associated with the active development of destruction processes and is confirmed by the lack of oxygen saturation of water in certain areas of the nearshore.
The non-Redfieldian stoichiometry is mentioned to be the character of many coastal and estuary water areas [49,50]. Under the eutrophic conditions, the biodiversity of the pelagic assemblages decreases; the input of diatoms lowers and is replaced by dinoflagellates and cyanobacteria [16]. The processes of enrichment with nutrients are accompanied by the development of small-cells green algae (in our case Monoraphidium sp.) and bigger cells of Miozoa (Tripos sp.). Such processes, associated with the long-term enrichment of the aquatic ecosystem with nutrients, lead to changes in biodiversity in the Danube River basin, the largest waterway in Europe [51].
In addition, the ratio of nitrogen and phosphorus is unbalanced, which was noted earlier in this region [2]. For example, in 2004-2010, the ratio of dissolved mineral and organic forms of nitrogen and phosphorus NT: PT in the photic layer in the nearshore area made 50:1; this is due to the input of nitrogen compounds from the dredging area being larger than that of phosphorus. Nitrogen compounds are mostly represented by their organic component-the share of DON in the total balance of nitrogen forms is more than 80% for the surface layer and about 90% for the bottom layer. The destruction of dissolved forms of organic nitrogen compounds, which leads to the formation of mineral forms (reduced, i.e., ammonia and oxidized, i.e., nitrites and nitrates), causes the further eutrophication of nearshore waters because the newly formed mineral nitrogen compounds enter the biotic turnover [1,2].
Calculated indices of saprobity and the state of the ecosystem, as well as the list of revealed species indicators with its ecological preferences, clearly demonstrate a dynamic and high trophicity of the studied area. According to the indicator species, the studied waters are classified as medium eutrophic (beta-mesosaprobes), but the values of the index of saprobity show the water quality class 3 and 4 describing not only organic pollution but also the intensive self-purification processes. The WESI indices are significantly higher than 1 that also reflects the active work of the ecosystem both on the surface and in the bottom layers (Figure 6c,d). The same assessment of the state of the environment in the Danube Delta was obtained by other researchers [19]. However, in our case, a search and replenishment of species-indicators in the studied area should be the focus of future research.
As to the impact of the dumping and dredging effect, we may indicate that there is no clear understanding of this effect on the phytoplankton in the studied region, which is supported by other studies for zooplankton [52]. Some authors represent increasing in the autotroph's production [53], whereas another, it is diminishing [54]. Together with this, it was shown that the release of sediments during aggradation might cause the intensification of microalgae production [55]. The main factor influencing the algal community is the river inflow and the catchment area [56]. In this regard, the social-economic activities in the Danube River basin result in deltas and nearshore marine ecosystems; thus, sustainable planning should improve the water quality in the lower regions. The indicators of salinity (mesohalobes) had better describe the character of the transportation of fresh riverine waters than salinity gradient. This fact also supports the efficiency of the bioindication approach to reveal the patterns of plankton functioning in transition water bodies.

Conclusions
As a result of the analysis in August 2018, the 66 species of microalgae from six divisions were revealed in the Ukrainian part of the Danube nearshore area. More than half of the species were bioindicators with several ecological preferences: salinity, habitat, level of saprobity. Based on the index of saprobity, classes 3 and 4 of water quality were identified for the studied area. The biomass of phytoplankton was presented for 63.2% by diatoms (523.1 mg m −3 ) and for 33.6% by Miozoa (303.2 mg m −3 ). The NT:PT ratio varied from 34:1 to 47:1 and from 27:1 to 72:1 in the surface and bottom layers, respectively, which indicate the imbalance of the ecosystem with respect to nitrogen.
Thus, we confirmed our hypothesis about the effectiveness of bioindication by phytoplankton communities and statistical mapping to identify the impact of the Danube Delta waters on the Black Sea's coastal ecosystems, which was impossible to do only using chemical data. The use of statistical maps as an assessment model helped to identify the paths and range of penetration of the Danube waters into the coastal sea masses, and also showed the main trends in the transformation of river waters on the coast and ranked environmental factors in terms of their importance for the plankton community and the direction of impact. Two pools of factors were identified with oppositely directed influence on the development of the phytoplankton community: salinity and oxygen, on one hand, and nutrients, on the other. The nearshore area is a eutrophic ecosystem where the river is the main source of nutrients. Moreover, our results demonstrate that works on dredging and dumping additional cause entrance of the nitrogen, phosphorus, and siliceous. The imbalance of nitrogen and phosphorus and the velocity of its recycling were the leading factor of the processes of production and destruction. The most indicative groups of microalgae for the salinity and oxygen gradients were Miozoa and Chlorophyta.
The interaction of the Danube River mouth that has a great social-economic impact from most European countries, with the Black Sea, challenges the resilient approach in coastal and land-based management. Further application of the ecological mapping can help to assess the temporal and spatial dynamic of trophic relationships in the ecosystem in detail in transition waters, such as the Danube nearshore area.