10urn:lsid:arphahub.com:pub:6BEA2F36-F41A-5C4C-8F71-33F322F204B5Nature ConservationNC1314-69471314-3301Pensoft Publishers10.3897/natureconservation.16.93129312Research ArticlePlantaeEcology & Environmental sciencesCenozoicAfricaAnalysing fragmentation in vulnerable biodiversity hotspots in Tanzania from 1975 to 2012 using remote sensing and fragstatsOjoyiM.M.1OdindiJohnodindi@ukzn.ac.za2MutangaO.1Abdel-RahmanE.M.1School of Agricultural, Earth and Environmental Sciences, Pietermaritzburg Campus, University of KwaZulu-Natal, Scottsville P/Bag X01, Pietermaritzburg 3209, South AfricaUniversity of KwaZulu-NatalPietermaritzburgSouth AfricaDepartment of Agronomy, Faculty of Agriculture, University of Khartoum, SudanUniversity of KwaZulu-NatalPIETERMARITZBURGSouth Africa
Corresponding author: John Odindi (Odindi@ukzn.ac.za)
Academic editor: M. Kleye
2016181120161619375B25FF86-FB40-384E-FF85-1A098133FFDD823FB494-A057-479B-9EFF-445F81F6E03A1674662305201630102016M.M. Ojoyi, John Odindi, O. Mutanga, E.M. Abdel-RahmanThis is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.http://zoobank.org/823FB494-A057-479B-9EFF-445F81F6E03A
Habitat fragmentation is a threat to conservation of biodiversity hotspots in the Morogoro region, Tanzania. However, on-going research on fragmentation has not kept pace with temporal lapses and how individual species respond to habitat transformation and heterogeneity. This study sought to model spatial and temporal fragmentation patterns. Cloud free multi-temporal Landsat imagery with similar spectral resolution were acquired in the same season in 1975, 1995 and 2012. The images were used to characterize the biophysical landscape characteristics and a range of metrics used to quantify the magnitude of fragmentation. Patches and classes in the landscape were assessed using Fragstats, a spatial statistics program useful in computing landscape metrics. Results show that patch number was higher in dense forest and woodland than in less dense forest and grassland in 1975, 1995 and 2012 while the interspersion Juxtaposition Index (IJI) ranged between 0 (for clumped patches) and 100 (for grassland). In 1975 and 1995, the grassland habitat had the highest IJI while in 2012 less dense forest had the highest IJI. The Games-Howell test showed a significant fragmentation trend in less dense forests class (p≤0.05). Generally, the study indicates a high fragmentation pattern in the vulnerable tropical eastern arc mountain region of East Africa. This finding demonstrates the value of remotely sensed data in understanding the impact of anthropogenic processes on natural landscape transformation. Furthermore, the study provides a basis for informed conservation policy design and implementation in the region.
Ojoyi MM, Odindi J, Mutanga O, Abdel-Rahman EM (2016) Analysing fragmentation in vulnerable biodiversity hotspots in Tanzania from 1975 to 2012 using remote sensing and fragstats. Nature Conservation 16: 19–37. doi: 10.3897/natureconservation.16.9312
Introduction
Habitat fragmentation, an indication of habitat transformation, degradation and loss is a great concern globally (Fahrig 2001, Fahrig 2003, Cushman et al. 2010, Forman and Godron 1986). It refers to habitat breakages or the degree of patchiness of a habitat, mainly as a result of anthropogenic impacts (Fahrig 2001, Fahrig 2003, McGarigal and Cushman 2002, Wiens 1995, Neel et al. 2004). Generally, habitat fragmentation interferes with the structural configuration of ecosystems and their ecological functioning (Abdullah and Nakagoshi 2007, Echeverria et al. 2006, Echeverria et al. 2008, Iida and Nakashizuka 1995, Forman and Godron 1986). Specifically, fragmentation reduces total habitat area, making species highly vulnerable to endemism and extinction (Yen et al. 2005, Murcia 1995, Aguilar et al. 2008, Yen et al. 2005). Hence fragmentation has long term impacts on species numbers (Aguilar et al. 2008) and species abundance (Fahrig 2003, Debinski and Holt 2001, Yen et al. 2005, Forman and Godron 1986) as it exposes natural ecosystems to external risks such as parasitism and dominance of invasive species (Wiens 1995).
Habitat fragmentation is an explicit challenge to conservation in the tropics (Vogelmann 1995). It is considered a major cause of species loss (Pelkey 2000, Adams et al. 2003, Bjørndalen 1992, Burgess et al. 2002, Burgess et al. 2001, Yen et al. 2005, Forman and Godron 1986). In Africa, approximately 310,000 hectares of forest is annually converted to agriculture, while 200,000 hectares is converted into woodlands, major causes of fragmentation (Achard 2002). Fragmentation acts synergistically with other factors like effects of solar radiation and open niches that lead to dominance of other invasive species. Consequently, native vegetation species are exposed to higher risks of extinction with a decline in the percentage area required for their survival (Rutledge 2003).
Ecosystems in Morogoro region, Tanzania contribute to the world’s climate regulation through large carbon stores (Burgess et al. 2007, Swetnam et al. 2011). These forests are also characterized by high levels of endemism and many species are vulnerable to extinction (Swetnam et al. 2011, Brooks et al. 2006, Myers et al. 2000). Increased anthropogenic disturbances in particular pose significant threats to their long term conservation (Hall et al. 2009, Hall 2009, Newmark 1998). Between 1955 and 2000 for instance, forest cover declined from 300 km2 to 220 km2 (Burgess et al. 2007). Despite the area’s global importance, few studies have been conducted with a focus on its spatial heterogeneity (Newmark 1998). Furthermore, mechanisms by which natural habitats respond to spatial heterogeneity across diverse fragmenting ecosystems remain largely unexplored (Swetnam et al. 2011, Yanda and Shishira 1999). Individual habitats may differ in their degree of response to fragmentation as the robustness of fragmentation may vary (Fahrig 2003, Neel et al. 2004, McGarigal 2006, Echeverría 2007). For instance, due to differences in their structural complexity and biological processes, what could be termed as fragmentation in homogeneous landscapes may be interpreted differently in a heterogeneous landscape (Murcia 1995, Fischer and Lindenmayer 2007, Wiens 2000). In this study, we tested the spatial extent and magnitude of fragmentation in four vulnerable habitats subjected to fragmentation in the region. Remote sensing was applied due to its increasing popularity in quantifying spatio-temporal patterns in diverse landscapes (Ojoyi et al. 2016, Nagendra et al. 2004, Lung and Schaab 2006, Southworth et al. 2002, Fjeldså 1999). Specifically, the study pursued the following objectives: (1) to investigate multi-temporal magnitude of fragmentation in diverse habitats; (2) to quantify the intensity of habitat fragmentation in each of the habitats.
Study area
Similar to this case study, most rich biodiversity hotspots in Tanzania are geographically located in the Eastern Arc Mountains (Burgess et al. 2007, Myers et al. 2000, Hall et al. 2009, Hall 2009, Newmark 1998, Olson and Dinerstein 1998). In this study, we selected a section of Morogoro region dominated by four major habitat types (Figure 1). The choice of the study location was based on previous ecological studies (e.g. Burgess et al. 2002, Burgess et al. 2001, Hall 2009, Luoga et al. 2000, Yanda and Shishira 1999) that attributed species losses to fragmentation. The study area is characterized by sub-montane (with trees 30-50m tall), montane (with trees 15-30m tall) and upper montane (with trees 15-20m tall) forest at 1200-1500, 1500-2100 and >2100m asl, respectively. Generally, forest density and height varies with elevation and aspect, with dense canopy dominating lower altitude and elfin forests dominating ridges above 1900m asl (Lovett et al. 1996). Stunted grass patches are also common at high altitudes. According to Burgess et al. (2002), the potential of closed natural forest cover is about 500km2, however, this has been reduced from 300km2 in 1955 to 230 km2 in 2001, with most decline recorded at 600-1600 m asl outside protected areas. The loss is mainly attributed to increasing population growth, estimated at about 2.5-3% per annum (Lovett 1996). The study area comprises four main habitats; woodland, dense forest, less dense forest and grassland. In this study, woodland is described as woody vegetation with scattered foliage cover (less than 30%) with mature stands of less than five meter tall while less dense forest consists of fields and patches with trees of more than six meters tall, with crown cover of less than 30%. The dominant tree species in the region include: Bersamaabyssinica, Cassipoureamalosana, Cornusvolkensii, Cussonialukwangulensis, C.spicata, Dombeyatorrida, Draceanaafromontana, Garciniavolkensii, and Xymalosmonospora. Bamboo thickets form dense stands of Sinarundinariaalpina 12-15 m tall and 15 cm diameter (Bjørndalen 1992, Shirima et al. 2011, Lovett 1993). The grassland habitat consists of Panicumlukwangulense and Andropogonthystinus with scattered trees of Agauriasaliciflora, Adenocarpusmannii, Myricasalicifolia and Berberis sp. thought to have replaced upper montane forest due to fires (Adams et al. 2003, Bjørndalen 1992). Kitulanghalo forest is located between Morogoro and Dar es Salaam within Morogoro region. The woodland forms part of Miombo woodland which covers 90% of the total forested ecosystem (Mugasha et al. 2013, Munishi et al. 2010). It is dominated by Brachystegia, Isoberlinia, Julbernardia, Pterocarpusangolensis, Afzeliaquanzesis and Albizia species (Chander 2009). The area falls within a semi-natural Miombo woodland which receives less than 1000 mm of rainfall per annum. Their proximity to Morogoro urban area increases their susceptibility to anthropogenic activities, altering their functioning and sustainable management (Mugasha et al. 2013).
88A95AB2-F2F2-53D4-907D-C5BFC47C784F
A 1975 Landsat composite overlaid on the study area.
Landsat MSS (20/08/1975), Landsat TM (30/09/1995) and Landsat ETM+ (20/07/2012) imagery with better visualization (less than 15% cloud cover) from the Global Landcover Facility (http://www.landcover.org) were selected for the study. Datum was set to WGS 84 and referenced to
Universal Transverse Mercator
(UTM) zone 37 South. All images were orthorectified using ground validation points,
Digital Elevation Model
(DEM) and aerial photos as a reference. Landsat images were resampled to a common resolution pixel (30 × 30 m) using bilinear resampling to ensure consistency in all image scenes. First order polynomial transformation was applied at image registration to correct for any shifts. It was deemed necessary to simulate atmospheric interactions between the sun and sensor pathways for the imagery used. Therefore, a radiative transfer model in
Atmospheric and Topographic Correction
(ATCOR) module in Erdas Imagine 2013 was used for atmospheric correction. ATCOR masks haze, cloud, water and enhances pixel visibility. In this study, we used the
MODerate resolution atmospheric TRANsmission
(MODTRAN) code to retrieve the atmospheric parameters for ATCOR from the look-up table as ground-based reflectance and atmospheric data were unavailable. Digital number values were then converted to reflectance based on metadata provided with the Landsat images (Chander et al. 2009, Guanter et al. 2009, ERDAS and Geosystems 2011).
Image classification
A supervised maximum likelihood classifier was adopted for classification (Liu et al. 2002, Manandhar et al. 2009, Tseng et al. 2008, Xi 2007). The technique is based on statistical probability that assigns pixel values to the category with the highest likelihood (Aldrich 1997, Dean and Smith 2003, Ince 1987). Spectral signatures were created and applied in categorizing similar pixels in the entire image using eight polygons representing training data sets for each habitat class. A color composite of 3, 4 and 5 bands were used to facilitate visual interpretation while the Gaussian distribution function was applied in the stretching process. The image was classified into four class categories namely: Woodland, Grassland, Dense forest and Less dense forest. A total of 82 field ground data points, archival high resolution aerial photographs, interviews and expert opinion were used to validate the classified images. Confusion matrices were then created to compare reference data with the maximum likelihood prediction and for calculation of the
overall accuracy
(OA),
producer’s accuracy
(PA) and
user’s accuracies
(UA). Overall accuracy is a percentage (%) between correctly classified classes and the total number of test reference data, while producer’s accuracy is the probability of a specific class being correctly classified. User Accuracy is the possibility that a sample of a specific class represents the category on the ground.
Modelling habitat fragmentation
Fragstats metrics were extracted from all processed Landsat images. Fragstat metrics offer a distinct capacity to determine a landscape’s spatial configuration, hence valuable in understanding landscape change arising from fragmentation (Cushman 2006, Jorge and Garcia 1997, Saikia et al. 2013, Millington et al. 2003). All classified images were converted to ASCII format in ArcGIS 10.2. A C-program, a raster version inbuilt within Fragstas that accepts ASCII image files was applied using the eight cell rule. The ASCII format scenes were imported into Fragstats and ASCII built-in-algorithm selected for running the Fragstats model. Three multi-level structure metrics were selected at patch, class and landscape level (McGarigal and Cushman 2002). Metrics relevant in explaining the magnitude and extent of fragmentation were then selected from the 1975, 1995 and 2012 image scenes. A total of 155 samples were randomly selected and extracted. As recommended by McGarigal and Cushman (2002) two metrics i.e. perimeter area relationship and patch area were statistically used to test the magnitude of fragmentation. Mann-Whitney U and Post hoc ANOVA tests were used to evaluate differences among patch areas in all the years. The Games-Howell was used to determine forest fragmentation. The indices used in this study are briefly described in Table 1.
Fragmentation Indices used in the current study.
Fragstats Metrics
Description
Patch Density
(PD)
Number of patches of the corresponding patch type.
Largest Patch Index
(LPI)
An index used to quantify the percentage of total landscape area characterized by the largest patch.
Edge density
(ED)
Used to assess edge length per unit area.
Patch Number
(NP)
A measure of the magnitude of fragmentation of patches
Interspersion Juxtaposition Index
(IJI)
A measure of adjacency of patches determined by dividing the length between patch edge by the number of patches within a landscape. Values approaching 0% indicate that a patch is adjacent to only one other patch and 100% indicate that a patch is in similar proximity to multiple patches within a landscape.
Patch Area
(MN)
The sum across all patches in the landscape of the corresponding patch metric values, divided by the total number of patches. Expressed in hectares.
Perimeter Area Ratio
(PARA)
Refers to the ratio of the patch perimeter (m) to area (m2).
Total Area
(CA)
Refers to the sum of areas (m2) of all patches for the patch type.
Percentage of Landscape
(PLAND)
Useful in computing the proportional abundance for each of the patch type across the landscape.
ResultsClassification and accuracy assessment
The overall accuracy for 1975, 1995 and 2012 image scenes was 78.26%, 84% and 76.54% respectively (Table 2). Changes in total area coverage were observed in all years (Figures 2a, b and c).
637526DD-C1B3-572B-9E33-19D4F5D2AE08
Land use land cover
(LULC) maps in 1975 (a), 1995 (b) and 2012 (c).
https://binary.pensoft.net/fig/111479
Accuracy assessment tests (
Producer’s Accuracy -
PA,
User’s Accuracy -
UA).
Habitat Class
1975
1995
2012
PA (%)
UA (%)
PA (%)
UA (%)
PA (%)
UA (%)
Dense Forest
100
75
100
100
80.77%
95.45%
Less Dense Forest
66.67
100
100
100
100.00%
60.00%
Woodland
66.67
100
66.67
66.67
75.00%
58.54%
Grassland
100
100
100
100
84.62%
91.67%
Overall Accuracy
78.26
84
76.54%
Kappa co-efficient
0.7416
0.812
0.7284
Change detection
The study findings showed substantial land modification in most of the cover types during the study period i.e. decline in dense forest (31, 675.70 hectares) and less dense forest (by 11, 267.38 hectares) and increase in grassland (21, 230.01 hectares). However, changes in areas covered by woodland were inconsistent, i.e. increase by 15,884.46 hectares between 1975 and 1985 and decline by 8, 182.03 between 1985 and 2012) – Figure 2.
Fragmentation trendsTemporal variability in fragmentation
Dynamic fragmentation trends were observed (Table 3). Patch number was relatively higher in dense forest and woodland in 1975, 1995 and 2012 than in less dense forest and grassland. The highest
percentage of landscape
(PLAND) were recorded in less dense forest than the rest of the habitats while woodland and less dense forest habitats had the highest edge density (Figure 3a–c). Furthermore, dense forest showed the most declining patch number during the study period. An analysis of the
largest patch index
(LPI) showed that less dense forest had the highest LPI, while woodland, dense forest and grassland had the least values, below five. Woodland had the highest PARA compared to the rest of the habitat types (Table 3).
9B55A2BF-BBC4-5A2D-A72E-66CBCB1FDA5B
Temporal patterns of total area coverage (A), percentage of landscape (B) and edge density (C).
https://binary.pensoft.net/fig/111480
Patch area compared by Mann-Whitney Tests.
Class
Year
z-value (1975–1995)
Prob > |z|
z-value (1995–2012)
Prob > |z|
Dense forest
1975
9.495***
0
-6.872
0.1895
1995
2012
Grassland
1975
13.680***
0
-7.441***
0
1995
2012
Less dense forest
1975
16.728***
0
-8.268***
0
1995
2012
Woodland
1975
-16.63***
0
2.461***
0
1995
2012
Spatial variation in fragmentation
Study findings indicated a higher probability of dispersion linked to woodland and less dense forest.
Interspersion Juxtaposition Index
(IJI) ranged between 0 (for clumped patches) and 100 (for grassland). In 1975 and 1995, the grassland habitat had the highest IJI while in 2012, less dense forest had the highest IJI. The
interspersion juxtaposition index
(IJI), was useful in characterizing the degree of adjacency for each patch type e.g. Burgess et al. (2007). Additionally, the largest patch number and mean patch area was evident in dense forest in 1975 and woodland in 2012 (Figure 4).
BC233461-9CA3-5EB8-82F8-AB37BD4C6D95
Spatial variability in number of patches (A), Interspersion Juxtaposition Index (B), Largest patch index (C), Patch density (D), Mean patch area (E) and Perimeter area ratio (F) in 1975, 1995, and 2012. DF
- dense forest
, LDF
- less dense forest
, WD
- woodland
, GR
- grassland
.
https://binary.pensoft.net/fig/111481Mann-Whitney test results
Mann-Whitney tests were applied to the data. Mann-Whitney test results showed distinct differences in patch area (p<0.01) as summarized in (Table 3). These results were strong indicator of a rapidly fragmenting landscape.
Games Howell test results for perimeter area relationship
Game-Howell test is ideal for unequal sample sizes charaterised by heterogeneity and has been widely used in vegetation mapping that include taxonomic profiles in the Atlantic and Caatinga biomes of northestern Brazil (Pacchioni et al. 2014), forest transformation in Uluguru mountains (Ojoyi et al. 2015), the effect of fire on Penderosa pine forest density, canopy cover, tree size and basal area (Stephens et al. 2015) and shrub density in Zegros forest, southwest Iran (Askari et al. 2013). Games Howell test results showed significant patterns of fragmentation between 1975 and 1995 in all habitats (p≤0.05). In 1975 and 2012, the trend was significant in less dense forest and woodland (p≤0.05), while in 1995 and 2012, the trend was significant in grassland, dense forest and less dense forest (p≤0.05) (Table 4). A highly significant trend with perimeter area relationship was evident with less dense forest across the years.
Games-Howell tests for the mean
parameter area ratio
(PARA) in 1975, 1995, 2012.
Class
Mean
p value
1975
1995
2012
1975 vs 1995
1975 vs 2012
1995 vs 2012
Grassland
565.28
606.21
560.00
0.0001
0.596
0.0001
Dense forest
498.12
549.14
483.7
0.0001
0.3
0.0001
Less dense forest
496.29
563.06
529.5
0.0001
0.0001
0.0001
Woodland
498.58
535.43
534.3
0.0001
0.0001
0.893
Discussion
This study showed a progressive fragmentation at both spatial and temporal domains. Variability in responses to fragmentation was also noted for different habitats. Fragmentation in the area is not only dependent on topography but also adjacency to land for agriculture, urbanization/settlement and infrastructure development, which are considered key drivers of landscape transformation in the region. All these anthropogenic activities contribute to habitat losses and species decline. Implications on the landscape are presented with a reflection on policy and future management.
Habitat modification
There was a transformation in habitat extents within the study area. Significant losses were recorded for dense forest (31, 675.70 hectares) and less dense forest (by 11, 267.38 hectares), however, there was a steady increase in areas covered by grassland. Based on field study observations, these changes can be attributed to expanding agricultural fields and increased exploitation of timber and non-timber products to meet the increasing urbanization demand in Morogoro district. This finding is in agreement with Burgess et al. (2001) and Burgess et al. (2002) who found a substantial decline of dense forest in the Uluguru mountains due to urbanization and agricultural (Burgess et al. 2002, Burgess et al. 2001). In other parts of Tanzania, related studies established effects of reduced tree density to land modification (Yanda and Shishira 1999, Munishi et al. 2010). Habitat modification could also be attributed to general population increase in non-urbanized areas, also known to influence its spatial configuration (Fischer and Lindenmayer 2007).
Spatial and temporal variation
As aforementioned, there was a general decrease in area covered by dense and less dense forest habitat. A decreasing trend in the extent of total habitat coverage relates to deleterious fragmentation as effects of habitat fragmentation are dependent on habitat size (Fahrig 2003). Furthermore, perimeter-area results in this study show distinct differences in woodland and grassland habitats. In most instances, high perimeter-area relationship characterizes rapid rate of fragmentation underlying the two landforms e.g. Jha et al. (2005) and McGarigal (2006). Woodland habitat displays a patchy type of deforestation, shown by an increased patch number between 1975 and 2012. The slight decline in patch number can be attributed to the strong traditional leadership forest maintance authority in the 1970s, a responsibility that has now been taken over by the Tanzanian Government that permits logging and farm allocations. Dynamics in mean patch area were observed in the woodland and less dense forest. Notable was the gradual decrease in patch size, while patch number increased by 412 and 391 in dense forest and woodland respectively, an indication of fragmentation patterns in the area earlier observed by Jha et al. (2005). On the other hand, patch area was ideal in characterizing distinct areas with analogous environmental conditions, where patch boundaries are distinguished by discontinuities in environmental character states relevant to the organism or ecological phenomenon under consideration. A combination of
patch density
(PD), PARA and mean nearest neighbor distance are considered profound in estimation of the extent of fragmentation in each of the habitats analyzed (Jha et al. 2005). Patch density and PARA are regarded as important in fragmentation assessments, particularly in natural ecosystems because they have a strong influence on ecosystem functioning and ecological processes (McGarigal 2006).
Similarly, a distinct variation in patch number was observed. Woodland and less dense forest had the highest patch number across the years. This can be attributed to the great extent of fragmentation resulting from natural resource exploitation. Furthermore, their vicinity to Morogoro town and management by local authorities may be possible drivers increasing their susceptibility to fragmentation (Fahrig 2001, Fahrig 2003, Wiens 1995, McGarigal 2006, Fischer and Lindenmayer 2007). The woodland habitat had a relatively greater patch density, signifying higher spatial heterogeneity. In addition, the largest patch index was associated with less dense forest while least values were associated with the grassland habitat. This provided information on least and most fragmented landscapes, a good indicator of minimum area requirements for species survival (McGarigal 2006). In addition, the largest patch index, another good indicator for species survival was significant in the less dense forest compared to the rest of the habitats (Rutledge 2003).
Dense forest and woodland had the greater edge density. This could be attributed to increased exposure to farmlands and settlements prevalent in the area. Edge effects characterize the biophysical state of ecosystems at the periphery or in the neighborhood. This is because increased habitat fragmentation exposes habitat to edge effects, compromising the ability of an ecosystem to provide relevant goods and services (Murcia 1995). This limits a habitat’s long-term ability to sustain a population as it intensifies species mortality rate (Fahrig 2003). It also influences occurrence of native species populations (Murcia 1995) and ensures that the interaction of species in disturbed environments remains restricted, advancing their mortality risk (Rutledge 2003). Related literature also found a high intensity of fragmentation associated with more edge effects through exposure of contiguous habitats to solar radiation and soil moisture to drier heat conditions (Rutledge 2003).
Games-Howell test results showed a significant level in the perimeter area relationship (p≤0.05). This could be explained by the fact that less dense forest adjoins dense forest, taking up regions dominated by woodland. It is also possible that the on-going fragmentation is a major driver of conversion of dense forest and woodland to less dense forest. Potential socio-economic drivers could be a result of the expanding Morogoro town and increasing agricultural fields in the adjacent local regions. Similarly, other studies showed how adjoining activities influence intact habitat ecosystems as a result of their structural configuration (Echeverría et al. 2007).
Drivers to habitat fragmentation and conservation implications
Anthropogenic activities significantly influence habitat fragmentation in the region. For instance, extensive farming and urban growth are possible drivers to habitat modification and fragmentation. The area has a conducive montane climate that supports subsistence farming, a prevalent socio-economic activity in the region (Burgess et al. 2007, Swetnam et al. 2011, Yanda and Shishira 1999). This seems to significantly influence all the four habitats. Increasing population growth and consequent increase in settlement and farmlands may have extirpated important fauna and flora in the Ulugurus (Bjørndalen 1992, Burgess et al. 2002, Burgess et al. 2001, Burgess et al. 2007, Swetnam et al. 2011, Hall 2009, Yanda and Shishira 1999) – Figure 5. Habitat fragmentation in the study area can also be attributed to a complex nexus of socio-economic processes (Kessy et al. 2016, Rosales 2008). These processes act at various scales i.e. international (global forest products market growth, commercialization and urbanization), national (changing population, growing local markets and national legislation and governance) and local conditions (livelihoods and levels of poverty) (Wehkamp et al. 2015, FAO 2007, Daly and Farley 2004, Czech 2013). Kessy et al. (2016) for instance notes that local and international demand for timber and agricultural commodities in a globalizing world are major drivers to forest fragmentation in the area. Globalization, with its characteristic scramble by the developing countries to increase their market share on the global marketplace has increased pressure on existing forests and forest land (Hecht and Saatchi 2007, Rosales 2008).
6CC23193-1F42-5C87-95DA-3648C38C24D1
Drivers to fragmentation, note the settlements in the valley and cleared forest in the background and fore ground for crop farming and grazing, respectively (A) and small scale maize and banana fields within the forest in (B).
https://binary.pensoft.net/fig/111482
To forestall some of the problems earlier highlighted, the study area, identified as biodiversity hotspots with important ecological functions such as groundwater recharge, surface flow and animal habitat need to be protected from the impacts of land modification and fragmentation. Implications of habitat modification and fragmentation in Morogoro region can be better deciphered through the impact on habitat structure and species losses. The increased habitat losses, mainly attributed to anthropogenic factors may negatively influence genetic diversity and lead to losses of potentially useful genes originally accommodated in intact areas (Ojoyi et al. 2015, Burgess et al. 2007, Swetnam et al. 2011, Hall 2009, Yanda and Shishira 1999, Shirima et al. 2011). Therefore, we recommend that mitigation measures should be adopted to ensure protection and management of these fragmenting habitat ecosystems. To optimize mitigation measures, the adverse effects of habitat modification and fragmentation need to be understood by all stakeholders. In addition, policy measures and sustainable bottom-up approaches to management and conservation of forest resources should be instituted in the region.
Conclusions
Distinct differences in magnitude of fragmentation were evident across the four habitat categories. The study findings show that fragmentation was highest in less dense forest. Subsistence farming, increasing human population and urban growth are thought to be key drivers to habitat modification and fragmentation, hence it is concluded that anthropogenic processes are the major drivers to habitat fragmentation in the area. The fragmenting landscape is expected to significantly influence floral and faunal vulnerability, likely to compromise the area’s ability to among others assimilate organic carbon and to supply socio-economic and environmental goods and services. It is therefore necessary that the study area, and indeed the entire eastern arc mountains region be protected from the impacts of land modification and fragmentation. The study further underscores the value of satellite imagery in concert with relevant reference data in understanding spatio-temporal transformation of vulnerable landscapes arising from anthropogenic processes.
Acknowledgements
The current study was financed by UNESCO L’Oreal Foundation for Women in Science, The International Development Research Centre (IDRC), Ontario Canada and University of KwaZulu-Natal Post-Graduate sponsorship. Training by the University of KwaZulu-Natal, South Africa and Faculty of Geo-Information Science and Earth Observation, University of Twente, Netherlands is greatly appreciated. Field support provided by Sokoine University of Agriculture, forestry department, Tanzania Government Forestry departments in Morogoro region and the Wami/Ruvu Basin Office is highly acknowledged.
ReferencesAbdullahSANakagoshiN (2007) Forest fragmentation and its correlation to human land use change in the state of Selangor, peninsular Malaysia.241(3): 9–48. doi: 10.1016/j.foreco.2006.12.016AchardFEvaHDStibigH-JMayauxPGallegoJRichardsTMalingreauJ-P (2002) Determination of deforestation rates of the World’s humid tropical forests.297(5583): 999–1002. doi: 10.1126/science.1070656AdamsMCooperJCollarN (2003) Extinct and endangered (E&E’) birds: a proposed list for collection catalogues.123: 338–54.AguilarRQuesadaMAshworthLHerrerias-DiegoYLoboJ (2008) Genetic consequences of habitat fragmentation in plant populations: susceptible signals in plant traits and methodological approaches.17(24): 5177–5188. doi: 10.1111/j.1365-294X.2008.03971.xAldrichJ (1997) R. A. Fisher and the making of maximum likelihood 1912–1922.12(3): 162–176. doi: 10.1214/ss/1030037906AskariYSaeiKERezaeiD (2013) Evaluation of different sampling methods to study of shrub density in Zagros forest.3(2): 121–128.BjørndalenJE (1992) Tanzania’s vanishing rain forests – assessment of nature conservation values, biodiversity and importance for water catchment.40(1): 313–334. doi: 10.1016/0167-8809(92)90100-PBrooksTMMittermeierRAda FonsecaGABGerlachJHoffmannMLamoreuxJFMittermeierCGPilgrimJDRodriguesASL (2006) Global biodiversity conservation priorities.313(5783): 58–61. doi: 10.1126/science.1127609BurgessNDButynskiTMCordeiroNJDoggartNHFjeldsåJHowellKMKilahamaFBLoaderSPLovettJCMbilinyiBMenegonM (2007) The biological importance of the Eastern Arc mountains of Tanzania and Kenya.134(2): 209–231. doi: 10.1016/j.biocon.2006.08.015BurgessNDoggartNLovettJ (2002) The Uluguru Mountains of eastern Tanzania: could we lose their globally significant biodiversity values?36: 140–152. doi: 10.1017/S0030605302000212BurgessNRomdalTSRahnerM (2001) Forest loss on the Ulugurus, Tanzania and the status of the Uluguru bush shrike Malconotus alius.8: 89–90.ChanderGMarkhamBLHelderDL (2009) Summary of current radiometric calibration coefficients for Landsat MSS, TM, ETM+, and EO-1 ALI sensors.113(5): 893–903. doi: 10.1016/j.rse.2009.01.007CushmanSA (2006) Effects of habitat loss and fragmentation on amphibians: A review and prospectus.128(2): 231–40. doi: 10.1016/j.biocon.2005.09.031CushmanSAComptonBMcGarigalK (2010) Habitat fragmentation effects depend on complex interactions between population size and dispersal ability: Modeling influences of roads, agriculture and residential development across a aange of life-history characteristics. In: CushmanSHuettmannF (Eds) . Springer, Tokyo, 350–369.CzechB (2013) New Society Publishers, Gabriola Island, 400 pp.DalyHEFarleyJ (2004) Island Press, Washington, 484 pp.DeanAMSmithGM (2003) An evaluation of per-parcel land cover mapping using maximum likelihood class probabilities.24(14): 2905–20. doi: 10.1080/0143116021055910DebinskiDMHoltRD (2001) A Survey and overview of habitat fragmentation experiments.14(2): 342–355. doi: 10.1046/j.1523-1739.2000.98081.xEcheverríaCCoomesDAHallMNewtonAC (2008) Spatially explicit models to analyze forest loss and fragmentation between 1976 and 2020 in southern Chile.212: 439–449. doi: 10.1016/j.ecolmodel.2007.10.045EcheverríaCCoomesDSalasJRey-BenayasJMLaraANewtonA (2006) Rapid deforestation and fragmentation of Chilean temperate forests.130: 481–494. doi: 10.1016/j.biocon.2006.01.017EcheverríaCNewtonACLaraABenayasJMRCoomesDA (2007) Impacts of forest fragmentation on species composition and forest structure in the temperate landscape of southern Chile.16(4): 426–439. doi: 10.1111/j.1466-8238.2007.00311.xERDAS and Geosystems (2011) . Leica Geosystems, Heerbrugg.FahrigL (2001) How much habitat is enough?100(1): 65–74. doi: 10.1016/S0006-3207(00)00208-1FahrigL (2003) Effects of habitat fragmentation on biodiversity.34: 487–515. doi: 10.1146/annurev.ecolsys.34.011802.132419FAO (2007) . FAO, Rome.FischerJLindenmayerDB (2007) Landscape modification and habitat fragmentation: a synthesis.16(3): 265–280. doi: 10.1111/j.1466-8238.2007.00287.xFjeldsåJ (1999) The impact of human forest disturbance on the endemic avifauna of the Udzungwa Mountains, Tanzania.9: 47–62. doi: 10.1017/S0959270900003348FormanRTGodronM (1986) Wiley, New York, 603 pp.GuanterLRichterRKaufmannH (2009) On the application of the MODTRAN4 atmospheric radiative transfer code to optical remote sensing.30(6): 1407–1427. doi: 10.1080/01431160802438555HallJM (2009) University of Florida, Florida, 131 pp.HallJBurgessNDLovettJMbilinyiBGereauRE (2009) Conservation implications of deforestation across an elevational gradient in the Eastern Arc Mountains, Tanzania.142(11): 2510–2521. doi: 10.1016/j.biocon.2009.05.028HechtSSaatchiS (2007) Globalization and Forest Resurgence: Changes in Forest Cover in El Salvador.57(8): 663–672.IidaSNakashizukaT (1995) Forest fragmentation and its effect on species diversity in sub-urban coppice forests in Japan.73(1–3): 197–210. doi: 10.1016/0378-1127(94)03484-EInceF (1987) Maximum likelihood classification, optimal or problematic? A comparison with the nearest neighbour classification.8(12): 1829–1838. doi: 10.1080/01431168708954819JhaCSGoparajuLTripathiAGharaiBRaghubanshiASSinghJS (2005) Forest fragmentation and its impact on species diversity: an analysis using remote sensing and GIS.14(7): 1681–1698. doi: 10.1007/s10531-004-0695-yJorgeLABGarciaGJ (1997) A study of habitat fragmentation in Southeastern Brazil using remote sensing and geographic information systems (GIS).98(1): 35–47. doi: 10.1016/S0378-1127(97)00072-8KessyJNsokkoEKaswamilaAKimaroF (2016) Analysis of drivers and agents of deforestation and forest degradation in Masito forests, Kigoma, Tanzania.6(2): 93–107. doi: 10.18488/journal.1/2016.6.2/1.2.93.107LiuX-HSkidmoreAVan OostenH (2002) Integration of classification methods for improvement of land-cover map accuracy.56(4): 257–68. doi: 10.1016/S0924-2716(02)00061-8LovettJ (1993) Eastern Arc moist forest flora. In: LovettCWasserSK (Ed.) Biogeography and ecology of the rain forests of eastern Africa., 33–57. doi: 10.1017/CBO9780511895692.004LovettJC (1996) Elevation and latitudinal changes in tree associations and diversity in the Eastern Arc Mountains of Tanzania.12: 629–650. doi: 10.1017/S0266467400009846LungTSchaabG (2006) Assessing fragmentation and disturbance of west Kenyan rainforests by means of remotely sensed time series data and landscape metrics.44(4): 491–506. doi: 10.1111/j.1365-2028.2006.00663.xLuogaEJWitkowskiEBalkwillK (2000) Differential utilization and ethnobotany of trees in Kitulanghalo forest reserve and surrounding communal lands, eastern Tanzania.54(3): 328–343. doi: 10.1007/BF02864785McGarigalK (2006) . John Wiley & Sons Ltd, Amherst. doi: 10.1002/9780470057339.val006.pub2McGarigalKCushmanSA (2002) Comparative evaluation of experimental approaches to the study of habitat fragmentation effects.12(2): 335–45. doi: 10.1890/1051-0761(2002)012[0335:CEOEAT]2.0.CO;2ManandharROdehIOAncevT (2009) Improving the accuracy of land use and land cover classification of Landsat data using post-classification enhancement.1(3): 330–44. doi: 10.3390/rs1030330MillingtonACVelez-LiendoXMBradleyAV (2003) Scale dependence in multitemporal mapping of forest fragmentation in Bolivia: implications for explaining temporal trends in landscape ecology and applications to biodiversity conservation.57(4): 289–299. doi: 10.1016/S0924-2716(02)00154-5MugashaWAEidTBollandsåsOMMalimbwiREChamshamaSAZahabuEKataniJZ (2013) Allometric models for prediction of above- and belowground biomass of trees in the miombo woodlands of Tanzania.310: 87–101. doi: 10.1016/j.foreco.2013.08.003MunishiPMringiSShirimaDLindaS (2010) The role of the Miombo woodlands of the southern highlands of Tanzania as carbon sinks.2(12): 261–269.MurciaC (1995) Edge effects in fragmented forests: implications for conservation.10(2): 58–62. doi: 10.1016/S0169-5347(00)88977-6MyersNMittermeierRAMittermeierCGda FonsecaGABKentJ (2000) Biodiversity hotspots for conservation priorities.403(6772): 853–858. doi: 10.1038/35002501NagendraHMunroeDKSouthworthJ (2004) From pattern to process: landscape fragmentation and the analysis of land use/land cover change.101(2–3): 111–115. doi: 10.1016/j.agee.2003.09.003NeelMMcGarigalKCushmanS (2004) Behavior of class-level landscape metrics across gradients of class aggregation and area.19(4): 435–55.NewmarkWD (1998) Forest Area, Fragmentation, and Loss in the Eastern Arc Mountains: Implications for the conservation of biological diversity.87(1): 29–36. doi: 10.2982/0012-8317(1998)87[29:FAFALI]2.0.CO;2OjoyiMMAntwi-AgyeiPMutangaOOdindiJAbdel-RahmanEM (2015) An analysis of ecosystem vulnerability and management interventions in the Morogoro region landscapes, Tanzania.8(3): 662–680. doi: 10.1177/194008291500800306OjoyiMMMutangaOOdindiOAbdel-RahmanEM (2016) Application of topo-edaphic factors and remotely sensed vegetation indices to enhance biomass estimation in a heterogeneous landscape in the Eastern Arc Mountains of Tanzania.31(1): 1–21. doi: 10.1080/10106049.2015.1041557OlsonDMDinersteinE (1998) The Global 200: A representation approach to conserving the Earth’s most biologically valuable ecoregions.12(3): 502–515. doi: 10.1046/j.1523-1739.1998.012003502.xPacchioniRGCarvalhoFMThompsonCEFaustinoALFNicoliniFPereiraTSSilvaRCCantaoMEGerberAVasconcelosATRAgnez-LimaLF (2014) Taxonomic and functional profiles of soil samples from Atlantic forest and Caatinga biomes in northeastern Brazil.3(3): 299–315 doi: 10.1002/mbo3.169PelkeyNWStonerCJCaroTM (2000) Vegetation in Tanzania: assessing long term trends and effects of protection using satellite imagery.94(3): 297–309. doi: 10.1016/S0006-3207(99)00195-0RosalesJ (2008) Economic growth, climate change, biodiversity loss: distributive justice for the global north and south.22(6): 1409–1417. doi: 10.1111/j.1523-1739.2008.01091.xRutledgeDT (2003) Landscape indices as measures of the effects of fragmentation: can pattern reflect process? In: O'LearyH (Ed.) New Zealand Department of Conservation, Wellington, 27 pp.SaikiaAHazarikaRSahariahD (2013) Land-use/land-cover change and fragmentation in the Nameri tiger reserve, India.113(1): 1–10. doi: 10.1080/00167223.2013.782991ShirimaDDMunishiPKLewisSLBurgessNDMarshallARBalmfordASwetnamRDZahabuEM (2011) Carbon storage, structure and composition of miombo woodlands in Tanzania’s Eastern Arc Mountains.49(3): 332–342. doi: 10.1111/j.1365-2028.2011.01269.xSouthworthJNagendraHTuckerC (2002) Fragmentation of a Landscape: Incorporating landscape metrics into satellite analyses of land-cover change.27(3): 253–269. doi: 10.1080/01426390220149511StephensSLydersenJMCollinsBMFryDLMeyerMD (2015) Historical and current landscape-scale ponderosa pine and mixed conifer forest structure in the Southern Sierra Nevada.6(5): 1–63. doi: 10.1890/ES14-00379.1SwetnamRDFisherBMbilinyiBPMunishiPKTWillcockSRickettsTMwakalilaSBalmfordABurgessNDMarshallARLewisSL (2011) Mapping socio-economic scenarios of land cover change: A GIS method to enable ecosystem service modelling.92(3): 563–574. doi: 10.1016/j.jenvman.2010.09.007TsengM-HChenS-JHwangG-HShenM-Y (2008) A genetic algorithm rule-based approach for land-cover classification.63(2): 202–212. doi: 10.1016/j.isprsjprs.2007.09.001VogelmannJE (1995) Assessment of forest fragmentation in Southern New England using remote sensing and geographic information systems technology.9(2): 439–449. doi: 10.1046/j.1523-1739.1995.9020439.xWehkampJAquinoAFussSReedE (2015) Analyzing the perception of deforestation drivers by African policy makers in light of possible REDD+ policy responses.59: 7–18. doi: 10.1016/j.forpol.2015.05.005WiensJA (1995) Habitat fragmentation: island v landscape perspectives on bird conservation.137: 97–104. doi: 10.1111/j.1474-919X.1995.tb08464.xWiensJ (2000) Ecological heterogeneity: an ontogeny of concepts and approaches. In: HutchingsMJJohnEAStewartAJA (Eds) The Ecological Consequences of Environmental Heterogeneity – 40th Symposium of the British Ecological Society, 23–25 March 1999. Blackwell Sciences, Falmer, 9–31.XiZ (2007) (MSc Dissertation). International Institute for Geoinformation Science and Earth Observation, Enschede.YandaPShishiraE (1999) Land cover changes and their Influencing factors in Tanzania. Regional land cover changes. In: RavichandranV (Ed.) Sustainable Agriculture and Their Interactions with Global Change – Proceedings of an International Workshop, 28–30 July 1997. Mozambique Universities Press (India) Limited (Maputo), 49–58.YenPZieglerSHuettmannFOnyeahialamAI (2005) Change detection of forest and habitat resources from 1973 to 2001 in Bach Ma National Park, Vietnam, using remote sensing imagery.7(1): 1–8. doi: 10.1505/ifor.7.1.1.64163