Geomorphic and Climatic Drivers Are Key Determinants of Structural Variability of Mangrove Forests along the Kenyan Coast

Mangrove forests occur across a diversity of coastal landforms that influence their structural development and productivity. Preliminary studies in Kenya indicate that mangroves growing in the region north and south of Tana River delta have different structural attributes. We hypothesise a close relationship between mangrove distribution, climate and landform types. Floristic composition of mangroves along the coast of Kenya was characterised and differences illustrated using non-metric multidimensional scaling (nMDS). Other structural properties of mangroves such as tree height, basal area, stand density and standing biomass were also assessed and their differences tested using analysis of variance (ANOVA). A hierarchical cluster analysis was used to compare mangrove species based on structural properties. Additionally, a regression fit model was used to investigate the relationship between mangrove standing biomass and possible drivers of variability. The study revealed significant differences in mangrove tree diameter, tree height, basal area, stand density and standing biomass across the sampled sites. High values of structural complexity were observed in estuarine and deltaic settings with high influence of freshwater input whereas low levels of structural complexity were observed for peri-urban with direct human influence. Our findings suggest that structural variability of mangroves in Kenya is highly influenced by geomorphological and climatic variability along the coast as well as the past and present management regimes of the forest.


Introduction
Mangrove forests grow in the intertidal areas of tropical and sub-tropical coasts [1]; occurring across a diversity of landforms such as protected bays, lagoons, estuaries, and small islands. These geomorphic features coupled with climatic and human drivers influence the structure and productivity of mangrove forests. As a result, the structural development of mangroves varies with the diversity of coastal geomorphology, climate and management regimes [1][2][3][4]. On a global scale, solar radiation controls the potential maximum development of a mangrove forest [2] while their distribution is limited by sea surface temperature, atmospheric temperature and precipitation [5]. On a regional scale, mangrove development and occurrence are controlled by geomorphological settings such as topographic expression and composition of landforms; climatic conditions such as temperature and precipitation; and oceanographic factors such as wave, tidal and river influences [6][7][8]. At the local scale, forest structure is a function of the tidal regime which is controlled by the topography of the intertidal area. At this scale, mangroves are classified into different forest types-fringing, riverine, over-wash, basin and dwarf forests [9,10]. There may exist all six types of mangrove forests within the regional boundaries of an environmental setting depending on the local effects of waves, tides and river flow [2].
The environmental factors that support and shape the ecological processes occurring within mangrove ecosystems are the forcing functions of that ecosystem [2]. These factors represent different types of energies such as solar radiation, river flows, tides, land formation and precipitation. The summation of energy sources to a mangrove system such as solar radiation or nutrient loading resulting from river flows minus the energy utilised for maintenance and overcoming stress such as high salinity or drought is the energy signature of that system [2,11]. An increase in the energy signature results in higher rates of energy flow that are evidenced by greater biomass productivity within the mangrove system [2]. This concept of energy signatures to describe mangrove ecosystem structure and function is similar to the concept by Thom [6,12] of using environmental settings to explain mangrove processes from a geomorphological perspective [2]. Thom [6,12] describes five environmental settings where terrigenous sediment inputs are dominant and three settings dominated by the accumulation of carbonate. This geomorphic approach of associating plants with diverse landforms and substrate conditions has been used in mangrove ecology studies across the world to classify different mangrove environmental settings [6,8,[12][13][14][15][16]. Recent work by Worthington et al. [16] classifies the settings into deltaic, estuarine, lagoonal, open coast and carbonate settings.
In Kenya, mangroves are a common feature in deltas, protected bays, creeks, river estuaries and lagoons distributed all along the 600 km coastline [17]. Earlier studies on mangrove forests in Kenya have observed that tree height, basal areas and biomass values vary between areas north and south of River Tana delta [18][19][20]. Classifications of environmental settings based on geomorphology and other forcing functions have been used extensively to describe variability in mangrove forest structure across the world [3,10,[14][15][16][21][22][23]. However, to our knowledge, only one study in Kenya has attempted its application [18]. Understanding the relationships between environmental settings and ecological functions is important for mangrove management because this ecosystem is often cited as being controlled by factors emanating from the land, sea and air [10,13,15,24].
The current study had two main objectives, (1) to characterise the structure and floristic composition of mangrove forests along the Kenyan coast, and (2) compare the structural variability of mangrove formations along the coast. The results of the study enable comparisons of ecological processes across different mangrove areas in Kenya. They also enhance our understanding of major environmental forcing functions influencing the growth and development of mangrove forests in Kenya. A major limitation of this study was limited geomorphic and oceanographic data for the Kenyan coast.

Description of the Study Area
The Kenyan coastline stretches over 600 km from Vanga at the Kenya-Tanzania border in the south to Ishakani at the Kenya-Somalia border in the north [25] (Figure 1). The coastline straddles five counties-Kwale, Mombasa, Kilifi, Tana River County and Lamu. Distinctive features along the coastline include a semi-continuous fringing reef system, sandy beaches, protected bays, estuaries and tidal creeks that support the natural growth of mangrove forests [17,[25][26][27].
All the nine mangrove species described in the Western Indian Ocean (WIO) region [28] occur in Kenya-Rhizophora mucronata, Bruguiera gymnorhiza, Ceriops tagal, Sonneratia alba, Avicennia marina, Lumnitzera racemosa, Xylocarpus granatum, Xylocarpus moluccensis and Heritiera littoralis [17]. These forests occur on reef platforms behind protective outcrops of coral limestone in Lamu; in the estuaries of Tana and Sabaki rivers; behind marine influenced barrier dunes in Ngomeni; in the drowned river valleys of Mombasa, Mtwapa, Kilifi and Mongoni-Dodori; and within the sheltered bays of Vanga and Gazi [17]. The total area of mangroves in Kenya is estimated at 61,271 ha; more than 60% of which occur in Lamu County [17]. Kenya's National Mangrove Ecosystem Management Plan (2017-2027) identifies 14 mangrove management units in Kenya-hereafter referred to as 'sites' (Figure 1). Five of these sites are in Lamu County: Northern Swamps, Northern Central Swamps, Mongoni and Dodori Creek Swamps, Pate Island Swamps, and Southern Swamps. Three of the sites are in Tana River County: Kipini, Mto Tana, and Ngomeni, while three sites are in Kilifi County: Mida, Kilifi, and Mtwapa. Mombasa County is classified as one unit, while Kwale County has two sites: Gazi and the Vanga-Funzi system [17]. This study has maintained this classification scheme for ease of reference and comparisons ( Figure 1).
The Vanga and Gazi sites are sheltered bays separated by the Shimoni Peninsula with the Vanga site located at the mouth of River Umba. The Mombasa, Mtwapa and Kilifi sites are peri-urban systems characterised by narrow creeks whose origins are drowned river valleys [29]. The Mida site has a similar morphology to Mombasa, Mtwapa and Kilifi but lacks a discharging river/stream [29]. The Ngomeni, Mto Tana and Kipini sites are located in Ungwana Bay, a wide bay in front of the River Tana Delta that is characterised by fringing dune complexes [29]. The Southern Swamps, Pate Island Swamp, Mongoni-Dodori Creeks, Northern Central Swamps and the Northern Swamps sites are in the coastal lagoons and multiple small islands that define the Lamu Archipelago.

Climate
The coast of Kenya is characterised by a hot and humid tropical climate and experiences a bimodal rainfall pattern influenced by monsoon winds. During March to May, the weather is dominated by the South East Monsoon (SEM) winds leading to heavy rainfall while the North East Monsoon (NEM) winds occur between October and December and are responsible for short rains [25]. The amount of rainfall varies along the coastline due to the inclination of the coast with the northern parts receiving between 500 and 900 mm yr −1 whereas areas south of Malindi receive rainfall ranging between 1000 and 1600 mm yr −1 [26]. Mean daily temperature ranges between 24 and 30 degrees Celsius (ºC) while humidity averages about 80% throughout the year [30].

Oceanography
Kenya's coast experience semi-diurnal tidal regimes with a maximum tide of 4 m [27]. Major currents influencing the coastal system in Kenya are the East Africa Coastal Current (EACC), the Somali Current (SC) and the Equatorial Counter Current (ECC) [31]. The EACC flows northward throughout the year while the SC is a seasonally reversing current. During the SEM, the SC flows northward as well but reverses during the NEM (Figure 2. The reversed SC meets the EACC somewhere between Malindi and Lamu depending on the strength of the monsoon winds and forms the ECC that flows eastward as an undercurrent [25,31].

Hydrology
Only two perennial rivers drain into the Indian Ocean along the coast of Kenya-Tana and Sabaki rivers. River Tana with a mean annual discharge of about 4000 million m 3 branches into a complex deltaic system that opens to the Indian Ocean at Kipini and Mto Tana sites [32]. South of the Tana delta, the Sabaki river drains into the Indian near Malindi town with an annual discharge of about 2000 million m 3 . The discharge of these two rivers is seasonal and characterised by a bimodal hydrological cycle which corresponds to rainfall patterns [25]. Seasonal rivers draining into the Indian Ocean are Ramisi that drains into the Indian Ocean between the Vanga and Gazi sites, Umba and Mwena that drain into the Vanga site, Mwache that drains in the Mombasa site, and Ndzovuni that drains in the Kilifi site [25,29].
Several freshwater aquifers are distributed along the coastal area of Kenya mainly within the sedimentary terrains. The largest aquifer system stretches from the north eastern parts of the country in Marsabit and terminates in Lamu, spanning about 250 km, part of which is the Merti aquifer [33]. Other aquifers are the Tiwi and Msambweni aquifers in the southern parts of the coast that stretch over 80 km along the coastline [33][34][35][36]. These aquifers are considered important in maintaining groundwater seepage into mangrove forests [35] though information on this remains scanty.

Geomorphology
The morphological structure of the coast of Kenya can be divided into three distinct sections. The southern section running from the Kenya-Tanzania border to Gazi is characterised by wide sheltered bays behind a broken chain of coral reef patches. The central section running from north of Gazi to the mouth of the Sabaki river is characterised by a straight fringing coral reef divided into several segments by narrow tidal outlets of branching bays and estuaries and bounded by steep cliffs on the landward side. The northern section running from the mouth of the Sabaki river to the Kenya-Somalia border is characterised by wide open bays in front of the Tana delta and near the mouth of the Sabaki river that are bordered by long beaches and high dune complexes. Reef patches and sheltered lagoons occur between these open bays [29].

Demography
The coastal counties of Kenya have about 3.9 million inhabitants which is about 8.4% of Kenya's total population [37]. More than 60% of this population resides in Mombasa and Kilifi counties. Kilifi County has the highest population along the coast with about 1.4 million inhabitants. Mombasa county has the highest population density with 5495 persons per km 2 . It is the largest urban centre in the coast with 97 % of the population dwelling within the city. Tana River County has the lowest population density with 8 persons per km 2 [38]. It is estimated that about 70% of communities living adjacent to mangroves derive their wood requirements from the forest [17].

Assessment of Vegetation Structure
The structural characteristics of mangroves in the 14 sites were assessed using a systematic random sampling design. Belt transects running perpendicular to the shoreline were randomly established. Plots measuring 20 m by 20 m were systematically established along these transects to capture the variability resulting from the zonation of mangroves. Across all the sites, a total of 372 plots were established representing an overall sampling intensity of 0.3%. Within each plot, all individual mangrove trees with a diameter at breast height (DBH) ≥ 2.5 cm were identified and counted. Data on tree height (m) and stem diameter (cm) were collected following the procedures outlined by Kauffman & Donato [39]. The ecological importance value (IV) of each species was calculated by summing its relative density, relative frequency, and relative dominance. The basal area (m 2 ha −1 ) and stand density (trees ha −1 ) were obtained following formulae outlined by Cintron & Schaeffer-Novelli [40]. The complexity indices (CI) of the study sites were calculated using Equation (1) below: Considering no robust biomass allometric equations have been developed for mangroves in Kenya, we used the generalised equation for estimating mangrove aboveground biomass (AGB) (Mg ha −1 ) developed by Komiyama et al. [41] (Equation (2)) but used localised species-specific wood density for mangroves in Kenya developed by Gillerot et al. [42].

Factors Affecting Mangrove Structural Variability
To compare the drivers of mangrove variability, mangrove sites were described based on environmental settings, average annual precipitation, population density and river influence. The environmental settings described in this study were modified to suit the various landforms of the Kenyan coast. The following categories were adopted: estuarine coasts are sheltered coasts with one or more rivers flowing into them and with a free connection to the open sea and include the Vanga and Mombasa sites; lagoonal coasts are shallow inland water bodies that are separated from the ocean by a barrier and include the Gazi, Southern Swamps, Pate Island Swamps and Northern Central sites; tidal creek coasts feature a narrow inlet or estuary that is affected by the flow and ebb of ocean tides and include Mtwapa, Kilifi, Mida, Mongoni-Dodori Creek and Mto Tana sites; deltaic coasts are areas of high river influence with sediment accumulating at the mouth of the river and include Ngomeni and Kipini sites; while open coasts are relatively exposed coasts that are only sheltered from the sea by minor reef segments and reef patches such as the Northern Swamps site. Table 1 below summarises key characteristics of the sampled sites. The average annual rainfall data  were accessed from the Kenya Meteorological Department website: http://kmddl.meteo.go.ke:8081/maproom/Climatology/ (accessed on 15 December 2021). The population density was the average population density per square kilometre in the administrative units where the mangrove sites occur sourced from the Kenya National Bureau of Statistics [38]. The influence of the river on the sites was described in an ordinal scale based on the freshwater discharge levels of the rivers draining into the site [27]. * River influence is based on the freshwater discharge levels of the rivers present in a site. High represents a discharge greater than 1500 million m 3 ; moderate (100-215 million m 3 ); low (<100 million m 3 ); no influence represents lack of a discharging river.

Data Analysis
Graphical presentation of data was used to describe the structure of mangroves in Kenya. Multivariate analyses were performed to examine the differences in species composition across the sites based on species abundance. An analysis of similarities (ANOSIM) using the Bray-Curtis similarity index was performed followed by a pairwise post-hoc test with a Bonferroni correction to examine variability in species composition across the sites. The data matrix consisted of the standardised abundance of each tree species at different size-class categories in each site. A percent similarities analysis (SIMPER) was used to determine which tree species contributed the most to the differences found.
These differences were then illustrated using a non-metric multidimensional scaling (nMDS) ordination plot. The nMDS was plotted in two dimensions (2D) using the Bray-Curtis similarity index. We subjected the structural data to normality tests before performing a Box-Cox transformation. A one-way analysis of variance (ANOVA) followed by a post-hoc Tukey pairwise test was performed (p < 0.05) to individually compare DBH, height, basal area and stand density and above-ground biomass across the different sites. A hierarchical cluster analysis (unweighted paired group mean average and squared Euclidian distances) was then performed to determine the degree of similarity of species from the 14 sites based on complexity index, biomass, DBH, and height. A multiple linear regression with a stepwise selection method was used to fit a model describing associations between standing biomass and the possible drivers of variability. The software, OriginPro 9.0, was used to develop the graphical presentations of the structural data while the multivariate analyses were performed using Paleontological Statistics software (PAST 4) and Minitab 18.

Species Composition
Across all the sites, a total of 34,050 individual mangrove trees were sampled. Based on the importance value (IV) index, Rhizophora mucronata was the most important species in 11 out of the 14 sites sampled-Vanga, Gazi, Mombasa, Mtwapa, Kilifi, Ngomeni, Southern Swamps, Pate Island Swamps, Mongoni-Dodori Creeks, Northern Central Swamps and Northern Swamps. Avicennia marina was the most important species in the estuarine area of Ungwana Bay, covering the Mto Tana and Kipini sites (Table A1). Mangrove forest species composition differed significantly across the sampled sites (ANOSIM R: 0.24, p = 0.001). The pairwise post-hoc test with a Bonferroni correction shows Kipini as being significantly different from all the other sampled sites (Table A2). The ordination plot from the nMDS shows little grouping of sites except Kipini that appears to form a cluster ( Figure 3). Based on the SIMPER analysis, Sonneratia alba (20.1% contribution), Ceriops tagal (19.64% contribution) and Avicennia marina (14.18% contribution) contributed more than half (53.87%) of the differences observed across the sites ( Table 2).   Table 3). A similar trend was observed with tree height. Northern Central Swamps in Lamu had the tallest trees overall with an average tree height of 11.54 ± 0.06 m while the lowest tree heights were recorded in Mtwapa (3.06 ± 0.03 m). There were significant differences in tree height across the sites [F (13, 34,050) =1827.28, p = 0.000]. The scatterplots of DBH against mangrove tree height suggest a positive, linear relationship between DBH and height across the sites. However, the regression analysis reveals this relationship is relatively weak, particularly in Mombasa and Mtwapa. (Figure 4).  There were significant differences in stand density across the sites [F (13, 358) = 8.68, p = 0.000]. The lowest stand density was recorded for the mangroves of Northern Swamps (1607 ± 129 trees ha −1 ) while the highest stand density was for the mangroves of Mtwapa (7856 ± 2094 trees ha −1 ). The pattern of stand density across the sampled sites showed an inverse relationship with mean DBH with the results of the Pearson correlation analysis indicating a negative correlation between stand density and mean DBH (r = −0.688, p = 0.0006, Table A3). The mangrove forest of Kipini had the highest mean basal area (29.78 ± 6.46 m 2 ha −1 ) while Mombasa had the lowest basal area (8.70 ± 0.831 m 2 ha −1 ). There were significant differences in mean basal area across the sites [F (13, 358) =5.45, p = 0.000].
The mangrove forest of Pate Island Swamps was the most structurally complex (CI = 35.79) while the mangrove forest of Mto Tana had the lowest complexity index (CI = 5.15). Patterns in complexity index along the sites were similar to the pattern observed for mangrove tree height (m) and basal area (m 2 ha −1 ). In forest analysis, complexity indices can be used to define if a forest is under stress [43]. The low values of structural complexity, tree DBH (cm), and tree height (m) observed for the mangroves of Mto Tana are a strong indication of stress within the forest. The mangroves of Mombasa had the lowest mean above-ground biomass (AGB) at 76.77 ± 9.61 Mg ha −1 while those in Kipini had the highest standing biomass at 358 ± 106 Mg ha −1 ( Figure 5). There were significant differences in AGB across the sites [F (13, 358) =15.36, p = 0.000].

Multivariate Analysis
Hierarchical cluster analysis using unweighted paired group mean average and squared Euclidian distances on mean AGB (Mg/ha) revealed distinct groupings across the mangrove sites ( Figure 6). Kipini has a unique mangrove assemblage with the highest standing biomass observed in the study. Mombasa, Mto Tana and Mtwapa, are also shown to be significantly different from the rest of the sites. These sites have the lowest biomass and complexity index values recorded in the study. Kilifi and Mida sites, both located in Kilifi County in the central region of the coastline are shown to be similar. These sites share some similarities with Northern Central swamps, Vanga, Gazi and Mongoni-Dodori creeks, located in the northern and southern regions of the coastline. Ngomeni, Pate Island Swamps, Southern Swamps and Northern Swamps are also shown to be similar. The multiple linear regression analysis indicated environmental settings and population density best explained the variability in standing biomass. The fit regression model explained 86.74% of the variance in standing biomass.

Discussion
This study revealed significant differences in mangrove forest structure along the coast of Kenya. High values of standing biomass and structural complexity were observed in the riverine mangroves of the Tana Delta as well as mangroves found in the protected islands of Lamu Archipelago in the northern parts of the coastline. High values of standing biomass were also observed within the sheltered bays of the southern parts of the coastline. In the central region characterised by drowned river valleys, a fringing coral reef and peri-urban settings, relatively lower levels of standing biomass and structural complexity were observed. These patterns of mangrove structural variability closely follow the patterns of geomorphological variability along the coastline [29].
Mangroves in the northern region of Kenya are influenced by the interplay between run-off from River Tana, hydrodynamics and air-sea interactions [44]. During the south east monsoon (SEM), high levels of sediment and nutrients are deposited from the hinterland into the ocean by River Tana and River [32,45]. Nutrient availability is key in enhancing mangrove growth and productivity [4]. However, it is probably the maintenance of balance between mineral nutrients and substrate salinity that is relevant, rather than the absolute nutrient levels [4]. For instance, the Heriteria littoralis dominated forest of the Kipini site at the mouth of River Tana receives copious amounts of freshwater, sediment and nutrients from the hinterland leading to higher productivity and overall biomass [32,45]. Conversely, the scrub mangroves growing in the Mto Tana site are majorly composed of dwarf trees due to low freshwater and sediment supply. Currently, the River Tana flows directly into an estuary at Kipini rather than in the complex system of channels and distributaries leading to its old mouth at Mto Tana. The little freshwater that still flowed into the old delta through the Kalota brook was blocked through the construction of a multi-purpose community dam in 1988 [29]. This, together with the presence of solar salt works around the area has ultimately limited mangrove growth in the site.
In the sites north of the Tana river (Southern swamps, Pate Island Swamps, Mongoni-Dodori Creek, Northern Central Swamps and Northern Swamps), the northern flowing EACC moves nutrient-rich sediment from the estuaries of the Tana and Sabaki rivers though long-shore transport and is responsible for productivity in the area [44]. These nutrients are later resuspended during the northeast monsoon (NEM) when the north flowing EACC meets the south flowing Somali Current (SC) causing upwelling and nutrient enrichment [44]. The presence of groundwater seepage in the area is another factor that could be driving the productivity of mangroves in the area by providing freshwater into the system [33,34,36]. The management system could also be responsible for the luxuriant growth of mangroves in some parts of Lamu. The Northern Swamps mangroves and some parts of Northern Central Swamps fall under the Kiunga Marine National Reserve (KMNR), where commercial exploitation of mangroves is prohibited [46]. This could explain the high DBH values observed in KMNR.
Compared to the northern parts of the coastline, nutrient levels in Ramisi and Umba river systems in the southern region are relatively low [29]. The differences in the sediment loading, nutrient levels and freshwater input between the river systems of the northern and southern regions explains the different levels of productivity between the northern and southern parts of the coastline of Kenya [29].
In the Vanga area, small, persistent, localised upwelling events occur around the narrow zonal strip that extends between Northern Pemba Island and the mainland, right at the border between Kenya and Tanzania during the NEM. These upwelling events are as a result of the instabilities of the EACC around the chain of islands (Mafia, Unguja and Pemba) and along the continent's lateral boundary [47]. Halo et al. [47] suggest that these upwelling events occur on-shore throughout the annual cycle evident by the development of near shore negative wind stress curls and consequent positive Ekman vertical velocities during the NEM. The presence of an intermittent stream (Mkurumudji river), ground water aquifers (the Tiwi and Msambweni aquifers) [35], as well as intensive communitybased mangrove conservation activities could explain the high structural complexity of mangroves in the Gazi site. The carbon offset project, Mikoko Pamoja in Gazi has played a role in the restoration and protection of mangroves in the area [48].
On the other hand, the mangroves in the central region are influenced by the geomorphology of the area as well as human influences. In this region, the transition of the littoral to the paralic zone is marked by a broken chain of prominent hills of the coastal range. The elevated areas of the coastal range and the Giriama hills lands have protected the littoral zone against strong fluvial erosion and sedimentation. This is evident by the development of a continuous reef complex in the littoral zone of the central region [29]. This coupled with the absence of perennial rivers partly explains the low biomass observed in the mangroves of this region comprising of Mtwapa, Kilifi and Mida. The availability of terrigenous sediments plays a crucial role not only in providing nutrients necessary for mangrove growth but also in providing the necessary accommodation space for mangrove colonisation. For instance, the mangroves of Mtwapa and Kilifi occur in an area characterised by steep cliffs on the margin of the land and the ocean with a narrow tidal zone. Hence, there is limited accommodation space available for mangrove colonisation. The dominance of Ceriops tagal in Mida resulting from past selective logging of desirable Rhizophora mucronata trees [18] could possibly explain the lower values of standing biomass. Ceriops tagal has a relatively lower DBH and mean height compared to the Rhizophora mucronata and Avicennia marina that dominate the rest of the sites (Table A1). The peri-urban mangroves in Mombasa are mainly threatened by human stressors such as overharvesting, habitat conversion, pollution and sedimentation [49].
High values of mangrove standing biomass were observed in estuarine and deltaic sites as well as in sites located in Lamu County (Table 4). Mangroves are facultative halophytes and are more productive in riverine systems with high inputs of freshwater and nutrients [4]. There is evidence of longshore transport of nutrients to the mangroves of Lamu area deposited from the Tana and Sabaki rivers [44]. This demonstrates linkages between mangrove structural properties and geomorphologically defined habitats as described by the concept of environmental settings [8,12,21]. However, relatively lower values of mangrove standing biomass were also recorded for some sites within similar estuarine settings such as Mombasa and Mto Tana. While these areas should have supported high levels of mangrove productivity, the lower biomass values observed could be explained by the presence of human-induced stressors that act to limit mangrove growth and productivity [29,49]. This shows that in addition to the interactions between ecological and geomorphological processes, human influence is also a key factor in driving mangrove forest structure in Kenya. A fitting explanation of mangrove structural variability would have to include the feedback processes resulting from the interactions between the forest and the communities living around them. This includes influences that improve the conservation of mangroves (positive feedback) as well as influence the degradation of the mangrove ecosystem (negative feedback). Table 4. Summary attribute table describing the environmental settings, forcing functions and structural properties of mangroves along the Kenyan coast. (1) denotes the number of species encountered in each site, (2) is basal area (BA) (m 2 ha −1 ), (3) is canopy height (m), (4) is stand density (SD) (trees ha −1 ), (5) is standing/aboveground biomass (AGB) (Mg ha −1 ), (6) is complexity index (CI). Basal area, canopy height, stand density and standing biomass are reported as mean ± standard error. Southern is Southern swamps, MD Creek is Mongoni Dodori Creek, NC is Northern Central swamps, Northern is Northern swamps.

Site Environmental Setting
Mangrove Twilley [2,11] explains that the amount of structure that develops in a mangrove system will be determined by the net energy available to the system-the difference between the energy available to the system-solar, chemical (organic matter and nutrient input) and mechanical energy (wind, tides and waves) and the energy required for maintenance or for overcoming stress-hyper salinity, drought, and actual biomass removal through harvesting. All external factors of coastal systems could be evaluated in the form of energy and measured in units of energy such as kilocalories or joules and this could be useful in explaining and predicting patterns of mangrove structural development in Kenya.

Conclusions
Mangrove forests in Kenya exhibit differences in floristic composition and structural development that closely follows the patterns of geomorphological and climate variability along the coastline. Effective mangrove management requires an adequate understanding of the forcing functions or potential energies that are present in an area. The management of mangrove ecosystems is most efficient when it prioritises the protection of energy flows and the processes that drive the system. This is because the system retains most of its natural capacity and requires minimal resources and effort to manage [10]. Our findings will be useful to mangrove managers and policy makers to tailor their strategies to ensure mangrove management works with the ecosystem processes in order to achieve more successful outcomes. Further investigation of the physical and chemical parameters within mangrove habitats is required to perform statistical tests and make meaningful comparisons of mangrove ecosystem properties and processes.

Acknowledgments:
The authors recognise the technical support received during the entire duration of the study by the team at KMFRI, KFS, WWF and TNC.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.