# Carbon Stocks in the Small Estuarine Mangroves of Geza and Mtimbwani, Tanga, Tanzania

Authors: Edmond Alavaisha; Mwita M. Mangora
Journal: International Journal of Forestry Research
 (2016)

---

## Abstract

Mangrove forests offer important ecosystem services, including their high capacity for carbon sequestration and stocking. However, they face rapid degradation and loss of ecological resilience particularly at local scales due to human pressure. We conducted inventory of mangrove forests to characterise forest stand structure and estimate carbon stocks in the small estuarine mangroves of Geza and Mtimbwani in Tanga, Tanzania. Forest structure, above-ground carbon (AGC), and below-ground carbon (BGC) were characterised. Soil carbon was estimated to 1 m depth using loss on ignition procedure. Six common mangrove species were identified dominated byAvicennia marina (Forsk.) Vierh. andRhizophora mucronata Lamarck. Forest stand density and basal area were 1740 stems ha−1 and 17.2 m2 ha−1 for Geza and 2334 stems ha−1 and 30.3 m2 ha−1 for Mtimbwani. Total ecosystem carbon stocks were 414.6 Mg C ha−1 for Geza and 684.9 Mg C ha−1 for Mtimbwani. Soil carbon contributed over 65% of these stocks, decreasing with depth. Mid zones of the mangrove stands had highest carbon stocks. These data demonstrate that studied mangroves are potential for carbon projects and provide the baseline for monitoring, reporting, and verification (MRV) to support the projects.


---

## Body

## 1. Introduction

Mangrove forests occur in fragmented stands along almost the entire coastline of Tanzania mainland. Mangroves flourish in river estuaries and deltas and in enclosed bays, lagoons, and tidal creeks. Despite representing only about 0.3% (108,000–115,000 ha) of the total forest area in the country [1, 2], these forests provide numerous ecosystem services including wood and nonwood products; buffering lands and coastal properties [1, 3]; and supporting fisheries as they serve as nursery, feeding, and shelter grounds for numerous commercially important fishes and invertebrates [4–6]. In global terms, mangroves are reported to be up to five times efficient carbon sinks compared to other forms of terrestrial forests [7–9]. However, there exists large nonlinearity in ecosystem services provided by different mangrove formations [8, 10] that hinder generalised estimates of mangrove forest ecosystems’ carbon sequestration and stocks which is complicated by local conditions such as climate and soil factors; forest age, growth, and structure; utilisation and management regime of respective forests [3, 7, 11–13]. As such, there is a considerable knowledge gap and uncertainty at local levels regarding the carbon pool size and variability of carbon sequestration within mangrove forests [14, 15]. Nonetheless, the reported high rate of mangrove carbon sequestration provides global benefits in mitigating the effects of climate change [9, 13, 14] and demonstrates potential for livelihood enhancement through community carbon market schemes for sustainable conservation [16, 17].Like in many other countries where mangroves occur, overexploitation and conversions to other land uses like solar salt pans, aquaculture, agriculture, and urban expansion are major threats to mangroves of Tanzania [1, 18]. While the global rate of mangrove loss is approximated at 1-2% per year [19–21], there are variations at regional, national, and local scales. In Tanzania, FAO [22] estimated the rate of mangrove loss at about 0.7% per year. In most cases there are no reliable records at local scales, and where they exist, they are potentially out dated [2]. In general, mangrove degradation jeopardizes their ability to provide ecosystem goods and services [7], but the full implication of the loss is not well understood and appreciated. This calls for urgent actions to design and implement new effective conservation strategies to protect and restore mangroves to sustain the ecosystem services that they provide [21, 23, 24].Over the recent past, international climate change agreements have highlighted Reduced Emissions from Deforestation and Degradation (REDD) as one of the possible cost-effective strategies for mitigating the effects of climate change [7, 14, 25]. Reducing emissions aim at maintaining forest carbon stores through financial incentives for forest conservation, for example, carbon credits [26]. With the emerging market-based conservation strategies involving carbon credits in schemes like REDD+ (Reducing Emissions from Deforestation and Forest Degradation, Sustainable Management of Forests, and Enhancement of Forest Carbon Stocks in Developing Countries) and PES (Payment for Ecosystem Services), mangroves are now increasingly recognized as promising natural solutions for carbon capture and long-term storage putting them on top of the agenda in the debates on climate change adaptation mechanisms due to their high rates of carbon sequestration [7, 8, 14]. This warrants for site-specific estimations of the carbon stocks in mangrove forests. In this paper, we report on results of the inventory in the mangrove forests of Geza and Mtimbwani villages in the north-eastern region of Tanga, Tanzania, that aimed to (i) characterise the tree composition and structure; (ii) estimate the carbon stocks; and (iii) demonstrate the potential for community based carbon market schemes.


## 2. Materials and Methods

### 2.1. Description of Study Sites

Geza village is about 26 km south of the city of Tanga whereas Mtimbwani village is 16 km north of the city on the north coast of Tanzania (Figure1). Mangroves of Geza are approximately 84 ha and are fed seasonally by Bongoa River while those of Mtimbwani are about 326 ha and receive freshwater discharge throughout the year from Ngole River. Common species in both forests areAvicennia marina (Forsk.) Vierh.,Rhizophora mucronata Lamarck,Ceriops tagal (Perr.) C. B. Robinson,Sonneratia alba J. Smith,Bruguiera gymnorrhiza (L.) Lamarck, andXylocarpus granatum König. The forests receive normal tidal inundation twice a day though the extent of tidal influence in upstream varies during rain and dry season. The average annual precipitation is between 1100 mm and 1900 mm, relative humidity of about 76% during daytime and 96% at night, and the average range of atmospheric temperature is from 18°C to 35°C.Figure 1
Map and Google Earth snapshots showing study sites and sampling plot locations.


### 2.2. Sampling Plots

Each mangrove site was divided into three zones along the river gradient from the sea landward using tidal range as criteria: lower, mid, and upper zone. In each zone, four circular nested plots of 10 m radius [27] spaced at approximately 200 m were established, giving a total of twelve sampling plots per each site (Figure 1). Plot locations were predetermined on Google Earth (GE), coordinates recorded and reached during field work by aid of the handheld GPS receiver (Garmin GPSMAP 64s).


### 2.3. Vegetation Inventory

Mangrove forest inventory protocols described by Kauffman and Donato [27] were adopted and modified accordingly to suite local conditions. In each plot, all trees with <2.5 cm stem diameter were counted and grouped as seedlings by species. Trees with ≥2.5 cm stem diameter at breast height (DBH, 1.3 m above the ground) were identified by species and measured for DBH (cm) and height (m). For the stilt rootedR.mucronata, DBH were measured at 30 cm above the highest root [28]. Standing dead wood trees were practically defined according to their death status. Trees with branches and twigs, which resemble a live tree, were categorized as “status 1,” trees with large branches only as “status 2,” and trees with bole (trunk) only, no branches, as “status 3.” For dead and downed wood, a line-intersect method was used for counting and measuring intersections of woody pieces with diameter >7.6 cm along a vertical sampling plane (transect). For the purpose of degradation, stumps were counted regardless of their species [27].


### 2.4. Soil Sampling for Carbon

Soil cores were retrieved from the centre of each sampling plot up to 1 m deep using a multistage sediment core sampler (AMS, Inc., American Falls, USA). The full length of core sampler was steadily inserted vertically into the soil at each depth stage. Soil samples were retrieved in divided segment depths of 0–30 cm, 30–60 cm, and 60–100 cm. Each retrieved core was subdivided in 5 cm and mean values of each segment depth were used for analysis.


### 2.5. Data Analyses

#### 2.5.1. Forest Structure and Composition

