Study area
The study area comprises two large estuarine systems on the Brazilian equatorial margin with recognized importance for the region: the São Marcos estuarine complex (SMEC) and the São José complex (SJEC).
The Itaqui Port Complex (IPC) is located at the SMEC, which has been internationally recognized for its operational capacity. The IPC comprises the Itaqui Organized Terminal (IOT), Ponta da Espera Maritime Terminal (PEMT), and Alumar Private Use Terminal (PUT Alumar). Intensive movement of solid mineral bulk consisting of iron and manganese ores occurs in the PEMT (MTPA, 2018). PUT Alumar is responsible for handling bauxite and alumina loads (MTPA, 2018). The processing of bauxite generates as a by-product the so-called red mud consisting of iron, aluminum, sodium, and calcium oxides, as well as titanium dioxide and metals such as gallium, vanadium, and rare earth.
SJEC was marked by intense urbanization, mainly at points 09 and 08. However, there was little to no industrial activity, especially in the absence of port activity.
Due to the strong fluvial input in both estuarine complexes, the waters are characteristically turbid, with high concentrations of particulate matter, mainly clay and silt, originating from continental and mangrove areas (Silva, 2011).
The region's climate is tropical, marked by well-defined periods of rain and drought, with an average annual rainfall of 1,896 mm/year and an average temperature of 27 °C (https://www.climatempo.com.br). Thus, it is part of the context of humid tropical regions between 15°N and 15°S latitude, whose characteristics are high and constant precipitation (> 1,500 mm/year), temperatures > 20°, and low thermal variation (Nittrouer et al., 1995).
Because of these characteristics, the study area is characterized by the presence of structurally complex mangroves (Rebelo-Mochel 1997), the most abundant of which are Rhizophora mangle Linnaeus (1753), Avicennia germinans (L.) Stearn (1958), and Laguncularia racemosa (L.) CF Gaertn.
Sediment and macroalgae sampling
Sampling was conducted between August 2020 and February 2021 at nine sampling points: P01 to P07 in SMEC and P08 and P09 in SJEC (Fig. 1).
At all points, the in-situ water temperature, pH, salinity, total dissolved solids (TDS), and dissolved oxygen (DO) were determined using a multiparametric probe model HANNA HI 98194.
The intertidal surface sediments (top 5 cm) were collected using a plastic scoop at each sampling point. Three replicates of each sample were conducted in each plot of the study area (n = 54; considering two sample surveys). The sediment samples were dried in the laboratory at 60oC for 48 h and ground for further treatment.
Macroalgae were collected from the trunks and pneumatophores in the intertidal zone at each sampling point. Five pneumatophores were randomly harvested from a single sampling plot and cut in the mudline using a clipper (Billah et al. 2017). Therefore, 90 pneumatophores (45 pneumatophores × 2 samples) were harvested. In the laboratory, the samples were carefully washed to remove adhered particles and identified based on marine macroalgae literature (JHA et al. 2009).
Determination of metal concentrations in macroalgae
In the laboratory, the algae samples were lyophilized at -20 °C for 48 h, crushed, and homogenized for the additional procedures. Extraction of inorganic elements from macroalgae was performed using 0.750 g (dry weight) of the crushed sample. Initially, the lyophilized samples were placed in Teflon tubes (X-press) to which concentrated HNO3 (8 ml) and concentrated HF (2 mL) were added.
The extracts were allowed to rest overnight at room temperature and then placed in a microwave, model Mars X-press (CEM), for 40 min (15 min –ramp and 25 min–wait) at a temperature of 175 °C and power of 1600 W, adapted from EPA 3052. After cooling (30 min), H3BO3 (12 mL) was added to neutralize the HF, and then the tubes were taken back to the microwave for 25 min (15 min, ramp and 10 min, hold) at 170 °C.
The final cooled extract was filtered (Whatman No. 40) and measured at a final volume of 25 mL with 0.5 N HNO3 in a volumetric flask. The elements were determined using inductively coupled plasma optical emission spectrometry (IPC-OES 720). Precision measurements were obtained from the digestion of macroalgae and sediment in triplicate every 10 samples, and accuracy was determined from certified macroalgae and sediment standards Apple Leaves NIST SRM 1515. Recovery rates for metals in the macroalgae were 106% for Zn, 92.7% for Cu, 99.6% for Cr and 104.5% for Fe.
Elementary extraction and determination of metals in sediment
Initially, dried and powdered sediment (0.5 g) was placed in Teflon tubes (X-press) to which HNO3 (9 ml), HF (4 ml), and HCl (2 ml) were added. The extracts were allowed at rest overnight at room temperature and then placed in a microwave, model Mars X-press (CEM), for 40 min (15 min - Ramp and 25 min - waiting) at a temperature of 175 °C and power of 1600 W. These adjustments were adapted from a method described by the United States Environmental Protection Agency (EPA 3052). After cooling (30 min), H3BO3 (25 mL) was added to neutralize the HF, and the tubes returned to the microwave where they remained for 25 min (15 min, ramp and 10 min, hold) at 170 °C. After cooling for 30 min, the final extract was filtered (Whatman No. 40) and measured at a final volume of 50 mL with 0.5 N HNO3 in a volumetric flask. The determination of metals was performed later in the IPC-OES 720 equipment, and NIST 1646a was used as the standard. Recovery rates for metals in the sediment were 98% for Zn, 92% for Cu, 111% for Pb, 95% for Cr, 94% for Mn, and 95% for Fe. All analyses were performed in triplicate.
Granulometric analysis and Pejrup diagram
An interval-programmed Sald-3101 Shimadzu laser diffraction particle analyzer (Table 1) was used to determine the particle size using an aliquot of 2 g of the crude sample.
Table 1 Ranges of granulometric values according to the Wentworth classification (1922)
Phi
|
Wentworth Classification
|
Diameter
|
0
|
very coarse sand
|
2mm - 1mm
|
1
|
coarse sand
|
1mm - 500µm
|
2
|
medium sand
|
500µm - 250µm
|
3
|
thin sand
|
250 µm - 125µm
|
4
|
very thin sand
|
125µm-63µm
|
5
|
coarse silt
|
63µm - 31µm
|
6
|
medium silt
|
31µm - 16µm
|
7
|
fine silt
|
16µm - 8µm
|
8
|
very fine silt
|
8µm - 4µm
|
9
|
clay
|
< 4µm
|
The particle size distribution was treated using the SYSGRAN program (version 3.0), producing the following parameters: mean, median, standard deviation, asymmetry, kurtosis, normalized kurtosis, Wentworth classification, selection degree, and kurtosis classification. The statistical parameters of the particle size distributions were calculated from the particle size diameters (fi, mean, median, standard deviation, asymmetry, and kurtosis) using the equations proposed by Folk and Ward (1957).
The values obtained from the particle size analysis were used to construct a Pejrup diagram (1988). This diagram indicates the hydrodynamic conditions of the environment based on the sand, silt, and clay content present in the sampled sediment. A higher percentage of grains in the fine fraction indicates calmer hydrodynamic conditions.
The triangular diagram has four hydrodynamic sections (I to IV) that indicate hydrodynamic conditions: low (I), moderate (II), high (III), and very high (IV). In addition, constant sand content lines are suitable for textural classification. Thus, the triangle was divided into four texture classes, indicating a decrease in the sand content in the sediment. The sand percentages ranged from A to D, with A (90–100%), B (50–90%), C (10–50%), and D (0–100%).
Geo-accumulation index (Igeo), contamination factor (CF), enrichment factor (EF), and pollution load index (PLI)
The geochemical background values for Zn, Cu, Pb, Cr, Mn, and Fe proposed by De Paula Filho et al. (2015) in a study carried out in the Parnaíba Delta were used to calculate the indices. The São Luís Island and Parnaíba sedimentary basins are classified in the same morphological unit type, known as the Barreiras do Tertiary Group (CPRM 2000).
Fe was adopted as a reference element because of its similarity to other trace metals, uniform natural concentration, and its characteristic association with thin solid surfaces (Varol 2011). Thus, it was considered suitable as a normalizing element for the evaluated indices.
Geo-accumulation index (Igeo)
This index quantifies the metal pollution in aquatic sediments (Muller, 1979) using the following formula:
Igeo = Log2 CN/1,5BN
where CN is the concentration of the “N” element measured in the sediment, 1.5 is a correction factor for the background matrix that includes any possible variations in values due to lithogenic effects (Muller 1979), and BN is the background value of the element (Zhiyuan et al. 2011; Ozkan 2012). Below, seven descriptive classes were proposed for the values obtained from the geo-accumulation index (Table 2).
Table 2 Geo-Accumulation Index (Igeo) classes and values for sediment quality
Igeo Value
|
Igeo Class
|
Sediment Quality
|
0
|
0
|
Uncontaminated
|
0 – 1
|
1
|
Uncontaminated to Moderately Contaminated
|
1 – 2
|
2
|
Moderately Contaminated
|
2 – 3
|
3
|
Moderately to strong contaminated
|
3 – 4
|
4
|
Strongly contaminated
|
4 – 5
|
5
|
Strongly to Extremely contaminated
|
> 5
|
6
|
Extremely contaminated
|
Contamination factor (CF)
Hakanson (1980) proposed contamination factor (CF) as a marker to assess the level of metal concentrations in soils. CF corresponds to the ratio between the concentration of metals measured in the sediments and the background value of the metal at a given location (Turekian and Wedepohl, 1961).

