Next Article in Journal
How Does Capital Endowment Impact Farmers’ Green Production Behavior? Perspectives on Ecological Cognition and Environmental Regulation
Next Article in Special Issue
Land Use/Cover Change Prediction Based on a New Hybrid Logistic-Multicriteria Evaluation-Cellular Automata-Markov Model Taking Hefei, China as an Example
Previous Article in Journal
The Effects of Urban Sprawl on Electricity Consumption: Empirical Evidence from 283 Prefecture-Level Cities in China
Previous Article in Special Issue
U-Net-STN: A Novel End-to-End Lake Boundary Prediction Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Assessment of Landscape and Land Use/Cover Change and Its Implications for Sustainable Landscape Management in the Chittagong Hill Tracts, Bangladesh

1
Research Center of Forest Management Engineering of National Forestry and Grassland Administration, Beijing Forestry University, Beijing 100083, China
2
Sino-France Joint Laboratory for Invasive Forest Pests in Eurasia, Department of Forest Protection, School of Forestry, Beijing Forestry University, Beijing 100083, China
3
Department of Forestry Sciences, University of Zalingei, Zalingei P.O. Box 6, Sudan
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Land 2023, 12(8), 1610; https://doi.org/10.3390/land12081610
Submission received: 16 July 2023 / Revised: 10 August 2023 / Accepted: 12 August 2023 / Published: 15 August 2023
(This article belongs to the Special Issue Assessment of Land Use/Cover Change Using Geospatial Technology)

Abstract

:
Human-caused environmental change has profoundly impacted resource management and land use patterns in Bangladesh’s Chittagong Hill Tracts. This study used multi-temporal Landsat images from 1998, 2008, and 2018 to analyze land use and land cover changes, particularly those associated with forest cover changes, in Bangladesh’s Chittagong Hill Tracts. Using object-based image classification, Landsat images from 1998, 2008, and 2018 were separated into four categories based on their dominant land use and land cover features: forest, grassland, water bodies, and bare land. Post-classification comparison was used to assess the degree and frequency of change, and this method was further developed to evaluate the balance, fluctuation, and adaptation of forests. In addition, the spatial structure of land cover and temporal trajectories related to changes in forest cover were studied. The CA–Markov chain model was also used to anticipate the 2048 LULC map. The image classification of the years 1998, 2008, and 2018 showed that the overall accuracy was 89.65%, 84.44%, and 86.26%; producer accuracy was 90.00%, 68.75%, and 72.22%; and the Kappa coefficient was 85.68, 82.84, and 76.36, respectively. The results showed that between 1998 and 2018, forest cover increased by 58.03%, transforming grassland to forest; grassland increased by 29.50%, converting bare land to grassland; and forest conversion to grassland was 13.34%. In addition, the result of the landscape metric revealed that during the whole study period, class level indicated a fragmentation of forest, bare land, grassland, and water in the CHT, and landscape level indicated by Shannon’s Diversity Index and Shannon’s Evenness Index showed a slight decrease in the land. Based on the CA–Markov model, forest area is predicted to expand to 9129 Km2 in 2048; however, other land uses (bare land and grassland) continue to decrease. This substantial increase in forest cover results from effective forest management based on community forestry practices and the successful execution of Bangladesh’s national forest strategy. However, as Bangladesh’s population rises, so does the country’s need for lumber/timber. Bangladesh’s government should revise its forest policy to meet the local community’s needs without endangering the forest, and policymakers must take climate change seriously. Our strategy for evaluating the critical indicators of changes in forest cover and pathways of change will aid in connecting these patterns to the dynamics of change, such as deforestation and reforestation. It would therefore serve as a framework for developing effective conservation and management plans for the Chittagong Hill regions in Bangladesh.

1. Introduction

The recognition of the invaluable services provided by tropical forests to both host nations and the global community is widely acknowledged [1]. These unique ecosystems, characterized by their unparalleled array of plant and animal species, possess irreplaceable biodiversity and genetic resources [2,3]. Since the Earth Summit, forest policies in numerous countries have prioritized the objective of sustainable forest management, regardless of the extent of human interventions within forested areas [4]. Forests and woodlands are vital to keeping the environment carbon-free [5,6,7]. However, ecological systems, such as forests, experience constant changes due to natural biological processes, resulting in continued instability [4].
Bangladesh has approximately 2.52 million hectares of tropical forest, accounting for 10% of its area [8]. Rural areas house over two-thirds of Bangladesh’s population, and their livelihoods are closely linked to the forests, either directly or indirectly [9]. For several decades, the national forests in Bangladesh have faced significant and rapid depletion, reaching a critical stage of concern [1]. The natural forests of Bangladesh experienced a consistent annual decline of 2.1% over 20 years until the early 1980s, which further accelerated to a rate of 2.7% between 1984 and 1990 [10]. Between 1964 and 1985, there was a decline in the growing stock within Chittagong Hill Tracts reserve forests, decreasing from 23.8 million m3 to below 19.8 million m3 [11]. Between 2000 and 2005, around 2000 hectares of forest cover were lost annually [12].
As a result of degradation, forest communities are motivated to utilize their traditional knowledge and practices to engage in activities such as conservation, reforestation, bushfire control, and the prevention of illegal forest exploitation and encroachment [13]. The combined efforts of these local communities have led to the regeneration of forested lands and the enhancement of biodiversity levels [14]. Several policies and practices in South and Southeast Asia have emerged centered around community-based forest management [15,16]. The global community places significant emphasis on preserving biodiversity, promoting forest health, ensuring sufficient forest productivity, and safeguarding the socio-economic functions associated with forest resources [17].
In Bangladesh, traditional forest management techniques have historically aimed to achieve economic and ecological objectives [18]. However, rapid deforestation occurred due to various socio-economic and socio-political factors [19,20], diminishing the effectiveness of traditional forest planning and management approaches. Unplanned human activities and unforeseen pressures exceeded planned conservation efforts, resulting in extensive deforestation and fragmentation of forest resources [19]. Given the country’s dense population and limited land area, policymakers had to explore alternative management practices. In the late 1970s, social forestry was introduced as a successful alternative that transitioned the Forest Department’s role from a custodial to a participatory model, involving local communities in forest protection, reforestation activities, and benefit-sharing arrangements [21].
During Bangladesh’s current period of sovereignty, the Forest Act underwent its first amendment in 1989, which aimed to enhance forest protection by imposing stricter penalties and limiting the discretionary powers of forest officials and local magistrates [18]. Although this amendment primarily focused on strengthening traditional forest protection measures, it did not introduce the concept of social forestry until 2000, when another amendment was introduced, leading to the emergence of social forestry in Bangladesh [22]. The Forest (Amendment) Act of 2000 marked a significant milestone as it paved the way for the formulation of the groundbreaking 2004 Social Forestry Rules (SFR) by the government [18]. Bangladesh’s fundamental principle of social forestry revolves around integrating local communities in reforestation activities, aiming to achieve multiple objectives encompassing ecological, economic, and social benefits [23].
One of the main requirements for global change research is to evaluate and track the condition of the earth’s surface [24,25]. As the foundation for all living things and a key factor in global climate change, vegetation classification and mapping are crucial technological undertakings for managing natural resources [25,26]. Land use/cover change (LUCC) is most frequently associated with logging, globalization, and agricultural expansion that alter the natural vegetation [27,28]. At both the local and global levels, LUCC causes several environmental issues, such as biodiversity loss brought on by greenhouse gas emissions [28,29], variations in land surface temperature (LST), and shifts in precipitation [30]. Urbanization’s adverse environmental effects, which include population increase, extensive infrastructure development, and constantly shifting landscapes, are a worldwide issue [31].
Bangladesh, known for its high population density, faces increasing land pressure, particularly in forested areas, due to food production, urban settlements, and industrial development [1]. Approximately 60% of the total land is utilized for agriculture, which serves as a primary source of livelihood for over two-thirds of the rural population [32]. However, limited land availability per person, constrained by geographical factors and inadequate farming practices, hinders the country’s food production capacity [9]. The forests and agricultural lands play a crucial role in the lives of people residing in Bangladesh’s Chittagong Hill Tracts (CHT) region [33]. In the past, these forested landscapes provided various local and regional benefits, including food, energy, timber, water, and healthcare, while also contributing to national revenue generation. Nonetheless, the exploitation and degradation of forests, which started in the previous century [34] and continue to this day, have significant implications for the sustainable livelihoods of forest-dependent communities in terms of both direct and indirect resources [35].
Bangladesh is not exceptional in these environmental changes. Urbanization has also affected the local environment to certain degrees and made it more susceptible to land degradation. Thus, timely and accurate information about local spatial coverage, distributions of LULC categories, and their dynamics are prerequisites for the country’s planning, socio-economic development, and sustainable land management. No evidence-based studies have been conducted in the Chittagong region to understand how implementing a community-based forestry policy affects forest cover change and regional land dynamics. This is the first study of its kind to examine the impact of community-based forestry on land use and land cover change in the Chittagong Hill Tracks and whether or not there are factors that influence this change. This study will highlight the aspects responsible for a land cover change in the Chittagong Hill Tracks. The study aimed to access the dynamic LULC change detection in the Chittagong Hill Tracks using remote sensing data and Geographic Information Systems technologies for 1998, 2008, and 2018. The study also used the cellular automata–Markov model (CA–Markov) to predict future land use changes under a simulated 2048 scenario. The main objectives of this study were: (i) to detect and identify the dynamics change in LULC from 1998 to 2018 in the Chittagong Hill Tracts, Bangladesh; (ii) to identify the driving forces behind these changes; and (iii) to predict the LULC map for the year 2048 using CA–Markov.

