1. Introduction
In recent years, the ecological and environmental impacts caused by land transformation have been reported to grow at a very rapid rate and are considered to be a major problem around the globe [
1,
2]. These impacts are connected with the exploitation of natural resources, destruction of biodiversity as well as disturbing the ecosystem structure that offers essential services to mankind including nutrient cycling, pollination, predator control, carbon sequestration, and soil fertility [
3]. According to [
4], Africa and Asia are the areas with the highest ecological problems caused by land use and land cover (LULC) changes. The drivers of LULC changes such as overpopulation, improper control of land resources, together with poor economic situations have made these regions highly affected by ecological problems [
5]. For instance, [
6] showed that between 1975 to 2000, more than 16% of forest cover and 5% of woodland and grassland had been destroyed on the African continent. In East Africa including Zanzibar, 48% of the tree cover was destroyed from 1986 to 2015 due to the rapid rate of deforestation and expansion of casual agricultural activities [
7]. Most of this destruction was attributed to the threat of an extremely limited number of natural resources, impairing agricultural lands, rising in irregular and unreliable urban areas as well as increasing in landscape fragility and ecological shock [
8]. In this sense, appropriate methods in assessing both ecological health together with environmental quality are essential to achieve policy making objectives through land use planning [
8,
9]. Additionally, policies required to prevent or mitigate ecological impacts of environmental changes call for robust scientific evidence to provide effective solutions to regional ecological problems and promote the useful interaction of socio–economic–ecological development [
10].
Assessing ecological risk conditions in the landscape through the analysis of LULC changes can be used to quantify the environmental effects of immediate and long-term damage or harm of certain stressors to an ecosystem at a regional scale [
11]. The concepts of ecological risk assessment (ERA) began back to the 1980s as an assessment tool to evaluate the likelihood of the potential adverse effects of one or more stress factors on ecosystems [
12,
13,
14]. These techniques can have an important role in formulating procedures to reduce negative shocks on ecosystem destruction as well as ensure proper management and precise spatial allocation of the natural resources [
15]. These can also be used in providing spatially explicit assessment in several important areas including land cover, biodiversity, carbon sequestration, and climatic regulation as well as water and soil-related ecosystems [
16,
17,
18].
ERA studies using LULC change maps have been applied in different areas around the globe. These have been developed to analyze the spatial-temporal distribution of landscape ERA in an watershed area [
19], to identify the landscape fragmentation level of green spaces as a result of built-up expansion [
20], to identify ecological networks and the most valued lands for protection [
21], among many other examples that can be found in the literature. In all these cases, remote sensing techniques have been used as a preliminary stage to acquire the land information before performing further analysis on calculating ecological indices to assess the likelihood of ecosystems to collapse on a given landscape. The ultimate goal in performing ERA is to generate relevant planning information so that climatic and other ecological damages can be reduced to a minimum scale [
22]. It seeks to use historical land information with that obtained from remote sensing to provide the scientific foundation for risk management in natural protection and biodiversity conservation, whose importance can serve in both the current and future sustainability of natural resources [
23]. Ecosystem functioning indicators such as the landscape ecological risk index (ERI) can be calculated from LULC maps and provide an estimation of the environmental effects as well as long-term damage to ecosystems of a particular region [
11]. It can also facilitate in maintaining sustainable development in natural resources and ensure the preservation of the integrity and stability of ecosystems, which is important for improving ecological security in fragile and intense ecosystem areas [
24]. Different management alternatives are needed to incorporate scenarios in a science-policy interface for environmental assessments [
25]. Thus, the use of plausible scenarios at regional scales may support sustainable development in a context of environmental changes on land systems [
26,
27,
28].
Studies for assessing ecological problems in Africa are still very rare [
5]. This is also influenced by insufficient data as well as the heterogeneity of the influencing factors. According to the review carried out by [
5], Tanzania (including Zanzibar) is among the countries with fewer ecological assessment studies in Africa having found only eight studies in a period of twelve years (2004–2016) in the Web of Science indexed journals. We updated and extended this assessment for Zanzibar for the period of 2000–2019 and found only seven articles, and none of these assessed the ecological risk conditions of Zanzibar’s landscape. This is problematic since the ecological health condition of the Zanzibar landscape is brutally affected by a combination of human and natural pressures [
29]. If not promptly checked, the complete extinction of endemic animals or plants together with intense changes in climatic conditions may occur in the near future [
29]. Therefore, the aim of this study was to assess the current trends and future projection of the ecological risk of Zanzibar as a consequence of thee LULC changes. Specifically, we aim to:
Produce LULC maps of Zanzibar Island for the years 2003, 2009, and 2018;
Identify LULC changes and project them according to scenarios; and
Estimate the ERI under different scenarios.
The outcome of this study may help in informing authorities to understand the ecological processes and land use dynamics of various land cover classes, along with preventing unmanaged growth and haphazard development of informal housing and infrastructure, which may result in a gradual deterioration of the environment and a decline in the quality of life and ecosystem services now and for future generations.
2. Study Area
Zanzibar (
Figure 1) is a small coastal island of about 2461 km² in total, located at 39°05’E to 39°55’E and 4°45’S to 6°30’S in the Western Indian Ocean, where tourism and agricultural activities are the major economic sectors contributing about 28% and 40% of its GDP, respectively [
30]. Zanzibar is part of the United Republic of Tanzania, and consists of two main islands—Unguja and Pemba—and almost 50 other small islets [
31,
32]. These two islands are located 40 km off the mainland coast of East Africa in the Indian Ocean. The two main islands are 50 km apart and are separated by the 700-meter-deep Pemba channel. However, the current study uses the term Zanzibar referring to the island of Unguja, as most of the people live on this island (Unguja) and because Zanzibar city is located on it.
According to the 2018 National Bureau of Statistics (NBS) census report, the overall population of Zanzibar Island (including Pemba) is about 1,348,776 people, and is growing at an annual rate of 3.4% [
33]. Natural population growth is considered as the main agent for this high rapid rate, but other factors including the growth of the tourism industry and economic development have also attracted a significant number of migrants from the mainland of Tanzania and other parts of East Africa, making this island very densely populated (546 hab./km
2).
Administratively, the Zanzibar Island (Ungula) is divided into three regions and six districts (Urban, West, Central, North A, North B, and Central). One-third of the population is concentrated within the small urban district covering the central parts of Zanzibar Town. Another one-third of the population is located within the West district. The remaining are distributed among the four other districts, occupying 86% of the land area. The Urban and West districts are the most developed districts, and have problems related to urbanization. Most agricultural and forest products are brought from other areas where the direct dependency on primary production is lower, though a gradual change towards rustic can already be seen in the West district. The four other districts can be characterized as more rural and primary production based, with less population pressure and lower development status. The lower incomes are compensated with domestic food and natural resource production, but still poverty and incapability are more present here than in the urban areas. Fishing is an important source of livelihoods in North A, South and Urban districts, while the people in South are less engaged to agricultural actions.
Zanzibar has a tropical climate with four distinct seasons: “Kaskazi” (the hot season, which is between December and February and associated with either little or no rains), “Masika” which is the long rainy season from March–May, “Kipupwe” (the cold season, with strong winds between June and September and “Vuli”, which is very short rainy season from October–December). The annual rainfall of the island ranges from 1600 mm for Unguja and 1900 mm for Pemba, respectively. The temperature ranges between 29 °C and 32 °C on average. Zanzibar Island was originally forested, but pressures like population increase, human habitation, and climate change have resulted in a widespread clearing of the forest and vegetation cover [
34].
The land cover patterns around Zanzibar are generally distributed under two different soil classes, namely the deep soil and the coral rag areas. The deep soil areas are mainly attributed by permanent cultivation and forest, while the coral rag is characterized by cultivation and scrublands [
31]. However, these areas are not stable in terms of land cover and biotope characteristics. The hilly deep soil areas are reported to be more vulnerable to soil erosion, high land use demand, and, simultaneously, face population pressure, which all influence the rapid and constant changes in land cover patterns and natural resources [
32,
34]. Alternatively, shifting cultivation in coral rag areas creates a constant element of change, which under the pressure of diminishing area due to commercial and conservation land use leads toward the deterioration of ecosystems and valuable natural resources [
32].
5. Discussion
The findings in this research support the evidence of a recent study that found that the rapid land use changes in Zanzibar have imposed an intensive degradation on the ecosystems of Zanzibar [
29]. The main drivers (i.e., informal housing, high rate of casual farming, and improvised planning in the tourism sector) behind this intensive degradation are increasing, uncontrolled, and poorly managed [
30]. During the first stages of analyzing land use dynamics, we have realized that Zanzibar City has experienced a trend of rapid land use changes and extensive built-up expansion, mostly in areas near existing buildings and some of the coastal zones. These changes are easily observable by visual comparison of the generated LULC maps for the different years. The distance to the buildings and distance to coast variables support this observation as they provided higher Crammer’s values, ensuring that these variables have a great influence on land cover dynamics in Zanzibar. Evidence has been clearly exposed by increasing the fragmentation index (Ci), isolation index (Ni), loss index (Ri) as well as overall ERI of the entire landscape. In order to reduce such types of threats, more systematic approaches in ecological assessment are required. Cooperative research among research groups, land planners, and all other stakeholders including the tourism sector should be highly enforced as a tool to monitor and recommend a long-term plan for sustainable ecosystems in Zanzibar.
In future scenarios, the findings have proven that disobeying land policies and regulations can negatively influence the ecological conditions of Zanzibar in the future. This has been shown by the results of the extreme scenario, which was the worst when compared to the other scenarios. The conservation scenario has shown little influence on directing future land cover changes and ecological risk control. These results were only moderate and, therefore, it is important for the government to take a special effort on setting up new land use approaches such as promoting the densification of vertical buildings rather than relying on one-floor horizontal buildings. If the current trend continues, the built-up area of Zanzibar City will be expanded more than 40% by 2027. The larger part of the agricultural and forested land has already been fragmented to housing parcels. The negative environmental effects relevant to land resources such as loss of natural forests, open cropland, and wetlands have been observed in Zanzibar [
52]. The government should enforce restrictions and impose strong land policies, leaving the remaining land cover classes to stop being further destroyed. These land use changes and extensive built-up expansion should slow or decline in the future. However, the situation is difficult to deal with because of the poor quality of housing and minimum government efforts to enforce the optimal usage of land resources on Zanzibar Island [
64]. Other important aspects such as population increase, the pressure of economic development and, particularly of the tourism sector, also need to be dealt with by the government to ensure sustainable land use policies for Zanzibar.
The scenarios used in this study were designed to show how different policies may impact the ecological risk in Zanzibar. These scenarios, designed by the authors of this study, aimed primarily at stressing differences and were supported by a coherent narrative [
65] aiming to facilitate initial discussions with stakeholders regarding land use planning. Thus, this study contributes with preliminary results so that an informed discussion with stakeholders can take place more effectively in Zanzibar or in other islands/coastal areas with similar contexts. Future versions of this study would benefit from the inclusion of stakeholders in an early stage to incorporate information that would make the scenarios more realistic (e.g., identification of other drivers not used in this study) [
66,
67]. The results obtained in this study are not easily comparable with studies of other islands since very few studies have addressed the impact of landcover changes in islands [
68]. However, the methodology described in this study, which is based on free data and tools, can easily be extended and adapted to other locations with similar characteristics.
It is important to highlight that in this study, we assessed the ERIs of Zanzibar Island with land cover change using freely available satellite and ancillary data. This process requires consistent and accurate land cover data for the purpose of providing reliable and accurate information. However, acquiring free satellite images with the desired standard and at the desired time was challenging. This is due to the geographical region of Zanzibar Island, making the most of their satellite images being covered with clouds. To have a more interesting time span for the study, we used an ETM+ image for 2018. The procedure followed [
46] for filling the data gaps of the SLC-off image, caused by a hardware error, with data from a SLC-on image was not completely successful. This may have impacted some of the results obtained for the land cover changes as well as for modeling future scenarios. However, we believe that this impact was not so important because only the eastern part of Zanzibar Island—with low levels of urbanization—was affected by this problem and the accuracy of assessment results were acceptable. In future works, these results can be improved by using more accurate and updated data such as Sentinel−2 as well as performing additional analysis such as comparing the evolution of the calculated indices inside and outside the protected areas to understand how effective these environmental protection instruments are [
69].