Aboveground Biomass Distribution in a Multi-Use Savannah Landscape in Southeastern Kenya: Impact of Land Use and Fences

: Savannahs provide valuable ecosystem services and contribute to continental and global carbon budgets. In addition, savannahs exhibit multiple land uses, e.g., wildlife conservation, pastoralism, and crop farming. Despite their importance, the e ﬀ ect of land use on woody aboveground biomass (AGB) in savannahs is understudied. Furthermore, fences used to reduce human–wildlife conﬂicts may a ﬀ ect AGB patterns. We assessed AGB densities and patterns, and the e ﬀ ect of land use and fences on AGB in a multi-use savannah landscape in southeastern Kenya. AGB was assessed with ﬁeld survey and airborne laser scanning (ALS) data, and a land cover map was developed using Sentinel-2 satellite images in Google Earth Engine. The highest woody AGB was found in riverine forest in a conservation area and in bushland outside the conservation area. The highest mean AGB density occurred in the non-conservation area with mixed bushland and cropland (8.9 Mg · ha − 1 ), while the lowest AGB density (2.6 Mg · ha − 1 ) occurred in overgrazed grassland in the conservation area. The largest di ﬀ erences in AGB distributions were observed in the fenced boundaries between the conservation and other land-use types. Our results provide evidence that conservation and fences can create sharp AGB transitions and lead to reduced AGB stocks, which is a vital role of savannahs as part of carbon sequestration.


Introduction
Savannahs are characterized by scattered tree cover and continuous coverage of grass-dominated herbaceous plants [1,2]. On the African continent, savannahs and woodlands play a particularly large role in the carbon cycle, and wildlife and biodiversity conservation, while providing livelihoods for a huge human population [3]. The area covered by savannahs is roughly three times larger than that of forests, corresponding to approximately 50% of the total area of the African continent. Savannahs therefore represent a major carbon stock in Africa despite having a lower carbon density compared to forests [4][5][6]. Another significant feature of the African carbon cycle is that emissions caused by land-use change are greater than fossil fuel emissions [7,8]. A large part of these emissions originates

Field Data
The field data were collected between 15 and 22 August 2018 to estimate the AGB of woody plants (trees and shrubs). The sample plots were selected subjectively to cover variation in land-use and land cover type based on high resolution satellite imagery in Google Earth, and tree cover and tree height based on ALS data ( Figure 1). In total, 49 sample plots were surveyed. The field plots were positioned using a Trimble GeoXH GNSS receiver with differential correction.
The sample plot design consisted of circular plots of different sizes. The main plot was 0.1 ha in size (radius 17.84 m) and was used for inventorying all the trees with a diameter at breast height (DBH, 1.3 m height from the ground) of more than 5 cm. Height (H) for the highest, median, and shortest tree were also measured at each plot using a hypsometer (Suunto). Tree species was identified for all of these trees. Furthermore, four "subplots" of 0.01 ha (radius 5.64 m) located within the main plot were used for inventorying shrubs with DBHs of 1-5 cm (see [56] for subplot locations), and four "micro plots" of 0.001 ha (radius 1.78 m) in the central points of the subplots for measuring shrubs with DBHs < 1 cm. Shrub measurements included count, DBH, basal diameter (BD), crown diameter (CD), and height for a median-sized shrub. The dominant woody species of each plot was also recorded.

Field Data
The field data were collected between 15 and 22 August 2018 to estimate the AGB of woody plants (trees and shrubs). The sample plots were selected subjectively to cover variation in land-use and land cover type based on high resolution satellite imagery in Google Earth, and tree cover and tree height based on ALS data ( Figure 1). In total, 49 sample plots were surveyed. The field plots were positioned using a Trimble GeoXH GNSS receiver with differential correction. The sample plot design consisted of circular plots of different sizes. The main plot was 0.1 ha in size (radius 17.84 m) and was used for inventorying all the trees with a diameter at breast height (DBH, 1.3 m height from the ground) of more than 5 cm. Height (H) for the highest, median, and shortest tree were also measured at each plot using a hypsometer (Suunto). Tree species was Aboveground biomass of trees with DBH > 5 cm (AGB trees ) was computed using pan-tropical biomass model [57] due to the absence of local, species-specific allometric equations. The model (Equation (1)) is based on DBH (cm), H (m), and wood-specific gravity (ρ, g/m 3 ). Wood densities were obtained from a species-specific list in the BIOMASS package [58] in the R software environment [59].
AGB shrubs = exp(−2.281 + 1.525 ln(BD) + 0.831 ln(CD) + 0.523 ln(H)) (2) Finally, we normalized AGB values per hectare and calculated the plot-level AGB as a sum of the tree and shrub AGB. Hereafter, by AGB, we refer to this aboveground biomass of woody plants unless specified otherwise.

