Impact of the Construction of New Port Facilities on Primary Production of Plankton in the Neva Estuary (Baltic Sea)

A significant increase in the human population on marine coast and steady growth of maritime water transport causes the construction of port infrastructure and the creation of new lands, which affects the ecosystems of coastal waters. Despite the widespread occurrence of such large-scale engineering projects in coastal areas, their impact on various components of aquatic ecosystems, including phytoplankton, is still poorly understood. The aim of the study was to assess the effect of the construction of ports and the alluvium of new lands in the Neva Estuary in 2000s on the productivity of phytoplankton. Digging and dredging of bottom sediments results in one order of magnitude elevation of suspended particulate matter (SM), which mostly consisted of sand and clayed deposits and in significant decrease water transparency, as compared to the average long-term values. Concentrations of total phosphorus in the estuarine waters during the works significantly positively correlated with the concentrations of SM. However, the multiple increase in nutrients was less important for phytoplankton development than expected. Analysis of variance and stepwise multiple regression analyses showed that the main predictor of the primary production of plankton in the periods of construction was water transparency. Gross primary production decreased significantly. In contrast to short-term effects caused by wind-induced events, which often stimulated phytoplankton development, long-term construction works of new port facilities negatively influenced phytoplankton productivity. Apart from pristine conditions when the phosphorus concentration was the main factor limiting the primary production in the estuary, the main limiting factor during long-term engineering projects became water transparency. Taking into account plans for further development of ports in coastal areas around the world, the influence of the large-scale engineering projects on the conditions for the development of phytoplankton may provide a new aspect of long-term regulation of algal blooms and ecosystem functioning in the coastal and estuarine zones.


INTRODUCTION
The most significant growth in human population and urbanization is predicted for coastal marine areas in the near future (Neumman et al., 2015). This causes the active development of coastal infrastructure, the construction of dams, the alluvium of new territories and digging fairways, which can affect the ecosystems of coastal zones (Bianchi, 2009;Wu et al., 2018). Coastal environments, in particular heavily populated semi-enclosed marginal seas and coasts like the Baltic Sea region, are strongly affected by human activities. A multitude of human impacts, including climate change, affects the different compartments of the environment, and these effects interact with each other (Reckermann et al., 2022). In recent decades, new lands, ports and fairways have been created at the mouth of the Neva River, which have become the infrastructure of the five million megapolis -St. Petersburg (Golubkov and Alimov, 2010).
Modern climate changes are causing a rise in the water level in the Baltic Sea, especially in its eastern part (Pindsoo and Soomere, 2020). Moreover, it has been shown that over the last 50 years the absolute (geocentric) mean sea level in the Baltic Sea has rising at a rate slightly higher than the world average, and this trend is likely to increase in the future (Weisse et al., 2021). The Flood Protective Facility that consists of 11 dams, has been protecting St. Petersburg from floods since the late 1980s (Ryabchuk et al., 2017), but future sea level rise is likely to lead to the need for their expansion and modernization. In the past, such works have already led to the appearance of a large amount of suspended matter in the waters of the Neva Estuary and affected it's ecosystem (Golubkov and Alimov, 2010).
The effect of suspended matter on phytoplankton has been studied in some detail in natural shallow water. In these ecosystems, significant amounts of suspended matter can rise from the bottom into the water column as a result of wind mixing (e.g.; Su et al., 2015;Kang et al., 2019) or turbulent tidal currents (e.g., de Jorge and van Beusekom, 1995;Guinder et al., 2009). In addition, according to climate models, coastal ecosystems will face an increased frequency of turbidity due to terrestrial sediment run-off induced by future increasing precipitation (Remy et al., 2017;Hoshiba et al., 2021).
A significant increase in the content of suspended particulate matter (SM) reduces light penetration in water and worsens conditions for the development of phytoplankton (Su et al., 2015;Horemans et al., 2020;Kumar et al., 2022). In opposite, the decrease in amount of SM with a concomitant increase in the penetration of solar radiation may be one of the main causes for the development of phytoplankton blooms (Guinder et al., 2009). Moreover, nutrients that have accumulated in the sediment over time can be recycled back into the water because of the intrusion of nutrient-rich pore water from sediments or bottom water layers into the whole water column (Su et al., 2015). This may release phytoplankton from nutrient limitation and stimulate plankton primary production (Su et al., 2015;He et al., 2022), but its response depends on the season and the species composition of the algae (Lagus et al., 2007;Su et al., 2015).
In some cases, the combined effect of both factors is important (Marzetz et al., 2020). Phytoplankton is likely to be light-limited at times of extreme turbidity because at lower turbidity light intensity is lower the threshold of lightlimitation of phytoplankton species (Wang et al., 2019). Thus, the response of phytoplankton to resuspension of bottom sediments is rather complicated. In addition, it should be taken into account that resuspension events resulting from strong winds or tides usually last for several hours (e.g. de Jorge and van Beusekom, 1995;Martyanov and Ryabchenko, 2016). On the contrary, resuspension of bottom sediments can last for months or even years with the creation of new alluvial lands or port infrastructure. For example, the creation of new lands and the construction of ports at the mouth of the Neva River took several years (Ryabchuk et al., 2017).
Despite the widespread occurrence of mega-engineering works in coastal areas, their impact on various components of near-shore ecosystems, including phytoplankton, is still poorly understood. The aim of this study was to assess the effect of increased concentration of suspended particulate matter resulting from new port facilities construction on the primary production of plankton in the Neva Estuary. We tested the hypothesis that an increase in the concentration of nutrients in the water due to the resuspension of bottom sediments can compensate for the deterioration of light conditions for the development of phytoplankton in the Neva Estuary during the construction of port infrastructure.

