Analysis of the Impact of Land-Use/Land-Cover Change on Land-Surface Temperature in the Villages within the Luki Biosphere Reserve

: Major land-use/land-cover change due to rapid urbanization has been known to increase the land-surface temperature around the world. Consequently, examining the variation of land-surface temperatures and mitigating the related impacts remain a challenge. The present study employed remote-sensing and geoinformational techniques to examine land-use/land-cover change and its effects on land-surface temperature variations in the villages within the Luki Biosphere Reserve, Democratic Republic of Congo. Land-use/land-cover change for the year 2038 was predicted by using the CA–Markov chain. Additionally, focus-group discussions (FGDs) with local communities from different villages were applied to better understand the impact of climate change, considering the increase of land-surface temperature. The results revealed major changes in land-use/land-cover in the four villages from 2002 to 2020, principally the expansion of fallow land and built-up areas, as well as the decline in forest land, and the complex of young secondary and degraded forest. There was an increase in mean LST values over all villages between 2002 and 2020. The highest value was observed in Tsumba kituti (25.12 ◦ C), followed by Kisavu (24.87 ◦ C), Kibuya (23.31 ◦ C) and Kiobo (21.82 ◦ C). Between 2002 and 2020, the mean LST of built-up areas increased from 23.18 to 25.12 ◦ C, 21.55 to 23.38 ◦ C, 21.4 to 25.78 ◦ C and 22.31 to 25.62 ◦ C in Tsumba kituti, Kiobo, Kisavu and Kibuya, respectively. Moreover, the mean LST of fallow land increased from 20.8 to 23.2 ◦ C, 21.13 to 22.12 ◦ C, 21.89 to 23.12 ◦ C and 20.31 to 23.47 ◦ C in Tsumba, Kiobo, Kibuya and Kisavu, respectively. This indicates that built-up and fallow land experienced the highest land-surface temperature compared to other land-use/land-cover categories. Meanwhile, the conversion of all land-use/land-cover categories into built-up areas in all the villages resulted in the increase of the land-surface temperature. FGDs results recognize the recurrent land-use/land-cover change as the major driver of the increase in LST (86%). However, it was predicted that farmland and built-up area will still increase within all the villages, while the forest land will decline. As for the complex of secondary and degraded forest, it will decrease in Tsumba kituti, while, in Kiobo and Kisavu, it is expected to increase. Through a combination of remote-sensing and primary data, this study provides accurate information that will beneﬁt decision-makers to implement appropriate landscape-planning techniques to mitigate the effect of the increased land-surface temperature in the villages. savannah; FF, fallow land and fields; BA, built ‐ up area).