Forest inventory data were processed using standard analysis procedures as described by Cintron and Novelli [29] to derive forest stand characteristics: stand frequency distribution, density (stems ha−1), basal area (m2 ha−1), relative frequency (1), relative density (2), and relative dominance (3). Ecological importance values (IV) of each species were determined by summing the respective relative frequency, relative density, and relative dominance. Importance value measures relative dominance of species by criteria of how often it occurred, number of species, and area it occupies in a community. The species that attained the highest IV was considered the principal species: (1)
Relative frequency
=
100
F
i
∑
i
=
1
m
​
F
i
,
(2)
Relative density
=
100
n
i
∑
i
=
1
m
​
n
i
,
(3)
Relative dominance
=
100
B
a
i
B
a
,where n
i is the number of trees sampled for species i; m is the number of species; F
i is number of plots in which species i occurred, multiplied by 100.


#### 2.5.2. Above- and below-Ground Biomass and Carbon Pools

Estimation of above-ground biomass (AGB) and below-ground biomass (BGB) for roots in live trees used general allometric equations (4) and (5), respectively, developed by Komiyama et al. [28, 30]:(4)
AGB
=
0.251
×
ρ
×
DBH
2.46
,
(5)
BGB
=
0.199
×
ρ
0.899
×
DBH
2.22
,where AGB and BGB are biomass (kg), DBH is diameter at breast height (cm), and ρ is average general wood density (0.752 g cm−3). These general equations and wood density were used because local or regional species-specific models and respective wood density for all the species have not been developed and established.Biomass of standing dead wood was obtained by subtraction of a percentage factor from the supposed biomass of live tree derived by a general formula to account for the loss of leaf and twigs biomass. For standing dead trees of status 1 and status 2, factors of 2.5% and 10% biomass of live tree were subtracted, respectively. For status 3, standing dead trees and downed dead wood biomass were estimated by factoring average wood density (0.752 g cm−3) to the volume of dead wood. Volume of downed dead wood was calculated using (6) by Brown [31]:(6)
Volume
m
3
h
-
1
=
π
2
d
1
2
+
d
2
2
+
d
3
2
+
⋯
+
d
n
2
8
L
.
π is constant; d
1
,
d
2
,
…
,
d
n are diameters of intersecting pieces of large dead wood (cm), and L is the length of the transect line for large size class (m).To obtain above-ground carbon (AGC) and below-ground carbon (BGC) density, default values of carbon concentration were used: 0.48 for above-ground live trees, 0.5 for dead and downed trees, and 0.39 for below-ground roots as suggested by Kauffman and Donato [27] and Howard et al. [32].


#### 2.5.3. Soil Carbon Pool

Dry combustion method was used to analyse soil carbon. Soil samples were dried at 60°C to constant weight for at least 72 hours and then allowed to cool to room temperature before weighing. Cooled samples were weighed and used for bulk density determination. Bulk density (g cm−3) was determined by dividing the oven-dry soil sample mass (g) by the volume (cm−3) of the sample: (7)
Soil bulk density
g
cm
-
3
=
oven-dry sample mass
g
sample volume
cm
3
.Twenty grams of oven dried subsample from each sample was ignited in a muffle furnace (AAF 11/7, Wolf Laboratories Limited, UK) at 540°C for 5 hours. The differences obtained between the dry weight and ignited weight (loss on ignition, LOI) indicated organic matter content expressed in percentage (% LOI) as(8)
%
LOI
=
dry mass before combustion
-
dry mass after combustion
dry mass before combustion
×
100
.Soil carbon density was determined as(9)
Soil carbon
Mg
ha
-
1
=
bulk density
g
cm
-
3
×
soil depth interval
cm
×
%
C
.To obtain organic carbon concentration (% C), an equation suggested by Kauffman and Donato [27] was adopted:(10)
Organic carbon concentration
%
C
=
0.415
×
%
LOI
+
2.89
.


#### 2.5.4. Total Ecosystem Carbon Stock

Estimation of ecosystem carbon density was done by summing all of the component carbon stocks of above-ground, below-ground, and soil.


#### 2.5.5. Statistical Tests

Statistical tests were performed for descriptive and inferential statistics on forest inventory variables. Mean values were subjected to one-way Analysis of Variance (ANOVA) at significance level of 0.05 to find the differences of mean forest density, basal area, stumps, regeneration, and carbon stocks within and between sites.


## 2.1. Description of Study Sites

Geza village is about 26 km south of the city of Tanga whereas Mtimbwani village is 16 km north of the city on the north coast of Tanzania (Figure1). Mangroves of Geza are approximately 84 ha and are fed seasonally by Bongoa River while those of Mtimbwani are about 326 ha and receive freshwater discharge throughout the year from Ngole River. Common species in both forests areAvicennia marina (Forsk.) Vierh.,Rhizophora mucronata Lamarck,Ceriops tagal (Perr.) C. B. Robinson,Sonneratia alba J. Smith,Bruguiera gymnorrhiza (L.) Lamarck, andXylocarpus granatum König. The forests receive normal tidal inundation twice a day though the extent of tidal influence in upstream varies during rain and dry season. The average annual precipitation is between 1100 mm and 1900 mm, relative humidity of about 76% during daytime and 96% at night, and the average range of atmospheric temperature is from 18°C to 35°C.Figure 1
Map and Google Earth snapshots showing study sites and sampling plot locations.


## 2.2. Sampling Plots

Each mangrove site was divided into three zones along the river gradient from the sea landward using tidal range as criteria: lower, mid, and upper zone. In each zone, four circular nested plots of 10 m radius [27] spaced at approximately 200 m were established, giving a total of twelve sampling plots per each site (Figure 1). Plot locations were predetermined on Google Earth (GE), coordinates recorded and reached during field work by aid of the handheld GPS receiver (Garmin GPSMAP 64s).


## 2.3. Vegetation Inventory

Mangrove forest inventory protocols described by Kauffman and Donato [27] were adopted and modified accordingly to suite local conditions. In each plot, all trees with <2.5 cm stem diameter were counted and grouped as seedlings by species. Trees with ≥2.5 cm stem diameter at breast height (DBH, 1.3 m above the ground) were identified by species and measured for DBH (cm) and height (m). For the stilt rootedR.mucronata, DBH were measured at 30 cm above the highest root [28]. Standing dead wood trees were practically defined according to their death status. Trees with branches and twigs, which resemble a live tree, were categorized as “status 1,” trees with large branches only as “status 2,” and trees with bole (trunk) only, no branches, as “status 3.” For dead and downed wood, a line-intersect method was used for counting and measuring intersections of woody pieces with diameter >7.6 cm along a vertical sampling plane (transect). For the purpose of degradation, stumps were counted regardless of their species [27].


## 2.4. Soil Sampling for Carbon

Soil cores were retrieved from the centre of each sampling plot up to 1 m deep using a multistage sediment core sampler (AMS, Inc., American Falls, USA). The full length of core sampler was steadily inserted vertically into the soil at each depth stage. Soil samples were retrieved in divided segment depths of 0–30 cm, 30–60 cm, and 60–100 cm. Each retrieved core was subdivided in 5 cm and mean values of each segment depth were used for analysis.


## 2.5. Data Analyses

### 2.5.1. Forest Structure and Composition

Forest inventory data were processed using standard analysis procedures as described by Cintron and Novelli [29] to derive forest stand characteristics: stand frequency distribution, density (stems ha−1), basal area (m2 ha−1), relative frequency (1), relative density (2), and relative dominance (3). Ecological importance values (IV) of each species were determined by summing the respective relative frequency, relative density, and relative dominance. Importance value measures relative dominance of species by criteria of how often it occurred, number of species, and area it occupies in a community. The species that attained the highest IV was considered the principal species: (1)
Relative frequency
=
100
F
i
∑
i
=
1
m
​
F
i
,
(2)
Relative density
=
100
n
i
∑
i
=
1
m
​
n
i
,
(3)
Relative dominance
=
100
B
a
i
B
a
,where n
i is the number of trees sampled for species i; m is the number of species; F
i is number of plots in which species i occurred, multiplied by 100.