2. Materials and Methods

2.1. Study Area

The Chittagong Hill Tracts are part of the Chittagong Division, located in the southeastern hilly area of Bangladesh, and play a vital role (Figure 1) [9]. It geographically lies between 21.025′ N to 23.045′ N latitude and 91.045′ E to 92.050′ E longitude and covers an area of approximately 13,198.16 Km2, almost 40% of the total forest land area, and 8% of the total land area in Bangladesh [36]. The study area borders Myanmar in the southeast, the Indian state of Tripura to the north, Mizoram to the east, and the Chittagong district in the west. The Chittagong Hill Tracts are divided into three hill districts: Khagrachhari, Bandarban, and Rangamati. The area comprises forest, grassland, watershed, bare land, and agricultural land [9]. Paths and rows cover this area: 136/44, 136/45, and 135/45. In this region, almost twelve ethnic groups are the principal inhabitants, highly dependent on natural resources for their livelihood [37]. The major sources of income are agriculture, livestock, and the harvesting of forest products. The landscape and local people are unique parts of the Chittagong Hill Tracts. The traditional landscape of Chittagong Hill Tracts combines a mix of land uses: natural forests and plantation forests cover more than 70% of the land area in three districts. The primary plant species in the region are tropical wet evergreen/semi-evergreen and deciduous, classified as ‘hill forests’. The large area covers mixed natural and planted forests [34]. This area is characterized by a tropical monsoon climate with three dominant seasons: the dry winter season (November to March), the pre-monsoon season (April to May), and the monsoon season (June to October). The rainfall usually starts to increase from March to June. Annual rainfall ranges from 2540 mm in the north to 2450 mm in the east, reaching nearly 3810 mm in the south and west.

2.2. Remote Sensing Data