Introduction
In recent years, most regions across the world are experiencing major land-use/landcover changes [1][2][3]. Numerous studies have reported that the observed changes in landuse/land-cover are mainly driven by human activities including agriculture [1], and rapid urbanization [2,4]. Urbanization plays a pivotal role in the transition of land use/land cover all over the globe [5][6][7]. The increase in population demands in order to sustain their livelihoods leads to land-use/land-cover change [8,9]. Recently, numerous natural and anthropogenic factors have extensively transformed natural forest land into builtup and farm land [10][11][12]. These transformations of natural lands to man-made land use have been negatively impacting the eco-environment [13], as well as the sustainable development of local communities [14,15]. Previous studies have reported the relationship between land-use/land-cover change caused by urbanization, and the rise in land-surface temperature [10,16,17]. It is worth noting that each land-use/land-cover category has its own properties that are related to the absorption and radiation of solar energy [12].
Other studies have pointed out that changes in land use/land cover due to urbanization, deforestation, and land degradation significantly affected the urban ecosystem [10,11,18]. These artificial changes in land use/land cover can lead to the increase in LST, due to the expansion of impervious surface [19][20][21]. Globally, change in land-surface temperature is one of the major issues that affect human populations, and often results from natural disasters and human induced land-use/land-cover changes. The extreme change in land use/land cover increases the related effects of land-surface temperature. Therefore, for a convenient life in the proximate environment, the implementation of alternative and mitigation solutions is necessary. Between 2010 and 2050, it is estimated that the number of people living in cities will increase by around~80% [22]. Democratic Republic of Congo is among African countries experiencing high rate of anarchic urbanization due to a lack of land planning and policies, and major changes in land use/land cover. Furthermore, in urban areas, the increased human population results in the increase of the demand of various infrastructures to support their livelihoods, such as the construction of commercial and residential areas, public utilities, and road infrastructures which ultimately leave footprints on the environment. Recently, villages within and around the Luki Biosphere Reserve are facing an expansion due to the increase in local populations. This region is particularly faced with the migration of people from different parts of the province due to the high forest dependency and high rate of poverty. This rapid increase of population in the villages within and around the reserve has negatively impacted the natural environment due to the conversion of forest land into farmland and built-up areas, leading to the increase of impervious surfaces and warming of the villages. Therefore, accurate and updated information on the trend of the expansion of villages is necessary to allow proper implementation of strategies related to sustainable management of the natural environment. Over years, mapping land-use/land-cover change represents an essential tool for land use planning [16,23] that can benefit decision-makers to reinforce their understanding on the factors that drive land-use/land-cover change, and to predict the future trends of the environmental changes [13,24]. Geoinformational and remote-sensing technologies present the most suitable tools for monitoring and planning land use [25,26], as well as examining land-use/land-cover change from a local to a global scale [16,17,20]. There are numerous methods for monitoring land-use/land-cover change. However, for several decades, remote sensing has been the most successfully used technique. This technology has been proved to be an indispensable tool, because it covers a wide area, provides spatial information and repeated coverage of large areas and reduces the cost and time for data acquisition. Thus, remotely sensed data, especially Landsat images, have been successfully used to model land-use/land-cover change [27][28][29][30][31][32] and assess its impact on the variation of LST [17,33]. Although numerous studies have been conducted across the world, the impact of land-use/land-cover change on the land-surface-temperature variation in the villages within the Luki Biosphere Reserve is yet to be investigated. As such, the present study was motivated by the need to provide the decision-makers with valuable, up-to-date and accurate information regarding land-use/land-cover change and its impact on land-surface temperature. This study is the first of its kind in Democratic Republic of Congo's rural zone. It is aimed at examining the past trend of land-use/land-cover change on the variation of the land-surface temperature over the last two decades (2000-2020), and evaluating the variation of land-surface temperature across different land-use/land-cover categories, using novel approaches such as Random Forest and CA-Markov chain. The local communities were also engaged through focus-group discussions (FGDs) to find the impact of land-use/land-cover change, land-surface temperature, and climate change in the study area. This is important in developing and implementing strategies of sustainable management of land use, temperature increase and climate-change-mitigating strategies. The specific objectives of the present research were to: examine land-use/land-cover change over villages within the Luki Biosphere Reserve between the years 2000 and 2020; assess the impact of land-use/land-cover change on the variation of land-surface temperature; predict the future trends of land-use/land-cover change. The findings from this study will enable the decision-makers and policy-makers to better understand the past and future trends of land-use/land-cover change, and its impact on land-surface temperature. As such, this will aid the implementation of appropriate strategies that could minimize problems related to the increased land-surface temperature.

Study Area
The present study was carried out in the Luki Biosphere Reserve, in the Democratic Republic of Congo (DRC) (Figure 1). The Luki Biosphere Reserve extends over an area of about 33,000 hectares [34]. It was created in 1937 and was recognized as a biosphere reserve by the UNESCO in 1979. The reserve is located at the eastern boundary of the Mayombe forest, and remains a sample relic of the Mayombe vegetation [35]. The annual average rainfall ranges between 1150 and 1500 mm. The annual average temperature ranges between 25 and 30 • C [36]. The region is located within a humid tropical climate Aw5 according to the Köppen classification [37]. The soils in the reserve are generally ferrallitic and acidic [38]. Additionally, these soils are characterized by a low content of cations. Although Luki Reserve is conserved with the support of the DRC government institutions and the WWF, natural resources are nevertheless prone to tremendous anthropogenic pressures, mainly due to local population activities, such as slash-and-burn agriculture, wood energy, bush fires, fuelwood, illegal logging, etc. This compromises the sustainability of natural resources in the Reserve. date and accurate information regarding land-use/land-cover change and its impact on land-surface temperature. This study is the first of its kind in Democratic Republic of Congo's rural zone. It is aimed at examining the past trend of land-use/land-cover change on the variation of the land-surface temperature over the last two decades (2000-2020), and evaluating the variation of land-surface temperature across different land-use/landcover categories, using novel approaches such as Random Forest and CA-Markov chain.
The local communities were also engaged through focus-group discussions (FGDs) to find the impact of land-use/land-cover change, land-surface temperature, and climate change in the study area. This is important in developing and implementing strategies of sustainable management of land use, temperature increase and climate-change-mitigating strategies. The specific objectives of the present research were to: examine land-use/land-cover change over villages within the Luki Biosphere Reserve between the years 2000 and 2020; assess the impact of land-use/land-cover change on the variation of land-surface temperature; predict the future trends of land-use/land-cover change. The findings from this study will enable the decision-makers and policy-makers to better understand the past and future trends of land-use/land-cover change, and its impact on land-surface temperature. As such, this will aid the implementation of appropriate strategies that could minimize problems related to the increased land-surface temperature.