### 2.5.2. Above- and below-Ground Biomass and Carbon Pools

Estimation of above-ground biomass (AGB) and below-ground biomass (BGB) for roots in live trees used general allometric equations (4) and (5), respectively, developed by Komiyama et al. [28, 30]:(4)
AGB
=
0.251
×
ρ
×
DBH
2.46
,
(5)
BGB
=
0.199
×
ρ
0.899
×
DBH
2.22
,where AGB and BGB are biomass (kg), DBH is diameter at breast height (cm), and ρ is average general wood density (0.752 g cm−3). These general equations and wood density were used because local or regional species-specific models and respective wood density for all the species have not been developed and established.Biomass of standing dead wood was obtained by subtraction of a percentage factor from the supposed biomass of live tree derived by a general formula to account for the loss of leaf and twigs biomass. For standing dead trees of status 1 and status 2, factors of 2.5% and 10% biomass of live tree were subtracted, respectively. For status 3, standing dead trees and downed dead wood biomass were estimated by factoring average wood density (0.752 g cm−3) to the volume of dead wood. Volume of downed dead wood was calculated using (6) by Brown [31]:(6)
Volume
m
3
h
-
1
=
π
2
d
1
2
+
d
2
2
+
d
3
2
+
⋯
+
d
n
2
8
L
.
π is constant; d
1
,
d
2
,
…
,
d
n are diameters of intersecting pieces of large dead wood (cm), and L is the length of the transect line for large size class (m).To obtain above-ground carbon (AGC) and below-ground carbon (BGC) density, default values of carbon concentration were used: 0.48 for above-ground live trees, 0.5 for dead and downed trees, and 0.39 for below-ground roots as suggested by Kauffman and Donato [27] and Howard et al. [32].


### 2.5.3. Soil Carbon Pool

Dry combustion method was used to analyse soil carbon. Soil samples were dried at 60°C to constant weight for at least 72 hours and then allowed to cool to room temperature before weighing. Cooled samples were weighed and used for bulk density determination. Bulk density (g cm−3) was determined by dividing the oven-dry soil sample mass (g) by the volume (cm−3) of the sample: (7)
Soil bulk density
g
cm
-
3
=
oven-dry sample mass
g
sample volume
cm
3
.Twenty grams of oven dried subsample from each sample was ignited in a muffle furnace (AAF 11/7, Wolf Laboratories Limited, UK) at 540°C for 5 hours. The differences obtained between the dry weight and ignited weight (loss on ignition, LOI) indicated organic matter content expressed in percentage (% LOI) as(8)
%
LOI
=
dry mass before combustion
-
dry mass after combustion
dry mass before combustion
×
100
.Soil carbon density was determined as(9)
Soil carbon
Mg
ha
-
1
=
bulk density
g
cm
-
3
×
soil depth interval
cm
×
%
C
.To obtain organic carbon concentration (% C), an equation suggested by Kauffman and Donato [27] was adopted:(10)
Organic carbon concentration
%
C
=
0.415
×
%
LOI
+
2.89
.


### 2.5.4. Total Ecosystem Carbon Stock

Estimation of ecosystem carbon density was done by summing all of the component carbon stocks of above-ground, below-ground, and soil.


### 2.5.5. Statistical Tests

Statistical tests were performed for descriptive and inferential statistics on forest inventory variables. Mean values were subjected to one-way Analysis of Variance (ANOVA) at significance level of 0.05 to find the differences of mean forest density, basal area, stumps, regeneration, and carbon stocks within and between sites.


## 2.5.1. Forest Structure and Composition

Forest inventory data were processed using standard analysis procedures as described by Cintron and Novelli [29] to derive forest stand characteristics: stand frequency distribution, density (stems ha−1), basal area (m2 ha−1), relative frequency (1), relative density (2), and relative dominance (3). Ecological importance values (IV) of each species were determined by summing the respective relative frequency, relative density, and relative dominance. Importance value measures relative dominance of species by criteria of how often it occurred, number of species, and area it occupies in a community. The species that attained the highest IV was considered the principal species: (1)
Relative frequency
=
100
F
i
∑
i
=
1
m
​
F
i
,
(2)
Relative density
=
100
n
i
∑
i
=
1
m
​
n
i
,
(3)
Relative dominance
=
100
B
a
i
B
a
,where n
i is the number of trees sampled for species i; m is the number of species; F
i is number of plots in which species i occurred, multiplied by 100.


## 2.5.2. Above- and below-Ground Biomass and Carbon Pools

Estimation of above-ground biomass (AGB) and below-ground biomass (BGB) for roots in live trees used general allometric equations (4) and (5), respectively, developed by Komiyama et al. [28, 30]:(4)
AGB
=
0.251
×
ρ
×
DBH
2.46
,
(5)
BGB
=
0.199
×
ρ
0.899
×
DBH
2.22
,where AGB and BGB are biomass (kg), DBH is diameter at breast height (cm), and ρ is average general wood density (0.752 g cm−3). These general equations and wood density were used because local or regional species-specific models and respective wood density for all the species have not been developed and established.Biomass of standing dead wood was obtained by subtraction of a percentage factor from the supposed biomass of live tree derived by a general formula to account for the loss of leaf and twigs biomass. For standing dead trees of status 1 and status 2, factors of 2.5% and 10% biomass of live tree were subtracted, respectively. For status 3, standing dead trees and downed dead wood biomass were estimated by factoring average wood density (0.752 g cm−3) to the volume of dead wood. Volume of downed dead wood was calculated using (6) by Brown [31]:(6)
Volume
m
3
h
-
1
=
π
2
d
1
2
+
d
2
2
+
d
3
2
+
⋯
+
d
n
2
8
L
.
π is constant; d
1
,
d
2
,
…
,
d
n are diameters of intersecting pieces of large dead wood (cm), and L is the length of the transect line for large size class (m).To obtain above-ground carbon (AGC) and below-ground carbon (BGC) density, default values of carbon concentration were used: 0.48 for above-ground live trees, 0.5 for dead and downed trees, and 0.39 for below-ground roots as suggested by Kauffman and Donato [27] and Howard et al. [32].


## 2.5.3. Soil Carbon Pool

Dry combustion method was used to analyse soil carbon. Soil samples were dried at 60°C to constant weight for at least 72 hours and then allowed to cool to room temperature before weighing. Cooled samples were weighed and used for bulk density determination. Bulk density (g cm−3) was determined by dividing the oven-dry soil sample mass (g) by the volume (cm−3) of the sample: (7)
Soil bulk density
g
cm
-
3
=
oven-dry sample mass
g
sample volume
cm
3
.Twenty grams of oven dried subsample from each sample was ignited in a muffle furnace (AAF 11/7, Wolf Laboratories Limited, UK) at 540°C for 5 hours. The differences obtained between the dry weight and ignited weight (loss on ignition, LOI) indicated organic matter content expressed in percentage (% LOI) as(8)
%
LOI
=
dry mass before combustion
-
dry mass after combustion
dry mass before combustion
×
100
.Soil carbon density was determined as(9)
Soil carbon
Mg
ha
-
1
=
bulk density
g
cm
-
3
×
soil depth interval
cm
×
%
C
.To obtain organic carbon concentration (% C), an equation suggested by Kauffman and Donato [27] was adopted:(10)
Organic carbon concentration
%
C
=
0.415
×
%
LOI
+
2.89
.


## 2.5.4. Total Ecosystem Carbon Stock

Estimation of ecosystem carbon density was done by summing all of the component carbon stocks of above-ground, below-ground, and soil.


## 2.5.5. Statistical Tests

Statistical tests were performed for descriptive and inferential statistics on forest inventory variables. Mean values were subjected to one-way Analysis of Variance (ANOVA) at significance level of 0.05 to find the differences of mean forest density, basal area, stumps, regeneration, and carbon stocks within and between sites.