The data used in this study were obtained from various Landsat sensors, including Thematic Mapper (TM), Enhanced Thematic Mapper Plus (ETM+), and Landsat 8 OLI for 1998, 2008, and 2018 (Table 1). A total of 3 senses were acquired between November and December. All images were obtained from the United States Geological Survey (USGS—https://earthexplorer.usgs.gov/ (accessed on 12 January 2023)). All raw satellite images had a cloud cover of less than 10%. Each image cloud was removed by supervised classification extraction by the mask tool in ENVI 5.1. For the Thematic Mapper (TM), the Enhanced Thematic Mapper Plus (ETM+) image band combination was red-green-blue (RGB) [4-3-2], and the Landsat 8 OLI image band combination was red-green-blue (RGB) [5-3-2].

2.3. Satellite Image Pre-Processing

Geometric rectification is essential for producing corrected LULC maps (Figure 2). For change detection, various pre-processing requirements, for example, geometric correction, radiometric collaboration, and atmospheric corrections, are the most important to avoid errors in the results. In this study, Landsat Collection Level-1 downloaded images were rectified and corrected—geometrically and topographically. Moreover, one of the more significant preconditions for remote sensing data analysis is an atmospheric correction, which has been done by applying the Model FLAASH (Fast Line-of-Sight Atmospheric Analysis of Spectral Hypercubes) to improve image information by transforming radiance as a sensor into surface reflectance values. Mosaic was applied for the paths and rows: 136/44, 136/45, and 135/45. The image pre-processing was done using ENVI 5.1 and ArcGIS 10.5 software.

2.4. Image Classification

Image classification was performed using a supervised classification method based on the Maximum Likelihood Classifier of given classes (Table 2). Training samples were selected with sufficient homogeneity to symbolize the spectral and spatial characteristics of each LULC class and maximize classification accuracy. This step is the most significant element of supervised classification because spectral signatures extracted from training samples will define the overall accuracy of the classification and, therefore, the final LULC map.

2.5. Accuracy Assessment

Accuracy assessment is a post-classification step accomplished through the correspondence evaluation of classified LULC maps against assumed factual geographical reference data [38]. The image classification accuracy was assessed with ground truth data in the form of reference data points obtained from Google Earth images of 1998, 2008, and 2018, used to obtain overall accuracy, user accuracy, producer accuracy, and Kappa coefficient. The classifier, carried out using 58, 90, and 211 points, was generated randomly for 1998, 2008, and 2018 supervised images, respectively. Each point has a specific pixel value based on ground-truth data and visual interpretation. Identified points were considered classified values, and pixel values were considered reference values.
Classification results were carried out on thematic maps using different truth reference data, and accuracy assessments were presented for all classifications. Overall accuracy, user accuracy, producer’s accuracy, and Kappa coefficient accuracy derived from the error matrices were used for the accuracy assessment of the final maps. An error matrix and Kappa statistics were generated from reference and classified data.
In this study, satellite images were accurately assessed with the Landsat images Landsat 5-TM 1988, Landsat 7-ETM+ 2008, and Landsat 8 OLI 2018, which the ground truth data likely equates to. Overall accuracy was calculated by dividing the sum of correctly classified pixels (diagonal) [39] by the total number of reference pixels.
Overall accuracy = (∑_(e = 1)^Lncc)/N
The overall classification accuracy was calculated for all classifications, as well as class-specific accuracy, which can be created on the user and producer levels. Users’ accuracy, i.e., the error of commission [40] and user accuracy assessment, was calculated by dividing the number of correctly classified pixels in each category by the total number of classified pixels in that category (row total).
Users’ accuracy = ncc/(n + c)
Producer accuracy is the measure of the error of omission [40], and producer accuracy was calculated by dividing the number of correctly classified pixels in each category by the total number of reference pixels in that category (total column).
Producer accuracy = ncc/(nc+)
The Kappa coefficient was calculated by using the following equations:
k^ = (∑_(1-1)^r〖P_ii -∑_(1-1)^r〖P_(i+) × P_(+i)〗〗)/(1-∑_(1-1)^rP_(i+) × P_(+i))
where r = the number of rows in the error matrix;
  • P_ii = the proportion of pixels in row ‘r’ and column ‘i’;
  • P_(i+) = the proportion of the marginal total of row ‘i’;
  • P_(+i) = the proportion of the marginal total of column ‘i’.

2.6. Cellular Automata Model

The cellular automata model mainly consists of cell, time, and rule neighbor. The filter of the cellular model determines the neighbors. Cellular automata (CA) are a spatiotemporal calculation of a dynamic process model used for LULC change. The commonly used neighborhoods are Moore, the extended Moore, and von Neumann. The expression of the Cellular Automata equation is:
S (t, t + 1) = f (St, N)
where S is the set of states of the finite cells, the (t, t + 1) are different times, N is the neighborhood of cells, and f is the transformation rule of local space [41].

2.7. CA–Markov Model

The CA–Markov model is valuable for modeling LULC changes and can simulate prediction changes. The spatial prediction accuracy can be effectively simulated simultaneously [42]. CA–Markov model effectively combines the Markov model and the CA model. This approach is based on a Markov stochastic probability matrix for predicting the transition from one condition to another [30]. The Markov chain model is most commonly used to simulate transitions, parameters, and trends. It created probability transition matrices to anticipate and categorize probable land use/cover change (LUCC) and urban development scenarios and investigate land change simulation trends [30]. Equations ((6)–(8)) based on the conditional probability formula were used to estimate trends:
t + 1 = P i j × S T ,
P i j = P i j = P 11 P 12 P 1 n P 21 P 22 P 2 n P n 1 P n 2 P n 3 ,
( 0 P i j < 1   a n d   j = 1 N P i j = 1 ,   i , j = 1 , 2 , n ,
where S(t) is the state of the system at time t, S (t + 1) is the state of the system at time (t + 1), and Pij is the matrix of the transition probability in a state. By projecting 2018 and 2048, the cellular automata (CA) and Markov chain models are employed to generate the LUCC future scenario. The land use change modeler (LCM) in TerrSet (Clark Labs TerrSet 18.31) was used to forecast LUCC for the projected time using the CA–Markov model [30].

2.8. Landscape Metrics

The FRAGSTATES software version 4.2.1 [43] was applied for computed landscape metrics in this study area. Landscape metrics characterized the spatial fragmentation and heterogeneity for 1998, 2008, and 2018, covering three decades. The selected landscape patterns were divided into three levels, i.e., patch, class, and landscape. FRAGSTATS is a spatial pattern analysis program implemented by the decision maker, forest manager, and ecologists to analyze landscape fragmentation or describe characteristics of the landscape and elements of those landscapes [44]. FRAGSTATS version 4.2.1 was used to extract the landscape metrics from each 1998, 2008, and 2018 classified map. In this study, a total of twelve landscape metrics were examined, including:
-
Class level: Class Area (CA), Total Landscape Area (TLA), Number of Patches (NP), Mean Patch Size (MPS), Patch Size Coefficients of Varian (PSCOV), Mean Shape Index (MSI), Edge Density (ED), Mean Nearest Neighbor Distance (MNN), and Interspersion Juxtaposition Index (IJI);
-
Landscape level: Total Landscape Area (TLA), Number of Patches (NP), Mean Patch Size (MPS), Mean Shape Index (MSI), Mean Nearest Neighbor Distance (MNN), Interspersion Juxtaposition Index (IJI), Mean Proximity Index (MPI), Shannon’s Diversity Index (SDI), and Shannon’s Evenness Index (SHEI).

2.9. LULC Change-Driving Forces Model

Previous studies have suggested that the main LULC of the Chittagong Hill Tracts area driving factor is climate change, while anthropogenic contributed to LULC slightly. Therefore, this study only considered climate change as a leading driving factor.
The monthly meteorological data from 1998 to 2018 were collected from the meteorological station of Bangladesh. Each year of the CHT area was obtained by substituting time into the linear regression model. The annual average temperature (AAT), annual average relative humidity (AARH), and annual average rainfall (AAR) were used in the analysis. Climate change factors such as the average yearly temperature and annual average relative humidity were used to establish corresponding data sets. All the meteorological data processing steps were carried out in linear regression.
Climate change factors such as the annual average temperature and annual average relative humidity were used to establish corresponding data sets. The relationship between dependent variables and independent variables is given below:
Y = a + b 1 x 1 + b 2 x 2 + b 3 x 3 + . + b n x n ,
where Y is a dependent variable, a is a constant turn, x1, x2, and x3 are independent variables, and b1, b2, and b3 are the coefficients of independent variables.
The constant coefficient illustrates the value of Y in all dependent variables at a zero time zone. In addition, the parameter coefficients state a change in Y for one unit increase in the dependent variable [45]. T-tests were applied to verify the statistical regression coefficients’ significance and constant in the linear regression model. The determination sample coefficient (R2) was applied to explain the contribution to the dependent variable in the linear regression model. R2 means how closely dependent variables are related to independent variables [46].

3. Results

3.1. LULC Classes and Their Distribution

Land use cover changes (LULCCs) were computed for 1998, 2008, and 2018, focusing on vegetation, grassland, water bodies, and the bare land area in the study area (Figure 3). In 1998, the maximum land area covered by grassland accounted for 44.71% (5901.11 Km2), which gradually decreased to 43.90% (5794.51 Km2) in 2008 and 24.10% (3180.34 Km2) in 2018. Whereas the forest area was calculated to be approximately (5576.26 Km2) in 1998, 5048.22 Km2 (2008), and 8284.62 Km2 (2018), with a slight decrease of −4% from 1998 to 2008 and periodic increments of 24.52% from 2008 to 2018. The bare land area was calculated to be 1188.75 Km2 in 1998, 1853.79 Km2 in 2008, and 1161.72 Km2 in 2018, with an increase of 5.04% from 1998 to 2008 and then a decrease of −5.25% from 2008 to 2018 (Table 3). The cumulative change calculated in water area was approximately 39.44 Km2 (1998–2018), 532.04 Km2 in 1998, and 571.48 Km2 in 2018. An increase of 0.3% was observed in water from 1998 to 2018 (Table 3, Figure 4).
Kappa coefficients and user accuracy for supervised classification (LULCC maps) were calculated with TerrSet IDRISI. The overall accuracy of the classification was observed at 89.65%, 84.44%, and 86.26%, while producer accuracies were 90.00%, 68.75%, and 72.22%, and the Kappa coefficient was 85.68, 82.84, and 76.36 in three different periods (Table 4), respectively.

3.2. LULC Change Detection

The post-classification comparison of LULC change for each class within the study period from 1998 to 2018 is shown in Figure 5, whereas the result analysis of LULC change is indicated in Table 5. From 1998 to 2008, only bare land increased by 665.04 Km2. In contrast, water, forest, and grassland decreased significantly by −30.40 Km2, −528.04 Km2, and −106.60 Km2, respectively, with a significant decrease in bare land and grassland, i.e., −692.07 Km2 and −2614.17 Km2, between 2008 and 2018 (Table 5). During the three decades, the grassland decreased by −2720.77 Km2, while the forest increased by an area of about 2708.36 Km2.
Regarding LULC conversions and transformations, the most dynamic change in the study area is the conversion of grassland to forest by 13.34% (Table 6), followed by forest into grassland by nearly 58.03%. The transformation of bare land to forest was 1.68%, and the transformation of bare land into grassland was calculated to be 6.28%. However, forest land conversion to bare land in the region was estimated to be 12.10%, while the transformation of grassland into bare land was observed to be 29.50%. Therefore, between 1998 and 2018, the region experienced many changes in its LULC.

3.3. Prediction of LULC Change Based on the Markov Model

The state transition area map was created according to LULC maps from 1998 to 2018, which can be used to predict, using the CA–Markov model in IDIRISI software version 17.0, the land requirements for the different LULC types in 2048. The predictive results map for 2018 is obtained with a 5 × 5 contiguity filter, whose running cycle is 30 years. The combination of cellular automata (CA) and the stochastic transition matrix of the Markov chain model resulted in LUCC for the projected period of 2048 (Figure 6). Map accuracy for the projected land use/cover change for predictive years was classified by a sufficient Kappa coefficient value of 0.97. A 2048 map predicted that the maximum area covered by forest accounts for 9129 Km2 (69.17%) and the minimum area covered by water accounts for 665 Km2 (5.05%). A decrease of 1.51% of bare land and 5.61% of grassland cover areas were estimated during 2018–2048, respectively (Table 7). However, forest and water area cover will expand by 6.40% and 0.72%, respectively, during 2018–2048.

3.4. Landscape Metrics Analysis of Land Use and Land Cover Structure

From LULC maps from 1998 to 2018, the most changing classes, such as water, bare land, forest, and grassland, were chosen to analyze spatial landscape metrics at class and landscape levels. The most significant change among land use and land cover classes is increased forest and decreased grassland. The statistical results of the landscape metrics in the Chittagong Hill Tracts area are shown below (Table 8).
The statistics of water showed that Class Area (CA) indices increased from 51,075.18 to 57,432.96 ha, the Number of Patches (NP) also increased from 217 to 256, the Patch Size Coefficient of Varian (PSCV) increased from 831.16 to 925.63 ha, and Mean Patch Size (MPS) decreased from 235.37 to 224.35 ha during the whole period from 1998 to 2018.
Regarding bare land area during the period 1998 to 2018, CA indices decreased from 119,352.9 to 116,615.1 ha, NP increased from 2030 to 2431, Edge Density (ED) increased from 5.5 to 6.3 m/ha, Interspersion Juxtaposition Index (IJI) increased from 60.38 to 79.38%, while the PSCV decreased from 669.55 to 473.06 ha, and MPS decreased from 58.79 to 47.97 ha.
The spatial analysis of forest areas showed that CA indices increased from 557,963.6 to 828,093.2 ha, NP decreased from 1639 to 744, ED increased from 13.7 to 16.09 m/ha, IJI increased from 36.37 to 59.96%, PSCOV increased from 2388.84 to 2548.7 ha, and MPS decreased from 58.79 to 47.97 ha.
Grassland areas showed that CA indices decreased from 591,653.8 to 317,904.2 ha, NP increased from 1704 to 3268, ED decreased from 17.07 to 15.17 m/ha, IJI decreased from 63.85 to 52.48%, PSCOV decreased from 2723.98 to 1057.64 ha, and MPS decreased hugely from 347.21 to 97.28 ha.
The analysis of the landscape level showed that the fragmentation of the landscape increased with the number of patches (NP) from 5590 to 6699, the Mean Proximity Index (MPI) increased from 2296.4 to 2676.39 m/ha, the Mean Patch Size (MPS) decreased from 236.14 to 197.05, Mean Shape Index (MSI) values are identical, Mean Nearest Neighbor Distance (MNND) decreased from 606.1 to 567.5 m, Interspersion Juxtaposition Index (IJI) increased from 56.95 to 62.83%, Shannon’s Diversity Index (SEI) decreased from 0.77 to 0.71, and Shannon’s Evenness Index (SDI) decreased from 1.07 to 0.99 (Table 9). There is no significant change in the heterogeneity of the landscape.

3.5. Analysis of the Driving Forces behind LULC Change

The linear regression model result indicated that all driving factors had significant values in adjusted R square for the model of change in the forest, bare land, grassland, and water of 0.619, 0.559, 0.718, and 0.752, respectively (Table 10). The study’s results explained that AT was highly significant and AR substantially impacted forests (p < 0.01, p < 0.05; Table 8 and Table 11). AR was highly influential, and ARH significantly impacted bare land (p < 0.01, p < 0.05; Table 8 and Table 11). AT had a highly significant impact on grassland (p < 0.01; Table 8 and Table 11). AT had a highly substantial effect, and ARH, AR had a considerable effect on water (p < 0.01, p < 0.05; Table 8 and Table 11) and the change of CHT in Bangladesh.

4. Discussion

4.1. Accuracy Assessment

After calibration, land cover data were extracted from photographs and corrected for atmospheric effects. Landsat 5 TM 1998 had six spectral bands plotted upon them, but Landsat 7 ETM+ 2008 and Landsat 8 OLI 2018 had seven bands, respectively. In three different periods, a summary of supervised classification accuracy was recorded as the highest in 1998 (89.65%) and 2008 (84.44%). This result is consistent with the study findings by the authors [44]. The study indicated that all classes of producers were classified over 85% accurately. Similarly, a study by Dewan and Yamaguchi [47] using MSS found that the lowest overall accuracy was 85.6%. Another study by Kayiranga et al. [48] revealed that from 1986 to 2015, all images’ overall accuracy and Kappa coefficient results were greater than 75%. In this study, among the three periods of 1998, 2008, and 2018, Kappa coefficients were above 0.8% in all classes except 2018. Very similar results were found in the study conducted by del Castillo [49], who investigated that all forest classes’ Kappa coefficients were above 0.8 except Q. pyrenaica in 2010.

4.2. Trends in Forest Cover Changes

According to estimates of forest cover change (loss and gain) derived from Landsat images during the past twenty years, forest and tree cover increased from 1998 to 2018. The national assessment led by FAO [50] showed that the net growth of tree cover has remained stable since 1990. Different spatial and temporal extents may, in part, explain these discrepancies. Furthermore, while our results show a net gain in forest cover over time, details on the visibility and psychological impact [51] on forest loss vs. gain may also shape perspectives on forest dynamics [52].
The most notable LULC changes were the transformation of grassland to forest. The analysis of the Landsat images reveals some trends that complement the LULC from 1998 to 2018. From 1998 to 2018, the change in LULC classes shows that the Chittagong Hill Tracts region has been impacted and degraded with forest (2708.36 Km2) and bare land (−27.03 Km2), respectively. A similar result was contrary to the range reported elsewhere [53,54]. A similar result was proposed by Baral et al. [55], who reported that forest cover increased with a similar loss in the extent of all other land cover types. The human population was the most important factor in increasing or decreasing land use and land cover change.
Here is the evidence that variations in forest cover changes may be associated with the types of management regimes [56]. Community-based management or secure tenure can result in positive social interactions, better forest conditions, and opportunities to increase tree cover and distribute economic benefits [57]. State-led plantation initiatives and the National Forest Policy of 1994 have been primarily responsible for expanding forests in the CHT region, focusing on plantations made possible by the relocation of native peoples [58]. Bangladesh has demonstrated encouraging progress in curbing the annual deforestation rate. It has reduced deforestation from 2.1% from 1960 to 1980 [59] to approximately 0.2% between 1990 and 2010 [18].

4.3. Future LULC Change Prediction

In contrast, between 2018 and 2048, the predicted LULC in Chittagong Hill Tracts grassland significantly decreased after 1998, the forest increased, and the remaining area increased by a certain amount in 2048. Similar predictive results were proposed by Liping et al. [41]; they reported that water and bare land increased and the woodland area decreased significantly. Predictive results were revealed by [60]: bare land will increase and potentially lead to the loss of ecosystem services. The Markov predictive results have shown that the probability of each class changing in LULC in the future is high.

4.4. Landscape Analysis

During the study period (1998–2018) in Chittagong Hill Tracts, changes occurred in four major land classes: water, bare land, forest, and grassland. The results of the landscape indices (forest) revealed that the values of NP decreased proportionally with increased CA, ED, IJI, PSCOV, and MPS. Similar results were found in the study conducted by Kayiranga et al. [46], who reported that during the whole study period (1986–2015), the values of MPS increased with a decrease of NP inside the Nyungwe–Kibira Park. del Castillo et al. [49] also reported that the forest inside the Moncayo Natural Park (Spain) was slightly fragmented. This study claims landscape pattern metrics at heterogeneity with a similar identic landscape configuration. According to research correlations on landscape metrics, Shannon’s Evenness Index (SEI) and Shannon’s Diversity Index (SDI) were also selected as indices for the landscape pattern analysis [61]. During the study period, Shannon’s Diversity Index and Shannon’s Evenness Index showed a slightly decreased landscape. The landscape metrics analysis provided valuable information about land use and land cover change, especially fragmentation and heterogeneity weighed at the landscape class and landscape level. The landscape structure became more fragmented and heterogeneous [62]. Moreover, during the study period, landscape fragmentation increased and forests became a major land area [63,64]. Landscape (fragmentation and heterogeneity) and land use/landcover results indicated that forest management significantly increases, decreases, and maintains landscape patterns [49].

4.5. LULC Change–Climate Factors

The linear regression results assumed that the climate change factors model variables indicated that three impact factors were selected for the forest: bare land, grassland, and water. These three factors are less than significantly related to forest, bare land, grassland, and water. According to Salman-Mahini and Kamyab [65], there was a linear relationship between dependent and independent variables in multiple linear regression. The result of linear regression indicated that the climate change factors (average rainfall, average humidity, and average temperature) model as independent variables estimated the influence of drivers behind changes in four kinds of LULC change: forest, bare land, grassland, and water. Grassland was most affected by rainfall, the forest was significantly affected by rainfall and temperature, bare land was significantly affected by rainfall and humidity, and water was significantly affected by rainfall, humidity, and temperature. Rainfall was the main factor driving LULC changes. A similar result was found in a Northern Iranian study by Jahanifar et al. [46]. The linear alternative of income per capita, rain, and temperature with a definite coefficient of 0.4 as independent variables qualified for estimating forest area reduction.
Although research into the spatial-temporal trends and connections between forest cover and local land use practices will aid in identifying regions where changes have taken place and making predictions about future changes, this research alone will not be sufficient to explain the neighborhood’s driving factors of change, which are necessary for making trained, evidence-based decisions. Understanding the mechanisms of forest cover changes, like deforestation and forest regeneration, in the Chittagong Hill Tracts at various geographical scales requires knowledge of the local driving variables at work. Therefore, future research should investigate the causes and consequences of land use and forest cover changes within this landscape’s socioeconomic, policy, institutional, and past histories and how these shifts relate to larger shifts happening at the regional and global levels.

5. Conclusions

Knowledge of past land use trends and land cover change is vital to comprehending the connection between landscape dynamics and ecological responses. Our method, which integrates multi-temporal remote sensing data with GIS techniques, has enabled us to quantify and characterize the spatial-temporal pattern of LULC changes, notably those associated with forest cover changes in mountainous regions. The findings revealed that the terrain of the Chittagong Hill Tracts underwent LULC shifts between 1998 and 2018. Between those years, forest cover and water body area expanded while grassland and bare land shrank. The spatial pattern shift demonstrated that between 1998 and 2008 and between 2008 and 2018, the forest cover gained and lost area with varying annual intensities and dynamics.
The annual increase was significant during both periods, while deforestation was prominent during the first period (1998–2008) but inactive during the second period (2008–2018). Between 1998 and 2018, there was a net change of 2708.36 Km2 in forest cover and a total change of 8284.62 Km2. The LULC structure, predicted based on the CA–Markov model that in 2048, the forest area would increase drastically with a consistent decrease in grassland and bare land areas. The linear regression results for the climate change factors model variables indicated that grassland was most affected by rainfall, the forest was significantly affected by rainfall and temperature, bare land was significantly affected by rainfall and humidity, and water was significantly affected by rainfall, humidity, and temperature. Rainfall was the main factor driving LULC changes. Currently, the area faces different environmental issues that threaten these resources, such as climate change, LULC change, disturbance of species diversity, ecosystem fragmentation, and flooding. It was found that massive population growth in settlements is an important factor influencing LULC changes and the implications of sustainable landscape management. Furthermore, the findings of this study have some important suggestions that a scientific forest management system should have a strict process. Additionally, efficient programs should be run to educate local communities about sustainable land management.

Author Contributions

Conceptualization, M.C., U.H. and J.M.; methodology, M.C. and U.H.; software, M.C., U.H. and M.A.H.; validation, M.C., U.H. and J.M.; formal analysis, M.C. and U.H.; investigation, M.C.; resources, U.H. and J.M.; data curation, M.A.H.; writing—original draft preparation, M.C. and U.H.; writing—review and editing, U.H. and J.M.; visualization, U.H.; supervision, J.M.; project administration, J.M.; funding acquisition, J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (32271871).

Data Availability Statement

Data will be available on demand.

Acknowledgments

We would like to thank all the lab members of the Research Center of Forest Management Engineering of the National Grassland Administration, Beijing Forestry University, for their technical and moral support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Salam, M.A.; Noguchi, T.; Koike, M. The causes of forest cover loss in the hill forests in Bangladesh. GeoJournal 1999, 47, 539–549. [Google Scholar] [CrossRef]
  2. WCMC. Global Biodiversity: Status of The Earth’S Living Resources; WCMC: Cambridge, UK, 1992. [Google Scholar]
  3. Bawa, K.S.; Dayanandan, S. Socioeconomic factors and tropical deforestation. Nature 1997, 386, 562–563. [Google Scholar] [CrossRef]
  4. Biswas, S.R.; Choudhury, J.K. Forests and forest management practices in Bangladesh: The question of sustainability. Int. For. Rev. 2007, 9, 627–640. [Google Scholar] [CrossRef]
  5. Khan, I.; Hayat, U.; Mujahid, A.; Majid, A.; Chaudhary, A.; Badshah, M.T.; Huang, J. Evaluation of growing stock, biomass and soil carbons and their association with a diameter: A case study from planted chir pine (Pinus roxburghii) forest. Appl. Ecol. Environ. Res. 2021, 19, 1457–1472. [Google Scholar] [CrossRef]
  6. Hayat, U.; Khan, I.; Roman, M. A response of tropical tree Terminalia arjuna (Arjun) [Combretaceae-Myrtales] seeds and seedlings towards the presowing seed treatment in forest nursery of Sargodha, Pakistan. Int. J. Biosci. 2020, 16, 495–505. [Google Scholar] [CrossRef]
  7. Farooqi, T.J.A.; Hayat, U.; Roman, M.; Abbas, H.; Hussain, S. Comparative study determining the impacts of broadleaved and Needle leaved forest harvesting on hydrology and water yield: State of knowledge and research outlook. Int. J. Biosci. 2020, 16, 231–240. [Google Scholar] [CrossRef]
  8. BFD. Bangladesh Forest Department. Forest Cover of Bangladesh, 2015. Available online: http://fd.portal.gov.bd/ (accessed on 28 April 2023).
  9. Ahammad, R.; Stacey, N. Forest and agrarian change in the Chittagong Hill Tracts region of Bangladesh. In Agrarian Change in Tropical Landscapes; Center for International Forestry Research (CIFOR): Bogor, Indonesia, 2016; p. 190. [Google Scholar]
  10. FMP. Forest Production. Forestry Master Plan 1992. ADB TA 1355-BAN UNDP/FAO BGD/88/025.
  11. FMP. Main plan—1993/2012. Forestry Master Plan 1993, Vol. 1. ADB TA 1355-BAN UNDP/FAO BGD/88/025.
  12. FAO. Food and Agriculture Organization of the United Nations. Global Forest Resources Assessment Country Report: Bangladesh; FAO: Rome, Italy, 2010. [Google Scholar]
  13. CIFOR. Community-Based Forest Management Key to Success of REDD+, Center for International Forestry Research (CIFOR). 2011. Available online: https://www.cifor.cgiar.org/mediamultimedia/newsroom/pressreleases/press-releases-detail-view/article/238/community-basedforest-management-key-to-success-of-redd.html (accessed on 13 May 2023).
  14. Panigrahi, R. Democratization of forest governance: Myths and realities. In Proceedings of the 11th Biennial Conference of the International Association for the Study of Common Property, Bali, Indonesia, 19–23 June 2006; Available online: https://hdl.handle.net/10535/2167 (accessed on 1 May 2023).
  15. Balooni, K.; Inoue, M. Decentralized Forest management in South and Southeast Asia. J. For. 2007, 105, 414–420. [Google Scholar] [CrossRef]
  16. Pulhin, J.M.; Inoue, M.; Enters, T. Three decades of community-based forest management in the Philippines: Emerging lessons for sustainable and equitable forest management. Int. For. Rev. 2007, 9, 865–883. [Google Scholar] [CrossRef]
  17. Muhammed, N.; Koike, M.; Haque, F. Forest policy and sustainable forest management in Bangladesh: An analysis from national and international perspectives. New For. 2008, 36, 201–216. [Google Scholar] [CrossRef] [Green Version]
  18. Jashimuddin, M.; Inoue, M. Community forestry for sustainable forest management: Experiences from Bangladesh and policy recommendations. Formath 2012, 11, 133–166. [Google Scholar] [CrossRef] [Green Version]
  19. FAO. Forest Resources of Bangladesh: Country Report; Working Paper 15; The Forest Resources Assessment Programme; Forestry Department, Food and Agriculture Organization of the United Nations: Rome, Italy, 2000. [Google Scholar]
  20. Muhammed, N.; Koike, M.; Sajjaduzzaman, M.D.; Sophanarith, K. Reckoning social forestry in Bangladesh: Policy and plan versus implementation. Forestry 2005, 78, 373–383. [Google Scholar] [CrossRef] [Green Version]
  21. BFD, Homepage, Bangladesh Forest Department. 2011. Available online: https://www.bforest.gov.bd (accessed on 15 March 2023).
  22. Alam, M. Evolution of forest policies in Bangladesh: A critical analysis. Int. J. Soc. For. 2009, 2, 149–166. [Google Scholar]
  23. Ahmed, F.U.; Akhtaruzzaman, A.F.M. A Vision for Agricultural Research in Bangladesh: Sub-Sector—Forestry 2030, Bangladesh Agricultural Research Council (BARC), Dhaka, Bangladesh, 2010. Available online: http://www.barc.gov.bd/documents/Final-Mr.%20Farid.pdf (accessed on 10 April 2023).
  24. Jung, M.; Henkel, K.; Herold, M.; Churkina, G. Exploiting synergies of global land cover products for carbon cycle modeling. Remote Sens. Environ. 2006, 101, 534–553. [Google Scholar] [CrossRef]
  25. Xie, Y.; Sha, Z.; Yu, M. Remote sensing imagery in vegetation mapping: A review. J. Plant Ecol. 2008, 1, 9–23. [Google Scholar] [CrossRef]
  26. Xiao, X.; Zhang, Q.; Braswell, B.; Urbanski, S.; Boles, S.; Wofsy, S.; Moore, B., III; Ojima, D. Modeling gross primary production of temperate deciduous broadleaf forest using satellite images and climate data. Remote Sens. Environ. 2004, 91, 256–270. [Google Scholar] [CrossRef]
  27. Foley, J.A.; DeFries, R.; Asner, G.P.; Barford, C.; Bonan, G.; Carpenter, S.R.; Chapin, F.S.; Coe, M.T.; Daily, G.C.; Gibbs, H.K. Global consequences of land use. Science 2005, 309, 570–574. [Google Scholar] [CrossRef] [Green Version]
  28. Beckline, M.; Yujun, S.; Yvette, B.; John, A.B.; Mor-Achankap, B.; Saeed, S.; Richard, T.; Wose, J.; Paul, C. Perspectives of remote sensing and GIS applications in tropical forest management. Am. J. Agric. For. 2017, 5, 33–39. [Google Scholar] [CrossRef] [Green Version]
  29. Hersperger, A.M.; Gennaio, M.P.; Verburg, P.H.; Bürgi, M. Linking land change with driving forces and actors: Four conceptual models. Ecol. Soc. 2010, 15, 1–17. Available online: https://www.jstor.org/stable/26268195 (accessed on 18 April 2023). [CrossRef] [Green Version]
  30. Siddique, M.A.; Wang, Y.; Xu, N.; Ullah, N.; Zeng, P. The Spatiotemporal Implications of Urbanization for Urban Heat Islands in Beijing: A Predictive Approach Based on CA–Markov Modeling (2004–2050). Remote Sens. 2021, 13, 4697. [Google Scholar] [CrossRef]
  31. Huang, Y.; Qiu, Q.; Sheng, Y.; Min, X.; Cao, Y. Exploring the Relationship between Urbanization and the Eco-Environment: A Case Study of Beijing. Sustainability 2019, 11, 6298. [Google Scholar] [CrossRef] [Green Version]
  32. Ghose, B.; Razib, B.; Sharmistha, G. Reviewing the status of agricultural production in Bangladesh from a food security perspective. Russ. J. Agric. Socio-Econ. Sci. 2014, 25, 19–27. [Google Scholar]
  33. Miah, M.D.; Chakma, S.; Koike, M.; Muhammed, N. Contribution of forests to the livelihood of the Chakma community in the Chittagong Hill Tracts of Bangladesh. J. For. Res. 2012, 17, 449–457. [Google Scholar] [CrossRef]
  34. Rasul, G. Political ecology of the degradation of forest commons in the Chittagong Hill Tracts of Bangladesh. Environ. Conserv. 2007, 34, 153–163. [Google Scholar] [CrossRef] [Green Version]
  35. Rasul, G. Ecosystem services and agricultural land-use practices: A case study of the Chittagong Hill Tracts of Bangladesh. Sustain. Sci. Pract. Policy 2009, 5, 15–27. [Google Scholar] [CrossRef]
  36. BFD. District Wise Forest Land; Bangladesh Forest Department: Dhaka, Bengal, 2017. Available online: http://fd.portal.gov.bd/ (accessed on 30 April 2023).
  37. [UNDP] United Nations Development Programme. Socio-Economic Baseline Survey of Chittagong Hill Tracts; Chittagong Hill Tracts Development Facility (CHTDF): Dhaka, Bengal, 2009. [Google Scholar]
  38. Congalton, R.G. A review of assessing the accuracy of classifications of remotely sensed data. Remote Sens. Environ. 1991, 37, 35–46. [Google Scholar] [CrossRef]
  39. Congalton, R.G.; Green, K. Assessing of Remotely Sensed Data: Principles and Practices; Lewis Publishers: Boca Rotan, FL, USA, 1999; pp. 34–46. [Google Scholar]
  40. Story, M.; Congalton, R.G. Accuracy assessment: A user’s perspective. Photogramm. Eng. Remote Sens. 1986, 52, 397–399. [Google Scholar]
  41. Liping, C.; Yujun, S.; Saeed, S. Monitoring and predicting land use and land cover changes using remote sensing and GIS techniques—A case study of a hilly area, Jiangle, China. PLoS ONE 2018, 13, e0200493. [Google Scholar] [CrossRef]
  42. Liao, F.; Ds, Z. Forestry landscape patterns changes and dynamic simulation of Nanling National Nature Reserve, Guangdong. Sci. Geogr. Sin. 2014, 34, 1099–1107. [Google Scholar]
  43. McGarigal, K.; Marks, B.J. FRAGSTATS: Spatial Pattern Analysis Program for Quantifying Landscape Structure; US Department of Agriculture, Forest Service, Pacific Northwest Research Station: Washington, DC, USA, 1995; Volume 351. [Google Scholar]
  44. Keleş, S.; Sivrikaya, F.; Çakir, G.; Köse, S. Urbanization and forest cover change in regional directorate of Trabzon forestry from 1975 to 2000 using landsat data. Environ. Monit. Assess. 2008, 140, 1–14. [Google Scholar] [CrossRef]
  45. Estman, J.R. IDRISI Selva Tutorial; Clark Labs, Clark University: Worcester, MA, USA, 2012; p. 354. [Google Scholar]
  46. Jahanifar, K.; Amirnejad, H.; Mojaverian, M.; Azadi, H. Land change Detection and Effective Factors on Forest Land Use Changes: Application of Land Change Modeler and Multiple Linear Regression. J. Appl. Sci. Environ. Manag. 2018, 22, 1269–1275. [Google Scholar] [CrossRef]
  47. Dewan, A.M.; Yamaguchi, Y. Land use and land cover change in Greater Dhaka, Bangladesh: Using remote sensing to promote sustainable urbanization. Appl. Geogr. 2009, 29, 390–401. [Google Scholar] [CrossRef]
  48. Kayiranga, A.; Kurban, A.; Ndayisaba, F.; Nahayo, L.; Karamage, F.; Ablekim, A.; Li, H.; Ilniyaz, O. Monitoring Forest cover change and fragmentation using remote sensing and landscape metrics in Nyungwe-Kibira park. J. Geosci. Environ. Prot. 2016, 4, 13–33. [Google Scholar] [CrossRef] [Green Version]
  49. del Castillo, E.M.; García-Martin, A.; Aladrén, L.A.L.; de Luis, M. Evaluation of forest cover change using remote sensing techniques and landscape metrics in Moncayo Natural Park (Spain). Appl. Geogr. 2015, 62, 247–255. [Google Scholar] [CrossRef]
  50. FAO. Global Forest Resources Assessment: Country Reports Bangladesh; Food and Agriculture Organisation of the United Nations: Rome, Italy, 2015. [Google Scholar]
  51. Alfonso, A.; Zorondo-Rodríguez, F.; Simonetti, J. Perceived changes in environmental degradation and loss of ecosystem services, and their implications in human wellbeing. Int. J. Sustain. Dev. World Ecol. 2016, 24, 1–14. [Google Scholar]
  52. Twongyirwe, R.; Bithell, M.; Richards, K.; Rees, W. Do livelihood typologies influence local perceptions of forest cover change? Evidence from a tropical forested and non-forested rural landscape in western Uganda. J. Rural. Stud. 2017, 50, 12–29. [Google Scholar] [CrossRef] [Green Version]
  53. Mengistu, D.A.; Salami, A.T. Application of remote sensing and GIS inland use/land cover mapping and change detection in a part of south western Nigeria. Afr. J. Environ. Sci. Technol. 2007, 1, 99–109. [Google Scholar]
  54. Diallo, Y.; Hu, G.; Wen, X. Applications of remote sensing in land use/land cover change detection in Puer and Simao Counties, Yunnan Province. J. Am. Sci. 2009, 5, 157–166. [Google Scholar]
  55. Baral, P.; Wen, Y.; Urriola, N.N. Forest cover changes and trajectories in a typical middle mountain watershed of western Nepal. Land 2018, 7, 72. [Google Scholar] [CrossRef] [Green Version]
  56. Ahammad, R.; Stacey, N.; Eddy, I.M.S.; Tomscha, S.A. Sunderland TCH. Recent trends of forest cover change and ecosystem services in eastern upland region of Bangladesh. Sci. Total Environ. 2018, 647, 379–389. [Google Scholar] [CrossRef]
  57. Robinson, B.; Holland, M.; Naughton-Treves, L. Does secure land tenure save forests? A meta-analysis of the relationship between land tenure and tropical deforestation. Glob. Environ. Chang. 2014, 29, 281–293. [Google Scholar] [CrossRef]
  58. Nath, T.K.; Inoue, M. The upland settlement project of Bangladesh as a means of reducing land degradation and improving rural livelihoods. Small-Scale For. 2008, 7, 163–182. [Google Scholar] [CrossRef]
  59. Chowdhury, J.K. National forest policy review: Bangladesh. In An Overview of Forest Policies in Asia; Enters, T., Qiang, M., Leslie, R.N., Eds.; EC-FAO Partnership Programme (2000–2002); European Commission and Food and Agriculture Organization: Bangkok, Thailand, 2003; Available online: http://www.fao.org/DOCREP/005/AC921E/AC921E00.HTM (accessed on 14 May 2023).
  60. Padmanaban, R.; Bhowmik, A.K.; Cabral, P.; Zamyatin, A.; Almegdadi, O.; Wang, S. Modelling urban sprawl using remotely sensed data: A case study of Chennai city, Tamilnadu. Entropy 2017, 19, 163. [Google Scholar] [CrossRef] [Green Version]
  61. Riitters, K.H.; O’neill, R.V.; Hunsaker, C.T.; Wickham, J.D.; Yankee, D.H.; Timmins, S.P.; Jones, K.B.; Jackson, B.L. A factor analysis of landscape pattern and structure metrics. Landsc. Ecol. 1995, 10, 23–39. [Google Scholar] [CrossRef]
  62. Linh, N.H.K.; Erasmi, S.; Kappas, M. Quantifying land use/cover change and landscape fragmentation in Danang City, Vietnam: 1979–2009. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2012, 39, B8. [Google Scholar]
  63. Baskent, E.M.İ.N.; Kadiogullari, A.L.İ. Spatial and temporal dynamics of land use pattern in Turkey: A case study in Inegol. Landsc. Urban Plan. 2007, 81, 316–327. [Google Scholar] [CrossRef]
  64. Teixido, A.L.; Quintanilla, L.G.; Carreno, F.; Gutierrez, D. Impacts of changes in land use and fragmentation patterns on Atlantic coastal forests in northern Spain. J. Environ. Manag. 2010, 91, 879–886. [Google Scholar] [CrossRef]
  65. Salman Mahini, A.; Kamyab, H. Applied Remote Sensing and GIS with Idrisi; Mehrmahdis: Tehran, Iran, 2012; p. 596. [Google Scholar]
Figure 1. Location map of Chittagong Hill Tracts showing the study area in Bangladesh.
Figure 1. Location map of Chittagong Hill Tracts showing the study area in Bangladesh.
Land 12 01610 g001
Figure 2. Flow chart of the remote sensing methodology.
Figure 2. Flow chart of the remote sensing methodology.
Land 12 01610 g002
Figure 3. LULC thematic maps of Chittagong Hill Tracts: (a) 1998, (b) 2008, and (c) 2018.
Figure 3. LULC thematic maps of Chittagong Hill Tracts: (a) 1998, (b) 2008, and (c) 2018.
Land 12 01610 g003
Figure 4. Spatial area coverage % of LULC in Chittagong Hill Tracts during the study interval (1998–2018).
Figure 4. Spatial area coverage % of LULC in Chittagong Hill Tracts during the study interval (1998–2018).
Land 12 01610 g004
Figure 5. LULC changes in Chittagong Hill Tracts from 1998 to 2018.
Figure 5. LULC changes in Chittagong Hill Tracts from 1998 to 2018.
Land 12 01610 g005
Figure 6. Predicted LULC classification maps of Chittagong Hill Tracts (a) predicted map of 2018 and (b) predicted map of 2048.
Figure 6. Predicted LULC classification maps of Chittagong Hill Tracts (a) predicted map of 2018 and (b) predicted map of 2048.
Land 12 01610 g006
Table 1. Characteristics of Landsat data used in the study.
Table 1. Characteristics of Landsat data used in the study.
YearsSatelliteSensor TypeSpatial ResolutionDate of AcquisitionImage
1998Landsat 5TM30 mDecember 1998GeoTiff
2008Landsat 7ETM+30 mNovember 2008GeoTiff
2018Landsat 8OLI30 mNovember 2018GeoTiff
Table 2. Classification scheme/categories.
Table 2. Classification scheme/categories.
CodeLULC ClassesDescription
1Water Area covered with ponds, lakes, rivers, and seasonal streams.
2Bare land Areas described with non-vegetative cover, rocks, and some limited settlement areas.
3Forest Areas with natural or artificial woody vegetation and canopy cover more than <10%.
4Grassland Areas covered by grasses, agriculture, and some sparse shrubs in between.
Table 3. Depicting the percentage of the spatial coverage of LULC in Chittagong Hill Tracts from 1998 to 2018.
Table 3. Depicting the percentage of the spatial coverage of LULC in Chittagong Hill Tracts from 1998 to 2018.
YearLand Use Classes
WaterForestGrasslandBare landTotal
Area
Km2%Km2%Km2%Km2%Km2%
1998532.044.035576.2642.255901.1144.711188.759.0113,198.16100
2008501.643.805048.2238.255794.5143.901853.7914.0513,198.16100
2018571.484.338284.6262.773180.3424.101161.728.8013,198.16100
Table 4. Accuracy computation of land use and land cover change (LULC) maps between 1990 and 2018.
Table 4. Accuracy computation of land use and land cover change (LULC) maps between 1990 and 2018.
YearLULC ClassesWaterBare LandForestGrasslandKappa ValueOverall Accuracy
1998User’s81.8277.7895.0094.440.85789.65
Producer’s90.0077.7895.0089.47
2008User’s84.6164.2891.1786.210.82884.44
Producers68.7569.2391.1892.59
2018User’s86.6776.1992.0675.510.76486.26
Producer’s72.2266.6692.1086.04
Table 5. Analysis of LULC change in Chittagong Hill Tracts.
Table 5. Analysis of LULC change in Chittagong Hill Tracts.
ClassesNet of LULC Changed (Km2)
1998–20082008–20181998–2018
Water−30.4069.8439.44
Bare land665.04−692.07−27.03
Forest−528.043236.402708.36
Grassland−106.60−2614.17−2720.77
Table 6. Change matrices of LULC classes in Chittagong Hill Tracts between 1998 and 2018.
Table 6. Change matrices of LULC classes in Chittagong Hill Tracts between 1998 and 2018.
Year1998
2018LULC Classes
WaterBare landForestGrassland
Area%%%%
Water95.721.260.730.51
Bare land2.7357.591.686.28
Forest1.3612.1084.2558.03
Grassland0.1929.5013.3435.19
Table 7. Probability area percentage of LULC classes in Chittagong Hill Tracts in 2048.
Table 7. Probability area percentage of LULC classes in Chittagong Hill Tracts in 2048.
LULC ClassesPeriod
2048Change Detection
2018–2048
Km2%Km2%
Water6655.0593.520.72
Bare land9637.29−198.72−1.51
Forest912969.17844.386.40
Grassland244118.49−739.34−5.61
Total area13,198100
Table 8. Metrics of landscape patterns for 1998, 2008, and 2018.
Table 8. Metrics of landscape patterns for 1998, 2008, and 2018.
Indices
LULC ClassesYearCA (ha)NPED (m/ha)IJI (%)PSCOV (ha)MPS (ha)
Water199858,223.882171.292.32831.16235.37
200858,223.887392.1693.48991.1178.79
201857,432.962561.3395.07925.63224.35
Bare land1998119,352.920305.560.38669.5558.79
2008178,306.826718.5685.47665.7866.76
2018116,615.124316.379.38473.0647.97
Forest1998557,963.6163913.736.372388.84340.43
2008505,352.3204913.257.752969.51246.63
2018828,093.274416.0959.962548.71113.03
Grassland1998591,653.8170417.0763.852723.98347.21
2008578,162.5149215.773.153197.04387.51
2018317,904.2326815.1752.481057.6497.28
CA = Class Area (ha); NP = Number of Patches; ED = Edge Density (m/ha); IJI = Interspersion Juxtaposition Index (%); PSCV = Patch Size Coefficient of Varian (ha); MPS = Mean Patch Size (ha).
Table 9. Analysis matrices of landscape patterns with references to the studied years 1998, 2008, and 2018.
Table 9. Analysis matrices of landscape patterns with references to the studied years 1998, 2008, and 2018.
Year Indices
NPMPIMPSMSIMNNDIJISEISDI
199855902296.4236.141.27606.156.950.771.07
200869512231.18189.911.23583.372.680.821.14
201866992676.39197.051.27567.562.830.710.99
NP = Number of Patches; MPI = Mean Proximity Index; MPS = Mean Patch Size (ha); MSI = Mean Shape Index; MNND = Mean Nearest Neighbor Distance; IJI = Interspersion Juxtaposition Index (%); SEI = Shannon’s Diversity Index; SDI = Shannon’s Evenness Index.
Table 10. Climate change impact driving model variables included in the CHT.
Table 10. Climate change impact driving model variables included in the CHT.
ModelR2Std ErrorF ValueSig.
Y10.62611.379.120.008 **
Y20.56138.887.350.015 *
Y30.72462.7313.750.002 **
Y40.7510.4911.100.004 **
Note: Y1 = forest land, Y2 = bare land, Y3 = grassland, and Y4 = water. ** = highly significant; * = significant.
Table 11. Estimating the influence of drivers behind changes in LULC.
Table 11. Estimating the influence of drivers behind changes in LULC.
Model Std ErrorCoefficientsT ValueSig.
Intercept3883.059−6113.87−1.57450.154 NA
Y1AR0.644799−1.51007−2.341930.047 *
AT159.0128658.40764.1405960.003 **
Intercept1046.292−2152.03−2.056810.073 NA
Y2AR0.1446350.5103663.5286590.007 **
ARH13.3150932.691152.4551950.039 *
Intercept2938.97418,078.936.1514430.000 **
Y3AR0.4880291.0268912.1041580.068 NA
AT120.3521−629.188−5.227890.000 **
Intercept156.4302−159.272−1.018170.342 NA
Y4AR0.0111290.0263542.3679970.049 *
ARH1.3089583.4918472.6676540.032 *
AT3.55111616.202954.5627770.002 **
Note: Y1 = forest land, Y2 = bare land, Y3 = grassland, and Y4 = water. AR = average rainfall, AT = average temperature, and ARH = average relative humidity. ** = highly significant; * = significant; NA = non-significant.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chakma, M.; Hayat, U.; Meng, J.; Hassan, M.A. An Assessment of Landscape and Land Use/Cover Change and Its Implications for Sustainable Landscape Management in the Chittagong Hill Tracts, Bangladesh. Land 2023, 12, 1610. https://doi.org/10.3390/land12081610

AMA Style

Chakma M, Hayat U, Meng J, Hassan MA. An Assessment of Landscape and Land Use/Cover Change and Its Implications for Sustainable Landscape Management in the Chittagong Hill Tracts, Bangladesh. Land. 2023; 12(8):1610. https://doi.org/10.3390/land12081610

Chicago/Turabian Style

Chakma, Masheli, Umer Hayat, Jinghui Meng, and Mohammed A Hassan. 2023. "An Assessment of Landscape and Land Use/Cover Change and Its Implications for Sustainable Landscape Management in the Chittagong Hill Tracts, Bangladesh" Land 12, no. 8: 1610. https://doi.org/10.3390/land12081610

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop