# Structure and Composition of Mangrove Associations in Tubli Bay of Bahrain as Affected by Municipal Wastewater Discharge and Anthropogenic Sedimentation

Authors: Kholoud Abou Seedo; Mohammad S. Abido; Ahmed Salih; Asma Abahussain
Journal: International Journal of Biodiversity
 (2017)

---

## Abstract

The effects of municipal wastewater discharge and anthropogenic sedimentation on the structure and composition of gray mangrove (Avicennia marina (Forsk.) Vierh.) communities along Tubli Bay coastlines in Bahrain were investigated. Growth and regeneration of mangrove were measured, and its community was characterized. Sediment profile was analyzed for texture, pH, and salinity. Mangrove area covered by sand depositions was measured using Google Earth Pro. ANOVA and regression tests were employed in the analysis of the data. Results indicated that mangrove overwhelmingly dominated plant community in the study area, which was zoned by a community of other salt-tolerant species. Three main habitats exist in the study area with high similarity in their floristic composition. Species richness and the number of habitats were low due to the aridity and high sediment salinity. The dilution effect of the secondary treated wastewater had a favorable effect on height and diameters of mangrove trees. However, no differences were observed in leaf area index, basal area, and density of mangrove. The long-term accumulation of anthropogenic sedimentation had a detrimental effect on the mangrove community, expressed in swath death of mangrove trees due to root burials and formation of high topography within the community boundaries.


---

## Body

## 1. Introduction

Gray mangrove (Avicennia marina (Forsk.) Vierh.) grows naturally in the intertidal zones of the southern coasts of the Arabian Peninsula forming dwarf stands intermingled with salt-tolerant scrub communities [1–3]. Remnant stands of the species are found along the coastlines of Tubli Bay in Bahrain, occupying 6% of the Bay habitats and 1% of the total intertidal and shallow water habitat in the country [3]. These stands are fragmented and subject to various forms of anthropogenic pressures, including the discharge of secondary treated wastewater and the sediment effluents of nearby sand washing plants [4, 5].Growth and community development of mangroves vary between regions due to natural, anthropogenic, and historic factors [6–12]. While climatic factors are decisive in controlling the worldwide distribution of mangroves, local site factors, including topography, soil properties, and tide fluctuations, are significant elements affecting variations in structure and composition of mangrove communities [13–16]. In dry environments, aridity and high salinity, among other site factors, are the most critical factors that control growth, structure, and composition of mangroves [17, 18]. Studies indicated that gray mangrove survives a broad range of salinity and pH conditions. However, best growth is reported to be in the range of 5–75% seawater concentrations [19].A. marina may reach a height of 1.6 m and 2.0 m when grows in pure and mixed stands consecutively in some of the hypersalinity spots along the Red Sea coasts of Egypt [20]. Nevertheless, increasing the salinity of mangrove media is inversely correlated with the plant water relationship, whereas inflow of freshwater modifies salinity gradients, thus yielding a positive effect on growth rates and species diversity of mangrove communities [21–24]. Often, species richness of mangrove communities is linked to the degree of salinity. High species richness and diversity were associated with moderate salinity and high moisture content of soil [7, 25].The effect of anthropogenic sedimentation on growth and survival of mangroves is variable and depends on the magnitude, type of sediment, and the degree of species tolerance [22]. Nevertheless, the high rate of sediment accumulation reduces growth, buries the roots, and kills the seedlings [22, 26]. The mechanism of reducing the growth of mangrove and subsequent death of its trees was associated with the smothering effect of sediment on pneumatophores [22]. On the other hand, mangrove architecture and health were related to sediment mud content and deposition patterns of sediment [23].Characterization of the structure and composition of mangrove communities as affected by modified site conditions is essential for enhancing the conservation effectiveness and implementation of the management plans in light of increasing threats to these forests [27–29]. In the Gulf region, the scientific knowledge on the effect of treated wastewater disposal and sedimentation burials on mangrove is scarce if not lacking. This study aims to investigate the effects of urban wastewater discharge and anthropogenic sedimentation on the structure and composition of mangrove communities along the coastlines of Tubli Bay in Bahrain.


## 2. Materials and Methods

### 2.1. Study Area

Tubli Bay is an enclosed sea area of 10 Km2 in size with an average depth of 2.5 m. It is bounded by urban and agricultural setting as well as causeways which limit water exchange with the open sea. The average salinity of the Bay is 17.5 dsm−1 compared to 59.7 dsm−1 of the outer sea. Tide occurs twice a day and ranges in heights between 0.20 and 2.62 m. The slope of the intertidal zones of the Bay is 2.2%.Three mangrove sites were selected along the intertidal zones of Tubli Bay coastlines based on accessibility (Figure1) [30]. The first site (Tubli) is directly affected by the discharges of approximately 1.2 m3 Sec−1 of secondary treated municipal wastewater [5, 6]. The second site (Sitra-1) had experienced a sedimentation discharge from a nearby sand washing plant for the past 15 years, whereas the third site (Sitra-2) is considered pollution-free [4]. The climate of the study sites is arid, mild in winters, and very hot and humid in summers. Rainfall is scanty and sporadic averaging 80 mm y−1. The yearly average temperature is 26.8°C and the annual evaporation rate reaches 2000 mm y−1 (Figure 2) [31].Figure 1
Avicennia marina study sites.Figure 2
Monthly average rainfall (mm) and temperature (°C).