## 3. Results

### 3.1. Forests Structure and Composition

Summary forest structural stand characteristics are presented in Table1. Six species of mangroves were found in both sites.A.marina was the most frequently observed and dominant species in both sites. The least observed species wasB.gymnorrhiza, whileS.alba was the least dominant species. In terms of density,R.mucronata had high relative stem density in Geza (58.9%) and Mtimbwani (41.8%) whereasS.alba was the least contributing species to stem density in both Geza (0.8%) and Mtimbwani (1.1%). Based on their IVs,R.mucronata was the principal species in Geza andA.marina in Mtimbwani. Seedlings and saplings count indicatedR.mucronata andA.marina as the most abundant in both sites.Table 1
Forest stand characteristics for Geza and Mtimbwani mangroves.
Site
Species
Frequency
Density (stems ha-
1)
Basal area(m2 ha-
1)
Seedlings/saplings (counts ha-
1)
Relative
Importance value
Frequency
Density
Dominance
Geza
A. marina
12
406
10.22
778
28.57
23.32
59.37
111.26
B. gymnorrhiza
2
27
0.24
17
4.76
1.52
1.39
7.67
C. tagal
10
231
0.43
531
26.19
13.26
2.49
41.94
R. mucronata
11
1026
5.82
1219
26.19
58.99
33.78
118.97
S. alba
3
13
0.11
0
7.14
0.76
0.67
8.57
X. granatum
3
37
0.4
24
7.14
2.13
2.31
11.59
Total
1740
17.21
2569
100
100
100
300
Mtimbwani
A. marina
12
621
20.78
326
26.09
26.59
68.64
121.31
B. gymnorrhiza
4
80
0.49
35
8.7
3.41
1.63
13.73
C. tagal
8
536
3.51
238
17.39
22.95
11.59
51.93
R. mucronata
11
976
4.49
483
23.91
41.82
14.83
80.56
S. alba
6
27
0.18
18
13.04
1.14
0.58
14.76
X. granatum
5
95
0.83
44
10.87
4.09
2.74
17.7
Total
2334
30.27
1144
100
100
100
300Mangroves of Mtimbwani had high mean stem density across the three zones compared to Geza mangroves (Figure2(a)). In both sites, mid zones had the highest mean stem density (Geza = 2252
±
810 stems ha−1 and Mtimbwani = 2944
±
634 stems ha−1) than other zones but not signifcantly different both between the sites (p
=
0.075) and within zones (p
=
0.233 for Geza and p
=
0.115 for Mtimbwani). For basal areas, differences were significant (p
=
0.041) between sites where Mtimbwani site had large mean basal area than that of Geza (Figure 2(b)). The mid zones represented large basal areas in both Geza (22.25
±
7.1 m2 ha−1) and Mtimbwani (46.10
±
12.1 m2 ha−1) although the differences were not significant among zones in Geza (p
=
0.397) and Mtimbwani (p
=
0.114).Figure 2
Forest structure indicating (a) stems density (b), basal area, (c) number of stumps, and (d) potential regeneration in different zones (mean ± SD).
(a)
(b)
(c)
(d)The mean number of stumps was not significantly different (p
=
0.112) between Geza and Mtimbwani sites (Figure 2(c)). Comparison within zones revealed a significant difference (p
=
0.031) in mean stumps count between the zones of Geza mangroves with highest number in the upper zone (49
±
14 stumps ha−1). In Mtimbwani, the difference in mean stumps count was not significant (p
=
0.539) between the zones although the lower zone had higher count (56
±
22 stumps ha−1). Fewer stumps were counted in mid zones for both sites. Overall, Geza had significantly higher mean number of seedlings than Mtimbwani (p
=
0.003; Figure 2(d)). Mean seedlings and saplings count was not significantly different between zones in both Geza (p
=
0.183) and Mtimbwani (p
=
0.168) although counts were higher in the mid zone of Geza (1107
±
437) and lower zone of Mtimbwani (512
±
91).


### 3.2. Above- and below-Ground Carbon Stocks

The difference in AGC and BGC density was significantly different between sites (p
=
0.041 for Geza and p
=
0.046 for Mtimbwani). There was no significant difference in AGC (p
=
0.214) and BGC (p
=
0.158) between zones of Geza but the difference was significant for Mtimbwani AGC (p
=
0.021) and BGC (p
=
0.038) where mid zone had higher AGC and BGC compared to lower and upper zones (Figure 3). Standing live trees contributed more than 70% of the AGC (Table 2). On average, BGC in Geza and Mtimbwani represented 24.26% and 21.93%, respectively, of vegetative carbon.Table 2
Mean AGC density (Mg C ha−1) for standing live, dead and downed trees (mean ± SD).
Site
Component
Zone
Low
Mid
Upper
Geza
Standing live tree
54.88± 19.16
102.63± 30.67
65.66± 14.45
Standing dead trees
9.02± 3.06
1.32± 0.48
1.16± 0.53
Downed trees
0.23± 0.08
0.03± 0.02
0.22± 0.18
Mtimbwani
Standing live tree
103.44± 39.53
167.84± 57.11
99.50± 45.15
Standing dead trees
7.69± 4.90
133.59± 126.43
10.12± 4.05
Downed trees
0.24± 0.20
0.10± 0.01
0.07± 0.03Figure 3
Mean AGC and BGC density (Mg C ha−1) in different zones for (a) Geza and (b) Mtimbwani (mean ± SD).
(a)
(b)


### 3.3. Mangrove Soil Carbon Stocks

#### 3.3.1. Vertical Profile of Carbon Stocks

The general profile in mean soil bulk density and carbon stock in different depths is summarized in Table3. The mean bulk density and soil carbon density significantly decreased with depth in both sites, Geza (p
=
0.038, p
=
0.024) and Mtimbwani (p
=
0.007, p
=
0.040). The top 30 cm had higher carbon density compared to 30–60 cm and 60–100 cm, demonstrating a higher C% in the top layers of the soil.Table 3
Vertical profile of mean soil bulk density and carbon stocks (mean± SD).
Site
Depth (cm)
Mean soil bulk density (g cm-
3)
Mean soil carbon (Mg C ha-
1)
Geza
0–30
0.62± 0.20a
106.17± 7.29a
30–60
0.65± 0.25a
104.80± 16.14b
60–100
0.70± 0.16b
100.36± 8.80c
Mtimbwani
0–30
0.83± 0.17a
160.00± 8.65a
30–60
0.92± 0.29b
152.59± 15.42a
b
60–100
0.94± 0.12c
149.32± 5.82b
Corresponding letters a, b, and c denote significant differences atp
<
0.05.


#### 3.3.2. Horizontal Profile of Soil Carbon Stocks

Generally, mangroves of Mtimbwani had significantly (p
=
0.002) higher mean soil carbon density in all zones than Geza. The variations of soil carbon density were significantly different within zones in Geza (p
=
0.037) and Mtimbwani (p
=
0.042). Mid zones had highest mean soil carbon density than lower and upper zones in both sites (Table 4).Table 4
Horizontal profile of soil carbon stocks (mean ± SD).
Zones
Mean carbon (Mg C ha-
1)
Geza
Mtimbwani
Lower
349.75 ± 11.66a
470.50 ± 21.29a
Mid
400.50 ± 31.90a
497.00 ± 48.91b
Upper
183.75 ± 15.26b
418.25 ± 56.51a
Corresponding letters a, b, and c denote significant differences atp
<
0.05.


#### 3.3.3. Total Carbon Stocks

The total ecosystem carbon stock for Mtimbwani was 684.99 Mg C ha−1 and for Geza was 414.64 Mg C ha−1 (Figure 4). Soil carbon pool was the main contributor to total carbon stock, accounting for 67.43% and 75.08% for Mtimbwani and Geza, respectively. In both sites BGC (roots) was less than 10% of the total carbon stock.Figure 4
Contribution of different carbon pools to the total ecosystem carbon stocks.