Airborne Laser Scanning Data (ALS) and Biomass Mapping
Airborne laser scanning data were used to generate a reference canopy height model and to predict a high-resolution wall-to-wall AGB map for the study area. The scanning was conducted in late March 2014 and covered an area of 433 km 2 . The sensor was a Leica ALS60 and a maximum of four returns per pulse were recorded. The pulse density was 1.04 pulses/m 2 .
The data vendor (Ramani Geosystems, Kenya) pre-processed the ALS data, including filtering of the ground returns using Terrascan software (Terrasolid Oy, Finland). The data were delivered as georeferenced point clouds in the UTM/WGS84 coordinate system with ellipsoidal heights. The ground-classified returns were used for generating digital elevation models (DEM) at a 1-m cell size. The ALS point cloud elevations were normalized to height from the ground levels using DEM. Furthermore, buildings, power lines, and outliers (high points) were filtered using Terrascan, LAStools (Rapidlasso GmbH), and manual editing.
A 3.5-m height threshold provided the best model between ALS metrics and field biomass and was used to separate understory and ground returns from the canopy returns. Height metrics were calculated separately using first and last returns and canopy cover metrics using all returns (single, first, and last) ( Table A1). The variables included all the variables available in the FUSION software [62] and ones used in our earlier study [63]. Square root transformation was applied to AGB, as it was found to improve the linear relationship between AGB and explanatory variables. The "regsubset" function in the "leaps" package [64] was used to fit multiple linear regression models between the ALS metrics calculated from the ALS point density clipped over the field plot and the AGB calculated from that field plot. The leave-one-out cross-validation root mean square error (RMSE) and the coefficient of determination (R 2 ) were used to select the best AGB model. The predictions were back-transformed (squared), and the square of the residual standard error was added to the predicted values to avoid back-transformation bias [45,65]. For AGB prediction at wall-to-wall level, spatial grids of ALS metrics were generated at a spatial resolution of 30 m × 30 m. Mean densities of AGB in each land-use and land cover class was calculated from the AGB map.

Satellite Imagery and Land Cover Mapping
We collected Sentinel-2 images (top of atmosphere reflectance) with cloud cover less than 20% in the images during the dry seasons [short dry season (January 1 to February 28) and long dry season (July 1 to September 30)] in 2017 and 2018, and pre-processed them in the GEE platform. In total, 103 Sentinel-2 images (bands with a resolution of 10 m and 20 m only) were used to calculate the median dry season image. Median dry season images were calculated for all bands in the blue to the shortwave infrared spectral range based on all available cloud-free pixels ( Figure 1B). In addition, a normalized difference vegetation index (NDVI) [66], an enhanced vegetation index (EVI) [67], EVI2 [68], two variants of normalized difference infrared index (NDII-1 and NDII-2) [69], and an optimized soil-adjusted vegetation index (OSAVI) [70] were calculated from the median image.
Additionally, land cover classification was performed in the GEE platform. In addition to median dry season Sentinel-2 composite and vegetation indices, input data included an ALS-based canopy height model (CHM). The land cover in the landscape was classified into four land cover types (cropland, grassland, forest and bushland) according to the Land Degradation Surveillance Framework [71]. Cropland is cultivated land with annual or perennial crops, while grassland contains grasses and other herbs with less than 10% woody cover. Forest in our classification is made up of a continuous stand of trees with partly interlocking crowns, typically along the riverbeds. Bushland is made up of mixed trees and shrubs with a canopy cover of 40% or more, while thickets are closed stands of bushes and climbers usually between 2 m and 7 m tall and shrubland are open or closed stands up to 3 m tall. For this study, thickets and shrubland were incorporated into bushland because we had few field plots for those classes and the classes were similar in reflectance and vegetation characteristics.
In the first step, training data were collected through visual interpretation using ArcGIS 10.3 for the four land cover classes. The pixels for training the classifier were selected based on image interpretation and CHM. Classification and regression trees (CART) [72] were observed to obtain the highest overall accuracy among the classifiers in GEE and was thus selected for the classification. The reference data set for accuracy assessment included the 49 points surveyed in 2018 in the field, which were not used as training points in the classification. Finally, manual editing was performed in ArcGIS to address some of the apparent misclassification in the land cover map.

Wildlife and Livestock Data
Elephant, buffalo, and cattle data were taken from the Tsavo-Mkomazi large mammal census of 2014 to be comparable with the 2014 ALS data used. The wildlife census is conducted by the Kenya Wildlife Service (KWS) every three years to establish the status of key species in the Tsavo ecosystem. The census is carried out from fixed-wing aircrafts and the data collection procedure is described in detail in [73]. The animal spatial distribution and densities were further compared with AGB in the studied landscape ( Figure 3).

Statistical Analyses of AGB Data
The plot-level AGB values were used for computing descriptive statistics (range, mean, median, and standard deviation) for the field data. The Kruskal-Wallis test was conducted to study whether differences in AGB were statistically significant between the land-use regions and land cover classes. Furthermore, median and mean values of the AGB per class were illustrated with a box plot for the different land-use regions and land cover classes. We also estimated the percentage area covered by each land cover in the respective land-use region. Finally, 500-m wide buffers were set in 11 segments of land-use region boundaries to assess local AGB differences ( Figure 1A). The buffers were categorized into fenced and non-fenced segments to determine the effect of the fence on AGB. Pixel values were studied separately for two sides of the boundary by calculating the percentage of zero AGB pixels. Furthermore, medians of the non-zero AGB values were studied using the Wilcoxon test. All analyses were performed in the R statistical environment [74].

Aboveground Biomass Estimates and Map
Woody AGB estimates based on the field plot measurements are summarized in Table 1. The maximum plot-level values are nearly 365 megagrams per hectare (Mg/ha) and were observed in the riverine forest. The plots with the lowest AGB had very little woody biomass and were located in the grassland areas. Table 1. Summary of the aboveground biomass (AGB) values based on the field data according to the diameter at breast height (DBH) class (n = 49). AGB was estimated based on diameter at ground level for shrubs with a DBH < 1 cm. SD = standard deviation, IQR = interquartile range. The final modeling results for mapping AGB using ALS data are shown in Figure 3. The model was based on two variables: CC.all (percentage of all returns above 3.5 m; p < 0.001) and Elev.min.fr (minimum elevation of the first returns above 3.5 m; p < 0.001). The model performed well in terms of model fit (R 2 = 0.88) although RMSE based on leave-one-out cross-validation was relatively large (26 Mg/ha, 75.6% of the mean AGB). However, the model did not show any signs of systematic over-or under-estimation ( Figure 3).

Aboveground Biomass Estimates and Map
Woody AGB estimates based on the field plot measurements are summarized in Table 1. The maximum plot-level values are nearly 365 megagrams per hectare (Mg/ha) and were observed in the riverine forest. The plots with the lowest AGB had very little woody biomass and were located in the grassland areas. The final modeling results for mapping AGB using ALS data are shown in Figure 3. The model was based on two variables: CC.all (percentage of all returns above 3.5 m; p < 0.001) and Elev.min.fr (minimum elevation of the first returns above 3.5 m; p < 0.001). The model performed well in terms of model fit (R 2 = 0.88) although RMSE based on leave-one-out cross-validation was relatively large (26 Mg/ha, 75.6% of the mean AGB). However, the model did not show any signs of systematic overor under-estimation ( Figure 3). The AGB map shows predicted biomass density patterns at 30 m × 30 m resolution ( Figure 4). The mean AGB in the study area was 5.9 Mg/ha. The Riverine forests in the southern and southeastern parts of the landscape within THWS had the largest AGB densities. We also observed relatively large AGB densities outside the protected areas towards the foothills of the Taita Hills, in the northeastern part of the landscape. Aboveground biomass spatial variations were also relatively large in the matrix The AGB map shows predicted biomass density patterns at 30 m × 30 m resolution ( Figure 4). The mean AGB in the study area was 5.9 Mg/ha. The Riverine forests in the southern and southeastern parts of the landscape within THWS had the largest AGB densities. We also observed relatively large AGB densities outside the protected areas towards the foothills of the Taita Hills, in the northeastern part of the landscape. Aboveground biomass spatial variations were also relatively large in the matrix and in LUMO Oza. On the other hand, the lowest AGB values were found in the nearly treeless grassland of THWS, LUMO East Mramba, LUMO West Mramba, and TWNP.
Wildlife (elephant and buffalos) and livestock (cattle) were highly evident in the conservation areas based on the 2014 KWS wildlife census. Elephants were present in LUMO Mramba East and THWS, and were absent in the matrix ( Figure 4, Table 2). Cattle were found in all the land-use regions, except in the small portion of TWNP captured during the ALS campaign ( Figure 4). Their density was highest in LUMO East Mramba (11.43 animals/km 2 ), a portion of the landscape secured for livestock grazing and was second highest in the matrix (4.40 animals/km 2 ), where agriculture is the most common land use. Buffalos were found in the conservation areas, showing the highest number per unit area in THWS (Table 2). We categorized the animals into three herd sizes, in which the number of animals per herd differed per animal species (Figure 4). We saw no elephants or buffaloes in the matrix during the 2014 wildlife census. Furthermore, no animals were visible in LUMO Oza.
Wildlife (elephant and buffalos) and livestock (cattle) were highly evident in the conservation areas based on the 2014 KWS wildlife census. Elephants were present in LUMO Mramba East and THWS, and were absent in the matrix ( Figure 4, Table 2). Cattle were found in all the land-use regions, except in the small portion of TWNP captured during the ALS campaign ( Figure 4). Their density was highest in LUMO East Mramba (11.43 animals/km 2 ), a portion of the landscape secured for livestock grazing and was second highest in the matrix (4.40 animals/km 2 ), where agriculture is the most common land use. Buffalos were found in the conservation areas, showing the highest number per unit area in THWS (Table 2). We categorized the animals into three herd sizes, in which the number of animals per herd differed per animal species (Figure 4). We saw no elephants or buffaloes in the matrix during the 2014 wildlife census. Furthermore, no animals were visible in LUMO Oza.

Land Cover Classification
The overall land cover classification accuracy was 88.78%. The producer's and user's accuracy are shown in Table A2. The land cover map shows the distribution of the land cover classes in the landscape ( Figure 5). Bushland and cropland dominate the matrix in northern and northeastern parts of the landscape, while grassland that is representative of the savannah biome dominates the southern and southeastern parts (THWS, LUMO Mramba, TWNP). LUMO Oza is almost completely bushland as there is less agriculture and livestock management. Forest is the land cover type with the Land 2020, 9, 381 11 of 24 smallest area, located mostly in THWS along the Bura River. THWS also has relatively large patches of bushland in its eastern parts bordering the matrix. Cropland is also present in the eastern part of THWS, while it is not observed in the other conservation areas. The THWS wardens consider it a form of informal encroachment.
The overall land cover classification accuracy was 88.78%. The producer's and user's accuracy are shown in Table A2. The land cover map shows the distribution of the land cover classes in the landscape ( Figure 5). Bushland and cropland dominate the matrix in northern and northeastern parts of the landscape, while grassland that is representative of the savannah biome dominates the southern and southeastern parts (THWS, LUMO Mramba, TWNP). LUMO Oza is almost completely bushland as there is less agriculture and livestock management. Forest is the land cover type with the smallest area, located mostly in THWS along the Bura River. THWS also has relatively large patches of bushland in its eastern parts bordering the matrix. Cropland is also present in the eastern part of THWS, while it is not observed in the other conservation areas. The THWS wardens consider it a form of informal encroachment.

Effect of Land Cover and Land Use on Aboveground Biomass
Aboveground biomass values for the land cover types are shown in Figure 6 and Table 3. The forest had the highest mean AGB (75.5 Mg/ha) followed by bushland (9.0 Mg/ha) and cropland (5.8 Mg/ha). Grassland had clearly the lowest mean AGB, as it is mostly treeless (mean 1.8 Mg/ha, median 0 Mg/ha). However, bushland, cropland, and grassland also had very high AGBs at certain locations (maximum values in Table 3). These areas correspond to "forest-like" bushland with trees and large shrubs. The highest values in the cropland were found in the fallowed fields and patches of bush and in the tree-covered areas next to the fields. In addition, certain farmers practice agroforestry, meaning that they grow trees for fruit and timber production and for providing shade for crops. Furthermore, the grasslands also have scattered large trees, e.g., in Figure 2C. We observed significant AGB