### 2.2. Methods

#### 2.2.1. Vegetation Analysis

Three replicated plots of 10 m × 10 m were established randomly in the seaward landward direction at each study site. Within these plots, heights and diameters of trees (30 cm above ground level) were measured with a graduated wooden staff and diameter tape accordingly [32, 33]. Heights and diameters of pneumatophores as well as seedling heights were measured in a 1 × 1 m subplots located in the middle of each larger plot. Leaf area index (LAI) was estimated in each plot using CI-110 Plant Canopy Imager. Species cover was measured using a 60-meter long transects stretched along the middle of each larger plot with a sea to landward direction; subsequently, species coverage, density, frequency, and importance value index (IVI) were computed for each plot [32, 33]. Species similarity and diversity were calculated using Sorensen-Dice coefficient [34] and Simpson’s index [35], respectively. Life forms and habitats of encountered species were categorized according to [36, 37], respectively. Identification of associated species was based on [38].


#### 2.2.2. Media Analysis

A sediment core (10 cm × 2 cm) was taken from the center of each plot in the three sites at three depths (0–5, 5–15, and 15–30 cm). After processing the samples, soil texture by depth was determined using a soil texture triangle calculator [39]. Sediment electrical conductivity (EC) and pH were measured according to [40]. For each plot in the three sites, water samples were taken, 30–60 minutes of low tide time. A subsample, 50 ml of seawater sample, was used to measure EC and pH using conductivity and pH meters, respectively.


#### 2.2.3. Statistical Analysis

Measured parameters were analyzed atp≤0.05 using ANOVA and regression analysis employing the JMP11 statistical software [41]. Accumulative sediment spread effects were measured as a percentage of mangrove area covered by sand washing plant depositions using Google Earth Pro and ground truthing. Fieldwork was conducted at times of low tide during the period of September 2014 to June 2015.


## 2.1. Study Area

Tubli Bay is an enclosed sea area of 10 Km2 in size with an average depth of 2.5 m. It is bounded by urban and agricultural setting as well as causeways which limit water exchange with the open sea. The average salinity of the Bay is 17.5 dsm−1 compared to 59.7 dsm−1 of the outer sea. Tide occurs twice a day and ranges in heights between 0.20 and 2.62 m. The slope of the intertidal zones of the Bay is 2.2%.Three mangrove sites were selected along the intertidal zones of Tubli Bay coastlines based on accessibility (Figure1) [30]. The first site (Tubli) is directly affected by the discharges of approximately 1.2 m3 Sec−1 of secondary treated municipal wastewater [5, 6]. The second site (Sitra-1) had experienced a sedimentation discharge from a nearby sand washing plant for the past 15 years, whereas the third site (Sitra-2) is considered pollution-free [4]. The climate of the study sites is arid, mild in winters, and very hot and humid in summers. Rainfall is scanty and sporadic averaging 80 mm y−1. The yearly average temperature is 26.8°C and the annual evaporation rate reaches 2000 mm y−1 (Figure 2) [31].Figure 1
Avicennia marina study sites.Figure 2
Monthly average rainfall (mm) and temperature (°C).


## 2.2. Methods

### 2.2.1. Vegetation Analysis

Three replicated plots of 10 m × 10 m were established randomly in the seaward landward direction at each study site. Within these plots, heights and diameters of trees (30 cm above ground level) were measured with a graduated wooden staff and diameter tape accordingly [32, 33]. Heights and diameters of pneumatophores as well as seedling heights were measured in a 1 × 1 m subplots located in the middle of each larger plot. Leaf area index (LAI) was estimated in each plot using CI-110 Plant Canopy Imager. Species cover was measured using a 60-meter long transects stretched along the middle of each larger plot with a sea to landward direction; subsequently, species coverage, density, frequency, and importance value index (IVI) were computed for each plot [32, 33]. Species similarity and diversity were calculated using Sorensen-Dice coefficient [34] and Simpson’s index [35], respectively. Life forms and habitats of encountered species were categorized according to [36, 37], respectively. Identification of associated species was based on [38].


### 2.2.2. Media Analysis

A sediment core (10 cm × 2 cm) was taken from the center of each plot in the three sites at three depths (0–5, 5–15, and 15–30 cm). After processing the samples, soil texture by depth was determined using a soil texture triangle calculator [39]. Sediment electrical conductivity (EC) and pH were measured according to [40]. For each plot in the three sites, water samples were taken, 30–60 minutes of low tide time. A subsample, 50 ml of seawater sample, was used to measure EC and pH using conductivity and pH meters, respectively.


### 2.2.3. Statistical Analysis