Study Sites
The Neva Estuary ( Figure 1) is located in the northeastern part of the Baltic Sea. The Neva River, with an average discharge of 2,490 m 3 s −1 (78.6 km 3 yr −1 ), flows out of Lake Ladoga and goes into the Gulf of Finland. The catchment area of the estuary exceeds 280,000 km 2 (Golubkov and Golubkov, 2020). It is located at the northern boundaries of the temperate zone and at the southern boundary of the subpolar zone (Meteoblue, 2021). Climate type in the region according to Köppen-Geiger climate classification (Kottek et al., 2006) is Dfc -Snowy climate, fully humid, with cool summers.
The Neva Estuary is shallow, non-tidal and with a smooth salinity gradient from fresh water in the upper part to slightly saline in the lower part (Golubkov and Golubkov, 2020). Since the five millionth megalopolis of St. Petersburg is located on the coast of the estuary, its ecosystem is subject to strong anthropogenic impact (Golubkov et al., 2019). The Flood Protective Facility that consists of 11 dams have separated the upper reach (UR) of the estuary from its middle reach (MR) since the end of 1980s ( Figure 1). Temperature stratification is absent in the UR of the estuary, and low water transparency and rather strong currents prevent the development of benthic aquatic vegetation. The MR of the estuary is located between Kotlin Island and the conditional boundary of 29°10`E. In summer, there is temperature stratification there. A detailed description of both parts of the estuary was given earlier in Golubkov, 2020, 2021). The main abiotic and biotic characteristics of the estuary are presented in Table 1.
In [2006][2007]476.7 hectares of new territory were created in the UR of the estuary for the new Passenger Harbor of St. Petersburg (Figure 1). At the same time, the shipping channel was deepened to 14 meters for the passage of large ferries. The construction of the new multifunctional maritime shipping complex «Bronka» was carried out in 2014-2015 in the southwestern part of the UR (Figure 1). During these mega engineering works, large amounts of bottom sediment were dredged, removed and transferred from one part of the estuary to another (Ryabchuk et al., 2017). This led to the appearance of large amounts of suspended matter in the water.

Sampling
Samples were taken at 17 stations ( Figure 1) in mid-summer (late Julyearly August) 2003-2020. Secchi-disk depth (Sec), salinity (S) and temperature (T) were measured at each station. T and S were measured by the CTD90m probe (Sea&Sun Tech., Germany) every 20 cm from the surface to the bottom in the water column. Since, according to these measurements, the whole water column in the shallow UR was mixed, samples were collected from the entire water column. Five water samples (2 l each) were taken at each station: from the surface, half a meter from the bottom and from three equal depths between them. Samples from different depths were taken in order to avoid errors associated with the vertical distribution of different phytoplankton species in the water column. These samples were composited and mixed to make up a pooled sample (10 l). Samples of total phosphorus, chlorophyll a, suspended particulate matter and suspended particulate organic matter (three replicates of water collection) were taken from these pooled samples.
In the MR, pooled water samples were taken from the layer above the thermocline (LAT). Five water samples (2 l each) were taken from LAT: from the surface, the thermocline and from three equal depths between them. These samples were mixed to create a pooled sample (10 l). The samples for chlorophyll a, total phosphorus, suspended particulate matter and suspended particulate organic matter analysis (three replicates) were taken from these pooled samples.

Sample Analysis
Three hundred millilitres of water were filtered through 0.85 µm membrane filters (Millipore AAWP) to determine the chlorophyll a (CHL) concentration, which was followed by 90% acetone extraction and spectrophotometric determination (Grasshoff et al., 1999). Total phosphorus (TP) was determined after acid hydrolysis with the molybdate blue method (Grasshoff et al., 1999). Suspended particulate matter (SM) samples were obtained by filtering water over preweighed and precombusted filters (Whatman GF/F). SM concentration was determined with a standard gravimetric technique (Grasshoff et al., 1999). Suspended particulate organic matter (SOM) was determined after filtration (Whatman GF/F filters) with the dichromate acid oxidation (Grasshoff et al., 1999). Suspended particulate mineral matter (SMM) was calculated as a difference between SM and SOM. After that, percentage of SMM in SM was calculated.
The primary production of plankton (PP) in the water column were measured by the oxygen method of light and dark bottles (Hall et al., 2007;Vernet and Smith, 2007). Since the depth of the LAT practically coincided with the depth of euphotic zone in mid-summer, 900 ml of water from the LAT pooled samples in the MR, and from the whole water column in the UR were used to determine PP. Three 100 ml light and three dark bottles were filled with the water from each sampling station and exposed in an aquarium on the ship's deck in shadow during 6 h at a surface water temperature to estimate PP. Three 100 ml bottles (control bottles) were filled with the water from each sampling station to determine the oxygen contents in water at the beginning of the experiment. The Winkler method was used to determine the oxygen contents in the control, the light and the dark bottles (Hall et al., 2007). A more detailed description of the method and experimental design is given in Golubkov et al. (2017) and Golubkov and Golubkov (2020).

Statistical Analyses
Concentration of SM and percentage of SMM in SM was averaged for each station for periods 2006 -2007, 2014 -2015  and for 2003 -2020 excluding the previous periods, and visualized using SURFER 8.0.
The pairwise Pearson correlation coefficients between concentrations of SM, SOM, SMM, Sec, TP and phytoplankton productivity indicators were calculated using Microsoft Excel. Before calculating correlation coefficients, the data were logarithm. Linear regression analysis was produced using Microsoft Excel, the original data were logarithm.
Statistical analysis was performed using R software (version 4.0.5) (R Development Core Team, 2022). Analysis of variance (ANOVA) was performed using the "anova" command of the R Stats package (R Development Core Team, 2022). Stepwise multiple linear regression was performed to determine the set of factors most influencing plankton primary production. Models "model.null" and "model.full" were built. The "step" function (direction = "both") was then used to add and remove factors from the model and find the model with the lowest Akaike Information Criteria (AIC). AIC was used as an estimator of out-of-sample prediction error and thereby relative quality of statistical models for a given set of data. The lower the value, the better for the AIC (Mangiafico, 2022). When the model with lowest AIC was found it called final.model. In this way, combination of factors most influencing PP was found. The residual standard error (RSE) for the resulting model was calculated. RSE is a measure of the inconsistency of the model with the available data. If the RSE value is very close to the actual outcome value, then model fits the data well (James et al., 2013). After that, the adjusted R 2 (Adj R 2 ) was calculated to avoid the error of overestimating the quality of the models, arising from the addition of each additional variable to the models. Adj R 2 is a version of R 2 , which is calculated to take into account the number of variables in the multiple regressions model. Adj R 2 is always lower than the R-squared. The Fisher distribution for the model was also checked and reliable F-values were obtained, taking into account the degrees of freedom (df) of the model. In addition, the model was tested for unbiasedness and homoscedacity.

Abiotic Environmental Factors
The most intensive construction of new port facilities in the Neva Estuary was carried out in 2006-2007. Because of these works, the concentration of suspended particulate matter sharply increased in all parts of the estuary, especially in its upper reach ( Figure 2A). Suspended particulate matter concentrations during the works in the UR and MR of the estuary were about 20 times higher than the long-term average summer values (Figures 2A, C). The construction of the new multifunctional maritime shipping complex «Bronka» in 2014-2015 were carried out in the near-dam zone ( Figure 1) and increased turbidity only in the western part of the UR ( Figure 2B). The suspended particulate matter mainly came in the southern part of the MR, which led to an increase in its concentration (up to 30 g m -3 ) in the area of station 17 ( Figure 2B). Concentrations of SM in this part of the Neva Estuary in 2014-2015 were approximately 6-7 times higher than the average long-term values ( Figures 2B, C).
The analysis of SM showed that both in 2006-2007 and in 2014-2015 it mainly consisted of the mineral fraction ( Figures 2D, E). In 2006-2007, the share of SMM in SM exceeded 90% in almost the entire UR and in a significant part of the MR of the estuary ( Figure 2D). In 2014-2015, the area of increased concentration of suspended particulate matter was significantly smaller compared to 2006-2007 ( Figure 2B), but in the southwestern zone near the dam, the share of the mineral fraction was also higher than 90% ( Figure 2E). For comparison, when no mega-engineering works were carried out, the share of the mineral fraction in the SM was about 60% in the UR of the estuary, and less than 20% in its MR ( Figure 2F).
The analysis of Pearson`s pairwise correlations showed that the concentrations of SM and SMM during mega-engineering works exhibit a high statistically significant positive correlation between each other (r=0.95, р<0.001, n=62) ( Figure 3). The relationship between the SOM and SM was weaker (r = 0.55, p<0.001, n = 62) than between SM and SMM. Such a strong increase in SM led to a significant decrease in water transparency (Sec) during the years of construction of new port facilities, as evidenced by the analysis of Pearson's pairwise correlations ( Figure 3) and the regression dependence of Sec on SM ( Figure 4A). The concentration of SOM also significantly negatively correlated with Sec, but the strength and reliability of their relationship were weaker (Figure 3). Mega-engineering works also influenced the concentration of total phosphorus (TP) in the estuary. Concentrations of TP in the waters of the UR and MR of the estuary during the works was significantly positively correlated with the concentrations of SM and SMM (Figure 3). Although the Pearson`s correlation coefficient between TP and SM concentrations was not very high and there was a noticeable scatter of points, the relationship was statistically significant ( Figure 4B). TP concentrations also showed a weak but significant correlation with SOM concentrations (Figure 3).