Effect of Land Cover and Land Use on Aboveground Biomass
Aboveground biomass values for the land cover types are shown in Figure 6 and Table 3. The forest had the highest mean AGB (75.5 Mg/ha) followed by bushland (9.0 Mg/ha) and cropland (5.8 Mg/ha). Grassland had clearly the lowest mean AGB, as it is mostly treeless (mean 1.8 Mg/ha, median 0 Mg/ha). However, bushland, cropland, and grassland also had very high AGBs at certain locations (maximum values in Table 3). These areas correspond to "forest-like" bushland with trees and large shrubs. The highest values in the cropland were found in the fallowed fields and patches of bush and in the tree-covered areas next to the fields. In addition, certain farmers practice agroforestry, meaning that they grow trees for fruit and timber production and for providing shade for crops. Furthermore, the grasslands also have scattered large trees, e.g., in Figure 2C. We observed significant AGB differences among the land cover types (p < 0.001) according to the Kruskal-Wallis mean rank test. Furthermore, the Dunn test indicated a significant difference (p < 0.05) between all the land cover types ( Figure 6).
When comparing the land-use regions, the mean AGB values in descending order were 8.9 Mg/ha in the matrix 8 (Table 3). According to the Kruskal-Wallis test, the AGB differences among land-use regions were significant (p < 0.001). These differences are mainly explained by dissimilarities in the land cover class distributions (Figure 7). The matrix has very little grassland with low AGB, but a large fraction of bushland with a relatively high AGB. The area also has some forest and cropland with high maximum values, which increase the mean AGB. LUMO Oza also mainly consists of higher AGB bushland, while lower AGB regions have larger fractions of grassland. This includes both the West Mramba grazing area and various protected areas. We conducted pairwise comparisons between the classes using the Dunn test, which indicates a significant difference (p < 0.05) between all the classes (Figure 7). differences among the land cover types (p < 0.001) according to the Kruskal-Wallis mean rank test. Furthermore, the Dunn test indicated a significant difference (p < 0.05) between all the land cover types ( Figure 6).
When comparing the land-use regions, the mean AGB values in descending order were 8.9 Mg/ha in the matrix 8.8 Mg/ha in LUMO Oza, 4.8 Mg/ha in THWS, 3.8 Mg/ha in TWNP, 2.6 Mg/ha in LUMO West Mramba, and 2.4 Mg/ha in LUMO East Mramba (Table 3). According to the Kruskal-Wallis test, the AGB differences among land-use regions were significant (p < 0.001). These differences are mainly explained by dissimilarities in the land cover class distributions (Figure 7). The matrix has very little grassland with low AGB, but a large fraction of bushland with a relatively high AGB. The area also has some forest and cropland with high maximum values, which increase the mean AGB. LUMO Oza also mainly consists of higher AGB bushland, while lower AGB regions have larger fractions of grassland. This includes both the West Mramba grazing area and various protected areas. We conducted pairwise comparisons between the classes using the Dunn test, which indicates a significant difference (p < 0.05) between all the classes (Figure 7).

Effect of Fences on Aboveground Biomass
Lastly, we compared AGB values in the fenced and non-fenced boundaries of the land-use regions (see Figure 1A for buffer numbers). Table 4 reports the fraction of zero AGB pixels for two sides of the buffer and the Wilcoxon test results for the non-zero AGB values.
The largest differences in the percentage of zero AGB occurred in the fenced boundaries (buffers 3, 5, 7 and 8). Most of the non-fenced boundaries (buffers 1, 2, 4, 6 and 11) showed only small differences. However, a greater difference was observed in the non-fenced buffer 9, which corresponds to the boundary between bushland part of THWS and the cropland-dominated matrix. Buffer 10 showed relatively small difference in the presence of zero AGB although there is a fence. This boundary is between THWS and the matrix in the eastern part of the study area.
The medians of the non-zero AGB values differed most substantially in the fenced buffers 7, 8 and 10 (all differences highly significant according to the Wilcoxon test) ( Table 4). Although percentage zero AGB was substantially higher in LUMO West and East Mramba than in the matrix in the fenced buffers 3 and 5, median AGB did not differ significantly (P > 0.05). However, smaller but highly significant differences were also observed in the non-fenced buffers 1, 6 and 9. Buffer 1 is located in the non-fenced boundary between TWNP and LUMO West Mramba, where bushland in the northern part of TWNP has a relatively high AGB compared to grassland-dominated West Mramba. Buffer 6 corresponds to the boundary between two conservation areas, LUMO East Mramba and THWS. Table 4. Percentage of zero woody aboveground biomass (AGB) and median AGB for non-zero AGB pixels. P value refers to the Wilcoxon test results made for the non-zero AGB values. Numbers in the end of land-use region names refer to the numbers of buffers in Figure 1A.