Measured parameters were analyzed atp≤0.05 using ANOVA and regression analysis employing the JMP11 statistical software [41]. Accumulative sediment spread effects were measured as a percentage of mangrove area covered by sand washing plant depositions using Google Earth Pro and ground truthing. Fieldwork was conducted at times of low tide during the period of September 2014 to June 2015.


## 2.2.1. Vegetation Analysis

Three replicated plots of 10 m × 10 m were established randomly in the seaward landward direction at each study site. Within these plots, heights and diameters of trees (30 cm above ground level) were measured with a graduated wooden staff and diameter tape accordingly [32, 33]. Heights and diameters of pneumatophores as well as seedling heights were measured in a 1 × 1 m subplots located in the middle of each larger plot. Leaf area index (LAI) was estimated in each plot using CI-110 Plant Canopy Imager. Species cover was measured using a 60-meter long transects stretched along the middle of each larger plot with a sea to landward direction; subsequently, species coverage, density, frequency, and importance value index (IVI) were computed for each plot [32, 33]. Species similarity and diversity were calculated using Sorensen-Dice coefficient [34] and Simpson’s index [35], respectively. Life forms and habitats of encountered species were categorized according to [36, 37], respectively. Identification of associated species was based on [38].


## 2.2.2. Media Analysis

A sediment core (10 cm × 2 cm) was taken from the center of each plot in the three sites at three depths (0–5, 5–15, and 15–30 cm). After processing the samples, soil texture by depth was determined using a soil texture triangle calculator [39]. Sediment electrical conductivity (EC) and pH were measured according to [40]. For each plot in the three sites, water samples were taken, 30–60 minutes of low tide time. A subsample, 50 ml of seawater sample, was used to measure EC and pH using conductivity and pH meters, respectively.


## 2.2.3. Statistical Analysis

Measured parameters were analyzed atp≤0.05 using ANOVA and regression analysis employing the JMP11 statistical software [41]. Accumulative sediment spread effects were measured as a percentage of mangrove area covered by sand washing plant depositions using Google Earth Pro and ground truthing. Fieldwork was conducted at times of low tide during the period of September 2014 to June 2015.


## 3. Results and Discussion

### 3.1. Community Structure