Productivity of Phytoplankton
Analysis of pairwise correlations showed that plankton primary production was significantly negatively correlated with the concentrations of SM and SMM (Figure 3). Regression analysis showed a statistically significant dependence of PP on Sec during the periods of construction of ports and alluvium of new lands ( Figure 5A). In other words, PP in the water column decreased at stations where the highest SM were observed during the works, and PP increased where there was a low concentration of SM. At the same time, PP was statistically significantly negatively correlated with the concentration of TP ( Figure 5B), which, apparently, was a consequence of the positive correlation of TP with SM and SMM ( Figure 3). Analysis of variance showed that water transparency had the greatest impact on the primary production of phytoplankton in the Neva Estuary during the periods of port infrastructure construction ( Table 2). Salinity and temperature of the water, as well as the concentration of total phosphorus also affected the PP, but much weaker. Factors not taken into account in the study also had a significant influence (29.7%), but their role was less than that of water transparency ( Table 2).  Stepwise multiple regression analyses showed that of all the studied factors, the primary production of plankton was influenced by a combination of four factors: transparency, salinity, water temperature and total phosphorus concentration ( Table 3). The resulting multi-regression model has a high degree of reliability, because adjusted for the number of measurements, R 2 is high, p-value<0.001. In addition, Fischer's test for the model showed that F-value is high and statistically reliable at these degrees of freedom. Residual standard error also has a low value, which means that the simulation results describe the initial empirical data well ( Table 3). The resulting model is homoscedastic and unbiased.
The resulting multiregression equation shows that plankton primary production was positively related to water transparency, i.e. with its increase, the PP increased and vice versa. An increase in water salinity had a similar effect. In other words, PP was higher in clearer and more saline waters. However, this was most likely due to an increase in water salinity with distance from the epicenter of the port and alluvial works that took place in the mouth of the river and the upper reach of the estuary. Conversely, primary production was negatively related to phosphorus concentration and water temperature ( Table 3). This is consistent with the results of the regression analysis, which also showed that phosphorus concentration was strongly positively related to suspended matter concentration ( Figure 4B) and negatively to PP ( Figure 5B).