## 3.1. Forests Structure and Composition

Summary forest structural stand characteristics are presented in Table1. Six species of mangroves were found in both sites.A.marina was the most frequently observed and dominant species in both sites. The least observed species wasB.gymnorrhiza, whileS.alba was the least dominant species. In terms of density,R.mucronata had high relative stem density in Geza (58.9%) and Mtimbwani (41.8%) whereasS.alba was the least contributing species to stem density in both Geza (0.8%) and Mtimbwani (1.1%). Based on their IVs,R.mucronata was the principal species in Geza andA.marina in Mtimbwani. Seedlings and saplings count indicatedR.mucronata andA.marina as the most abundant in both sites.Table 1
Forest stand characteristics for Geza and Mtimbwani mangroves.
Site
Species
Frequency
Density (stems ha-
1)
Basal area(m2 ha-
1)
Seedlings/saplings (counts ha-
1)
Relative
Importance value
Frequency
Density
Dominance
Geza
A. marina
12
406
10.22
778
28.57
23.32
59.37
111.26
B. gymnorrhiza
2
27
0.24
17
4.76
1.52
1.39
7.67
C. tagal
10
231
0.43
531
26.19
13.26
2.49
41.94
R. mucronata
11
1026
5.82
1219
26.19
58.99
33.78
118.97
S. alba
3
13
0.11
0
7.14
0.76
0.67
8.57
X. granatum
3
37
0.4
24
7.14
2.13
2.31
11.59
Total
1740
17.21
2569
100
100
100
300
Mtimbwani
A. marina
12
621
20.78
326
26.09
26.59
68.64
121.31
B. gymnorrhiza
4
80
0.49
35
8.7
3.41
1.63
13.73
C. tagal
8
536
3.51
238
17.39
22.95
11.59
51.93
R. mucronata
11
976
4.49
483
23.91
41.82
14.83
80.56
S. alba
6
27
0.18
18
13.04
1.14
0.58
14.76
X. granatum
5
95
0.83
44
10.87
4.09
2.74
17.7
Total
2334
30.27
1144
100
100
100
300Mangroves of Mtimbwani had high mean stem density across the three zones compared to Geza mangroves (Figure2(a)). In both sites, mid zones had the highest mean stem density (Geza = 2252
±
810 stems ha−1 and Mtimbwani = 2944
±
634 stems ha−1) than other zones but not signifcantly different both between the sites (p
=
0.075) and within zones (p
=
0.233 for Geza and p
=
0.115 for Mtimbwani). For basal areas, differences were significant (p
=
0.041) between sites where Mtimbwani site had large mean basal area than that of Geza (Figure 2(b)). The mid zones represented large basal areas in both Geza (22.25
±
7.1 m2 ha−1) and Mtimbwani (46.10
±
12.1 m2 ha−1) although the differences were not significant among zones in Geza (p
=
0.397) and Mtimbwani (p
=
0.114).Figure 2
Forest structure indicating (a) stems density (b), basal area, (c) number of stumps, and (d) potential regeneration in different zones (mean ± SD).
(a)
(b)
(c)
(d)The mean number of stumps was not significantly different (p
=
0.112) between Geza and Mtimbwani sites (Figure 2(c)). Comparison within zones revealed a significant difference (p
=
0.031) in mean stumps count between the zones of Geza mangroves with highest number in the upper zone (49
±
14 stumps ha−1). In Mtimbwani, the difference in mean stumps count was not significant (p
=
0.539) between the zones although the lower zone had higher count (56
±
22 stumps ha−1). Fewer stumps were counted in mid zones for both sites. Overall, Geza had significantly higher mean number of seedlings than Mtimbwani (p
=
0.003; Figure 2(d)). Mean seedlings and saplings count was not significantly different between zones in both Geza (p
=
0.183) and Mtimbwani (p
=
0.168) although counts were higher in the mid zone of Geza (1107
±
437) and lower zone of Mtimbwani (512
±
91).


## 3.2. Above- and below-Ground Carbon Stocks

The difference in AGC and BGC density was significantly different between sites (p
=
0.041 for Geza and p
=
0.046 for Mtimbwani). There was no significant difference in AGC (p
=
0.214) and BGC (p
=
0.158) between zones of Geza but the difference was significant for Mtimbwani AGC (p
=
0.021) and BGC (p
=
0.038) where mid zone had higher AGC and BGC compared to lower and upper zones (Figure 3). Standing live trees contributed more than 70% of the AGC (Table 2). On average, BGC in Geza and Mtimbwani represented 24.26% and 21.93%, respectively, of vegetative carbon.Table 2
Mean AGC density (Mg C ha−1) for standing live, dead and downed trees (mean ± SD).
Site
Component
Zone
Low
Mid
Upper
Geza
Standing live tree
54.88± 19.16
102.63± 30.67
65.66± 14.45
Standing dead trees
9.02± 3.06
1.32± 0.48
1.16± 0.53
Downed trees
0.23± 0.08
0.03± 0.02
0.22± 0.18
Mtimbwani
Standing live tree
103.44± 39.53
167.84± 57.11
99.50± 45.15
Standing dead trees
7.69± 4.90
133.59± 126.43
10.12± 4.05
Downed trees
0.24± 0.20
0.10± 0.01
0.07± 0.03Figure 3
Mean AGC and BGC density (Mg C ha−1) in different zones for (a) Geza and (b) Mtimbwani (mean ± SD).
(a)
(b)


## 3.3. Mangrove Soil Carbon Stocks

### 3.3.1. Vertical Profile of Carbon Stocks

The general profile in mean soil bulk density and carbon stock in different depths is summarized in Table3. The mean bulk density and soil carbon density significantly decreased with depth in both sites, Geza (p
=
0.038, p
=
0.024) and Mtimbwani (p
=
0.007, p
=
0.040). The top 30 cm had higher carbon density compared to 30–60 cm and 60–100 cm, demonstrating a higher C% in the top layers of the soil.Table 3
Vertical profile of mean soil bulk density and carbon stocks (mean± SD).
Site
Depth (cm)
Mean soil bulk density (g cm-
3)
Mean soil carbon (Mg C ha-
1)
Geza
0–30
0.62± 0.20a
106.17± 7.29a
30–60
0.65± 0.25a
104.80± 16.14b
60–100
0.70± 0.16b
100.36± 8.80c
Mtimbwani
0–30
0.83± 0.17a
160.00± 8.65a
30–60
0.92± 0.29b
152.59± 15.42a
b
60–100
0.94± 0.12c
149.32± 5.82b
Corresponding letters a, b, and c denote significant differences atp
<
0.05.


### 3.3.2. Horizontal Profile of Soil Carbon Stocks

Generally, mangroves of Mtimbwani had significantly (p
=
0.002) higher mean soil carbon density in all zones than Geza. The variations of soil carbon density were significantly different within zones in Geza (p
=
0.037) and Mtimbwani (p
=
0.042). Mid zones had highest mean soil carbon density than lower and upper zones in both sites (Table 4).Table 4
Horizontal profile of soil carbon stocks (mean ± SD).
Zones
Mean carbon (Mg C ha-
1)
Geza
Mtimbwani
Lower
349.75 ± 11.66a
470.50 ± 21.29a
Mid
400.50 ± 31.90a
497.00 ± 48.91b
Upper
183.75 ± 15.26b
418.25 ± 56.51a
Corresponding letters a, b, and c denote significant differences atp
<
0.05.


### 3.3.3. Total Carbon Stocks

The total ecosystem carbon stock for Mtimbwani was 684.99 Mg C ha−1 and for Geza was 414.64 Mg C ha−1 (Figure 4). Soil carbon pool was the main contributor to total carbon stock, accounting for 67.43% and 75.08% for Mtimbwani and Geza, respectively. In both sites BGC (roots) was less than 10% of the total carbon stock.Figure 4
Contribution of different carbon pools to the total ecosystem carbon stocks.