Study Area
The present study was carried out in the Luki Biosphere Reserve, in the Democratic Republic of Congo (DRC) (Figure 1). The Luki Biosphere Reserve extends over an area of about 33,000 hectares [34]. It was created in 1937 and was recognized as a biosphere reserve by the UNESCO in 1979. The reserve is located at the eastern boundary of the Mayombe forest, and remains a sample relic of the Mayombe vegetation [35]. The annual average rainfall ranges between 1150 and 1500 mm. The annual average temperature ranges between 25 and 30 °C [36]. The region is located within a humid tropical climate Aw5 according to the Köppen classification [37]. The soils in the reserve are generally ferrallitic and acidic [38]. Additionally, these soils are characterized by a low content of cations. Although Luki Reserve is conserved with the support of the DRC government institutions and the WWF, natural resources are nevertheless prone to tremendous anthropogenic pressures, mainly due to local population activities, such as slash-and-burn agriculture, wood energy, bush fires, fuelwood, illegal logging, etc. This compromises the sustainability of natural resources in the Reserve.

Data Collection and Preprocessing
The present research examines the effect of the land-use/land-cover change on LST in the Luki Biosphere Reserve from 1987 to 2020. Satellite images were freely downloaded from the United States Geological Survey (USGS) earth explorer website (http: //earthexplorer.usgs.gov/ (accessed on 10 July 2021)) and the Google Earth Engine website (https://code.earthengine.google.com (accessed on 10 July 2021)). Landsat images (Landsat 4, 7 and 8) were collected for the year 1987, 2002 and 2020 (Table 1). Upon performing the radiometric and atmospheric corrections [12,17], the LST [39,40] was derived from spectral bands. To derive the LST values, the digital numbers of spectral bands were first converted into radiance. The radiance values were then converted into brightness temperature values. In addition, other spectral bands were used to provide land-use/land-cover maps by means of the maximum likelihood classifier. The ENVI 5.3 (L3Harris Geospatial, Boulder, CO, USA), ArcGIS 10.1 (ESRI, New York, NY, USA), QGIS 2.18 (Open Source Geospatial, Beaverton, OR, USA) and R 3.61 (R Foundation for Statistical Computing, Vienna, Austria) softwares were used for all image processing operations.

Land-Use/Land-Cover Classification
After the conversion of radiance into reflectance, the extraction of land-use/landcover categories was carried out for the years 2002 and 2020. The supervised classification was performed by using the maximum-likelihood algorithm to categorize the pixels into different land-use/land-cover types [41]. The landscape of the study area was classified into five categories: forest land, built-up area, savannah, fallow land and fields, complex of young secondary and degraded forest.

Estimation of Land Surface Temperature
The land surface temperature was estimated using the Landsat thermal bands. The following process was used to convert Landsat thermal bands as Digital Number (DN) to land surface temperature: (a) Conversion of digital number into radiance [12] where L λ is the spectral radiance; L MAXλ is the maximum spectral radiance; L MINλ is the minimum spectral radiance; Q Calmin is the minimum quantized calibrated pixel value in DN; Q Calmax is the maximum quantized calibrated pixel value in DN; and Q Cal is the quantized calibrated pixel value in DN.

Cellular Automata
In general, the Cellular Automata models predict the real nature regulations. Numerous studies have reported that predicting land-use/land-cover change based CA technique is the frequently used method, which is capable of providing explicit spatial modelling results based on defined transition rule [44][45][46]. Additionally, Cellular Automata models are appropriate for representing, analyzing and stimulating geographic processes due to the relationships among raster grids [47,48]. The model incorporates both static and dynamic aspects of LULC transformation, and has been commonly used to simulate LULC change [45,47,49]. Recently, CA models have been widely used in the study of landuse/land-cover change [2,50]. According to Abubakr et al. [51], the CA models consider time and space as discrete units, where the space is considered as a regular two-dimensional grid. The main characteristic of this model is that it is based on local interactions reflecting the dynamics of system evolution [52,53]. Batty et al. [50] noted that, CA models have the ability to simulate stochastic, nonlinear and spatial processes. Numerous studies have illustrated the capability of Cellular Automata models to model complex spatiotemporal processes of land-use change [46,[54][55][56][57]. The CA model comprises five components, namely space, state, time step, neighborhood and transitional rule [58]. The space in a CA is represented as a grid, and it comprises cells of equal size and shape. There are transitional rules that guide the mutual interaction between cells in a grid, and they are based on the spatial relationship between cells, including distance and neighborhood. The interactions between the cells take place at each instant and the cell state update takes place as progress is made. With the presence of multiple layers in CA, the interaction of cells in each grid is based on multiple predictors.
The transition rule that defines the state of each cell in the subsequent time step depends on the current state of that cell and its surrounding neighborhood cells. Afterwards, land-use/land-cover change suitability maps are required, and the dynamics should be defined into the system. Filters are used to determine the neighbors in CA model. The magnitude of the weight factor increases with the proximity of the nuclear cell to its neighbor. Therefore, land use/cover is not an entirely random decision, and the simulation of the adjacent cells in a grid should result from the combination of the weight factor and the probabilities of transition. Two types of neighborhood are widely used in CA, namely Von Neumann and Moore [59]. The present research applied a filter with 5-by-5 Moore neighborhood to create transition rules with the aim of simulating the future change in the land-use/cover. The rules represent the suitability maps, which illustrate the probability of each pixel (cell) to change its state during the period under study.
The basic expression of CA model as illustrated in Figure 2 can be expressed as follows: where S denotes a set of finite cell states. The instant and future times are represented by t and t + 1, respectively. N illustrates the cell neighborhood, and f represents the transition rule of cellular states in local space.
Sustainability 2021, 13, 11242 6 of 25 land-use/land-cover change suitability maps are required, and the dynamics should be defined into the system. Filters are used to determine the neighbors in CA model. The magnitude of the weight factor increases with the proximity of the nuclear cell to its neighbor. Therefore, land use/cover is not an entirely random decision, and the simulation of the adjacent cells in a grid should result from the combination of the weight factor and the probabilities of transition. Two types of neighborhood are widely used in CA, namely Von Neumann and Moore [59]. The present research applied a filter with 5-by-5 Moore neighborhood to create transition rules with the aim of simulating the future change in the landuse/cover. The rules represent the suitability maps, which illustrate the probability of each pixel (cell) to change its state during the period under study.
The basic expression of CA model as illustrated in Figure 2 can be expressed as follows: where S denotes a set of finite cell states. The instant and future times are represented by t and t+1, respectively. N illustrates the cell neighborhood, and f represents the transition rule of cellular states in local space.

Markov Chain Analysis.
Markov Chain approaches have been effectively used to simulate land use/land cover. The model fundamentally evaluates the transition probability trend of landuse/land-cover categories within a certain area during a certain period. The model is made up of an ensemble of alternative states over time. This assumes that the prediction of system state at time t + 1 is based on its previous state at time t. This study considered the

Markov Chain Analysis
Markov Chain approaches have been effectively used to simulate land use/land cover. The model fundamentally evaluates the transition probability trend of land-use/land-cover categories within a certain area during a certain period. The model is made up of an ensemble of alternative states over time. This assumes that the prediction of system state at time t + 1 is based on its previous state at time t. This study considered the period between 2002 and 2020 to generate the transition probability matrices [60]. The generated matrices offer the transition probabilities for all land-use/land-cover classes over the defined time, and were used to predict the land-use/land-cover change in 2020. The mathematical expression of the Markov chain is illustrated as follows [61,62]: where S (t) and S (t + 1) are states of the system at time t and t + 1, respectively; P ij represents the transition probability matrix for the land-use/land-cover category, and transition from the state i to the state j, and from the time t to the time t + 1. The mathematical expression of P ij is presented as follows [63,64]: ||P ij || = P 11 P 12 P 13 . . . P 1n P 21 P 22 P 23 . . . P 2n P 31 P 32 P 33 . . . P 3n . . . · · · · · · . . . · · · P n1 P n2 P n3 . . . P nn 0 ≤ P ij ≤ 1) and ∑ n j=1 P ij=1, i , j=1, 2,...,n where Pn represents the state probability at any time, n is the LUC type number, S is LUC status and t or t + 1 is the time point. Thus, any class with a low transition value will be expected to have a probability close to (0), while a high transition will have a probability close to (1) [64].

Implementation of CA-Markov Chain Model
The improvement of the precision of land-use/land-cover change prediction techniques can be achieved by integrating two or more techniques, considering that each technique has its advantages [65][66][67]. Several studies have recently used the CA-Markov model to predict dynamic spatial phenomena and forecasting future land-use changes [53]. Indeed, the CA-Markov chain model benefit from the integration of Markov chain and the CA model. As such, CA-Markov models have been found to be suitable for efficient prediction of change in land use/land cover [53,[67][68][69]. Additionally, the integration of Cellular Automata and Markov chain approaches offers accurate results in land-use/landcover change prediction and takes into account the absence of statistical, socioeconomic and historical data [70][71][72]. In the CA-Markov method, the component related to temporal change of land-use/land-cover categories is guided by Markov chain process based on generated transition matrices [73]. On the other hand, the spatial change component is directed by transition potential maps, neighborhood configuration and the rules of local transition throughout the process of CA model [68,[74][75][76]. In this study, prediction of land-use/land-cover change was performed by using IDRISI-Selva 17.0 (Clark University, Worcester, MA, USA) software environment. Statistics of change in each land-use/landcover class was computed by using cross tabulation analysis. The CA-Markov model was then run to predict the future change of each land-use/land-cover category for each village, as illustrated in Figure 3. Thus, the following three major steps were followed: 1.
The computation land-use transition probability matrix and transition rules, using Markov chain analysis. The precedent state of each land-use/land-cover category was used to predict its future state. To predict land-use/land-cover change in 2038, we used the transition probabilities map between 2002 and 2020. Moreover, the land-use/land-cover transformation rules were provided by the transition probability matrices. These matrices provide the description of transition probability of each land-use/land-cover category into others. However, quantities of change in each land-use/land-cover categories to other categories were provided by the transition area matrices.

2.
Determining the CA filter. Here, several options of standard contiguity kernels were applied as neighborhoods to identify the suitable contiguity filter for the purpose of land-use/land-cover prediction. These are 7 × 7, 5 × 5 and 3 × 3 contiguity kernels.
For this study, we selected the contiguity filter 5 × 5. This filter considers that the center of each pixel is surrounded by matrix space of 5 × 5 dimension in order to reflect the change in each pixel significantly.

3.
Determining the number of iterations and the time of starting point of CA. To identify the suitable iterations, we applied several iteration numbers, starting from 1 to 300. Ultimately, we selected four iterations to carry out the land-use/land-cover change prediction.

Model Validation
The present study used the three map comparison method (a reference map of time 1, a reference map of time 2 and a predicted map of time 2) to evaluate the model accuracy. This approach was implemented by Pontius et al. [77], and is composed of two components of agreement and three components of disagreement [78]. The components of agreement includes correctly simulated persistence (correct rejection) and correctly simulated The present study used the three map comparison method (a reference map of time 1, a reference map of time 2 and a predicted map of time 2) to evaluate the model accuracy. This approach was implemented by Pontius et al. [77], and is composed of two components of agreement and three components of disagreement [78]. The components of agreement includes correctly simulated persistence (correct rejection) and correctly simulated change (hits), while the components of disagreement comprise change predicted as persistence (misses), persistence predicted as change (false alarm) and change predicted as change to wrong category (wrong hits). We then used the reference map of 2002, the reference map of 2017 and the predicted map of 2017. The validation results showed a low percentage of disagreement components of 3.60%, 3.39%, 2.49% and 1.98% for Tsumba kituti, Kisavu, Kibuya and Kiobo, respectively; meanwhile, the percentages of agreement components were high for all villages ( Table 2).  The decrease in forest area was estimated at 5.4 hectares. Consequently, this research revealed tremendous land-use/land-cover transformation that varied with space and time in the four villages. These land-use/land-cover transformations were mainly driven by rapid urbanization. Our results are consistent with previous findings [16,17,21,33,79,80]. Major land-use/land-cover transformations were the main causes various environmental issues, including the increase in land-surface temperature.   49°C). The present research has revealed highest LST values in the center of built-up zones of all villages, and it decreased with the distance from the center. This increase was due to the presence of settlements and other impervious surface with high reflectance rate. These findings were consistent with previous research [16,19,20,21,80]. Additionally our findings for the four villages denote that highest LST values for 2002 and 2020 were related to dry objects such as buildings, rocks, and bare lands. According to previous studies, the aforementioned objects have the capability of raising the efficacy of the surface to absorb solar energy and convert it into heat energy  This increase was due to the presence of settlements and other impervious surface with high reflectance rate. These findings were consistent with previous research [16,[19][20][21]80]. Additionally, our findings for the four villages denote that highest LST values for 2002 and 2020 were related to dry objects such as buildings, rocks, and bare lands. According to previous studies, the aforementioned objects have the capability of raising the efficacy of the surface to absorb solar energy and convert it into heat energy [39,81]. This results in the       Figure 8 illustrates the mean LST for different land-use/land-cover types for 2002 and 2020. Given the variation in reflectance properties related to each land-use/land-cover type, there was observed variation in the mean LST value for different land cover types The spatial structure of each land-use/land-cover type was also not homogeneous. Thus Traore [85] and Ahmed [86], observed the highest mean LST values in built-up area, in Bangui and Dhaka cities, respectively. In addition, the low increase in vegetation coverage (forest land) can be related to the absorption of solar radiation by leaves and canopy, leading to the temperature regulation of the proximate environment. Our results were in line with previous findings [87,88]. It is worth noting that, the conversion of natural landscape to artificial land-use as result of rapid urbanization increases the LST [16,[19][20][21]80].
Traore [85] and Ahmed [86], observed the highest mean LST values in built-up in Bangui and Dhaka cities, respectively. In addition, the low increase in vegetatio erage (forest land) can be related to the absorption of solar radiation by leaves and ca leading to the temperature regulation of the proximate environment. Our results w line with previous findings [87,88]. It is worth noting that, the conversion of natura scape to artificial land-use as result of rapid urbanization increases the LST [16,19-2   Table 3 illustrates LST values in different land-use/land-cover classes newly converted during the period under study (2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016)(2017)(2018)(2019)(2020). It was revealed that LST values increased in all newly converted areas. In all the villages, areas newly converted into built-up and farmland exhibited high increase in LST values compared to other land-use/land-cover classes. In Tsumba kituti, the increase in mean LST values was estimated at 3.13, 3.62 and 3.12 • C, from the conversion into built-up areas from forest lands, complex of degraded and young secondary forest, and farmland, respectively. The conversion of forest lands and complex of degraded and young secondary forest into farmland resulted in an increase in LST values by 2.43 • C and 2.2 • C, respectively. In Kisavu, transformation of forest land into built-up area resulted to an increase in mean LST by 3.25 • C, while the transformation of complex of degraded and young secondary forest into the built-up area increase in mean LST by 2.91 • C. Transformation of complex of degraded and young secondary forest into farmland, forest land into built-up areas, and forest land into farmland increased the mean LST by 2.17, 3.86 and 3.07 • C, respectively. In Kiobo, mean LST increased by 2.23 • C by transforming forest land into built-up area, 1.99 • C by transforming farmland into the built-up area, 1.1 • C by transforming forest land into farmland, and 1.18 • C by transforming complex of degraded and young secondary forest into farmland. In Kibuya, transforming farmland into built-up areas, and forest land into built-up areas increased the mean LST by 2.74 and 3.6 • C, respectively. The increase in the mean LST values in areas newly converted to built-up and farmland confirms that the increase of LST in the villages of Luki results from the increase in bare land and imperviousness of soils.

Spatiotemporal Distribution of LST and Its Variations with LULC Types
The present study was aimed at examining the land-use/land-cover change and its effects on the variation of the land-surface temperature in the villages within the Luki Biosphere Reserve for the years 2002 and 2020. The increasing trend of built-up areas and decreasing trend of forest land were revealed, and this was in line with previous findings from the region [1]. Indeed, built-up area increased by 25.3, 3.7, 2.5 and 7.5 hectares for Tsumba kituti, Kiobo, Kibuya and Kisavu, respectively. On the other hand, forest land decreased by 49.9, 27, 5.4 and 157.7 hectares for Tsumba kituti, Kiobo, Kibuya and Kisavu, respectively. Fallow land and fields increased in all the villages during the period under study, while the complex degraded and young secondary forest decreased in Tsumba kituti and Kiobo, increasing in Kisavu. These dramatic changes experienced in the four villages were induced by anthropogenic pressures, especially urbanization, agriculture and forest resources utilization [1]. In each village, large amount of forest lands were converted into built-up areas or fallow land and fields to establish settlement, houses and farmland. Additionally, the expansion of built-up and farmland was attributed to population growth in the villages within and around the Luki Biosphere Reserve [1]. In order to meet the demand of housing for the growing population, extensive built-up expansion has happened with dramatic impacts on natural resources as a result of over-exploitation. This expansion of built-up area did not occur equally in the four villages, and the high expansion rate was observed in villages located near the national road (i.e., Tsumba kituti and Kisavu). The expansion of built-up area, converted from other land-use/land-cover classes, has led to the increase of the land-surface temperature in the four villages. This could be due to the increase of impervious areas, which are known for their capacity to store solar heat during the day, resulting in the overall increase of the land-surface temperature. A previous study conducted by Pal and Ziaul [84] reported similar trend in the change of land use/land cover on the shifting of heat zone from lower to a higher rate in English Bazar urban center, west Bengal, India. Indeed, the built-up expansion has been found to play a key role in the increase of the land-surface temperature in the four villages. Based on the spatial distribution of LST, it was found that high temperature corresponded to built-up areas. The increase of LST in newly grown built-up was higher than the other land use categories.
The effect of land-use/land-cover change on land-surface temperature has been addressed on the previous sections. However, the impact of climate change was also taken into account. The present research revealed that that the minimum, mean and maximum land-surface temperature progressively increased in the four villages, from 2002 to 2020 ( Figure 6). In Tsumba kituti, the minimum LST value decreased with 0.26 • C from 2002 to 2020, while the mean and maximum LST values increased by 4.84 and 4.03 • C, respectively. In Kisavu, the minimum, mean and maximum temperatures increased by 0.04, 4.85 and 4.74 • C, respectively, from 2002 to 2020. In Kiobo, the minimum LST value decreased by 0.31 • C from 2002 to 2020, while the mean and maximum temperatures increased by 2.81 and 1.49 • C, respectively. In Kibuya, the minimum LST value decreased by 0.63 • C from 2002 to 2020, while the mean and maximum temperatures increased by 2.2 and 3.3 • C, respectively. However, the increasing trend of the land-surface temperature revealed in the present research cannot only be linked to land-use/land-cover change, but could also be linked to climate change; this was, however, beyond the scope of this study. These findings show that the increase of land-surface temperature in these villages was mainly caused by the changes in land cover and possibly climate change.

Impact of Land-Use/Land-Cover Change, Climate Change and Land-Surface Temperature in the Villages
In the present section, the results of focus discussion groups are briefly presented in order to highlight the perception of local communities on the impacts of land-use/landcover change and temperature changes in the villages over the last 20 years.

Perception of Local Communities on the Impact of Land-Use/Land-Cover Change and Climate Change in the Region
In the present research, focus-group discussions were conducted to collect the perception of local communities regarding the impact of land-use/land-cover change on temperature, soil erosion, biodiversity, forest lands, and water resources (Figure 9). The 200 respondents made their choices from a wide range of possibilities, from very low to very high, for each question related to the potential impact of land-use/land-cover change as aforementioned. The results from the discussion groups reveal that land-use/land-cover change significantly increases temperature (86% very high), transforms forest lands into settlements and agriculture land (96% very high), reduces water-body surface (85% very high) and leads to the increase in surface runoff and soil erosion (71% very high). Additionally, land-use/land-cover change causes biodiversity loss (97% very high) due to the destruction or fragmentation of habitats. Based on these results, it is possible to validate observations acquired from satellite images that revealed dramatic land-use/land-cover change in the region. The conversion of forest lands to built-up and farmland leads to the increase of land-surface temperature due to the appearance of impervious surfaces.
On the other hand, a series of questions were asked to the local communities, regarding the impact of climate change on seasonal variations, agricultural production, frequency and intensity of precipitation, biodiversity loss and human health ( Figure 10). The respondents chose among the different options offered to them, ranging from very less to very high, for each question related to the impact of climate change as abovementioned. There was significant seasonal variation (91% very high) and the change in the frequency and intensity of precipitation (92%), leading to the reduction of crop production (91% very high). This situation compromises food security and causes an increase in malnutrition and outbreaks of epidemic diseases in the region. Moreover, climate change increases the risk of some infectious diseases (89% very high) in the Luki region, including diseases that occur in hot areas and are transmitted by mosquitoes and other insects. This includes malaria, a disease that is very common in the region. Finally, the majority of respondents reported a significant reduction of biodiversity (93% very high). These results were consistent with those obtained from satellite image processing, which revealed a significant increase in temperature and the reduction in forest land and biodiversity due to climate change, and the conversion of forest lands into built-up area and farmland.  In conclusion, the summary of the findings from the present study is outlined as follows: 1.
The present study revealed a major growth of built-up area in all the four villages, which was transformed from forest land, complex degraded and young secondary forest, fallow land and fields, and savannah. 2.
The highest and lowest land-surface temperature in built-up area and forest land, respectively, were determined for all the years under study. The spatial mean LST significantly increased by 2 Changes in land use/land cover will continue to be experienced in all the villages (Figures 11 and 12). In Tsumba kituti, the complex of degraded and young sec-ondary forest, forest land and savannah will continue to decrease with 40 Consequently, this will impact the variation of land-surface temperature in the future. 5.
Land-use/land-cover change and climate-change-impacted land-surface temperature in the four villages. According to the perception of local communities, these changes can cause different environmental issues such as biodiversity loss, soil erosion, change in frequency and intensity of precipitation, seasonal variation and increase in temperature. Additionally, this change negatively impacts human health and crop production in the region. In conclusion, the summary of the findings from the present study is outlined as follows: 1. The present study revealed a major growth of built-up area in all the four villages, which was transformed from forest land, complex degraded and young secondary forest, fallow land and fields, and savannah. 2. The highest and lowest land-surface temperature in built-up area and forest land, decrease by 4.2 hectares. Consequently, this will impact the variation of land surface temperature in the future.
(5) Land use/land cover change and climate change impacted land surface temperature in the four villages. According to the perception of local communities, these changes can cause different environmental issues such as biodiversity loss, soil erosion, and change in frequency and intensity of precipitation, and seasonal variation and increase in temperature. Additionally, this change negatively impacts human health and crop production in the region.

Conclusions
The degradation of natural land-cover types has been found to influence climate change in many cities and counties. Major changes in land-use/land-cover type impacts the land-surface temperature all over the world. Thus, monitoring land-use/land-cover change is important to better understand the land-surface temperature variation. Remotesensing and geoinformational techniques were used to examine the change in landuse/land-cover and its impact in the variation of land-surface temperature. Moreover, fo-

Conclusions
The degradation of natural land-cover types has been found to influence climate change in many cities and counties. Major changes in land-use/land-cover type impacts the land-surface temperature all over the world. Thus, monitoring land-use/land-cover change is important to better understand the land-surface temperature variation. Remote-sensing and geoinformational techniques were used to examine the change in land-use/landcover and its impact in the variation of land-surface temperature. Moreover, focus-group discussions with local communities were applied to understand the perception of local communities towards the drivers of climate change in the villages. The study results revealed tremendous changes in land-use/land-cover in the four villages from 2002 to 2020, principally the expansion of fallow land and fields and built-up areas, and the decline in forest land and complex of degraded and young secondary forest. The present study found an increase in the mean LST values for all villages from 2002 to 2020. Tsumba kituti presented the highest mean LST (25.12 • C), followed by Kisavu (24. Results provided by focus-group discussions with local communities showed that land-use/land-cover change is the main cause of the increase in LST. The majority of respondents (86%) revealed that heat waves are among the most important drivers of adverse climate change in the villages. However, the prediction results found that the farmland and built-up area will still increase within all the villages, while the forest land will decline. As for the complex of secondary and degraded forest, it will decrease in Tsumba kituti, while, in Kiobo and Kisavu, it is expected to increase. Therefore, to mitigate the effects of land-use/land-cover change on land-surface temperature, decision-makers should implement sustainable management strategies.