The results of the study revealed that the heights of mangroves in the study sites ranged between 1 and 5.5 m with an average of 2.7 m. Trees were mainly of multiple trunks, ranging from 2 to 5 stems per plant. Diameters of trees ranged in size from 2.2 to 12.5 cm with a mean of 5.3 cm. Mean density and basal area of mangrove trees were 4576.8 tree ha−1 and 11.4 m2 ha−1, respectively. The heights and diameters of trees in Tubli site were 3.25 m and 5.82 cm, respectively. Tree heights were significantly higher (p≤0.05) in the Tubli stand as compared to Sitra (2) but not Sitra (1) (Table 1). Moreover, no significant differences in the height of trees were observed between the Sitra (1) and Sitra (2) stands. High growth parameters could be attributed to the dilution effect of discharged wastewater where measured EC of seawater at Tubli site was 13.2 ± 2 compared to 18.8 ± 2 dS m−1 and 17.1 ± 2 dS m−1 in Sitra-1 and Sitra-2 sites, respectively. However, the insignificance in the height of trees between Tubli site and Sitra (1) could be related to the larger density of trees in the second site, though not significant per se between the two sites. Also, tree diameters were significantly higher (p≤0.05) in the Tubli stand as compared to the Sitra (1) and Sitra (2) stands with no observed significant differences between diameters of trees in the last two sites. On the other hand, no differences in density, basal area, and LAI were observed between the three sites (see Table 1).Table 1
Avicennia marinaattributes in the study sites∗.
Stand
Trees attributes
Stand attributes
Pneumatophores
Seedling
Height (m)
Diameter (cm)
Basal area (m2ha−1)
Leaf area index (m2m−2)
Density (tree ha−1)
Height (cm)
Diameter (cm)
Number m−2
Height (cm)
Number m−2
Tubli
3.25
±
0.3
a
5.82
±
0.1
a
12.99
±
3.1
a
1.28
±
0.05
a
4565.8
±
770
a
15.34
±
0.9
a
1.01
±
0.01
a
106.33
±
12.2
b
4.50
±
2.7
b
1.0
±
1.7
b
Sitra-1
2.58
±
0.3
a
b
5.07
±
0.2
b
11.95
±
1.2
a
1.18
±
0.0
4
a
5085.7
±
79
7
a
6.29
±
0.9
b
0.50
±
0.0
1
b
202.19
±
42.6
a
3.84
±
1.2
b
9.2
±
3.8
a
Sitra-2
2.15
±
0.3
b
5.04
±
0.2
b
9.15
±
1.2
a
1.24
±
0.05
a
4079.1
±
56
5
a
15.50
±
0.7
a
1.01
±
0.01
a
140.33
±
15.1
a
b
15.10
±
1.3
a
5.5
±
1.04
a
b
Mean
2.7
±
0.2
5.3
±
0.1
11.4
±
1.2
1.2
±
0.03
4576.85
±
406
12.4
±
0.97
0.83
±
0.05
149.6
±
16.9
8.4
±
1.96
6.9
±
2.05
p value
0.0200
0.0003
0.4463
0.3308
0.6178
<0.0001
<0.0001
0.0570
0.0012
0.0439
∗Levels not connected with the same letter are significantly different (p≤0.05).In the meantime, the sediment salinity level of soil profiles varied among the sites and ranged from a low value of 10.3 dS m−1 in Tubli to 23.2 dS m−1 in Sitra-2. The average salinity levels were 13.72 dS m−1, 19.57 dS m−1, and 22.27 dS m−1 in Tubli, Sitra-1, and Sitra-2, respectively. In this context, it was reported that wastewater discharge had a positive effect on seedlings of mangroves through enhancement of growth and reducing the salt-induced damage to the physiological processes of plants [42].Physiognomically, the lower heights and small diameters of Tubli Bay mangroves reflect low structural stand development according to [43]. The shrubby form of mangrove in this part of the world could be attributed to the aridity of the environment in the first place and secondly to the high salinity in the study area [44–46]. The effect of salinity on mangrove is a controversial issue due to variations in site salinities. Salinity was reported as a necessity for the better growth of mangrove [47] while other scholars considered mangrove insensitive to wide ranges of salinity fluctuations [10, 48, 49]. In our study, negative correlation existed between tree heights and EC (−0.9) and between trunk diameters and EC (−1.0), although not significant, probably due to a small sample size. Likewise, positive correlation (0.9) existed between the heights of trees and silt content of the soil profile. Also, the correlation coefficient was not significant between diameters of trees and the silt content of the soil profile (Table 2).Table 2
Correlation coefficient (r) between growth and soil variables at p<0.05.
Growth variables
Soil variables
Cor. coefficient (r)
Prob.
Tree height (m)
EC
-
0.90
0.284
Trunk diameter (cm)
EC
-
1.00
0.052
Tree height (m)
pH
0.95
0.198
Trunk diameter (cm)
pH
0.78
0.430
Tree height (m)
Sand
-
0.59
0.602
Trunk diameter (cm)
Sand
-
0.84
0.370
Tree height (m)
Silt
0.90
0.290
Trunk diameter (cm)
Silt
1.00
0.060
Tree height (m)
Clay
-
0.12
0.925
Trunk diameter (cm)
Clay
0.24
0.843The canopy of trees was relatively open as LAI ranged between 0.58 and 2.4 averaging 1.23 ± 0.03, which is within the range values reported by [50] for the species. Pneumatophores mean height (12.4 ± 1 cm) and density (149.6 ± 16.9 m−2) were significantly different between sites. Lowest number of pneumatophores was observed at Tubli site (106.33), which could be attributed to the effect of sewage discharge as indicated by [51]. The least heights of pneumatophores (6.29 cm) and their maximum number (202) per square meter in the Sitra-1 site could be attributed to sediment deposition where nearly 80% of pneumatophores heights were buried in a 30 cm white-layer depth of silt. Differences in pneumatophores density were ascribed to the microtopographic features [52] and the anaerobic circumstance of sites caused by mud content as well as waterlogging [53]. The silt percentage of the soil profile in the Sitra-1 site ranged from 38.65 at 0–5 cm depth to 43.31 at 15–30 cm, averaging 41.75%, which is considerably higher than the silt particle percentages in comparable depths of the other two sites due to the siltation from the nearby sand washing plant (Table 3). Average portions of silt particles in soil profile were 27.5% and 13.10% in Sitra-2 and Tubli site, respectively.Table 3
Sediment properties ofAvicennia marinastand in the study sites.
Site
Depth (cm)
Soil properties
Texture (%)
Soil texture
Salinity (dS m−1)
pH
Clay
Silt
Sand
Sandy loam
Tubli
0–5
17.99bcd
19.99cd
62.02ab
Sandy loam
16.2bc
8.02a
5–15
12.00d
9.32d
78.68a
Sandy loam
14.7cd
8.04a
15–30
19.99bc
10.00d
70.02a
Loam
10.3d
8.00ab
Sitra-1
0–5
18.66bcd
38.65a
42.68cd
Loam
20.9ab
7.93abc
5–15
23.32abc
43.30a
33.37cd
Clay loam
18.0abc
7.83bcd
15–30
29.32a
43.31a
27.36d
Sandy clay Loam
19.8abc
8.02a
Sitra-2
0–5
23.30abc
26.63bc
50.07bc
Loam
18.2abc
7.71d
5–15
24.62ab
35.92ab
39.46cd
Sandy loam
23.2a
7.71d
15–30
15.96cd
19.95cd
64.09ab
Sandy loam
16.4bc
7.78cd
Levels not connected with the same letter are significantly different (p≤0.05).Finally, the number of seedlings among sites varied from one seedling per square meter in Tubli site to more than nine seedlings in Sitra-1. The presence of microtopography created by the siltation process in Sitra-1 site compared to the other two sites may facilitate seedlings establishments and regeneration of trees [54–56]. The number of seedlings in the Sitra-1 site was significantly higher than those of the other two sites. In the meantime, the height of seedlings (3.84 cm) in Sitra-1 was the least compared to 4.5 cm and 15.1 cm in Tubli and Sitra-2 sites consecutively due to sediment deposition. Statistically, the height of seedlings in Sitra-2 differed significantly than those of Tubli and Sitra-1.


### 3.2. Community Composition

In Tubli site, the relative cover of mangrove averaged 41.4% compared to 52.1% and 54.7% in Sitra-1 and Sitra-2 consecutively. Relative density was 98.5% and 93.6% in Tubli and Sitra-1, respectively, compared to 100% in Sitra-2. Relative frequency was 75% in the site of Tubli, while it registered 92.4% and 91.7% in Sitra-1 and Sitra-2. High relative density figures reveal the purity of mangrove stands in all study sites. IVI for mangrove were relatively similar in the three sites ranging between 138.7 and 155.7 reflecting a clear dominance of the existing plant communities. In the meantime, IVI for the other associated species were of low values reaching, for instance, a maximum of 33 forArthrocnemum macrostachyum in both Tubli and Sitra (1) sites to 39 in Sitra (2) (Figures 3–5). On the other hand, infrequent tide flushing and hypersalinity formed mangrove-salt marsh ecotones where salt-tolerant species likeHalocnemum strobilaceum,Suaeda aegyptiaca,Suaeda maritima,Halopeplis perfoliata,Arthrocnemum macrostachyum, andSalicornia herbacea (Table 4) were zoned beyond the dense mangrove areas. The inadequate flushing and the hypersalinity conditions were highlighted by [17] as the two decisive factors affecting mangrove growth and distribution, which is the case in this region.Table 4
List of mangrove associated species, families, life forms, and habitats in the study sites.
Number
Plant species
Family
Life form
Habitat
Sites
Tubli
Sitra-1
Sitra-2
(
1
)
Arthrocnemum macrostachyum
Chenopodiaceae
Ch
LSM
-
√
√
(
2
)
Avicennia marina
Avicenniaceae
Ph
MSR
√
√
√
(
3
)
Cistanche tubulosa
Orobanchaceae
Cr/P
USM
-
√
√
(
4
)
Citrullus colocynthis
Cucurbitaceae
He
USM
-
√
-
(
5
)
Halocnemum strobilaceum
Amaranthaceae
Ch
LSM
√
√
√
(
6
)
Halopeplis perfoliata
Chenopodiaceae
Ch
LSM
-
√
-
(
7
)
Limonium axillare
Plumbaginaceae
Ch
USM
√
-
-
(
8
)
Phragmites australis
Poaceae
Cr
USM
-
√
-
(
9
)
Pluchea ovalis
Asteraceae
Ph
USM
√
-
-
(
10
)
Salicornia herbacea
Chenopodiaceae
Th
LSM
√
√
√
(
11
)
Sesuvium verrucosum
Aizoaceae
Ch
USM
√
-
-
(
12
)
Suaeda aegyptiaca
Chenopodiaceae
Ch
LSM
√
√
√
(
13
)
Suaeda maritima
Chenopodiaceae
Ch
LSM
-
√
-
(
14
)
Tamarix arabica
Tamaricaceae
Ph
USM
√
√
√
(
15
)
Zygophyllum qatarense
Zygophyllaceae
Ch
USM
√
-
-
Ch: Chamaephytes, Cr: Cryptophytes, He: Hemicryptophytes, P: parasite, Ph: Phanerophyte, Th: Therophytes, LSM: lower salt marsh, MSR: mangrove submerged root, and USM: upper salt marsh.Figure 3
Relative density, cover, frequency, and importance values for Tubli site mangrove stand.Figure 4
Relative density, cover, frequency, and importance values for Sitra-1 mangrove stand.Figure 5
Relative density, cover, frequency, and importance values for Sitra-2 mangrove stand.Three main habitats were identified in the study sites. These habitats were mangrove submerged roots (MSR), lower salt marshes (LSM), and the upper salt marshes (USM). Chamaephytes constitute more than 50% of life forms in the study sites. Species similarity between Tubli and Sitra-1 was 67%, compared to 71% between Tubli and Sitra-2 and 78% between Sitra-1 and Sitra-2 consecutively. Simpson diversity index of species was 0.72, 0.82, and 0.72 for Tubli, Sitra-1, and Sitra-2 sites sequentially.


### 3.3. Siltation Effects

Sediment texture in the study sites was sandy loam in Tubli site, loam to clay loam in Sitra-1, and sandy clay loam/loam to sandy loam in Sitra-2 (see Table3). Sand dominated sediment profile in Tubli site is due to reclamation activities, while in Sitra-2 the sediment profile was more balanced as no human intervention was observed on the site. On the other hand, silt dominated the main component of the soil profile in Sitra-1 plots (39–43%) due to prolonged siltation effects of nearby sand washing plants. Significant differences existed in soil texture between the three sites and between depths of each site.The effect of siltation on the structure and composition of mangroves was evident in elevating site topography which prevented frequent tide inundation of these areas. The occasional flooding permittedHalopeplis perfoliata, an annual halophyte species, to proliferate in these spots. Further, in May of 2015, sudden mangrove mortality was observed at the site along the path of sediment release as a swath measuring 40 m in width and 160 m in length which are equivalent to 22% of the stand area, possibly due to the effect of long-term silt accumulation (Figure 6). Various authors reported the detrimental effects of siltation on root aeration of mangroves, which eventually leads to the death of mangroves [22, 24, 57].Figure 6
Mangrove death due to silt accumulation in Sitra-1 site.