## 3.3.1. Vertical Profile of Carbon Stocks

The general profile in mean soil bulk density and carbon stock in different depths is summarized in Table3. The mean bulk density and soil carbon density significantly decreased with depth in both sites, Geza (p
=
0.038, p
=
0.024) and Mtimbwani (p
=
0.007, p
=
0.040). The top 30 cm had higher carbon density compared to 30–60 cm and 60–100 cm, demonstrating a higher C% in the top layers of the soil.Table 3
Vertical profile of mean soil bulk density and carbon stocks (mean± SD).
Site
Depth (cm)
Mean soil bulk density (g cm-
3)
Mean soil carbon (Mg C ha-
1)
Geza
0–30
0.62± 0.20a
106.17± 7.29a
30–60
0.65± 0.25a
104.80± 16.14b
60–100
0.70± 0.16b
100.36± 8.80c
Mtimbwani
0–30
0.83± 0.17a
160.00± 8.65a
30–60
0.92± 0.29b
152.59± 15.42a
b
60–100
0.94± 0.12c
149.32± 5.82b
Corresponding letters a, b, and c denote significant differences atp
<
0.05.


## 3.3.2. Horizontal Profile of Soil Carbon Stocks

Generally, mangroves of Mtimbwani had significantly (p
=
0.002) higher mean soil carbon density in all zones than Geza. The variations of soil carbon density were significantly different within zones in Geza (p
=
0.037) and Mtimbwani (p
=
0.042). Mid zones had highest mean soil carbon density than lower and upper zones in both sites (Table 4).Table 4
Horizontal profile of soil carbon stocks (mean ± SD).
Zones
Mean carbon (Mg C ha-
1)
Geza
Mtimbwani
Lower
349.75 ± 11.66a
470.50 ± 21.29a
Mid
400.50 ± 31.90a
497.00 ± 48.91b
Upper
183.75 ± 15.26b
418.25 ± 56.51a
Corresponding letters a, b, and c denote significant differences atp
<
0.05.


## 3.3.3. Total Carbon Stocks

The total ecosystem carbon stock for Mtimbwani was 684.99 Mg C ha−1 and for Geza was 414.64 Mg C ha−1 (Figure 4). Soil carbon pool was the main contributor to total carbon stock, accounting for 67.43% and 75.08% for Mtimbwani and Geza, respectively. In both sites BGC (roots) was less than 10% of the total carbon stock.Figure 4
Contribution of different carbon pools to the total ecosystem carbon stocks.


## 4. Discussion

### 4.1. Vegetation Structure and Composition

Six mangrove species were found in the two sites withA.marina andR.mucronata having high relative frequency compared to other species. According to Ksawani et al. [33] the mixture of different species influences the health of forest and enhancement of carbon storage. High frequency ofA.marina andR.mucronata might be attributed to their high regeneration capacity despite these species having high use preferences, particularly the latter [2]. Hamad et al. [34] reported high frequency ofA.marina andR.mucronata in the mangroves of Pemba despite their high use preference for building poles and firewood. The high seedlings and saplings (2569 counts ha−1) implying high regeneration in Geza forest depicted the mangroves to be under human pressure such as cutting for fuelwood, building poles, and conversion for salt pans construction [2]. The higher regeneration might be contributed by opening of canopy gaps that enhances light availability to support seed germination and seedling growth [34, 35].The species richness creates stable ecosystem that is more likely to self-sustain in the events of harvesting pressures [12]. Kirui et al. [36] reported that changes in species richness in Gazi Bay were likely to reduce resilience of mangrove ecosystem and make it vulnerable to natural and anthropogenic activities. Results revealed higher stem density and larger basal area in Mtimbwani compared to Geza (Figures 2(a) and 2(b)). But in general, Geza and Mtimbwani mangroves suffer similar pattern of exploitation, largely the conversion for solar salt production. The effect of human pressure is indicated by the low basal area and stem density. Although the observed stem densities and basal areas were low by standards of a healthy forest according to Bundotich et al. [37], they were high compared to reports from other areas in the region. Kairo et al. [38] reporting for Mida Creek indicated tree density to vary from 1197 to 1585 trees ha−1. Recently Trettin et al. [39] reported that overstory trees in the large Zambezi deltaic mangroves averaged to 2000 trees ha−1. Comparing the two study sites, the observed higher stem density and basal area in Mtimbwani might be attributed to the influence of freshwater from the Ngole River, unlike in Geza where Bongoa River is only seasonal. The freshwater inflow brings nutrient from upstream and creates brackish condition, which is an important aspect for mangroves growth.Mid zones had higher stem density and basal area (Figures2(a) and 2(b)) which may be attributed to the accessibility as it was similarly reported by Mangora and Shalli [40]. Both ends of the forests are easily accessible from land and sea, with relatively less effort required to harvest the products. Large number of stumps observed in the lower and upper zones (Figure 2(c)) demonstrates the exploitation pressures in these zones. Nonetheless, the observed higher basal areas similar to those reported by Trettin et al. [39] in the well-developed mangroves of Zambezi Delta indicate that mangroves in the study sites are still viable. Mangrove species dominance values indicatedA.marina cover large areas in both sites, despite low stem density (Table 1). This might be attributed to the fact that most ofA.marina species were large in size, an indication that the species is less preferred for cutting as compared to the species of the family Rhizophoraceae and therefore has opportunities to grow into large trees. Similar explanation was given forR.mucronata of Mida Creek [38] and Ngomeni forest [37], whereR.mucronata was not a dominant species but had high stem density. Importance values indicatedA.marina andR.mucronata as the most common species in both sites, suggesting their potential for restoration to maximise carbon stocks and sustain wood stocks.


### 4.2. Above- and below-Ground (Roots) Carbon Pools

Standing live trees were the major contributor to AGC followed by standing dead trees and downed trees (Table2). The present reported AGC for Geza and Mtimbwani were comparable to those reported by Kirui et al. [41] and Cohen et al. [42] for Gazi Bay and Vanga which are a little further north of the study sites across the border in Kenya, where carbon density estimate was between 125 and 226 Mg C ha−1, and to those recently reported by Trettin et al. [39] for the large deltaic mangroves of Zambezi in Mozambique which ranged from 111 to 483 Mg ha−1 but exceeded those reported by Cohen et al. [42] for Mtwapa Creek (36.4 Mg C ha−1) and Mida Creek (38.5 Mg C ha−1) in Kenya and Sitoe et al. [43] for Sofala Bay (33.3 Mg C ha−1) in Mozambique. The possible explanation for the variation may be the structural characteristics of the forests [33], whereby mangroves of Mtimbwani and Geza had comparably large basal area and high stem density, which also suggest the low level of degradation. The BGC (roots) was the least pool of carbon after soil and AGC components (Figure 4). Nonetheless, the obtained mean estimates of BGC density were above those reported by Sitoe et al. [43] and Fatoyinbo et al. [44] in Mozambique. These differences in the amount of carbon stored in roots between sites and zones might be explained by differences in nutrient cycling, sediment characteristics, and related growth rates [9, 45].


### 4.3. Soil Carbon Pools