Concentrations of Suspended Matter and Total Phosphorus, Water Transparency
Periodic high increases in the concentration of suspended matter are common for shallow water, including estuaries (de Jorge and van Beusekom, 1995;Guinder et al., 2009;Su et al., 2015). Periodic rather high concentrations of SM and low Sec are also common for the Neva Estuary, especially for its shallow upper reach, even in the absence of mega-engineering works (Martyanov and Ryabchenko, 2016). Hydrobiological studies in the upper reach of the Neva Estuary, carried out at the beginning of the 20th century (1911)(1912), showed that the concentration of suspended particulate matter in the summer period was within 0.9-75 g m -3 , with an average concentration about 5 g m -3 and average Secchi-disk depth was 1.4 m (Zalessky and Wulf, 1913). The concentration of SM in 1959-1961 was 5-8 g m -3 in calm weather, but with a wind of 9-11 m s -1 , it rapidly increased to 20-34 g m -3 (Ryabinina, 1965). The maximum concentration of SM (89 g m -3 ) in those years was recorded along the southern coast of UR with a long east wind and a wave height of 1 meter. At the same time, the Secchidisk depth in summer was in the range of 0.3-2.1 m (Ryabinina, 1965). In the 2000s, when no mega-engineering works were carried out, the concentration of SM in the upper reach was in the range of 5-20 g m -3 , and in the middle reach, it decreased below 5 g m -3 (Figure 2). The average Secchi-disk depth in the 2000s, when no works were carried out, was 1.3 meters. Mega-engineering works also led to an increase in turbidity of water in various regions of the Baltic Sea. For instance, the SM concentration reached 126 g m -3 , and the water transparency decreased to 0.2 meters along the Estonian coast during the dredging period (Liblik and Lips, 2011). The Neva Estuary has been under strong anthropogenic impact for over a hundred years (Golubkov et al., 2003;Golubkov et al., 2019). Dredging and digging of bottom sediments in the Neva Estuary were carried out periodically and before 2000s (Ryabchuk et al., 2017), which led to increased turbidity. For example, the concentration of SM at some stations in the upper reach of the estuary exceeded 110 g m -3 during the construction of the Flood Protective Facility of St. Petersburg in 1982-1984(Umnova, 1987. However, the maximum concentration of suspended particulate matter recorded at station 7 (186 g m -3 ) in 2006-2007 was one and a half times higher than its maximum value in the 1980s.
Concentration of suspended particulate matter can affect phytoplankton in different ways (de Jorge and van Beusekom, 1995). For instance, SM resuspended from the bottom was an important source of nitrogen and phosphorus in the shallow estuary in the Gulf of Bothnia (Nilsson and Jansson, 2002). At the same time, an increase in turbidity has a negative effect on water transparency, the second most important factor in the development of phytoplankton (Kirk, 2011). An increase in the SM in the Neva Estuary during the works related to the creation of a new port infrastructure and the alluvium of new lands also correlated with these two environmental factors. On the one hand, this led to an increase in the concentration of TP (Figures 3, 4B), the main element limiting the development of algae in the estuary (Golubkov and Golubkov, 2020). It was also shown that the bulk of nutrients entering the Neva Estuary through various mechanisms are associated with inorganic substances in the suspended particulate matter (Lehtoranta et al., 2004). On the other hand, the increase in concentration of suspended matter due to these works negatively affected the transparency of water in the estuary (Figures 3, 4A).

Primary Production of Plankton
The large-scale construction of mega-engineering facilities can have a greater impact on estuarine ecosystems than the facility itself (He et al., 2022). This usually leads to an increase in turbidity and nutrient concentration, which affects phytoplankton (Wang et al., 2019).
In the Nevа Estuary, an increase in the concentration of phosphorus as a result of the resuspension of a large amount of suspended particulate matter into the water during the periods of mega-engineering construction should have resulted in the intensive development of phytoplankton, but this did not happen. The concentration of chlorophyll a did not increase, and the primary production of plankton during the years of construction was the lowest in the entire 2000s . This was obviously a consequence of a decrease in water transparency due to the high concentration of SM ( Figure 4A), which, in turn, led to a decrease in phytoplankton productivity because of a lack of light ( Figure 5A and Table 3). Light limitation of plankton primary production in the Neva Estuary during the periods of construction is confirmed by the negative correlation between PP and SM in those years ( Figure 3).
The composition of substances resuspended from the bottom can also affect the productivity of phytoplankton. For example, the use of active clay as a technique to combat the mass development of toxic algal species in coastal waters is now gaining popularity (de Magalhães et al., 2017;Yu et al., 2017). This method is quite simple and cheap (Song et al., 2021). The addition of clay particles to the water causes sedimentation of joint aggregates of clay and algae cells, resulting in a decrease in the number of algae in the water column (de Magalhães et al., 2019;Cao et al., 2022). Moreover, clay particles absorb nutrients from the water and thus become inaccessible to algae (de Magalhães et al., 2019). Suspended material during the period of mega-engineering projects in the Neva Estuary mainly contained fine sand, consisting of silica, and clayey Quaternary deposits, which make up the underlying bedrock of the bottom of the upper reach of the estuary (Ryabchuk et al., 2017). Sandy material was excavated from sand pits in the lower part of the estuary and delivered by barges to its upper part, where this material was poured out to create new territories, and the clayey deposits fell into the water column when new deep fairways were dug. It is likely that a significant positive relationship  between total phosphorus and the concentration of suspended matter ( Figure 4B) and a negative relationship with primary production (Figure 3) observed during construction periods is due to the appearance of a large amount of clay particles in estuary water. Salinity and water temperature also influenced the primary production of plankton in the estuary (Table 2, 3). However, the influence of these predictors was more likely to be associated with the action of other factors. Previously, it was shown that in the middle reaches of the Neva Estuary, primary production is higher than in the upper reaches (Golubkov et al., 2017). In the upper part of the estuary, the runoff is faster, and the residence time of water is therefore shorter than in the middle part of the estuary (Golubkov and Golubkov, 2020). Due to the short residence time of water, developing algae are washed out from the upper to the middle reaches of the estuary, where they accumulate in large biomasses. The salinity of water in the middle reaches is higher than in the upper reaches, and therefore, salinity appears in the multi-regression as a positive predictor ( Table 3). In addition, it has been shown that the biomass of algae is higher in less flowing and saline parts of the estuary . In addition, the upper shallow part of the estuary is usually warmer in summer compared to its middle part. As a result, water temperature is also a negative predictor in the multi-regression equation for the dependence of PP on environmental factors (Table 3).
Thus, our study showed that light was the main factor limiting the primary production of plankton in the Neva Estuary during the periods of mega-engineering construction. This factor limited the primary production despite the increase in the phosphorus concentration in the water during these periods. Similar limitation of primary production with an excess of nutrients by light was observed in some other estuaries under natural conditions as a results of shallowness and high resuspension of bottom sediments or overeutrophication (Yoshiyama and Sharp, 2006;Guinder et al., 2009;Horemans et al., 2020). The Neva Estuary is one of the most eutrophic water areas in the Baltic Sea; the primary production and biomass of phytoplankton are quite high at the peak of summer phytoplankton development (Golubkov et al., 2017;. However, under natural conditions, the light intensity is not a factor limiting the development of phytoplankton in the estuary. This factor became limiting during the construction of port infrastructure due to the resuspension of a large amount of suspended matter into the water column. An increase in the concentration of nutrients during these periods did not have a stimulating effect on the development of phytoplankton in the estuary under conditions of light limitation.

CONCLUSIONS
Large-scale engineering projects carried out in the Neva Estuary are a reflection of the general global trend towards the rapid development of coastal territories and port infrastructure at the mouths of large rivers. This trend is closely related to the steady growth of international trade and the large role of maritime water transport in it. An example of the construction of a new port infrastructure and the creation of new lands in the Neva Estuary shows that the implementation of such mega-projects is associated with an intense impact on the estuary ecosystem, including its primary productivity. The conducted study shows that construction work in the e st ua ry water a re a c au sed incre ased concentrations of suspended particulate matter in the estuary water, which significantly exceeded the natural background. This led, on the one hand, to a significant increase in the content of nutrients in the water, and, on the other hand, to a multiple increase in turbidity. As a result, there was a change in the factor limiting the productivity of phytoplankton. The multiple increase in nutrients turned out to be less important than expected. Against the background of high concentrations of phosphorus, the main factor limiting the productivity of phytoplankton became transparency of water. In the case of the Neva Estuary, this led to a decrease in plankton primary production. Taking into account plans for further development of port areas in the Neva Estuary and other coastal areas around the world, the influence of the large-scale engineering projects on the conditions for the development of phytoplankton may provide a new aspect of long-term regulation of algal booms and ecosystem functioning in the coastal and estuarine zones.