## 3.1. Community Structure

The results of the study revealed that the heights of mangroves in the study sites ranged between 1 and 5.5 m with an average of 2.7 m. Trees were mainly of multiple trunks, ranging from 2 to 5 stems per plant. Diameters of trees ranged in size from 2.2 to 12.5 cm with a mean of 5.3 cm. Mean density and basal area of mangrove trees were 4576.8 tree ha−1 and 11.4 m2 ha−1, respectively. The heights and diameters of trees in Tubli site were 3.25 m and 5.82 cm, respectively. Tree heights were significantly higher (p≤0.05) in the Tubli stand as compared to Sitra (2) but not Sitra (1) (Table 1). Moreover, no significant differences in the height of trees were observed between the Sitra (1) and Sitra (2) stands. High growth parameters could be attributed to the dilution effect of discharged wastewater where measured EC of seawater at Tubli site was 13.2 ± 2 compared to 18.8 ± 2 dS m−1 and 17.1 ± 2 dS m−1 in Sitra-1 and Sitra-2 sites, respectively. However, the insignificance in the height of trees between Tubli site and Sitra (1) could be related to the larger density of trees in the second site, though not significant per se between the two sites. Also, tree diameters were significantly higher (p≤0.05) in the Tubli stand as compared to the Sitra (1) and Sitra (2) stands with no observed significant differences between diameters of trees in the last two sites. On the other hand, no differences in density, basal area, and LAI were observed between the three sites (see Table 1).Table 1
Avicennia marinaattributes in the study sites∗.
Stand
Trees attributes
Stand attributes
Pneumatophores
Seedling
Height (m)
Diameter (cm)
Basal area (m2ha−1)
Leaf area index (m2m−2)
Density (tree ha−1)
Height (cm)
Diameter (cm)
Number m−2
Height (cm)
Number m−2
Tubli
3.25
±
0.3
a
5.82
±
0.1
a
12.99
±
3.1
a
1.28
±
0.05
a
4565.8
±
770
a
15.34
±
0.9
a
1.01
±
0.01
a
106.33
±
12.2
b
4.50
±
2.7
b
1.0
±
1.7
b
Sitra-1
2.58
±
0.3
a
b
5.07
±
0.2
b
11.95
±
1.2
a
1.18
±
0.0
4
a
5085.7
±
79
7
a
6.29
±
0.9
b
0.50
±
0.0
1
b
202.19
±
42.6
a
3.84
±
1.2
b
9.2
±
3.8
a
Sitra-2
2.15
±
0.3
b
5.04
±
0.2
b
9.15
±
1.2
a
1.24
±
0.05
a
4079.1
±
56
5
a
15.50
±
0.7
a
1.01
±
0.01
a
140.33
±
15.1
a
b
15.10
±
1.3
a
5.5
±
1.04
a
b
Mean
2.7
±
0.2
5.3
±
0.1
11.4
±
1.2
1.2
±
0.03
4576.85
±
406
12.4
±
0.97
0.83
±
0.05
149.6
±
16.9
8.4
±
1.96
6.9
±
2.05
p value
0.0200
0.0003
0.4463
0.3308
0.6178
<0.0001
<0.0001
0.0570
0.0012
0.0439
∗Levels not connected with the same letter are significantly different (p≤0.05).In the meantime, the sediment salinity level of soil profiles varied among the sites and ranged from a low value of 10.3 dS m−1 in Tubli to 23.2 dS m−1 in Sitra-2. The average salinity levels were 13.72 dS m−1, 19.57 dS m−1, and 22.27 dS m−1 in Tubli, Sitra-1, and Sitra-2, respectively. In this context, it was reported that wastewater discharge had a positive effect on seedlings of mangroves through enhancement of growth and reducing the salt-induced damage to the physiological processes of plants [42].Physiognomically, the lower heights and small diameters of Tubli Bay mangroves reflect low structural stand development according to [43]. The shrubby form of mangrove in this part of the world could be attributed to the aridity of the environment in the first place and secondly to the high salinity in the study area [44–46]. The effect of salinity on mangrove is a controversial issue due to variations in site salinities. Salinity was reported as a necessity for the better growth of mangrove [47] while other scholars considered mangrove insensitive to wide ranges of salinity fluctuations [10, 48, 49]. In our study, negative correlation existed between tree heights and EC (−0.9) and between trunk diameters and EC (−1.0), although not significant, probably due to a small sample size. Likewise, positive correlation (0.9) existed between the heights of trees and silt content of the soil profile. Also, the correlation coefficient was not significant between diameters of trees and the silt content of the soil profile (Table 2).Table 2
Correlation coefficient (r) between growth and soil variables at p<0.05.
Growth variables
Soil variables
Cor. coefficient (r)
Prob.
Tree height (m)
EC
-
0.90
0.284
Trunk diameter (cm)
EC
-
1.00
0.052
Tree height (m)
pH
0.95
0.198
Trunk diameter (cm)
pH
0.78
0.430
Tree height (m)
Sand
-
0.59
0.602
Trunk diameter (cm)
Sand
-
0.84
0.370
Tree height (m)
Silt
0.90
0.290
Trunk diameter (cm)
Silt
1.00
0.060
Tree height (m)
Clay
-
0.12
0.925
Trunk diameter (cm)
Clay
0.24
0.843The canopy of trees was relatively open as LAI ranged between 0.58 and 2.4 averaging 1.23 ± 0.03, which is within the range values reported by [50] for the species. Pneumatophores mean height (12.4 ± 1 cm) and density (149.6 ± 16.9 m−2) were significantly different between sites. Lowest number of pneumatophores was observed at Tubli site (106.33), which could be attributed to the effect of sewage discharge as indicated by [51]. The least heights of pneumatophores (6.29 cm) and their maximum number (202) per square meter in the Sitra-1 site could be attributed to sediment deposition where nearly 80% of pneumatophores heights were buried in a 30 cm white-layer depth of silt. Differences in pneumatophores density were ascribed to the microtopographic features [52] and the anaerobic circumstance of sites caused by mud content as well as waterlogging [53]. The silt percentage of the soil profile in the Sitra-1 site ranged from 38.65 at 0–5 cm depth to 43.31 at 15–30 cm, averaging 41.75%, which is considerably higher than the silt particle percentages in comparable depths of the other two sites due to the siltation from the nearby sand washing plant (Table 3). Average portions of silt particles in soil profile were 27.5% and 13.10% in Sitra-2 and Tubli site, respectively.Table 3
Sediment properties ofAvicennia marinastand in the study sites.
Site
Depth (cm)
Soil properties
Texture (%)
Soil texture
Salinity (dS m−1)
pH
Clay
Silt
Sand
Sandy loam
Tubli
0–5
17.99bcd
19.99cd
62.02ab
Sandy loam
16.2bc
8.02a
5–15
12.00d
9.32d
78.68a
Sandy loam
14.7cd
8.04a
15–30
19.99bc
10.00d
70.02a
Loam
10.3d
8.00ab
Sitra-1
0–5
18.66bcd
38.65a
42.68cd
Loam
20.9ab
7.93abc
5–15
23.32abc
43.30a
33.37cd
Clay loam
18.0abc
7.83bcd
15–30
29.32a
43.31a
27.36d
Sandy clay Loam
19.8abc
8.02a
Sitra-2
0–5
23.30abc
26.63bc
50.07bc
Loam
18.2abc
7.71d
5–15
24.62ab
35.92ab
39.46cd
Sandy loam
23.2a
7.71d
15–30
15.96cd
19.95cd
64.09ab
Sandy loam
16.4bc
7.78cd
Levels not connected with the same letter are significantly different (p≤0.05).Finally, the number of seedlings among sites varied from one seedling per square meter in Tubli site to more than nine seedlings in Sitra-1. The presence of microtopography created by the siltation process in Sitra-1 site compared to the other two sites may facilitate seedlings establishments and regeneration of trees [54–56]. The number of seedlings in the Sitra-1 site was significantly higher than those of the other two sites. In the meantime, the height of seedlings (3.84 cm) in Sitra-1 was the least compared to 4.5 cm and 15.1 cm in Tubli and Sitra-2 sites consecutively due to sediment deposition. Statistically, the height of seedlings in Sitra-2 differed significantly than those of Tubli and Sitra-1.