Soil bulk density increased while mean soil carbon decreased with depth (Table3). Similar trends were reported by Lupembe [46] in the Rufiji Delta (Tanzania), Rahman et al. [47] in the Sundarbans (Bangladesh), Sitoe et al. [43] in Sofala Bay (Mozambique), and Stringer et al. [13] in the Zambezi Delta (Mozambique). Increasing bulk density is associated with the increase of soil compactness as depth increases and decrease in concentration of organic matter [13, 32] which also explains high carbon in upper soil layer as concentration of organic matter triggered by biological activity, particularly litter deposition and decomposition on the soil surface, is high [7, 45].The average soil carbon density estimated for up to 1 m deep was311.3
±
87.2 Mg C ha−1 for Geza and 461.9
±
88.1 Mg C ha−1 for Mtimbwani. The estimated soil carbon was comparable to that reported from Madagascar by Jones et al. [16] who found 446.2 Mg C ha−1. The present values of soil carbon density were higher than those reported by Sitoe et al. [43] who found an average of 218.5 Mg C ha−1 in Sofala Bay (Mozambique) but low compared to estimates reported by Kauffman et al. [48] for the Mangroves of Yap Micronesia (724.8 Mg C ha−1) and Donato et al. [7] for the coastal fringe mangroves of Republic of Palau (527.7 Mg C ha−1). Again, site-specific difference in forest structure, composition, history of exploitation, and management regimes may explain the variations and warrant site-specific conservation prescriptions.


### 4.4. Total Ecosystem Carbon Stocks

Total ecosystem carbon stocks were different between sites whereby Geza had 414.6 Mg C ha−1 and Mtimbwani had 684.9 Mg C ha−1. These values were within the range of ecosystem carbon stocks reported from other sites in the region [13, 16, 17, 42, 43] and those reported from other regions [32, 49], which ranges from 55 to 800 Mg C ha−1. Higher ecosystem carbon stocks have been reported from the Indo-West Pacific with values ranging from 830 to 1131 Mg C ha−1 [7, 48]. In terms of the pattern of carbon density with tidal gradient, mid zones had large stocks. This may be attributed to forest density and basal area, similar to the observation by Stringer et al. [13] in the Zambezi Delta. The observed low density of stumps in the mid zones, implying low deforestation, further supports this. Other reports suggest that differences in mangrove carbon stocks between sites and zones may as well be due to differences in forest age, forest conservation and management status, soil depth, and soil water content [13, 16, 27], factors which need further research for the study sites. The influence of these local factors demonstrates that estimation of carbon stocks in mangroves cannot be generalised, especially for MRV and their potential in climate change mitigation and adaptation.Despite the differences in mangrove carbon stocks for different ecosystems across the world, the relative contribution of various ecosystem components shows a similar trend. Above-ground carbon density contributed 18.87% and 25.41% of total ecosystem carbon estimates in Geza and Mtimbwani, respectively. Kauffman et al. [48] reported that AGC accounts for about 18 to 25% of total ecosystem carbon. Observations on BGC density were also coherent with values reported from other places, which range from 6 to 23% [7]. Similar to reports from other assessments [7, 13, 43, 48], the soil carbon pool was dominant representing 75.08% and 67.43% of total ecosystem carbon stock for Geza and Mtimbwani mangroves, respectively. Soil carbon values of the present study and those from other areas demonstrate the role of mangrove soil as an important carbon sink. This carbon pool has a significant contribution to primary productivity and health of the forests [9] that can ensure mangroves are able to sustain in providing other ecosystem services that are valuable to local communities.


### 4.5. Potential for Carbon Markets

Detailed analysis of the socioeconomic factors is not presented here but indicated that communities in Geza and Mtimbwani are substantially dependent on ecosystem services provided by mangroves. Similar to other reports [40], low income and low level of environmental education contributed to rapid degradation of mangrove forests and eventually structure and carbon stocks. Most households in Geza and Mtimbwani are engaged in harvesting and selling of mangrove products to maximise their income in addition to household consumption [50]. Mangroves are widely used as cheap construction materials to build and maintain their homes and a source of household fuelwood. Mangrove areas are also converted to solar salt pans. The observed large amount of carbon stored in the two mangrove communities is vulnerable to loss from such human pressures that are often coupled with the natural phenomena such as sea level rise, disrupted water budgets, and increased salinity. Bhomia et al. [14] suggested that carbon stocks >500 Mg C ha−1 confirm the potential of mangroves as sinks and that any loss in such mangrove forests exposes them to significant carbon sources. These data sets provide a clear message that concerted efforts in conservation and management of mangrove areas are warranted to minimize the emissions. The quantification of fairly large carbon stocks in the mangroves of Mtimbwani and Geza provides key baselines for informed policy and planning decisions to safeguard mangroves and restore degraded areas.Opportunities exist because of the growing realization of the impacts of mangrove degradation coupled with those of climate change that shape the community willingness for conservation initiatives to address mangrove degradation. Emerging market-based conservation through carbon projects (carbon credits) is a potential opportunity to incentivize Geza and Mtimbwani communities, providing them with necessary resources needed to sustain conservation and the carbon sink function of mangroves. This is further supported by the observed vegetative structural indices that proved that mangroves of Geza and Mtimbwani are still viable for self-sustenance given appropriate conservation measures. Incentivized communities would conduct self-enforcement in restoring and protecting the mangroves as they would devise alternative sources of wood products for local consumption through tree planting of fast growing terrestrial based species.


## 4.1. Vegetation Structure and Composition

Six mangrove species were found in the two sites withA.marina andR.mucronata having high relative frequency compared to other species. According to Ksawani et al. [33] the mixture of different species influences the health of forest and enhancement of carbon storage. High frequency ofA.marina andR.mucronata might be attributed to their high regeneration capacity despite these species having high use preferences, particularly the latter [2]. Hamad et al. [34] reported high frequency ofA.marina andR.mucronata in the mangroves of Pemba despite their high use preference for building poles and firewood. The high seedlings and saplings (2569 counts ha−1) implying high regeneration in Geza forest depicted the mangroves to be under human pressure such as cutting for fuelwood, building poles, and conversion for salt pans construction [2]. The higher regeneration might be contributed by opening of canopy gaps that enhances light availability to support seed germination and seedling growth [34, 35].The species richness creates stable ecosystem that is more likely to self-sustain in the events of harvesting pressures [12]. Kirui et al. [36] reported that changes in species richness in Gazi Bay were likely to reduce resilience of mangrove ecosystem and make it vulnerable to natural and anthropogenic activities. Results revealed higher stem density and larger basal area in Mtimbwani compared to Geza (Figures 2(a) and 2(b)). But in general, Geza and Mtimbwani mangroves suffer similar pattern of exploitation, largely the conversion for solar salt production. The effect of human pressure is indicated by the low basal area and stem density. Although the observed stem densities and basal areas were low by standards of a healthy forest according to Bundotich et al. [37], they were high compared to reports from other areas in the region. Kairo et al. [38] reporting for Mida Creek indicated tree density to vary from 1197 to 1585 trees ha−1. Recently Trettin et al. [39] reported that overstory trees in the large Zambezi deltaic mangroves averaged to 2000 trees ha−1. Comparing the two study sites, the observed higher stem density and basal area in Mtimbwani might be attributed to the influence of freshwater from the Ngole River, unlike in Geza where Bongoa River is only seasonal. The freshwater inflow brings nutrient from upstream and creates brackish condition, which is an important aspect for mangroves growth.Mid zones had higher stem density and basal area (Figures2(a) and 2(b)) which may be attributed to the accessibility as it was similarly reported by Mangora and Shalli [40]. Both ends of the forests are easily accessible from land and sea, with relatively less effort required to harvest the products. Large number of stumps observed in the lower and upper zones (Figure 2(c)) demonstrates the exploitation pressures in these zones. Nonetheless, the observed higher basal areas similar to those reported by Trettin et al. [39] in the well-developed mangroves of Zambezi Delta indicate that mangroves in the study sites are still viable. Mangrove species dominance values indicatedA.marina cover large areas in both sites, despite low stem density (Table 1). This might be attributed to the fact that most ofA.marina species were large in size, an indication that the species is less preferred for cutting as compared to the species of the family Rhizophoraceae and therefore has opportunities to grow into large trees. Similar explanation was given forR.mucronata of Mida Creek [38] and Ngomeni forest [37], whereR.mucronata was not a dominant species but had high stem density. Importance values indicatedA.marina andR.mucronata as the most common species in both sites, suggesting their potential for restoration to maximise carbon stocks and sustain wood stocks.