Remote Sensing-Based Biomass and Land Cover Maps
We used field data and ALS metrics to create a wall-to-wall high-resolution AGB map. The model fit and accuracy were similar [75,76] or compared favorably with previous studies in temperate and tropical forests [65,[77][78][79]. Our model was based on two predictors: minimum elevation of the first returns above 3.5 m and percentage of all returns above 3.5 m. These variables characterize canopy height and cover, both of which are related to AGB. Similar combinations of height and cover variables have also been used in previous studies in sub-Saharan Africa [78,80,81]. Field-measured AGB included both shrubs and trees. According to the field data, shrubs (DBH 1-5 cm) can make an important contribution to woody AGB. However, as a height threshold of 3.5 m was used to separate canopy and ground returns, woody vegetation less than 3.5 m in height does not affect the ALS variables. Therefore, areas where shrubs are less than 3.5 m in height appear as zero AGB in the map. We selected the height threshold from the tested values, as it provided the most accurate predictions. Further research should be conducted to map AGB variations in the smallest shrubs and grasses.
We used Sentinel-2 satellite images and the CART algorithm in the GEE platform for creating the LULC map. We achieved a good overall accuracy of 88.78% when using dry season composites. Previous studies have shown that the dry season is best suited for separating variations in woody AGB [82,83]. One topic for further research would include classifying various grassland types within the study area.
Spatially explicit AGB and LULC maps offer additional knowledge of AGB variations across the savannah landscapes compared to spatially limited field inventories. In this study, maps demonstrated the link between LULC and AGB, and sharp AGB gradients in certain boundaries of the land-use regions. Furthermore, maps enable geospatial analyses of the AGB patterns, e.g., together with wildlife and livestock inventories, and can inform land management interventions [45]. In our study, maps showed that grassland concentrated in the wildlife conservation areas, where AGB was reduced due to the browsing effect on trees [83]. As there are fewer large mammals outside the conservation areas, their negative impact on woody vegetation is less in these areas. Therefore, wildlife and livestock frequency in the multi-use landscape contributes to the low biomass densities in the region.

Effect of Land Cover and Land Use on Aboveground Biomass
Our results reveal a significant difference (P < 0.05) in woody AGB among the land cover and land-use classes in the studied landscape. In general, AGB is concentrated in areas with larger tree densities. According to the field data, shrubs and smaller trees can also have a considerable effect on woody AGB density. We observed the highest AGB densities in the forest along the Bura River Valley and towards the foothills of the Taita Hills, while grassland had the lowest AGB densities. As the forest class only occupies a small area, other land cover classes contributed more to the total AGB stock at the landscape level. This emphasizes that a greater amount of AGB is stored in open savannah and bushland than in the forest. Bushland occupies more than half of the total area, and therefore contributes the most to the total AGB stock. The contribution of cropland to the total AGB in the landscape is due to agroforestry (i.e., trees growing on cropland). The mean AGB densities in the landscape were low compared to montane forest, exotic plantation, and woodland in the higher altitudes of the Taita Hills [65,84]. Furthermore, biomass in the bushland was comparable to the leaf biomass of sisal (Agave sisalana) in a commercially owned plantation established in the savannah landscape in Taita Taveta [85]. Low precipitation [39,86], small-scale farming by resource-poor farmers [87], low CO 2 concentrations in arid and semiarid regions [88], and disturbance from fire and herbivores [89,90] are among factors responsible for the generally low AGB in the savannah landscape.
We categorized the multi-use savannah landscape into conservation (TWNP, THWS, and LUMO) and non-conservation areas (matrix) based on land use. Furthermore, the conservation types were categorized based on ownership and management. The TWNP, LUMO, and THWS are government, community, and privately owned and managed, respectively. The AGB differences between land-use regions are driven by the land cover differences. We observed the highest woody AGB densities in non-conservation areas (matrix), which are mainly bushland and cropland, while LUMO West Mramba and LUMO East Mramba, community owned and managed wildlife sanctuaries that are mainly grassland, showed the lowest mean AGB densities. THWS and TWNP had similar mean AGB densities, while LUMO Oza had a much higher AGB compared to grassland-dominated regions because of its larger fraction of bushland.
Our results support the hypothesis that there is a link between the land use (conservation and non-conservation) and dominant land cover type, which affect the observed AGB patterns. Presence of wildlife is important for grassland to remain sparsely wooded, and hence, wildlife conservation contributes to open grassland with relatively low woody AGB. Furthermore, ranches for livestock contribute to the low AGB. THWS and LUMO West Mramba serve as a migratory corridor for elephants moving between Tsavo West and Tsavo East NPs in search of food and water [51]. Contrary to the February 2011 elephant census [73], the elephant density in the region increased from less than 0.5 elephants/km 2 in 2011 to > 2 elephants/km 2 in February 2014. This considerable increase in elephant population contributes to low AGB densities in the region. Habitat improvements through water supplementation in the protected areas also attract wildlife and further create pressure on the vegetation. Waterholes attract large congregations of herbivores particularly during the dry season [73]. Williams et al. [91] have also suggested the presence of surface water acts as a determinant of the distribution of water-dependent wildlife species. The wildlife and livestock census data also showed that private (THWS) and government (TWNP) owned conservation areas had more wildlife (elephants and buffalos) while the community owned conservation areas attract more livestock. This could be Land 2020, 9,381 16 of 24 associated with the management strategies employed by the respective agencies. Therefore, policies and management strategies geared towards woody vegetation protection should be introduced into wildlife conservation management plan in order to reduce AGB decline in conservation areas.
Recent studies in the same region show that conversion of bushland to treeless cropland [92] increases land surface temperatures and decreases evapotranspiration, and low tree canopy cover areas cause higher land surface temperatures and higher temperatures in general [93]. Together with increasing proportion of agricultural areas, conservation areas have a negative contribution to the local climate, and furthermore, to the regional climate. Furthermore, bushland protection is vital for the conservation of flora and fauna, and for habitat conservation [91,94]. Furthermore, high AGB bushland supports, for example, the mitigation of wildfire, poor water quality, soil erosion, soil PH, air temperature and other ecosystem services of importance to the ecology, climate and wildlife [91,92,95]. Restoration of degraded areas by fencing, enrichment planting of woody plants and translocation of wildlife (browsers) to high biomass areas, agroforestry, and sustainable environmental regulation are some ways to mitigate these effects. Therefore, the trade-offs between the wildlife conservation and benefits of woody vegetation should be considered carefully in the conservation area management and land-use planning.
Although not addressed in this study, in addition to land use, natural factors, such as soil type, ground water table level, and rainfall, may contribute to land cover and AGB patterns. The soil type is typically red laterite, but parts of the landscape are characterized by sedimentary carbonites, which are drier and less fertile soils, thus introducing sparser woody vegetation. The water table level is high, especially along the Bura River Valley, enabling better tree growth. Furthermore, rainfall and mist emergence in topographically higher areas, such as Maktau Hill in LUMO Oza, may increase tree cover and height. Further studies should aim to clarify the roles of land use and natural factors on land cover and AGB in the study area.

Effect of Wildlife Fences on Biomass Distribution and Density
Fencing conservation areas is primarily done to prevent wildlife from intruding into surrounding communities and farmlands, in other words, to reduce human-wildlife conflicts [96][97][98][99]. Fences additionally help minimize wildlife poaching and the illegal extraction of other vital resources from protected areas [33] and hinder the transmission of vector-borne diseases between livestock and wildlife, as production animals and wildlife are kept separate. In Kenya, 60% of all protected areas are fully or partially fenced [35].
The ecosystem in the Taita Hills lowlands faces challenges, including livestock incursion, poaching, drought, land-use change, human-wildlife conflict, unprescribed fires, invasive species, and vegetation damage by elephants [100]. Electric and non-electric fences have therefore been constructed on the borders of the protected areas to minimize some of these challenges. The fence from Maktau to Bura Village was built in 1999 [33]. It restricts the movement of wildlife from conservation areas and hinders unauthorized access to the areas [99]. Fences also protect degraded habitats and support forest regeneration trials. Furthermore, fences around farms restrict wildlife and livestock from entering the farms.
According to our analysis of the AGB variation in the boundaries of the land-use regions (buffers), the largest differences in the percentage of zero AGB and median AGB occurred in the fenced boundaries. In the buffers 3, 5 and 7, which correspond to the boundary between the conservation areas (LUMO West and East Mramba, THWS) and the matrix, the percentage of zero AGB was considerably higher in the conservation area sides of the fence. The zero AGB pixels refer to pixels without any woody AGB, which may indicate large pressure from herbivores on woody vegetation. This is supported by the high density of wildlife close to the fence in LUMO East Mramba and THWS (Figure 4). In the buffer 7, matrix side also had higher median AGB but buffers 3 and 5 did not show significant difference in median AGB. The latter suggest that although woody vegetation is considerably less in LUMO West and East Mramba sides of the buffers, woody vegetation in both sides has similar character and median AGB. Buffer 8 matches the fenced boundary in the northern part of THWS and it is associated with a sharp transition from grassland to relatively dense bushland. This explains both the larger fraction of zero AGB pixels and the lower median AGB in the grassland side. Furthermore, buffer 10 is located in the fenced boundary between THWS and the matrix. THWS side of this boundary had slightly higher percentage of zero AGB than matrix side, similar to other conservation areas but the difference was smaller. However, THWS side of the buffer had significantly higher median AGB. This can be explained by the presence of riverine forest in that side of the boundary with greater AGB. Furthermore, the matrix in this area lies on a flood plain dominated by cropland interspersed with bushland in contrast to grassland in the THWS side, which may explain this difference.
Among the non-fenced boundaries, buffer 9 had the most apparent difference between the two sides of the boundary. This buffer corresponds to the northern boundary of THWS with rapid change from bushland to cropland-dominated area within the matrix. Matrix-side had clearly more zero AGB pixels corresponding to cropland and lower median AGB. Although not fenced, this boundary follows a road, which makes it clearly visible. Furthermore, the fence south of the area protects it from herbivores in THWS. In addition, statistically significant differences in median AGB were observed in the non-fenced buffers 1 and 6. In the boundary between TWNP and LUMO West Mramba (buffer 1), bushland in the northern part of TWNP has a relatively high AGB compared with grassland-dominated West Mramba, which explains higher median AGB in the TWNP side. Buffer 6 corresponds to the boundary between two LUMO East Mramba and THWS. Slightly higher AGB in the THWS side could relate to higher grazing pressure in LUMO East Mramba. However, differences in these two unfenced boundaries are very small in comparison to the fenced boundaries with obvious differences.
Our results support our hypothesis that fences play a role in the distribution of wildlife and livestock, and woody AGB patterns in the landscape. This creates sharp land cover transitions to the fenced boundaries of the land-use regions. The conservation and grassland sides of the buffers 5, 7 and 8 experience high pressure from wildlife and cattle while pressure is particularly low in the matrix sides of buffers 7 and 8 with fewer cattle (Figure 4). In buffers 3 and 10, the difference in herbivore density between the conservation areas and the matrix were not as evident at the time of counting. However, free ranging wildlife are constantly moving based on resource conditions. In general, fencing can increase the wildlife population in the conservation areas and enhance biodiversity conservation [101][102][103]. However, an increased abundance of (mega)herbivores [104] reduce biomass densities due to tree mortality caused by browsing. The browsers suppress woody plant recruitment in the grassland and have a long-term impact on their growth and mortality rates [105]. This is particularly true for non-selective feeders, such as elephants, who debark trees and thus suppress recruitment and vegetation generation. According to Ogutu et al. [52], the landscape experienced a moderate growth in elephant density between 1977 and 2016. Similar pressure on woody plants was observed during 1970-1973, when the elephant population was large [55]. The problem is further aggravated by the fence, which restricts wildlife dispersal, and hence, reduces the ecosystem's resilience [98]. Thus, fencing combined with heavy browsing may reduce the biomass in conservation areas.

Conclusions
Taita Hills lowland savannah landscape, similar to other typical African savannah biomes, exhibits multi-use functionality, which results in heterogeneous land cover. We combined field data with ALS metrics to predict a woody AGB map in the study area and created a land cover map using Google Earth Engine. AGB densities in the region were comparatively low and influenced by wildlife conservation. The highest AGB densities were observed in the forest class (riverine forest) in THWS in the conservation area. Greater AGB densities were also found in the bushland in the matrix, LUMO Oza, and southern parts of THWS. The western parts of the landscape dominated by grassland and influenced by wildlife conservation and livestock grazing had a lower woody AGB density. Wildlife and livestock densities in the conservation area are high compared to the matrix. Bushland and cropland dominate the matrix, which support the livelihood of community members through farming and other livelihood options (fuelwood, etc.). The electric fence restricts the movement of wildlife, creating grassland within protected areas and contributing to the low densities of woody AGB. In addition to human-wildlife conflict mitigation, fencing also influences the spatial distribution and density of woody AGB in a multi-use savannah landscape. Further investigating the effect of wildlife and livestock fencing on land cover and biomass (including grass biomass) in multi-use savannah landscapes at various spatial and temporal scales is important. Furthermore, our results need to be scaled up and contributions of livestock management and conservation areas to climate change require investigation. The impact of wildlife conservation on land cover change, and plant species diversity and composition also deserve further investigation.

Predictor Description
CC.first First returns > 3.5 m/Total first returns * 100 CC.all All returns > 3.5 m/Total all returns * 100 CC.all.first All returns > 3.5 m/Total first returns * 100 CC.first.mean First returns above mean/Total first returns * 100 CC.all.mean All returns above mean/Total all returns * 100 CC.all.mean.first All returns above mean/Total first returns * 100 CC.first.mode First returns above mode/Total first returns * 100 CC.all.mode All returns above mode/Total all returns * 100 CC.all.mode.first All returns above mode/Total first returns * 100 All height variables (beginning with 'H') were calculated separately using first and last pulse returns, which are indicated by the prefix 'FR_' or 'LR_', respectively. All canopy variables (beginning with "CC") were calculated using all returns only.