Variation in Seagrass Carbon Stocks Between Tropical Estuarine and Marine Mangrove-Fringed Creeks

Seagrass and associated blue carbon ecosystems are important carbon sinks, and hence understanding their spatial and temporal variability is vital in appreciating their potential roles in climate change mitigation and adaptation. The Indo-Pacific region has the highest seagrass biodiversity, yet little focus has been made to compare seagrass habitat extent and carbon dynamics with their temperate counterparts. The present study assessed habitat characteristics and seagrass species distribution, diversity, and carbon storage in Eastern (marine) and Western (estuarine) mangrove-fringed creeks of Gazi Bay, Kenya. Data on species composition, canopy cover, biomass, and sediment organic carbon were collected in 80 plots of 0.25 × 0.25 m laid along transects established perpendicular to the waterline. Five species formation, viz., Thalassia hemprichii, Cymodocea rotundata, Cymodocea serrulata, Enhalus acoroides, and Thalassidendron ciliatum, were encountered as either single or mixed stands. There was a significant difference in total seagrass biomass between creeks (p < 0.01), with the Eastern creek recording a mean of 10.2 ± 0.6 Mg C ha–1 while the Western creek recording 4.3 ± 0.3 Mg C ha–1. In addition, sediment carbon to 1-m depth varied significantly (p < 0.01) between species in the two creeks and ranged from 98 to 302 Mg C ha–1, with the Eastern and Western creeks recording means of 258 ± 90 and 107 ± 21 Mg C ha–1, respectively. The total carbon stock from 50 ha of seagrasses in the Eastern creek was 13,420 Mg C, whereas in the 70 ha of the Western creek it was 7,769 Mg C. The study shows that seagrass community attributes such as species composition and productivity can vary dramatically over a small spatial extent due to differences in biophysical conditions and caution estimations of site-specific carbon stocks using generalized global values.