## 4.2. Above- and below-Ground (Roots) Carbon Pools

Standing live trees were the major contributor to AGC followed by standing dead trees and downed trees (Table2). The present reported AGC for Geza and Mtimbwani were comparable to those reported by Kirui et al. [41] and Cohen et al. [42] for Gazi Bay and Vanga which are a little further north of the study sites across the border in Kenya, where carbon density estimate was between 125 and 226 Mg C ha−1, and to those recently reported by Trettin et al. [39] for the large deltaic mangroves of Zambezi in Mozambique which ranged from 111 to 483 Mg ha−1 but exceeded those reported by Cohen et al. [42] for Mtwapa Creek (36.4 Mg C ha−1) and Mida Creek (38.5 Mg C ha−1) in Kenya and Sitoe et al. [43] for Sofala Bay (33.3 Mg C ha−1) in Mozambique. The possible explanation for the variation may be the structural characteristics of the forests [33], whereby mangroves of Mtimbwani and Geza had comparably large basal area and high stem density, which also suggest the low level of degradation. The BGC (roots) was the least pool of carbon after soil and AGC components (Figure 4). Nonetheless, the obtained mean estimates of BGC density were above those reported by Sitoe et al. [43] and Fatoyinbo et al. [44] in Mozambique. These differences in the amount of carbon stored in roots between sites and zones might be explained by differences in nutrient cycling, sediment characteristics, and related growth rates [9, 45].


## 4.3. Soil Carbon Pools

Soil bulk density increased while mean soil carbon decreased with depth (Table3). Similar trends were reported by Lupembe [46] in the Rufiji Delta (Tanzania), Rahman et al. [47] in the Sundarbans (Bangladesh), Sitoe et al. [43] in Sofala Bay (Mozambique), and Stringer et al. [13] in the Zambezi Delta (Mozambique). Increasing bulk density is associated with the increase of soil compactness as depth increases and decrease in concentration of organic matter [13, 32] which also explains high carbon in upper soil layer as concentration of organic matter triggered by biological activity, particularly litter deposition and decomposition on the soil surface, is high [7, 45].The average soil carbon density estimated for up to 1 m deep was311.3
±
87.2 Mg C ha−1 for Geza and 461.9
±
88.1 Mg C ha−1 for Mtimbwani. The estimated soil carbon was comparable to that reported from Madagascar by Jones et al. [16] who found 446.2 Mg C ha−1. The present values of soil carbon density were higher than those reported by Sitoe et al. [43] who found an average of 218.5 Mg C ha−1 in Sofala Bay (Mozambique) but low compared to estimates reported by Kauffman et al. [48] for the Mangroves of Yap Micronesia (724.8 Mg C ha−1) and Donato et al. [7] for the coastal fringe mangroves of Republic of Palau (527.7 Mg C ha−1). Again, site-specific difference in forest structure, composition, history of exploitation, and management regimes may explain the variations and warrant site-specific conservation prescriptions.


## 4.4. Total Ecosystem Carbon Stocks

Total ecosystem carbon stocks were different between sites whereby Geza had 414.6 Mg C ha−1 and Mtimbwani had 684.9 Mg C ha−1. These values were within the range of ecosystem carbon stocks reported from other sites in the region [13, 16, 17, 42, 43] and those reported from other regions [32, 49], which ranges from 55 to 800 Mg C ha−1. Higher ecosystem carbon stocks have been reported from the Indo-West Pacific with values ranging from 830 to 1131 Mg C ha−1 [7, 48]. In terms of the pattern of carbon density with tidal gradient, mid zones had large stocks. This may be attributed to forest density and basal area, similar to the observation by Stringer et al. [13] in the Zambezi Delta. The observed low density of stumps in the mid zones, implying low deforestation, further supports this. Other reports suggest that differences in mangrove carbon stocks between sites and zones may as well be due to differences in forest age, forest conservation and management status, soil depth, and soil water content [13, 16, 27], factors which need further research for the study sites. The influence of these local factors demonstrates that estimation of carbon stocks in mangroves cannot be generalised, especially for MRV and their potential in climate change mitigation and adaptation.Despite the differences in mangrove carbon stocks for different ecosystems across the world, the relative contribution of various ecosystem components shows a similar trend. Above-ground carbon density contributed 18.87% and 25.41% of total ecosystem carbon estimates in Geza and Mtimbwani, respectively. Kauffman et al. [48] reported that AGC accounts for about 18 to 25% of total ecosystem carbon. Observations on BGC density were also coherent with values reported from other places, which range from 6 to 23% [7]. Similar to reports from other assessments [7, 13, 43, 48], the soil carbon pool was dominant representing 75.08% and 67.43% of total ecosystem carbon stock for Geza and Mtimbwani mangroves, respectively. Soil carbon values of the present study and those from other areas demonstrate the role of mangrove soil as an important carbon sink. This carbon pool has a significant contribution to primary productivity and health of the forests [9] that can ensure mangroves are able to sustain in providing other ecosystem services that are valuable to local communities.


## 4.5. Potential for Carbon Markets

Detailed analysis of the socioeconomic factors is not presented here but indicated that communities in Geza and Mtimbwani are substantially dependent on ecosystem services provided by mangroves. Similar to other reports [40], low income and low level of environmental education contributed to rapid degradation of mangrove forests and eventually structure and carbon stocks. Most households in Geza and Mtimbwani are engaged in harvesting and selling of mangrove products to maximise their income in addition to household consumption [50]. Mangroves are widely used as cheap construction materials to build and maintain their homes and a source of household fuelwood. Mangrove areas are also converted to solar salt pans. The observed large amount of carbon stored in the two mangrove communities is vulnerable to loss from such human pressures that are often coupled with the natural phenomena such as sea level rise, disrupted water budgets, and increased salinity. Bhomia et al. [14] suggested that carbon stocks >500 Mg C ha−1 confirm the potential of mangroves as sinks and that any loss in such mangrove forests exposes them to significant carbon sources. These data sets provide a clear message that concerted efforts in conservation and management of mangrove areas are warranted to minimize the emissions. The quantification of fairly large carbon stocks in the mangroves of Mtimbwani and Geza provides key baselines for informed policy and planning decisions to safeguard mangroves and restore degraded areas.Opportunities exist because of the growing realization of the impacts of mangrove degradation coupled with those of climate change that shape the community willingness for conservation initiatives to address mangrove degradation. Emerging market-based conservation through carbon projects (carbon credits) is a potential opportunity to incentivize Geza and Mtimbwani communities, providing them with necessary resources needed to sustain conservation and the carbon sink function of mangroves. This is further supported by the observed vegetative structural indices that proved that mangroves of Geza and Mtimbwani are still viable for self-sustenance given appropriate conservation measures. Incentivized communities would conduct self-enforcement in restoring and protecting the mangroves as they would devise alternative sources of wood products for local consumption through tree planting of fast growing terrestrial based species.


## 5. Conclusion

This study presented findings that demonstrate that forest structure and ecosystem carbon densities are key elements to voluntary carbon market schemes. Inventory of both mangrove vegetation and soil carbon demonstrated comparatively higher carbon stocks in the study sites and therefore reflecting their potential for carbon markets. Geza and Mtimbwani communities have experience in collaborative arrangements in mangroves management from the former Tanga Coastal Zone Conservation and Development programme [51]. The estimated carbon stocks have provided baseline data for future MRV that is important in supporting carbon projects. Based on present findings and the potential for carbon projects, it is recommended that (a) permanent plots should be established for MRV of additionality in terms of health and carbon stocks; (b) vulnerability assessment and monitoring should be done to assess mangrove cover changes over time and predict extents of human impacts on mangroves; (c) collaborative management should be strengthened by harmonising rules and regulations across stakeholders; (d) detailed studies should be done to investigate whether existing governance mechanism can favour initiation of carbon market schemes.



---

*Source: 2068283-2016-09-28.xml*