Phytoplankton Productivity and Hydrology in an Impacted Estuarine Complex in Northeastern Brazil ()
1. Introduction
Estuaries are dynamic environments found in the transition zone between the continent and the ocean. In addition, they have high variability in terms of physical and hydrochemical parameters [1] .
They are also unique systems, seen as aquatic-life nurseries, since they provide ideal protection and abundant food, which enable different organisms to inhabit them and reproduce [2] .
Certain aspects of estuarine ecosystems, such as their high productivity, ease of navigation and the shelter they provide against waves and currents, make them a desirable place for the establishment of residential, recreational, industrial and port endeavors [3] .
The economic importance of ports to international trade has significantly increased, mainly in developing countries. Thus, ports have become a vital component of the global economy [4] . However, several anthropogenic activities, such as vegetation suppression, earthwork, embankments, dredging [5] , and changes in sediment distribution patterns and coastal currents [6] , are linked to port construction and expansion processes.
Changes in estuarine geomorphology, sedimentology and hydrodynamics directly affect phytoplankton physiology [7] , biomass and community structure [8] , in addition to influencing local primary production [9] .
Primary production always depends on the synergistic effects of meteorological, physical, chemical and biological factors that limit or stimulate the activity of chlorophyll-containing planktonic organisms [10] [11] .
The determination of the rate of primary productivity based on 14C has been used to characterize the trophic state of water bodies [12] , since phytoplankton reacts rapidly to physical and chemical changes taking place in this medium, and since it is a sensitive and precise technique [13] .
The Suape estuarine system is one of the regions of greater ecological and economic relevance in Pernambuco State. In recent years, the estuary has faced a continuous degradation process due to the construction and expansion of an industrial and port complex. Several studies have focused on evaluating the quality of the water ( [14] [15] [16] [17] ) and investigating the local biodiversity ( [18] [19] [20] [21] ).
Accordingly, the aim of the current study was to help improve our understanding of the primary productivity capacity of the phytoplankton community (based on the 14C technique) and of various environmental variables by considering the constant anthropic changes that have been taking place in the Suape coastal area.
2. Materials and Methods
2.1. Study Site and Sampling Protocol
The Suape estuarine system is located approximately 50 km away from the capital city of Pernambuco State (08˚23'45''S 34˚58'04''W) (Figure 1). According to
Figure 1. Location of sampling points and Suape Industrial and Port Complex, Pernambuco, Brazil. Source: Produced by the authors, 2016.
the Köppen classification, the climate in the region is hot and humid (As’) with well-defined rainy (March to August) and dry seasons (September to February). Annual rainfall ranges from 1850 to 2364 mm, mean air temperature is 24˚C, relative humidity is higher than 80% and southeast winds prevail in the region [22] [23] .
With respect to oceanographic features, the region has semidiurnal tides, which are classified as mesotides based on their amplitude [23] .
Prior to the construction of Suape Industrial and Port Complex (SIPC), four rivers (Massangana, Tatuoca, Ipojuca and Merepe) converged towards Suape Bay, where their water was channeled through a continuous line of sandstone reefs, which ended in northern Santo Agostinho Cape. The construction of the port blocked the flow of the Ipojuca and Merepe rivers, which resulted in water retention and in the accumulation of substantial amounts of suspended material. A partial opening was made in the reef line near the mouth of the Ipojuca Riverto minimize the problem [24] . Years later, a second opening was made to allow vessels to have access to the inner portion of the port, dividing the original reef into two parts: the Muro Alto Reef (south) and the Suape Reef (north).
Water samples for nutrients (inorganic nitrogen compounds, phosphate, silicate), dissolved oxygen, chlorophyll and phytoplanktonic productivity were collected in November 2015, January 2016 and January 2017 (dry season), as well as in April, July and August 2016 (rainy season), during low spring tides.
Sampling was carried out in the internal (P1, P2, P3 and P4) and external portions of the port (P5, P6 and P7); P5 was under the influence of the plume of the Ipojuca and Merepe rivers.
2.2. Climatology and Environmental Variables
Pluviometric data were provided by the National Institute of Meteorology (INMET) [25] .
The hydrological parameters determined herein comprised bathymetry, based on echo sounding (digital echo sounder, mark LCD resolution: 0.1 m); water transparency, based on a Secchi disk; suspended particulate matter (SPM), based on the gravimetric volatilization method adopted by [26] ; temperature and salinity, from CTD measurements; dissolved oxygen (DO) concentration, based on the modified method by Winkler, described by [26] (precision of ±1.3 μM); oxygen saturation rate, based on the table by [27] ; dissolved nutrients such as ammonia (
), nitrite (
) and nitrate (
), based on the method described by [26] ; phosphate (
) and silicate (
), based on the method described by [28] . Dissolved inorganic nitrogen (DIN) was calculated based on the sum of
+
+
, dissolved inorganic silicate (DIS) was analyzed via
, and dissolved inorganic phosphorus (DIP) was analyzed via
. Stoichiometric phytoplankton limitations were evaluated with respect to the Redfield ratio, based on [29] and [30] . The precision was ±0.02 μmol for
, ±0.02 μmol for
, ±0.02 μmol for
, and 0.01 μmol for
. The accuracy was ±2% for
, ±3% for
and
, and ±5% for
.
2.3. Chlorophyll a and Phytoplanktonic Productivity
Total and fractionated chlorophyll a were analyzed based on the spectrophotometry method described by [31] ; the calculations with the equation from [32] , and the results were expressed in mg∙m−3.
Primary productivity was determined based on the 14C method described by [33] : 120 ml aliquots of each sample were inoculated with 1 ml sodium bicarbonate solution (NaHC14O3), which had an equivalent of 10 μCi of radioactive tracer, and subjected to a 4-hour incubation in situ, which was followed by vacuum filtration using 0.45-μm porosity and 47-mm diameter filters. The results were calculated based on [34] and expressed in mg∙C∙m−3∙h−1.
2.4. Statistical Analyses
The normality and homogeneity of the assessed variables were analyzed using the Komogorov-Smirnov and Levene tests, respectively. After nonnormality of the data was identified, the nonparametric Mann-Whitney test was used to check for seasonal differences (2 samples), whereas the nonparametric Kruskal-Wallis test was used to check for spatial differences (>3 samples); variables with p < 0.05 were significant. All calculations were performed using Statistic 7® software. The multivariate analysis applied to the investigated parameters was based on principal component analysis (PCA) using the program XLStat 2019®.
3. Results
3.1. Climatology and Environmental Variables
The total monthly rainfall ranged from 292.4 mm in July 2017 to 28.7 mm in January 2017. The study period did not show a significant difference from the historical values (Mann-Whitney; p = 0.19 and p = 0.93 for the dry and rainy seasons, respectively; α = 0.05) (Table 1), although the months of sampling (except January 2016) showed monthly rainfall below the average for the last 30 years (1985-2014) (Figure 2).
The local depth ranged from 21.0 m for P6 (July 2016) to 1.30 m for P3 (April 2016). There was significant spatial variation (Kruskal-Wallis test; α: 0.05; p = 0.00). The internal points were the shallowest, except for P4 (located at the entrance and exit of the port), whereas the external points presented the greatest depths (Table 1).
Water transparency ranged from 4.40 m (P4) to 0.80 m (P1) in January 2017. There were significant spatial (Kruskal-Wallis test; α: 0.05; p = 0.01) and seasonal (p = 0.04) differences, and the values recorded for points internal to the reef line were significantly lower (Mann-Whitney test; α: 0.05; p: 0.04) during the rainy season (1.5 m; Table 1).
The SPM concentration ranged from 49.00 mg∙L−1 in P1 (November 2015) to 23.30 mg∙L−1 for P7 (January 2017). There was a significant seasonal difference, and the highest SPM values were recorded during the rainy season (Mann-Whitney test; α: 0.05; p = 0.01) (Table 1).
The water temperature ranged from 29.44˚C for P7 (April 2016) to 26.55˚C for P2 (July 2016), and the thermal amplitude was 2.89˚C. There was a significant seasonal difference (Mann-Whitney test; α: 0.05; p = 0.04), with the lowest temperatures occurring in the rainy season (Table 1).
Salinity hada maximum value of 37.10 for P6 and P7 (August 2016) and a minimum value of 32.15 for P5 (November 2015, January 2016) and P6 (April 2016) (Figure 3(a)). There were significant spatial (Kruskal-Wallis test; α: 0.05; p = 0.02) and seasonal (Mann-Whitney test; α: 0.05; p = 0.04) differences, and the lowest salinity values were recorded for P5 during the rainy season (Table 1).
Dissolved oxygen and oxygen saturation had a maximum concentration of 5.19 ml∙L−1 and a maximum rate of 114.46%, respectively, for P7 (August 2016), as well as a minimum concentration of 2.60 ml∙L−1 and a minimum rate of 58.53%, respectively, for P1 (January 2016). There was significant spatial and seasonal variation in the DO concentration (Kruskal-Wallis test; α: 0.05; p = 0.00; Mann-Whitney test; α: 0.05; p = 0.02, respectively) and in the oxygen
Table 1. Abiotic and biological variables in rainy and dry seasons (Min = minimum, Max = maximum, Med = median, and SD = standard deviation) and the non-parametric tests (MW = Mann-Whitney; KW = Kruskal-Wallis) with p-values for the various treatments (seasonal and spatial) *p < 0.05 (SPM = Suspended particulate matter; Chla = Chlorophyll a total; Chla < 20 μm = Chlorophyll a fractional).
Figure 2. Rainfall data from the Experimental Station of Curado-Pernaambuco state between 2015 and 2017 and the average historical monthly values (1985-2014). Source: INMET.
Figure 3. Seasonal and spatial variation of abiotic variables: Salinity (a), Ammonia (b), Nitrite (c), Nitrate (d), Phosphate (e) and Silicate (f), during the study period in the port region of Suape, Pernambuco, Brazil.
saturation rate (p = 0.00; p = 0.03) (Table 1). Both parameters had their highest values at points external to the reef line and during the rainy season.
Overall, the nitrogenous compound (
,
and
) concentrations were low.
concentrations had the highest values, which ranged from 4.40 μmol∙L−1 to 0.01 μmol∙L−1 (Figure 3(d)), in comparison with those of
and
, there was no significant spatial or seasonal difference.
concentrations ranged from 0.84 μmol∙L−1 to 0.01 μmol∙L−1 (Figure 3(b)), whereas
concentrations ranged from 0.14 μmol∙L−1 to 0.01 μmol∙L−1 (Figure 3(c)). The concentrations of both had significant seasonal variation during the dry (Mann-Whitney test; α: 0.05; p = 0.01) and rainy seasons (Mann-Whitney test; α: 0.05; p = 0.00), respectively (Table 1).
The
concentrations ranged from 0.35 μmol∙L−1 (P5) to 0.04 μmol∙L−1 (P3) in April 2016 (Figure 3(e)). There was only significant spatial variation (Kruskal-Wallis test; α: 0.05; p = 0.00), especially for P5 (Table 1).
ranged from 24.20 μmol∙L−1 for P1 (July 2016) to 1.90 μmol∙L−1 for P6 (January 2016) (Figure 3(f)). There were spatial (Kruskal-Wallis test; α: 0.05; p = 0.00) and seasonal variations (Mann-Whitney test; α: 0.05; p = 0.02) (Table 1). The highest silicate concentrations were recorded at points internal to the reef line and during the rainy season.
The DIN:DIP ratio ranged from 0.13 for P4 (January 2017) to 28.00 for P1 (July 2016); there was no spatial (Kruskal-Wallis test; α: 0.05; p = 0.08) or seasonal (Mann-Whitney test; α: 0.05; p = 0.10) variation. The DIN:DIS ratio ranged from 0.00 for P6 (January 2017) to 0.86 for P7 (January 2016); there was no spatial (Kruskal-Wallis test; α: 0.05; p = 0.21) or seasonal (Mann-Whitney test; α: 0.05; p = 0.30) variation. The DIS:DIP ratio ranged from 13.18 for P6 (November 2015) to 268.88 for P1 (July 2016); there were significant spatial (Kruskal-Wallis test; α: 0.05; p = 0.01) and seasonal (Mann-Whitney test; α: 0.05; p = 0.00) variations (Table 1). Thus, nitrogen was the limiting element regardless of spatial (internal and external) and seasonal variations.
3.2. Chlorophyll a and Phytoplanktonic Productivity
Total chlorophyll-a concentrations ranged from 0.02 mg∙m−3 for P5 and P6 (January 2017) to 2.45 mg∙m−3 for P3 (November 2015) (Figure 4(a)), while
Figure 4. Seasonal and spatial variation of biological variables: Chlorophyllaa total (a), Chlorophyllaa < 20 µm (b) and Primary productivity (c), during the study period in the port region of Suape, Pernambuco, Brazil.
fractionated chlorophyll ranged from 0.02 mg∙m−3 (for several months and sampled points) to 2.17 mg∙m−3 for P3 (November 2015) (Figure 4(b)). There were no statistically significant variations in the total and fractionated chlorophyll-a concentrations (seasonal; Mann-Whitney test; α: 0.05; p = 0.71; spatial; Kruskal-Wallis test; α: 0.05; p = 0.90; Table 1), although they were slightly higher in the dry season than in the rainy season and at the innermost points of the port compared to those farther outside the port. The fraction of plankton < 20 μm, corresponding to the peak and nanophytoplankton, contributed the most to the increase in the phytoplankton biomass (88.57%).
Primary productivity ranged from 0.34 mg∙C∙m−3∙h−1 for P7 (January 2017) to 4.32 mg∙C∙m−3∙h−1 for P3 (November 2015); the median value was 2.01 mg∙C∙m−3∙h−1 (Figure 4). The phytoplankton assimilation rate ranged from 1.14 mg∙C∙m−3∙h−1 for P7 (January 2017) to 23.93 mg∙C∙m−3∙h−1 for P3 (November 2015); the median value was 4.99 mg∙C∙m−3∙h−1 (Figure 4(c)). There was no significant variation in any of the parameters, although they presented a pattern similar to that of chlorophyll-a (Table 1).
3.3. Statistical Analyses
According to the PCA, the two components explained 52.32% of the data variance (Figure 5). For the first component, 28.66% of the variance was directly correlated with total chlorophyll, chlorophyll < 20 µm, primary productivity and temperature, indicating that the increase in temperature in the dry period promoted an increase in the metabolism of phytoplanktonic organisms. For the second component (23.66%), dissolved oxygen showed a direct correlation with salinity and an inverse correlation with silicate, indicating that in the area external to Suape port, hydrodynamics keep the water more oxygenated (Table 2).
4. Discussion
The incidence of rainfall is among the main factors responsible for modulating and conditioning changes in the hydrological and geomorphological parameters of estuaries and other coastal ecosystems [10] .
The rainfall rate recorded in the coastal zone of Suape during the months of sampling (except January 2016) was lower than that of the historical average over the last 30 years (1985 to 2014). Although precipitation was lower than expected, it was sufficient to increase nitrite, silicate, and SPM concentrations as well as decrease water temperature and transparency within the SIPC area.
Given the relevant morphological changes that have taken place in the area, saltwater intrusion has also been a determining factor for changes found in some of the hydrological parameters. Among them, one finds an increase in salinity, a dilution of nutrient salts and, consequently, a decrease in chlorophyll-a concentrations and primary productivity. Studies conducted in the past have reported that changes in the coastal dynamics in the Suape port area affected nutrient and algal biomass contents [15] , as well as the species composition of phytoplankton [35] and zooplankton populations [36] .
Figure 5. Principal Component Analysis (PCA) of the abiotic and biological variables of the port region of Suape, Pernambuco, Brazil. Analyzed for the two principal component. Subtitles: CHLA T = Chlorophylla a total; CHLA F = Chlorophylla a < 20 µm; DIN = Ammonia + Nitrite + Nitrate; DO = Oxygen Dissolved; PO4 = Phosphate; PP = Primary productivity; SIO2 = Silicate; TEMP = Temperature.
Table 2. Correlation analysis with factor loading of the variables (DO = Oxygen dissolved; TEMP= Temperature; SAL = Salinity; PO4 = Phosphate; SIO2 = Silicate; CHL T = Chlorophylla a total; CHL F = Chlorophylla a < 20 µm; PP = Primary productivity; DIN = Ammonia + Nitrite + Nitrate. The bold values indicate significant load factors.
Water transparency in estuarine ecosystems may vary depending on rainfall intensity and duration. Such transparency follows an increasing gradient towards the portion of the ecosystem with greater marine influence, in addition to its being a limiting factor for the development of phytoplanktonic organisms [37] .
In Suape, water transparency followed the spatial and seasonal pattern of estuarine environments in that the waters were less transparent in the inner portion of the bay due to increased concentrations of the allochthonous material arriving in the region during the rainy season. These trends in water transparency and particulate matter have been reported in other estuaries in Pernambuco State by [38] in the estuary of the Sirinhaém River, by [39] in the coastal zone of Ipojuca, and by [40] in the estuarine ecosystem of Barra das Jangadas.
The water temperature was always high and presented a low thermal amplitude; the lowest values were recorded during the rainy season. It is noteworthy that such low thermal variations occurred gradually throughout the months, which allowed planktonic organisms to adapt to new conditions [41] .
Low thermal amplitude and seasonal differences were also reported in other coastal regions of Pernambuco State, such as in the inner portion of Suape Port [42] , the Porto de Galinhas beach [43] and the estuary of Maracaípe River [39] . The aforementioned studies associated this pattern with efficient water renewal resulting from tidal actions.
Salinity is another important parameter for the coastal zone; its dynamics are controlled by the influx of freshwater from the rivers, by tides and by meteorological and geomorphological factors. In addition, salinity is able to strongly influence the distribution and succession of phytoplankton species [44] .
According to previous data recorded by [45] , the local ecosystem was divided into three areas before the implementation of the SIPC: the first area covered Suape Bay and varied from euhaline to polyhaline; the second area comprised the Massangana and Tatuoca rivers and was characterized as an estuarine zone with salinity ranging from polyhaline to mesohaline; and the third area comprised the estuary of the Ipojuca River, with salinity ranging from polyhaline to limnetic.
The opening made on the reef allowed access to the internal port and the loss of communication between the Ipojuca and Merepe rivers, causing changes in the current salinity regime in Suape Bay and in the mouth of the Massangana River to euhaline, a fact that was also reported by [18] [46] and [47] . According to [36] , Suape Bay had a typical estuarine zooplankton community prior to the implementation of the port; however, this community was replaced by a group of coastal-water copepod species, which accounts for approximately 73% of the total zooplankton population.
Contrary to our expectations, the lowest salinity value was recorded in the external portion of the port (P5 and P6), possibly due to the inflow of less saline water from the Ipojuca and Merepe rivers. It is worth emphasizing that [24] and [35] found over a 2-hour delay in tidal dynamics at the mouths of the Ipojuca and Merepe estuaries due to the previously mentioned changes in coastal dynamics and morphology; such a delay changed the saline regime in the region, which currently oscillates between polyhaline and euhaline. According to [35] , increased salinity levels have significantly increased the incidence of marine dinoflagellate species (from 3.33% to 27.41%) in the lower estuary of the Ipojuca River.
With respect to gases dissolved in water, oxygen is one of the most important in the dynamics of aquatic environments, and its concentration changes depending on biochemical and biological processes, such as respiration and production. Oxygen is also an important indicator of the quality of water bodies [48] .
Based on the classification system developed by [49] , the oxygen saturation rate in the waters adjacent to the SIPC varied from saturated to supersaturated, thus indicating that this environment receives low organic contributions and substantial amounts of seawater. This outcome corroborated results found by [50] on the northern coast of Bahia State, by [51] in Maracajaú (RN), and by [43] in Porto de Galinhas.
Nutrients such as nitrogen and phosphates play a key role in aquatic environments because, along with light, they are the main limiting factors for primary production and, consequently, they affect the entire trophic web [9] .
Nitrogenous compound concentrations were particularly low in Suape; nitrogen was the element limiting the development of phytoplanktonic organisms.
had the lowest concentrations, followed by
and
. [43] and [52] conducted studies in Porto de Galinhas and Serrambi, respectively, and they reported
values (0.01 - 0.17 μmol∙L−1 and DIN - 0.13 μmol∙L−1, respectively) very close to those found in the aforementioned study. The authors associated such outcomes with the good environmental quality of the investigated ecosystems.
had significant seasonal variation; the highest values were recorded during the dry season, when there was lower organic matter dilution. P1 had the highest spatial concentration of
, probably due to the influence of Barbosa stream, which crosses a village and disembogues shortly before the mouth of the Massangana River. Among the forms of dissolved organic nitrogen,
is preferred by phytoplanktonic organisms, which consume it rapidly, thus reducing its levels in aquatic ecosystems.
The maximum
value recorded in the present study was 4.40 μmol∙L−1, which is unlike what was observed by [35] in the estuary of the Ipojuca River, The Brazilian environmental legislation [53] indicates a limit value of 0.4 mg∙L−1 (23.5 μmol∙L−1).
PO4- only had spatial variation; the highest concentrations were found for P5, and they were possibly influenced by the Ipojuca and Merepe rivers. The values observed here were close to the values specified by the Brazilian environmental resolution [53] (Table 1).
According to [40] , the highest
concentrations were found on the continental shelf of Pernambuco State, which is strongly influenced by the Jaboatão River. On the other hand, [54] showed that the plumes of the Capibaribe and Beberibe rivers contributed to the coastal environment by increasing nutrient salts and chlorophyll-a contents in the surrounding seawater.
Silicate had the highest concentrations in the waters adjacent to the SIPC in comparison to the concentrations of other inorganic nutrients; however, since it is an estuarine region, silicate values were classified as low, which can be explained by the reduced silicate transport in the Massangana and Tatuoca rivers (coastal and lowland rivers with small hydrographic basins) and by the loss of silicate contribution from the Ipojuca River.
With respect to spatial-temporal variation, silicate had higher concentrations in the inner points of Suape and during the rainy period. The same pattern was found in the estuarine environment of the Ariquindá River [55] and in the estuary of the Timbó River [56] , which was explained by increased continental contributions and influenced by the incidence of rainfall.
Phytoplankton found in tropical coastal regions are susceptible to environmental changes to which they respond rapidly due to their short life cycle [10] . Consequently, they are globally used as bioindicators of water quality and trophic status, in addition to being valuable tools for the environmental management of coastal zones [41] .
In Suape, the chlorophyll-a concentration and phytoplankton productivity were moderately higher at inner points of the port during the dry period, which conformed to increases in the photic layer,
concentration and water temperature. According to [57] , the water temperature increase plays a fundamental role in changing photosynthetic and microalgae respiration rates.
The spatial-temporal pattern of chlorophyll-a observed in the present study was also observed by [58] in the estuarine zone of the Ilhetas and Mamucabas rivers, by [59] in Golfão Maranhese (Maranhão state, Brazil), and by [40] in the estuary of Barras das Jangadas. These authors attributed the observed outcome to an increase in suspended material concentrations and a decrease in the photic layer. These phenomena often occur during the rainy season when light is a limiting factor for phytoplankton development.
Chlorophyll-a < 20 μm was the one that most contributed to the phytoplankton biomass of the system (88.57%). A similar result was observed by [42] in Suape Bay and by [55] in the estuarine zone of the Ariquindá River. These outcomes can be explained by the better photosynthetic efficiency of this fraction of the phytoplankton population in regions facing nutrient scarcity and hydrological parameter variability, as well as by the ability of this fraction to tolerate variations in environmental conditions.
Based on the classification by [60] , the waters adjacent to the SIPC were typical of oligotrophic environments.
Other researchers also observed low chlorophyll-a levels (values < 6.00 mg∙m−3) along the Brazilian coast, among them were the following: [50] in Todosos Santos Bay on the northern coast of Bahia State, [61] in Bragantina Peninsula (Pará state, Brazil), and [62] in the estuary of Paraíba do Sul River (Rio de Janeiro state, Brazil).
Opposite results were recorded by [63] in the estuary of the Jaboatão River, by [64] in the port region of Recife, and by [65] in the Pina basin (values > 100 mg∙m−3). According to the aforementioned authors, these values indicated strong eutrophication, which was mainly influenced by untreated domestic and industrial effluents released in these ecosystems.
With respect to primary productivity and considering the classification of [66] , it is possible to state that the values recorded in Suape were typical of unpolluted coastal environments, a fact that was confirmed by the low nutrient salt and chlorophyll-a concentrations that is found in such environments.
The low primary phytoplankton productivity values observed in the present study were similar to those recorded by [67] in Tamandaré Bay, by [68] in a profile perpendicular to the coast at Piedade beach, by [51] in Maracajaú (Rio Grande do Norte state, Brazil), and by [43] in the reef environments of Porto de Galinhas. On the other hand, the present results were different from the high phytoplankton productivity levels observed by [69] in Pina basin, by [38] upstream in the Sirinhaém River, and by [70] in the estuary of the Timbó River, which are highlyeutrophic ecosystems.
The phytoplankton assimilation rate corresponds to the ratio between productivity and chlorophyll-a values. Based on [71] , it is possible to say that the water adjacent to the SIPC was mesotrophic, although it had a strong tendency towards oligotrophy.
The comparison of previous and current data allowed us to identify significant loss of phytoplankton production efficiency in the estuarine system of Suape. Before implementation of the SIPC, the total phytoplankton density in the estuary of the Ipojuca River varied from 416.00 to 5,748.00 cells∙L−1 [72] . After its implementation, cell density decreased by 70% [35] . With respect to zooplankton, [73] indicated an expressed cellular biovolume decrease in comparison to that recorded before such development occurred [74] .
[75] reported a significant decrease in the number of ichthyoplankton families and species in comparison to previous information recorded by [76] . It is worth emphasizing that Suape Bay was one of the largest fish larvae nurseries and distributing centers in the coastal region of Pernambuco State before the port was built [77] .
The local decrease in fishing resources has caused severe ecological damage. Shark species that use the Suape estuarine complex as a source of food and protection now swim towards highly urbanized adjacent coastal areas where they cause serious accidents involving human beings [20] [21] .
5. Conclusion
In light of the foregoing results, it is possible to conclude that decreased rainfall indices and constant changes in the coastal morphology and dynamics of the estuarine ecosystem investigated herein were the strongest factors affecting the environmental parameters analyzed in the current study. The loss of contribution from the Ipojuca River and the strong saltwater intrusion have changed the salinity regime of Suape Bay (currently classified as euhaline), leading to nutrient dilution and reduced primary productivity. Picoplankton and nanoplankton were the fractions of phytoplankton that most contributed to the system (88.6%), reflected the low nutrient content and indicated that the area is free from eutrophication processes. However, the decreased fertility of the waters adjacent to the SIPC has had negative impacts on other links in the trophic web.
Acknowledgements
The authors are grateful to LABFITO, LOQUIM-UFPE and LARADOB-USP for the analysis and data availability. C. Noriega acknowledges the National Council for Scientific and Technological Development (CNPq) for financial support for the project (No. 465634/2014-1; INCT-AmbiTropic; Process No. 370821/2019-0).
Financial Support
The current research is part of a larger project entitled “Studies on Dissolved Inorganic Carbon Cycle in Northeastern and Coastal Areas in Brazil and its Relation to Ocean Acidification Processes (DICAM)” funded by the Coordination for the Improvement of Higher Education Personnel (CAPES).