Where for the “N” element, CFN is the ratio between the concentration of the element in the sediment (CN) and the background value (BN). Using the classification of Hakanson (1980), we obtain:
FC < 1: Low Contamination
1 ≤ FC < 3: Moderate Contamination
3 ≤ FC < 6: Considerable Contamination
FC ≥ 6: High Contamination
Enrichment factor (EF)
The enrichment factor (EF) corresponds to the standardization of a tested element with a reference element (Fe) that has low occurrence variability. The most common reference elements are Sc, Mn, Ti, and Al (Reimann and De Caritat, 2000; Sutherland, 2000).

Ms and Mb are the concentration of a given metal in the sediment and its background value, respectively. Fes and Feb are the concentration in the sediment and the background value of the reference metal Fe, respectively. The Enrichment Factor has seven categories (Taylor 1964) which are interpreted as:
EF < 1: No Mineral Enrichment;
EF = 1 – 3: Minimum Enrichment
EF = 3 – 5: Moderate Enrichment
EF = 5 – 10: Significant Enrichment;
EF = 10 – 25: Severe Enrichment;
EF = 25 – 50: Very Severe Enrichment;
EF > 50: Extreme enrichment.
EF values ranging between 0.5 and 1.5 suggest the contribution of the crust as a metal source, while values > 1.5 indicate anthropogenic influence (Nowrouzi and Pourkhabbaz 2014).
Pollution load index (PLI)
The Pollution Load Index (PLI) is obtained by multiplying the previously calculated contamination factors (CFs) and subsequent extraction of the root, according to the following equation:

Where: PLI < 1 No Polluted
PLI = 0 Normal levels
PLI > 1 Polluted Environment
Bioconcentration factor (BCF)
Bioconcentration is the accumulation of contaminants in aquatic organisms through non-dietary absorption pathways, for example, from the soluble phase. The bioconcentration factor has been used to assess the potential for metal bioaccumulation (Conti and Cecchetti 2003; Akcali and Kucuksezgin 2011). The calculation corresponds to the ratio of the metal concentration in the macroalgae to that in the sediment (BCFs) according to the following formula:

Threshold effect level (TEL), probable effect level (PEL), effects range low (ERL), and effects range medium (ERM)
TEL and PEL are interpretive toxicity criteria for marine and estuarine sediments that were established by Canadian legislation based on compiled biological and chemical data from laboratory studies, field measurements, and numerous models. TEL indicates the level below which there are no adverse effects on the biological community, and PEL is the level above which adverse effects are expected. The range above the TEL and below the PEL represents a possible adverse effect on the community (Hortellani et al. 2008; De Paula Filho et al. 2021).
ERL and ERM are American criteria established by Long et al. (1995) using chemical and biological data from field studies in marine and estuarine sediments. The ERL is the concentration limit below which it is rarely toxic. The ERM, in turn, is the limit above which there is a high likelihood of toxicity, and the range between the ERL and ERM indicates possible toxicity (Hortellani et al. 2008; Long et al. 1995). The ERL and ERM limits were adopted as reference values by resolutions 344/2004 and 454/2012 of the National Environmental Council (CONAMA, Brazil) for estuarine and marine sediments. The TEL, PEL, ERL, and ERM limits are not applicable to Al, Fe, and Mn.
Table 3 TEL, PEL, ERL and ERM limits for the metals Cd, Cr, Cu, Hg, Ni, Pb and Zn in estuarine and marine sediments. 1Hortellani et al, 2008; 2Long et al. 1995
Elements
|
1TEL
(μg g-1)
|
1PEL
(μg g-1)
|
2ERL
(μg g-1)
|
2ERM
(μg g-1)
|
Cd
|
0.7
|
4.21
|
1.2
|
9.6
|
Cr
|
52.3
|
160
|
81
|
370
|
Cu
|
18.7
|
108
|
34
|
270
|
Hg
|
0.13
|
0.70
|
0.15
|
0.71
|
Ni
|
15.9
|
42.8
|
20.9
|
51.6
|
Pb
|
30.2
|
112
|
46.7
|
218
|
Zn
|
124
|
271
|
150
|
410
|
Statistical Analysis
The Shapiro-Wilk test was used to verify the normality of all data. ANOVA (One-way) and Kruskal-Wallis tests were used to test for significant differences. Spearman's correlation coefficient was used to understand the paired relationships between the measured variables. For all tests, values of p < 0.05 and p < 0.01 were considered to indicate statistical significance. Analyses were performed using SPSS v26 software.