## 3.2. Community Composition

In Tubli site, the relative cover of mangrove averaged 41.4% compared to 52.1% and 54.7% in Sitra-1 and Sitra-2 consecutively. Relative density was 98.5% and 93.6% in Tubli and Sitra-1, respectively, compared to 100% in Sitra-2. Relative frequency was 75% in the site of Tubli, while it registered 92.4% and 91.7% in Sitra-1 and Sitra-2. High relative density figures reveal the purity of mangrove stands in all study sites. IVI for mangrove were relatively similar in the three sites ranging between 138.7 and 155.7 reflecting a clear dominance of the existing plant communities. In the meantime, IVI for the other associated species were of low values reaching, for instance, a maximum of 33 forArthrocnemum macrostachyum in both Tubli and Sitra (1) sites to 39 in Sitra (2) (Figures 3–5). On the other hand, infrequent tide flushing and hypersalinity formed mangrove-salt marsh ecotones where salt-tolerant species likeHalocnemum strobilaceum,Suaeda aegyptiaca,Suaeda maritima,Halopeplis perfoliata,Arthrocnemum macrostachyum, andSalicornia herbacea (Table 4) were zoned beyond the dense mangrove areas. The inadequate flushing and the hypersalinity conditions were highlighted by [17] as the two decisive factors affecting mangrove growth and distribution, which is the case in this region.Table 4
List of mangrove associated species, families, life forms, and habitats in the study sites.
Number
Plant species
Family
Life form
Habitat
Sites
Tubli
Sitra-1
Sitra-2
(
1
)
Arthrocnemum macrostachyum
Chenopodiaceae
Ch
LSM
-
√
√
(
2
)
Avicennia marina
Avicenniaceae
Ph
MSR
√
√
√
(
3
)
Cistanche tubulosa
Orobanchaceae
Cr/P
USM
-
√
√
(
4
)
Citrullus colocynthis
Cucurbitaceae
He
USM
-
√
-
(
5
)
Halocnemum strobilaceum
Amaranthaceae
Ch
LSM
√
√
√
(
6
)
Halopeplis perfoliata
Chenopodiaceae
Ch
LSM
-
√
-
(
7
)
Limonium axillare
Plumbaginaceae
Ch
USM
√
-
-
(
8
)
Phragmites australis
Poaceae
Cr
USM
-
√
-
(
9
)
Pluchea ovalis
Asteraceae
Ph
USM
√
-
-
(
10
)
Salicornia herbacea
Chenopodiaceae
Th
LSM
√
√
√
(
11
)
Sesuvium verrucosum
Aizoaceae
Ch
USM
√
-
-
(
12
)
Suaeda aegyptiaca
Chenopodiaceae
Ch
LSM
√
√
√
(
13
)
Suaeda maritima
Chenopodiaceae
Ch
LSM
-
√
-
(
14
)
Tamarix arabica
Tamaricaceae
Ph
USM
√
√
√
(
15
)
Zygophyllum qatarense
Zygophyllaceae
Ch
USM
√
-
-
Ch: Chamaephytes, Cr: Cryptophytes, He: Hemicryptophytes, P: parasite, Ph: Phanerophyte, Th: Therophytes, LSM: lower salt marsh, MSR: mangrove submerged root, and USM: upper salt marsh.Figure 3
Relative density, cover, frequency, and importance values for Tubli site mangrove stand.Figure 4
Relative density, cover, frequency, and importance values for Sitra-1 mangrove stand.Figure 5
Relative density, cover, frequency, and importance values for Sitra-2 mangrove stand.Three main habitats were identified in the study sites. These habitats were mangrove submerged roots (MSR), lower salt marshes (LSM), and the upper salt marshes (USM). Chamaephytes constitute more than 50% of life forms in the study sites. Species similarity between Tubli and Sitra-1 was 67%, compared to 71% between Tubli and Sitra-2 and 78% between Sitra-1 and Sitra-2 consecutively. Simpson diversity index of species was 0.72, 0.82, and 0.72 for Tubli, Sitra-1, and Sitra-2 sites sequentially.