INTRODUCTION
Conservation of the blue carbon sink such as mangroves, tidal salt marshes, and seagrasses is important because of their contribution to climate change mitigation and adaptations. Despite their role in capturing and storing huge carbon stocks (Fourqurean et al., 2012;Rogelj et al., 2016), less information is available from blue carbon ecosystems when compared to their terrestrial counterparts (Huxham et al., 2018). The natural and human pressures on these habitats expose them to loss and degradation, leading to the release of CO 2 back into the atmosphere and thereby hindering their role in carbon storage Howard et al., 2018).
Seagrass meadows are keystone ecosystems along the coasts providing a range of benefits to humans and other organisms, among them regulatory, support, provisioning, and cultural services (Hejnowicz et al., 2015;Nordlund et al., 2018). They contribute to nutrient cycling (Costanza et al., 1997) and sediment stabilization and provide nursery to over one fifth of the world's 25 largest fisheries (Unsworth et al., 2019). They are efficient in carbon capture and storage (Duarte and Prairie, 2005;Gedan et al., 2009), hence necessitating their inclusion in carbon offsetting schemes (Howard et al., 2017;Herr et al., 2017). Despite their value, seagrasses across the world are threatened by a combination of both human and natural factors (Cabaço et al., 2008;Mazarrasa et al., 2017;Githaiga et al., 2019), contributing to the global estimated loss of 1.5% per year (Waycott et al., 2009;Pendleton et al., 2012).
Seagrass species thrive in varying habitats ranging from estuaries, mud flats, sheltered bays, and lagoons across wide latitudinal distributions ranging from temperate to tropical latitudes with varying geomorphological, hydrological, and ecological conditions (Lyimo et al., 2018;Bulmer et al., 2020). These affect the distribution and abundance of the seagrass species in their physiological and survival processes (Lyimo, 2016;Gullström et al., 2017;Mazarrasa et al., 2018). The integrity of these marine habitats is also likely to be affected by biophysical factors such as water currents, day length, epiphytes, wave action, temperature, light availability, salinity, and substrate type and depth (Hemminga and Duarte, 2000; Green and Short, 2003). Additionally, anthropogenic stressors upstream or along the coastline, such as effluent discharge and agricultural activities, facilitate sediment and nutrient loading, thereby affecting seagrass health status. Direct effects are on structural traits such as cover, density, and biomass accumulation that subsequently influence seagrass sediment organic carbon (C org ) sequestration and deposition in the soil through three major processes: biomass accumulation and productivity, allochthonous C org load in the water column, and the burial rates in sediment (Serrano et al., 2016a,b;Mazarrasa et al., 2018). Largest of the carbon pools is the sediment organic carbon comprising autochthonous and allochthonous carbon which may be deposited for millennia in sediment, though in some conditions get remineralized, leading to further CO 2 emissions (Kennedy et al., 2010). Seagrass meadows have a larger global spatial extent when compared to other blue carbon sinks (Mcleod et al., 2011). They occupy a global coverage of 3,000 km 2 with more than 70 species (Small and Nicholls, 2003;Short et al., 2011). Species vary in abundance and distribution across latitudes, with the tropical bioregions having the highest diversity at 24 out of the known 60 species globally (Short et al., 2007). Compared to other blue carbon sinks, there is less research focusing on seagrass carbon dynamics, leading to uncertainties in the seagrass habitat extent and contribution to the blue carbon budget. In the African coastline for example, there is still a huge paucity of information regarding seagrass habitat extent, quantitative carbon estimates, and habitat influence despite the extensive meadows (Githaiga et al., 2016(Githaiga et al., , 2017. This study assessed habitat characteristics and the carbon stocks of two seagrass meadows in Gazi Bay, Kenya, that lies within the East African coastline. We hypothesized that the seagrass distribution, abundance, and carbon stocks differ significantly between the estuarine and marine seagrass meadows. The bay is known to host a high seagrass diversity, with up to 12 species described in the Western Indian Ocean (WIO) region occurring here (Ochieng and Erftemeijer, 2003). Our study on habitat variations and influence on seagrass carbon stocks is of significance as it provides a better understanding of tropical seagrass ecosystems for comparison with their temperate and subtropical counterparts. Additionally, as the first study assessing habitat influence on carbon storage in seagrass within the African tropics, it provides insights for the enhanced management of these critical ecosystems.

The Study Site
The study was carried out in the two creeks of Gazi Bay -Eastern creek (04.41610 • S, 039.52610 • E) and Western creek (04.41661 • S, 039.51253 • E) -in Kenya between May and September 2018. The Eastern creek lies on the eastern part of the bay bordering Mikoko Pamoja, a community-managed managed mangrove carbon offsetting project area in Gazi. It is characterized by shallow subtidal areas and channels that intersect the intertidal flats (Hemminga et al., 1994). Some areas have rocky substrate and macroalgae lying in close association with the seagrass species. The area of the eastern creek is estimated at 50 ha and lacks direct influence of freshwater (Figure 1). Key mangrove species fringing the eastern creeks are the Sonneratia alba and Rhizophora mucronate (Supplementary Figure S1).
The Western creek, on the other hand, lies on the northwestern part of the bay outside the Mikoko Pamoja project area. This creek covers about 70 ha and is open to the seasonal river channel, River Kidogoweni, allowing inflow of freshwater during the rainy season. The creek is shallow with exposed sand banks during low tides (Figure 1). Due to proximity to human settlement in Gazi, the seagrasses in the western creek are highly influenced by fishing and tourism activities in the area (Supplementary Figure S2).

Sampling Design
The sampling approach adopted the protocol recommended by the Coastal Blue Carbon manual and the Intergovernmental Panel for Climate Change (IPCC, 2014). In each creek, eight transects at intervals of 100 m were initially laid perpendicular to the creek from the upper zone, adjacent to the mangrove ecosystem, all the way to the seagrass meadows near the open waters of the bay. Five quadrats measuring 0.5 × 0.5 m were laid along the transects at intervals of 25 m to determine the distribution and abundance of seagrass species in each creek. In assessing the spatial variability in sediment organic carbon and biomass between creeks, two sub-transects were laid within Thalassia hemprichii, Cymodocea rotundata, Cymodocea serrulata and mixed species stands in Western creek; and within Enhalus acoroides, Thalassodendron ciliatum, C. rotundata and mixed stands in the Eastern creek based on dominance. Five quadrates measuring 0.5 × 0.5 m were then laid along each sub-transect at intervals of 25 m. A total of 40 quadrats were established in each creek and 480 samples obtained for laboratory processing. Data collection was done during low spring tides, when the seagrass beds are easily accessible by foot.

Measurement of Physiochemical Parameters
Measurements for the analysis of the physiochemical parameters were done in situ, within the plots. Total dissolved solids (in milligrams per liter), conductivity (in microseconds per centimeter), water temperature (in degree Celsius), and pH were measured using the HANNA Combo PH and EC multimeter Hi-98129. The probe was immersed to 5 cm in water for every reading. Salinity (in parts per thousand) was measured in the creeks using a refractometer, while depth was measured using a tape.

Determination of the Species Type, Cover, Canopy Height, and the Shoot Density
Seagrass species were identified in situ using field manuals appropriate for the region included (Richmond, 2002). The percentage canopy cover was determined by visual estimate. Canopy height estimation involved measuring the total blade length of 10% of individual shoots randomly selected from the total number within the quadrat and calculating their mean heights. Shoot density was determined by counting the entire shoots within the quadrats and expressed in square meters (Howard et al., 2014).

Estimation of Above-and Belowground Biomass
Aboveground seagrass materials were harvested within the 0.25m 2 quadrats and cleaned using freshwater. Sorting was done followed by scraping using a razor blade to remove the epiphytes. The fronds were then washed in 10% HCl, rinsed to remove the calcareous materials, and then dried in an oven at 60 • C for 72 h.
Estimation of the belowground biomass was done by taking four cores from each of the 0.25 m 2 quadrats using a Russian peat sampler. The cores were washed through a 500 µm sieve. Cleaned material was sorted into necromass, rhizomes, and roots and then dried in an oven for 72 h at 60 • C. Total biomass was determined by multiplying the biomass with the carbon conversion factor of 0.34 and extrapolation done per hectare (IPCC, 2014).

Determination of Sediment C org
Two sediment cores, each having an extension of 50 cm, were obtained using a peat sampler within the quadrats. Since the relative content of the sedimentary carbon may be influenced by the sediment compaction, the difference in length from the upper part of the core to the surface of the sediment, outside and inside the corer, was assessed when passing it down into the sediment (IPCC, 2014). Since we intended to analyze for sediment C org , plant material, infauna, and larger shells were removed from the sediment samples. The samples were sliced into subsections of 5 cm and oven-dried at 60 • C for 72 h in the laboratory to attain a constant weight. The top 50 cm estimates of the sediment were extrapolated to 1 m in estimating the carbon stocks (Howard et al., 2014).

Measurement of Percentage Organic Matter
The organic matter content was determined using the loss of ignition (LOI) technique at 450 • C for 6 h. The percentage organic matter (OM) was calculated using the formula (1) where % LOI is the percentage loss of ignition.
Depending on the organic matter in each of the sample, the sediment C org values were calculated using the following relations (Howard et al., 2014): where % LOI is the percentage loss of ignition and % C org is the percentage sediment organic carbon.

Calculation of the Total Carbon Stocks
In order to determine the total carbon stocks in the two creeks, the soil dry bulk density (DBD) (sediment dry weight per unit volume) for each of the subsections per core was calculated using the formula: DBD (g/cm 3 ) = Dry weight/Original volume of sediment (4) where DBD is the dry bulk density. The soil carbon density was then calculated as follows: Soil carbon density (g/cm 3 ) = dry bulk density (g/cm 3 ) The amount of carbon in each core section was then calculated by multiplying the soil carbon density (in grams per cubic centimeter) by 5 cm (the thickness interval). The amount of carbon per core was then obtained by the addition of the carbon amount in all the sections of one core and the value extrapolated to 1 m (Howard et al., 2014;Githaiga et al., 2017). An average amount of carbon was then obtained for the two cores in each quadrat for sediment C org .
Total core carbon (Mg C/hectare − cm) = Summed corecarbon(g/cm 3 ) × (1 Mg/1, 000, 000 g)× (100, 000, 000 cm 2 /1 hectare) (6) To determine the total amount of carbon in the two meadows, the average carbon stocks from each pool was summed up for each creek and the sum multiplied by the area of the creek.
The variabilities and errors associated with the measurements were determined by calculating the standard deviation for each pool in each creek and multiplying by the area of the corresponding creek.

Data Analysis
Species richness, diversity, and abundance were determined by the Shannon-Wiener diversity index (H) and index of evenness.
To analyze for statistically significant variations between creeks and among species, data were tested for normality and homogeneity of variance and log-transformed, where assumptions of normality were not found. One-way ANOVA was used to test for variations in the aboveground, belowground, and sediment carbon among species. Where significant variations were detected, Tukey's post hoc test was used to compare the means. Two-sample t-test was used to test for significant variations between the creeks. In all the statistical tests, the significant level was set at α = 0.05. A generalized linear model (GLM) with Gaussian error distribution was fitted to establish the relationship between sediment carbon and the environmental variables in the Eastern and Western creeks.

Assessment of Physical-Chemical Parameters in Eastern and Western Creeks
Most of the physical-chemical parameters measured were found to vary significantly between the creeks. Depth ranged between 0.20 and 0.80 m in the Western creek (mean = 0.53 ± 0.31 m), while the values for the Eastern creek ranged between Turbidity values were higher in the Western creek (mean = 0.52 ± 0.07 mg/L, range = 0.35-0.68 mg/L) than in the Eastern creek (mean = 0.32 ± 0.06 mg/L, range = 0.2-0.41 mg/L). The means recorded for the two creeks revealed a significant difference (t = −13.6, df = 76, p < 0.05) at the 95% confidence level. The pH for the Western creek ranged between 7.4 and 8.3, while the pH values for the Eastern creek ranged between 7.8 and 8.0 (Table 1). There was significant variation in the water temperature (t = −15.66, df = 51, p < 0.05) at the 95% confidence level, with the Western creek having higher values (mean = 31.9 ± 0.83 • C, range = 30.3-33.2 • C) while the Eastern creek had a mean of 29.7 ± 0.40 • C, range = 29.0-30.2 • C.

Seagrass Species Diversity and Distribution Between the Creeks
Nine seagrass species out of the 12 species recorded in the entire bay were encountered in the creeks ( Table 2). Six species -C. rotundata, Halodule uninervis, T. hemprichii, C. serrulata, Syringodium isoetifolium, and E. acoroides -were observed in both creeks. However, Halophila stipulacea and Halophila ovalis were only observed in the Western creek, while T. ciliatum was observed in the Eastern creek. The most common species in the Eastern creek was T. ciliatum, while C. rotundata was the most common in the Western creek ( Table 2). The diversity of seagrass species (expressed as the Shannon-Wiener index, H) was higher in the Eastern creek (H = 1.71) than in the Western creek (H = 1.67). Similarly, the Eastern creek had a higher evenness value (E = 0.88) than the Western creek (E = 0.80). However, there was no significant difference in seagrass species diversity between the two creeks (t = 0.11, df = 13, p = 0.46).

Comparison of Sediment Organic Carbon Among the Species and in the Creeks
Soil properties (dry bulk density and percentage organic carbon) significantly varied between the creeks (p < 0.001). These differences strikingly persisted to 50 cm depth (Figure 3). The Western creek recorded a higher mean dry bulk density at 1.49 ± 0.38 g cm −3 , while the Eastern creek recorded 0.87 ± 0.32 g cm −3 . On the other hand, the Eastern creek recorded a higher organic carbon concentration at 3.56%, while the Western creek recorded an average of 0.76%.
Sediment organic carbon differed among the species, with the highest values recorded in the mixed stands of the Eastern creek at 302.45 ± 43.23 Mg C ha −1 while C. rotundata of the Western creek had the lowest at 97.57 ± 7.74 Mg C ha −1 (mean ± 95% CI). One-way ANOVA revealed a significant FIGURE 2 | Variations in biomass between the Eastern and Western creeks (mean ± 95% CI). Letters indicate significant differences in biomass (p < 0.01). difference among the species C org [F (7, 79) = 20.28, p < 0.001]. Sediment C org also varied between the creeks, with the Eastern creek recording higher values at 258.21 ± 90.12 Mg C ha −1 (range = 117.85-544.65 Mg C ha −1 ) while the Western creek recorded 106.66 ± 21.36 Mg C ha −1 (range = 67.25-160.48 Mg C ha −1 ) (Figure 4). The two-sample t tests revealed a significant difference in the sediment C org of the two creeks (t = −10.86, df = 44, p < 0.001).
FIGURE 4 | Variations in sediment C org to 1 m depth among species in the Eastern and Western creeks (mean ± 95% CI). Sed Carbon, sediment organic carbon; E, Eastern; W, Western. Letters a, b and c indicate significant difference in sediment carbon among species (p < 0.001).

Environmental Factors Explaining C org in Eastern and Western Creeks
A GLM, with Gaussian error distribution, was fitted separately for each creek with sediment carbon as the continuous response variable. In the Eastern creek, GLM was constructed for seagrass species, belowground biomass, aboveground biomass, temperature, salinity, depth, shoot density, and seagrass shoot height as the explanatory variable. Turbidity was found to correlate above 0.7 with depth while temperature correlated to pH. Using the Akaike's information criterion (AIC), species (T. ciliatum, p < 0.01) and temperature (p < 0.03) had a significant influence on sediment carbon. The residuals of the best model (AIC = 22) are presented in Supplementary Figure S3.
In the Western creek, a correlation above 0.7 was found between depth and turbidity. A GLM was fitted for aboveground biomass, belowground biomass, species, shoot density, shoot height, pH, salinity, temperature, and turbidity. Aboveground biomass (p < 0.03), turbidity (p < 0.05), total shoot density (0.008), and species (C. serrulate: p = 0.03; mixed species: 0.05; and T. hemprichii: p < 0.06) were found to significantly influence sediment carbon stocks in the Western creek. The residuals of the model with the lowest AIC of 336 are shown in Supplementary Figure S4.

Sediment C org Stock Variations Between the Creeks
The current study compared seagrass habitats in the Eastern and Western creeks of Gazi Bay, Kenya, in terms of their species composition, distribution and the carbon stocks in the above-and belowground components. The Eastern creek was dominated by T. ciliatum and E. acoroides that are characterized by long leaves and have higher shoot density, signifying more sediment trapping. Higher values for rhizomes, roots, and necromass were recorded in this creek. On the other hand, Western creek was dominated by small-leaved and shortlived species such as C. rotundata, which could be attributed to human perturbations and sediment erosion. It was also characterized by lower diversity, abundance, and shoot density of the seagrass communities. Such meadows are known to have lower dissolved organic carbon (DOC) supply, which is likely to affect the percentage of organic carbon through the reduced decomposition of the refractory organic compounds. This difference in species composition and shoot density between the creeks may therefore partly explain the differences in the concentrations of organic matter between them Enríquez et al., 2019), resulting in variations in sediment carbon between the creeks.
The study found that the Eastern creek had higher seagrass biomass and sediment C org values for the top 1 m when compared to the Western creek and the values obtained by Githaiga et al. (2017) in the open waters of the bay (236 ± 24 Mg C ha −1 ). Similarly, the average C org values obtained for the two creeks are above the mean derived from the global seagrass data set at 166 Mg C ha −1 , although they are within the global range of 115.5-829.2 Mg C ha −1 (Fourqurean et al., 2012). Recent studies have indicated variations in organic carbon stock estimates and other seagrass community attributes such as species over such smaller spatial scales (Lavery et al., 2013;Röhr et al., 2016).
In the two creeks, sediment C org stocks varied with species. This was the same case in biomass and other structural attributes such as shoot density and canopy heights, which could explain the variations in the carbon storage capacities. Large-leaved species are efficient in sequestering carbon as well as in trapping allochthonous carbon (Agawin and Duarte, 2002;Mellors et al., 2002). They also tend to establish more permanent and stable beds than the ephemeral and small species due to their higher resistance to hydrodynamic energy, hence the high carbon stocking (Ondiviela et al., 2014;Huxham et al., 2018). Additionally, rhizomes, roots, necromass, and DOC differ among species and contribute substantially to the primary production, resulting in more C org in complex meadows with high shoot density and belowground biomass than in sparse meadows (Kaldy et al., 2006;Armitage and Fourqurean, 2016). The burial rates of both autochthonous and allochthonous organic carbon are estimated at 50% originating from other ecosystems (Kennedy et al., 2010;Duarte and Krause-Jensen, 2017). Past studies have shown that C org from the adjacent mangrove ecosystems were transported to the seagrass meadows in both creeks (Bouillon et al., 2007;Huxham et al., 2018).
Nine seagrass species out of the 12 recorded in Gazi Bay were encountered in the two creeks during the study (Ochieng and Erftemeijer, 2003), with the Eastern creek having higher seagrass species diversity and abundance. Turbidity and depth were found to influence the species distribution and biomass between the creeks. The dominance of T. ciliatum in the Eastern creek can be attributed to the stable state of this creek. It is usually intolerant of any freshwater input and occupies hard and rocky substrates (Waycott et al., 2004). On the other hand, the pioneer and short-lived species H. ovalis and H. stipulacea were recorded only in the Western creek. The availability of an inflow channel (River Kidogoweni) in the Western creek and the shifting sand banks make this area highly disturbed, which could explain the existence of H. ovalis and H. stipulacea and the higher abundance of C. rotundata, which are usually pioneer species and occupy disturbed areas (Noel et al., 2012). Additionally, destructive fishing practices such as boat dredging were observed in this area and could contribute to the disturbance. The higher turbidity recorded in the Western creek may also account for the low seagrass diversity in this creek. Reduction of irradiance as a result of turbidity hinders the photosynthetic process, thereby affecting seagrass morphology and growth (Lee et al., 2007;Lavery et al., 2013).

Aboveground and Belowground Biomass
The total biomass estimates for the current study were within range when compared to other published data. The biomass values compare well with the 5.9 ± 0.9 Mg C ha −1 recorded by Githaiga et al. (2017) in the open water of the bay, although higher than the global mean of 2.51 ± 0.49 Mg C ha −1 (Fourqurean et al., 2012). The biomass values significantly varied between the creeks and were highly correlated to species and plant shoot density. Seagrass meadows with larger species have higher aboveground and belowground biomass compared to smaller species (Duarte and Chiscano, 1999). The variations in biomass among species and between creeks enhance the meadows' efficiency to sequester autochthonous C org . The estuarine environment of the Western creek is likely to subject this seagrass habitat to periodic freshwater inundation, leading to low light levels and burial by suspended sediments, runoff turbidity, and nutrient loading (McDonald et al., 2016). This may lower the plant growth rate and shoot density in this creek, hence lowering biomass (Peralta et al., 2008). Additionally, the proximity of this creek to human settlement increases the prevalence of anthropogenic impacts on seagrass through boat dredging and seine net fishing, hence affecting seagrass survival and the expansion of meadows. This was observed during the study and may explain the lower biomass values recorded in the Western creek. The Eastern creek had species with numerous stems that have the ability to withstand wave action when compared to the other species. Enhalus, which also dominated this area, has big rhizomes, extensive root system, and large fronds, accumulating greater biomass in the belowground components during growth and development (Waycott et al., 2004). In general, the variation is attributed to differences in the biophysical parameters leading to different physiological, phenotypic, and morphological growth patterns of seagrass.

Total Carbon Stocks
The current study focused on seagrass carbon stocks within the mangrove-fringed creeks of Gazi Bay. It compared the carbon stocks between the Eastern and Western creeks that have varying biophysical features. In the Eastern creek (50 ha), the total carbon stocks in the three carbon pools is 13,419.5 ± 4,547 Mg. On the other hand, the Western creek (70 ha) has a total stock of 7,769.3 ± 1,518 Mg of carbon. Although the sediment carbon values were extrapolated to 1-m depth following IPCC (2014) and Githaiga et al. (2017), this may be a limitation since sediment C org decreases with depth. Additionally, future studies may focus on carbon sources and the sediment accumulation rates. Based on the IPCC's tier 1 emission factor of 7.9 tons of C per hectare, protecting seagrass in these two creeks will prevent emissions of 948 Mg CO 2 eq. year-1 (IPCC, 2014). (Tier 1 emission factor uses a gain-loss concept to estimate carbon dioxide emissions, based on general conversion factors, obtained from globally combined databases; Howard et al., 2014). The sediment C org data contributes information on the seagrass carbon budget of the greater Gazi Bay region and adds to our understanding of the total C org of seagrass along the Kenyan coast and the African coastline.
This study shows that the community attributes (species) and carbon stocks (biomass and soil) of seagrass meadows can differ dramatically over small spatial scales and complements prior studies. Intertidal areas that are under the influence of riverine inputs (and greater human disturbance) can have lower species diversity and reduced carbon stocks when compared to those areas that do not experience these impacts. This type of spatial variability should be considered when designing blue carbon projects.

DATA AVAILABILITY STATEMENT
All datasets generated for this study are included in the article/Supplementary Material.

AUTHOR CONTRIBUTIONS
GJ performed conceptualization of the draft, data curation, worked on the methodology, research administration, research investigation, the formal data analysis, writing the original draft, and performed editing and review. AM carried out conceptualization of the research, methodology development, supervision, validation of methods, visualization, writing, review and final editing. GM performed conceptualization, helped in developing methodology, supervision, project administration, validation of methods, and writing of the manuscript (review and editing). JK is a Pew Fellow (2019). He carried out research conceptualization, aligning methodology, field supervision, visualization, and manuscript writing (review and editing).
All authors contributed to the article and approved the submitted version.

FUNDING
The research was funded by the Newton Fund from the British Council (Grant No. 275670159).

ACKNOWLEDGMENTS
The study was done under the scientific framework of Kenya Marine and Fisheries Research Institute under "Punguza Hewa Kaa Project, " with funding from the UK's Newton Institutional Link. We thank Chuka University for technical support. We acknowledge the contribution of the KMFRI's team, Derrick Omollo, Charles Muthama, Peter Musembi, Tom Kisiengo, George Onduso, and Laitani Suleiman, for their field and laboratory assistance.

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars. 2020.00696/full#supplementary-material FIGURE S1 | A section of seagrass meadows in Mangrove fringed Eastern creek.