## 3.3. Siltation Effects

Sediment texture in the study sites was sandy loam in Tubli site, loam to clay loam in Sitra-1, and sandy clay loam/loam to sandy loam in Sitra-2 (see Table3). Sand dominated sediment profile in Tubli site is due to reclamation activities, while in Sitra-2 the sediment profile was more balanced as no human intervention was observed on the site. On the other hand, silt dominated the main component of the soil profile in Sitra-1 plots (39–43%) due to prolonged siltation effects of nearby sand washing plants. Significant differences existed in soil texture between the three sites and between depths of each site.The effect of siltation on the structure and composition of mangroves was evident in elevating site topography which prevented frequent tide inundation of these areas. The occasional flooding permittedHalopeplis perfoliata, an annual halophyte species, to proliferate in these spots. Further, in May of 2015, sudden mangrove mortality was observed at the site along the path of sediment release as a swath measuring 40 m in width and 160 m in length which are equivalent to 22% of the stand area, possibly due to the effect of long-term silt accumulation (Figure 6). Various authors reported the detrimental effects of siltation on root aeration of mangroves, which eventually leads to the death of mangroves [22, 24, 57].Figure 6
Mangrove death due to silt accumulation in Sitra-1 site.


## 4. Conclusions

Mangrove formed pure plant communities along the coastlines of Tubli Bay, intermingled in some cases with other salt-tolerant species. The aridity of the environment, salinity, and the extent of the tide flooded area contributed to the formation of a few species salt-tolerant plant community in the study sites. The structure and composition of mangrove communities were similar in all of the study sites. Nevertheless, wastewater discharge had a positive effect on heights and diameters of mangrove trees in Tubli site due to the dilution effect of secondary treated wastewater. In the meantime, anthropogenic siltation had a detrimental impact on mangroves on the exposed site (Sitra-1) by excluding mangroves from some spots and burying roots of the species which led eventually to swath death in the stand. Further investigation is needed to reveal the exact cause and mechanism of positive mangrove response to wastewater discharge and the death mechanism caused by sedimentation in mangroves of Tubli Bay.



---

*Source: 2084256-2017-12-19.xml*