Next Article in Journal
Microplastic Pollution in the Surface Waters from Plain and Mountainous Lakes in Siberia, Russia
Previous Article in Journal
Understanding the Mechanical Biases of Tipping-Bucket Rain Gauges: A Semi-Analytical Calibration Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Monitoring Land Use Changes and Their Future Prospects Using GIS and ANN-CA for Perak River Basin, Malaysia

Department of Civil and Environmental Engineering, Universiti Teknologi PETRONAS, Seri Iskandar 32610, Perak, Malaysia
*
Authors to whom correspondence should be addressed.
Water 2021, 13(16), 2286; https://doi.org/10.3390/w13162286
Submission received: 12 July 2021 / Revised: 9 August 2021 / Accepted: 12 August 2021 / Published: 21 August 2021
(This article belongs to the Section Hydrology)

Abstract

:
Natural landscapes have changed significantly through anthropogenic activities, particularly in areas that are severely impacted by climate change and population expansion, such as countries in Southeast Asia. It is essential for sustainable development, particularly efficient water management practices, to know about the impact of land use and land cover (LULC) changes. Geographic information systems (GIS) and remote sensing were used for monitoring land use changes, whereas artificial neural network cellular automata (ANN-CA) modeling using quantum geographic information systems (QGIS) was performed for prediction of LULC changes. This study investigated the changes in LULC in the Perak River basin for the years 2000, 2010, and 2020. The study also provides predictions of future changes for the years 2030, 2040, and 2050. Landsat satellite images were utilized to monitor the land use changes. For the classification of Landsat images, maximum-likelihood supervised classification was implemented. The broad classification defines four main classes in the study area, including (i) waterbodies, (ii) agricultural lands, (iii) barren and urban lands, and (iv) dense forests. The outcomes revealed a considerable reduction in dense forests from the year 2000 to 2020, whereas a substantial increase in barren lands (up to 547.39 km2) had occurred by the year 2020, while urban land use has seen a rapid rise. The kappa coefficient was used to assess the validity of classified images, with an overall kappa coefficient of 0.86, 0.88, and 0.91 for the years 2000, 2010, and 2020, respectively. In addition, ANN-CA simulation results predicted that barren and urban lands will expand in the future at the expense of other classes in the years 2030, 2040, and 2050. However, a considerable decrease will occur in the area of dense forests in the simulated years. The study successfully presents LULC changes and future predictions highlighting significant pattern of land use change in the Perak River basin. This information could be helpful for land use administration and future planning in the region.

1. Introduction

The process of determining changes in any process or object through analysis at different time periods is known as change detection [1]. Over time, humans have brought significant changes to the earth’s surface to produce food through agricultural techniques. Nearly one-third of the surface of the earth is estimated to be agricultural, and over half of the earth’s surface has been altered during the past few years [2]. This transition from naturally arising farming land to agricultural land is still underway [3]. These significant changes have drawn the attention of land use administrators and researchers to the influence of land use changes on hydrological processes [4]. Land use managers and decision makers can better understand the interactions between human and natural activities by examining the trends in change detection.
The assessment of LULC changes at various levels, such as global and regional, includes the evaluation of changes within river basins [5]. The tremendous growth in the rate of population is the dominant factor at the global scale in the transition of land use, according to [6]. The drastic variations in land cover, especially in emerging countries, are mainly due to comprehensive urban development and the transformation of natural areas into industrial or agricultural lands [7]. The depletion of natural land, thick forests, and watersheds puts considerable pressure on the hydrological regimes and mechanisms of river basins [8]. It is important to provide multi-temporal sets of data for the evaluation of changes in the spatial characteristics of land [9]. The use of multi-temporal datasets makes it simple to explain the major LULC changes along with their pattern [10]. The advancement of computer technologies and the introduction of Landsat satellites have made it simpler to trace the changes and advancements that have occurred over the last several decades. Remote sensing technology linked with a GIS has been usefully implementated in the identification of multiple environmental characteristics, i.e., vegetation covering, urban sprawls, transition in forests, and particularly variations in LULC changes for certain time periods [11]. It has been observed that remote sensing and GIS techniques offer more accurate and cost-efficient data evaluation in comparison to other traditional approaches and surveys [7]. Remote sensing is defined as observing spatial variations in different objects without physical interaction with them. To classify the various characteristics of the earth, remote sensing uses space-borne satellites. Because of their frequent monitoring of the earth’s features, they have the ability to assist in tracking changes on the land surface [12]. The information gained on temporal and spatial levels allows scientists and researchers to recognize the large-scale changes in pattern of land and permits regional policymakers and authorities to make future decisions. The application of remote sensing in natural disaster management has been reported in many studies. It had been used for monitoring floods in India [13]. Geomorphic factors were examined for the establishment of a decision support system by [14] in order to avoid the occurrence of landslides using remote sensing. Furthermore, it has proved useful in identifying changes in agricultural patterns, land cover changes, and urban sprawls [15]. GIS are regarded as essential for identifying the changes in remotely sensed images [16]. It offers the ability to integrate data from multiple resources for the identification of change. The combined effect of the hydrological maps, soil and topographic maps, and classified images obtained through GIS may provide helpful information in extracting land use for a given area. In addition, it may show the tendencies of land use changes because of its ability to develop the model by using provided statistics and datasets. Moreover, GIS and remotely sensed images are widely utilized to recognize LULC changes [17,18,19]. For tracking land cover changes, the combined use of GIS and remote sensing has been demonstrated to be a reliable and profitable technique [20,21,22,23]. The efficacy of space-borne imaging in mapping LULC changes has been demonstrated in many studies [24,25]. Several Landsat images of various time spans (1972–2008) have been used by [26] to investigate the land use and land cover changes in Egypt. In order to examine LULC changes in Rwanda, [27] used Landsat images from various years (1987–2016). The recent study by [28] analyzed LULC changes in Kenya by using GIS and remote sensing techniques. Similarly, numerous studies have been reported to observe LULC changes in different catchments in Malaysia [29,30,31].
For the forecasting and simulation of LULC changes, a wide range of models have been developed so far [32,33,34]. For land use modeling, statistical-based methods and machine learning methods have also been used for decades [35,36]. These methods statistically describe the relationships between variables [33]. Markov chain analysis (MCA) or Markov models [37], cellular automata (CA) [38], cellular automata–Markov models (CA–Markov) [39], artificial neural networks (ANNs), binary logistic regression, and fractal models [40,41] are among the most common models used for simulation and prediction of LULC changes [42]. For LULC simulation, CA constitute an appropriate and commonly used model [32]. However, it appears from the literature that so far, no study has been conducted for prediction of LULC changes in the Perak River basin. The combination of ANN-CA for prediction of LULC changes has also not been utilized for this region. Performing such a study is imperative, and its outcomes could be helpful for regional policymakers and authorities in making decisions for future monitoring and planning. Therefore, in this study, ANN-CA simulation was used to forecast the land cover change scenarios for the Perak River basin until the year 2050.
The Perak River basin is an important asset to the people in the region and the primary source of raw water in the region. Hence, the sustainability of the environment in and around the Perak River basin is important, as urbanization is one of the possible challenges to the sustainable growth of the region. Therefore, the aim of this study was to examine the LULC changes in the Perak River basin of Malaysia from the years 2000, 2010, and 2020 with the specific objectives to measure the extent and causative factors of land use in the Perak River basin over the given time span. ANN-CA simulation was employed to determine the future aspects of LULC changes in the Perak River basin for the years 2030, 2040, and 2050.

2. Materials and Methods

2.1. Study Area

Perak is situated in the western part of peninsular Malaysia, bordered to the north side by Kedah and the south coast of Thailand, to the northwest by Palau Penang, to the east by Kelantan and Pahang, and to the south by Selangor, and occupies an area of 21,035 km2. It ranges geographically from 3°30′ N to 6°0′ N longitude and100°0′ E to 102°0′ E latitude, as shown in Figure 1 [31]. The shapefile of the study area was obtained from a website providing data-interpolating variational analysis (DIVA–GIS) (Retrieved on 27 February 2021, https://gadm.org/maps/MYS.html) from which study area map was clipped. The second-largest river basin in peninsular Malaysia is the Perak River basin. It occupies about 70% of the state of Perak and has a river basin area of approximately 14,908 km2. Its reservoirs occur between the Perak–Kelantan–Thailand mountainous areas, bordering the Belum forest in the north of the region. Sg. Pelus, Sg. Kinta, Sg. Batang Padang, and Sg. Bidor are the main tributaries in the surrounding area. River system communities inhabited along the main river include Lenggong, Kuala Kangsar, Gerik, and Parrit. The range of temperature in the Perak River basin is from 22 °C to 35 °C, with sunny days and cold nights all around the year. The rate of humidity varies from 80 to 82% annually, which is relatively high. The mean annual precipitation is recorded to be 3200 mm.
The United States Geological Survey (USGS) earth explorer website was used for acquiring the digital elevation model (DEM) from which the elevation data of this study area were extracted, as shown in Figure 2. The radar system collects information that results in the most precise and detailed topographical map of the surface of Earth ever produced. The Perak River basin area consists primarily of gentle, elongated highlands with mild to steep slopes. The study area has a range of elevation between 27 m and 2168 m above mean sea level (MSL). From Figure 2, it can be observed that areas towards the north side such as Hulu Perak, Kuala Kangsar, Kinta and Batang Padang have the highest elevation rates above mean sea level. Regulations for planning and development in hilly and highlands regions have been established by the government of Malaysia. The regulations were categorized into various elevation ranges: lowlands (less than 150 m), hilly lands (150–300 m), highlands (300–1000 m), and mountainous lands (more than 1000 m). The elevation of the study area is directly proportional to the steepness of slopes, which has a direct connection with soil erosion rate [43].
A slope map of the study area was made in order to estimate the steepness of the region. From the slope map, five classes of distinct steepness were generated. The variation in steepness was based on the topographical features of the study area. The main ranges of steepness in the five classes of the study area were from 0–3.48°, 3.48–7.67°, 7.67–11.83°, 11.83–16.26°, and 16.26–36.14°, as shown in Figure 3. The factor of slope steepness has a significant effect on the rate of soil erosion, soil stability, and sedimentation. The amount of soil erosion increases with the increase in slope gradient due to a rise in the velocity of surface runoff [44].

2.2. Data Acquisitions and Preparation

Landsat images acquired from Landsat 5 and Landsat 8 satellites were utilized to analyze the LULC changes, whereas dates were chosen focusing on the quality and accessibility of Landsat data and weather conditions, as shown in Table 1. The data for Landsat 5 and Landsat 8 showed a minimal consistency issue from the perspective of spatial composition. These issues were rectified during pre-classification phase of the analysis and have been successfully addressed by some previous studies [4,5,6] as well. The visible bands (red, green, and blue) were selected for land use classification. The bands selected for Landsat 5 classification were 2, 3, and 4, whereas 4, 5, and 6 bands were selected for Landsat 8 image classification [45]. Four Landsat images were obtained for the years 2000, 2010, and 2020 from the USGS earth explorer website (Retrieved on 27 February 2021, https://glovis.usgs.gov) using Path/Row 127/56, 127/57, 128/56, and 127/57 for the Perak River basin (Table 1). The datasets were added to ArcGIS to generate LULC maps. The ArcGIS 10.8 software package was used at the different stages of the analysis. The images obtained from the satellite were sensed at various times of the year and all had the same spatial resolution (30 m). The resulting four images were clipped to extract the study area, after checking the Landsat scene by date.

2.3. Classification

In order to achieve high precision of classification, false color composites and indexes of water and vegetation were produced with normalized differences. For improved visualization, the mixture of green, near-infrared (NIR), and red bands (R) was the false color composite that was utilized. Moreover, NDVI, i.e., normalized difference vegetation index, and NDWI, i.e., normalized difference water index, were created for all images. The usage of NDVI in estimation of vegetation covering has been effectively in past years for different purposes [46]. It can be described by Equation (1):
N D V I = N I R R N I R + R
where N I R and R indicate the reflectance of the surface over nearly 0.8 µm and visible band (0.6 µm) in the spectrum of light, respectively. Normalized difference water index ( N D W I ) as shown in Equation (2), on the other hand, demonstrates the maximum water reflectance around an area by using the visible green band (G) (0.5 µm) and the near-infrared band ( N I R ) (0.8 µm) [47,48,49].
N D W I = G N I R G + N I R
False color images helped in detection and visualization of different land characteristics, while N D V I and N D W I were used to distinguish the vegetative and water surfaces, respectively, in the study area. The supervised classification technique has been used effectively in the case of spectral variability in individual types of cover, and hence it was applied for the digital classification of Landsat images [50]. Supervised classification yields superior results in comparison to unsupervised classification, as shown by many recent studies [51,52]. All the satellite images were analyzed by applying per-pixel signature to the same digital number (DN) of various landscape features and used to distinguished the study area into four classes. The defined classes were dense forests, barren and urban lands, waterbodies, and agricultural lands, as shown in Table 2. Training samples were selected by defining polygons around the representative locations for each of the predetermined LULC types. An appropriate spectral signature is one that helps to ensure that there is negligible confusion between land covers that are to be mapped [53]. A total of 120 spectral signatures for each class were captured; therefore, an accumulated 480 training points from all four classes were selected. Then, the signature files were created, which are representative of a class or a cluster. After the creation of signature files, they were added with the base map, which provides the background of the geographical context of the study area.
Later, for the supervised classification of the Landsat images, the maximum-likelihood classification technique was applied. With this classification method, selection of pixel values related to different classes can be controlled by the user [54]. The outcomes obtained from the supervised classification and images of higher resolution show the classification of the study area. Then, kappa coefficient was calculated for the accuracy assessment of classified maps. The kappa coefficient is an accuracy indicator for measuring the overall coherence between the classified images and reference data. Multi-separate variables technique was used to calculate kappa coefficient, which was used to measure the accuracy of classified maps. It was derived from the confusion matrix that relates the classification map and referenced data [55]. The kappa coefficient was calculated by using Equation (3) [56]:
K = N i = 1 m C M i i i = 1 m C i c o r r C i p r e d N 2 i = 1 m C i c o r r C i p r e d
where K represents Cohen’s kappa coefficient, C i c o r r represents the corrected random samples in rows of confusion matrix, C i p r e d represents the predicted random samples in columns of confusion matrix of each class, C M i i represents the diagonal elements of the confusion matrix, and n represents the total number of random samples. The flow chart of all the processes involved in the methodology is shown in Figure 4.

2.4. Change Detection

Change detection is considered as a primary method for identifying the variations in land use pattern of different classes (i.e., waterbodies, barren and urban lands, agricultural lands, and dense forests) in distinct periods of time. The essential characteristic of this method is the capacity to identify the changes in data from various sources and time spans [57]. Change detection was calculated by finding the difference in the area covered by different land use classes between the two periods of 2000 to 2010 and 2010 to 2020. For the calculation of area covered by different land use classes, the “attribute table” was obtained after the land use classification was analyzed. The pixels covered by each class were divided by the total pixels of the study area to obtain percentage change in cover for each class, as shown in Equation (4) [58].
P e r c e n t a g e   c o v e r   ( % ) = p i x e l s   c o v e r e d   b y   l a n d   u s e   c l a s s t o t a l   p i x e l   o f   t h e   s t u d y   a r e a × 100

2.5. Artificial Neural Network–Cellular Automata Modeling

During the last two decades, several models based on simulations have been developed that are used for the modeling of land cover changes around the globe [59]. In this study, a combination of ANN and CA was used to simulate and evaluate the LULC trends of the Perak River basin up to the year 2050 using open source QGIS software version 2.18.25 [60]. The CA feature in QGIS is based on the Markov chain algorithm; i.e., it relies on the present state of land use rather than the previous state [15]. This model generates the output data in the form of tables and maps by combining previous and current land use maps with spatial input parameters [61]. Based on that data, MOLUSCE uses algorithms to train the model. Since ANN algorithms are more accurate than other algorithms, it was used to train the transition potential model of LULC [60]. In recent decades, ANNs have become most common in remote sensing for proper LULC modeling and classification [62]. An ANN is composed of neurons, which are the same as those found in the brains of humans, and it uses them to recognize the trends in data [63]. The most popular form of ANN is the multilayer perceptron (MLP) [64]. The MLP-ANN preprocesses the provided data from land use groups such as barren and urban lands, water bodies, and agricultural lands by using dummy coding of different groupings into a collection of independent variables such as 0 and 1.
The transition potential model used in this study was trained with a momentum of 0.050 and a learning rate of 0.100 for the stabilization of learning graph. Furthermore, the number of iterations was set to 100 to prevent the issue of overfitting in the model. ANN-CA simulation was utilized to simulate the land use changes. The state of the new cell was determined by the existing state of a current cell and changes in the neighborhood cells in CA [65,66]. The ANN-CA simulation selects raster data, such as classes of LULC, raster of spatial parameters, and transition potential model, based on ANN algorithm [67]. Potential changes are determined for each class, and the simulation creates a raster of the most likely transitions. The simulation examines a fixed number of pixels, with the greatest certainty for each transition corresponding to the most likely transitions, and then adjusts the class of the pixel [68]. Multiple iterations of the simulation were performed to achieve subsequent prediction maps of 2030, 2040, and 2050.
Validation was conducted after the ANN-CA simulation, which allows for verifying, comparing, and validating the outcomes achieved. The method of validation was carried out by comparing simulated outcomes to the reference data [68]. Calibration and validation processes are critical factors in validating a simulation model, and [69] stated that there is a clear distinction between these two terms. The predicted map of 2020 was obtained by inserting classified maps of 2000 and 2010 as input data. The predicted map of 2020 was compared with the observed classified map of 2020 to assess the degree of agreement between the pixels of the two maps. The overall kappa coefficient was calculated, the value of which ranges from 0 to 1 [15]. The high degree of agreement and satisfactory value of kappa coefficient indicated the validation of the simulation model. This validated simulation model formed the basis for future predicted maps of 2030, 2040, and 2050. The detailed methodology flow chart of prediction maps is shown in Figure 5.

3. Results and Discussions

3.1. Land Use Change and Accuracy Assessment

The spatial representation of LULC types of four major classes of land cover, namely (i) waterbodies, (ii) agricultural lands, (iii) barren and urban lands, and (iv) dense forests, from the year 2000 to 2020 is shown in Figure 6a–c. These figures illustrate the land use condition in the Perak River basin. Firstly, for the composition of the band, Landsat images were inputted into an image analysis tool. Then, by using Landsat images, the signatures of all the spectral classes were identified. For the purpose of supervised classification, the classes with the same spectral signatures were merged. This allows all the pixels that are included in an image to be automatically assigned to the land cover classes [70]. The maximum-likelihood classifier was employed by applying training sets developed from Landsat images.
The accuracy assessment of the classified images was the next stage after classification. The stratified samplings were generated by using Google Earth images, and then accuracy assessment of LULC maps was performed. Confusion matrix was created between testing samples and classified images for the measurement of accuracy. It is a simple way to evaluate how frequently a pixel misidentifies class [15]. For the formation of confusion matrix table, 189, 202, and 177 random points were selected on the classified maps for years 2000, 2010, and 2020 respectively (Table 3). These random points represented different land use classes that were compared with Google Earth images. These random points represented the producer’s values and user’s corrected values. The dataset table was transformed into a confusion matrix by comparing the total producer’s values and the user’s corrected values in the respective classes. According to Anderson’s classification scheme, the kappa values of classified images are considered satisfactory if they reach the minimum accuracy of 0.85 [45]. Several attempts were made by increasing training sample size to enhance accuracy until reaching kappa coefficient of more than 0.85 as the desired accuracy. The overall accuracies of classified images for the years 2000, 2010, and 2020 were 0.86, 0.88, and 0.91 respectively, as shown in Table 3. Additionally, the table also describes the accuracy of each class for all the images. Based on the outcomes shown in the table, it is apparent that the accuracy of classification of waterbodies (96%) is relatively higher as compared to agricultural and barren land in the year 2000. This is because the forests and waterbodies represent clear pixels and make the process of classification simple [71]. In comparison to this, the accuracy of the barren and agricultural classes was comparatively poor due to the mixed environmental characteristics. The accuracy assessment for all images indicates that it is difficult to determine agricultural lands and areas with fewer forests with high precision due to the spatial resolution restrictions of the Landsat images. The existence of waterbodies and farmlands in barren areas contributes to a rise in the amount of mixed pixels. For appropriate mapping of agricultural and barren lands, ancillary data and visual representation of images was used. For the classified images of the years 2010 and 2020, the confusion matrix in Table 3 demonstrates the increase in accuracy of classifications. The accessibility of ancillary data for the classification was the reason for the enhancement of accuracy. The accuracy for agricultural and barren land classification also improved due to the interconnection of various data sources, i.e., the data from Google Earth and Landsat images data. These findings offer a significant foundation for the future study of LULC changes.

3.2. Land Use Changes in Perak River Basin

The details of land use changes in the Perak River basin for the years 2000, 2010, and 2020 are given in Table 4 and Table 5, as well as Figure 7. In the year 2000, the major land cover in the study area was occupied by dense forests comprising 56% of the total area, followed by agricultural lands (27.73%), barren and urban lands (14.63%), and waterbodies (1.62%). Compared to the year 2000, in 2010, the land cover experienced a tremendous decrease in dense forests (43.69%), followed by a considerable increase in barren and urban lands (15.71%) as well as a slight increase in the area of waterbodies (1.68%) and agricultural lands (39.05%). However, for the year 2020, the observed LULC changes were dominated by a sizeable increase in the area of barren and urban lands (20.39%), followed by a dramatic decrease in dense forests (35.47%) and a slight increase in waterbodies (1.75%) and agricultural lands (42.46%), as shown in Table 4.
The main factors contributing to the increase in barren and urban lands in the Perak River basin were rapid expansion of urban development, the increase in the growth of agricultural cultivation that eventually contributed to the degradation of forest, and the variation in the amount of annual rainfall [72,73,74]. Such LULC change analysis is important for a deeper understanding of how these developments were changing the hydrological process of the Perak River basin. Remotely sensed image integration with GIS has the ability to provide a good basis for comparing the influencing factors to the dynamics of the river basin. The outcomes obtained from the analysis have shown that in the last two decades, the land use in Perak, Malaysia has been changed. The main components of the river basin have been decreased dramatically, such as dense forests. Barren land, by comparison, has seen more changes and increased more considerably than any other class of land use. Due to the increase in barren and urban land, the huge number of light forests have been converted into development zones, and thus the area of impermeable surfaces has been enhanced [75]. The rapid urbanization was mainly due to industrial development and growth in the rate of population. A slight positive increase was detected in the class of waterbodies in the year 2010, which can be attributed to the rise in water bodies, such as reservoirs that perform the function of retention ponds in the season of floods. Increases in barren lands and rainfall played significant roles in changing the hydrological system of the river basin [76,77,78]. Precipitation over recently built areas results in increased rates of runoff, which creates the issue of floods in the low-lying region, according to various studies [79]. Additionally, in urban areas where the capacity for flood retention and drainage is not adequate, heavy rainfall may also lead to an increase in stagnant water. Factors influencing the river basin, i.e., change in climate, land use, and changes in soil infiltration rate, may also result in significant deterioration of the river basin hydrological process [80]. Therefore, the analysis of land use changes can also provide key details on potential improvements to the hydrology of the river basins.

3.3. Prediction of Land Use Changes Using ANN-CA

3.3.1. Transition Potential Modeling Using ANN

The transition potential model was trained using an ANN for identifying LULC changes between the years 2000 and 2010. The changes in the area of different classes were used to create the LULC transition matrix, which is an input for the ANN to obtain the transition probability. The ANN was run with the learning rate of 0.09 and a maximum iteration of 100, and momentum was set to 0.050. These inputs were found to be optimal to train the ANN. ANN-CA simulated the prediction map of 2020 based on the transition potential model and spatial parameters.

3.3.2. Validation of ANN-CA Simulation

The simulation results were statistically validated by kappa coefficient and percentage of correctness. The kappa coefficient was calculated using the guidelines in QGIS [60]. The same method of validation has been successfully used in the previous studies concerning land use prediction using ANN-CA simulation [15,81]. The parameters used for validation were obtained from the simulated and classified reference maps of 2020. The comparison of these maps produced the kappa coefficient and percentage of correctness. The value obtained for the kappa coefficient was 0.83. The graph shown in Figure 8 represents the comparison and correlation between the predicted and reference maps of 2020. The validation graph tracks the agreement between data points of the two maps, which suggests an accurate prediction has been conducted, whereas a divergence from the reference line signifies the inaccurately simulated land cover. The minimal errors related to the simulation can be attributed to the inaccurate identification of land use patterns and misinterpretation of spatial parameters (i.e., network of roads, elevation, and slope maps). Similar trends have been reported during the validation of simulated maps in previous studies [61,81,82].
The predicted maps of 2030, 2040, and 2050 were obtained by implementing multiple iterations of the validated simulation results of 2020. The percentage correctness of the simulation was 73.41%. The percentage of correctness suggests that there would be a 73.41% chance of correctness in predicted outcomes of 2030, 2040, and 2050 simulations when compared with the actual LULC maps of 2030, 2040, and 2050. It also indicates that if the current LULC changes would proceed in the same trend, then the future land use pattern for the Perak River basin would be comparable to the maps shown in Figure 9.

3.3.3. Prediction Maps of Year 2030, 2040, and 2050

The outcomes of ANN-CA simulation as shown in Figure 9 illustrate that the Perak River basin will see a slight increase in areas of barren lands and agricultural lands in the years leading to 2030, 2040, and 2050. Based on the simulation results, there will be a minor decrease in the area of waterbodies of (−1.56 km2) between the years 2030 to 2040 and (−2.53 km2) from the year 2040 to 2050 due to the expected increases in bedload and sedimentation in the flood retention ponds and rainfed catchments, as shown in Table 6. Similarly, an increase in the area of agricultural lands of (60.6 km2) was predicted from the year 2030 to 2040 and (46.9 km2) from the year 2040 to 2050. However, dense forests will decrease (−335.91 km2) from the year 2030 to 2040 and (−94.77 km2) from the year 2040 to 2050 due to the expected increase in rates of deforestation and conversion of forests into developed areas. Similarly, there will be an increase in the area of barren and urban lands (65.95 km2) from the year 2030 to 2040 and (51.24 km2) from the year 2040 to 2050 due to an expected increase in commercial, residential, and industrial areas.

3.3.4. Land Use Transition Matrix of Simulated Outcomes

A transition matrix table illustrates how future land use classes may change during the periods of 2020–2030, 2030–2040, and 2040–2050. The transition matrix table was formed by comparing the observed classified map of 2020 with the simulated map of 2030, and similarly comparing the simulated map of 2030 with 2040 and 2040 with 2050. The transition matrix indicates the changes in land use patterns in the simulated years. The observed changes in land use classes are shown in Table 7. From the table, it is evident that from 2020–2030, the major transition occurs in dense forest to agricultural lands. Similarly, there is a shift of agricultural lands to barren and urban lands. However, minor transformations occur in the classes of waterbodies, barren and urban lands. This trend of changes in the land use patterns remains constant in the years leading up to 2040 and 2050.

4. Conclusions

The present study has successfully provided an application of ANN-CA for the monitoring and prediction of LULC changes and spatial distribution patterns. The growth in LULC changes and future prediction of the Perak River basin were analyzed by using multi-temporal data from the year 2000 to 2020. The LULC changes maps classified for the years 2000, 2010, and 2020 demonstrated the major changes in the study area. Furthermore, areas of various LULC classes were calculated to monitor decadal changes in each class. According to the findings, barren and urban lands in the Perak River basin grew rapidly between 2000 and 2020. In 2000, the total area of barren and urban lands in the region was 3270.01 km2, which had increased up to 3817.07 km2 by the year 2020, leading to extreme soil erosion risk in the catchment. This shows that there was an approximately five-fold increase in barren and urban lands during this period of time. The decadal change in LULC mapping of the Perak River basin indicated that the dense forests significantly decreased from 12,512.30 km2 in the year 2000 to only 7925.04 km2 in 2020. However, a slight increases were observed in the areas of waterbodies, from 363.18 km2 in 2000 to 391.89 km2 in 2020, and agricultural lands, from 6195.98 km2 in 2000 to 9603.39 km2 in 2020. Dense forests near the Perak River basin were transformed into non-agricultural lands initially and later converted into a developed built-up area.
The study suggested that an ANN-CA simulation model is one of the most appropriate and effective approaches for simulating complex types of LULC changes, and it can be used for other regions with a similar level of complexity. From the outcomes of the ANN-CA simulation, barren and urban lands will expand by many folds in commercial, industrial, and residential areas over the next decades up to 2050. The outcomes of this study on LULC changes and its future prediction might be useful for regional policymakers and authorities in developing sustainable urban planning and improving living standards.

Author Contributions

Conceptualization, M.R.U.M.; data curation, M.T.Z. and M.F.B.; formal analysis, M.T.Z. and M.R.U.M.; funding acquisition, M.R.U.M.; investigation, M.T.Z. and M.R.U.M.; methodology, M.T.Z. and M.F.B.; project administration, M.R.U.M.; resources, M.R.U.M.; supervision, M.R.U.M.; validation, M.T.Z.; writing—original draft, M.T.Z.; writing—review and editing, M.R.U.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by YUTP research project (cost center; 015LC0-044).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Further information on data presented in this work can be obtained from the authors.

Acknowledgments

The authors would like to gratefully acknowledge the financial support provided under the YUTP grant (cost center; 015LC0-044). The first author would like to acknowledge Universiti Teknologi PETRONAS for providing support under Graduate Assistantship scheme. Free data availability of Landsat images at the website of United States Geological Survey (USGS) Department is also greatly acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Singh, A. Review Article Digital change detection techniques using remotely-sensed data. Int. J. Remote. Sens. 1989, 10, 989–1003. [Google Scholar] [CrossRef] [Green Version]
  2. Houghton, R.A. The worldwide extent of land-use change. Bioscience 1994, 44, 305–313. [Google Scholar] [CrossRef]
  3. Hathout, S. The use of GIS for monitoring and predicting urban growth in East and West St Paul, Winnipeg, Manitoba, Canada. J. Environ. Manag. 2002, 66, 229–238. [Google Scholar] [CrossRef]
  4. Fei, L.; Shuwen, Z.; Jiuchun, Y.; Liping, C.; Haijuan, Y.; Kun, B. Effects of land use change on ecosystem services value in West Jilin since the reform and opening of China. Ecosyst. Serv. 2018, 31, 12–20. [Google Scholar] [CrossRef]
  5. Guzha, A.; Rufino, M.C. Impacts of land use and land cover change on surface runoff, discharge and low flows: Evidence from East Africa. J. Hydrol. Reg. Stud. 2018, 15, 49–67. [Google Scholar] [CrossRef]
  6. López, E.; Bocco, G.; Mendoza, M.; Duhau, E. Predicting land-cover and land-use change in the urban fringe: A case in Morelia city, Mexico. Landsc. Urban Plan. 2001, 55, 271–285. [Google Scholar] [CrossRef]
  7. Jat, M.K.; Garg, P.; Khare, D. Monitoring and modelling of urban sprawl using remote sensing and GIS techniques. Int. J. Appl. Earth Obs. Geoinf. 2008, 10, 26–43. [Google Scholar] [CrossRef]
  8. Guerra, F.; Puig, H.; Chaume, R. The forest-savanna dynamics from multi-date Landsat-TM data in Sierra Parima, Venezuela. Int. J. Remote. Sens. 1998, 19, 2061–2075. [Google Scholar] [CrossRef]
  9. Lu, D.; Mausel, P.; Brondizio, E.; Moran, E. Change detection techniques. Int. J. Remote. Sensing 2004, 25, 2365–2401. [Google Scholar] [CrossRef]
  10. Roy, A.; Inamdar, A.B. Multi-temporal Land Use Land Cover (LULC) change analysis of a dry semi-arid river basin in western India following a robust multi-sensor satellite image calibration strategy. Heliyon 2019, 5, e01478. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Helmer, E.H.; Brown, S.; Cohen, W.B. Mapping montane tropical forest successional stage and land use with multi-date Landsat imagery. Int. J. Remote. Sens. 2000, 21, 2163–2183. [Google Scholar] [CrossRef]
  12. Hussain, M.; Chen, D.; Cheng, A.; Wei, H.; Stanley, D. Change detection from remotely sensed images: From pixel-based to object-based approaches. ISPRS J. Photogramm. Remote. Sens. 2013, 80, 91–106. [Google Scholar] [CrossRef]
  13. Singh, Y.; Ferrazzoli, P.; Rahmoune, R. Flood monitoring using microwave passive remote sensing (AMSR-E) in part of the Brahmaputra basin, India. Int. J. Remote. Sens. 2013, 34, 4967–4985. [Google Scholar] [CrossRef]
  14. Alexakis, D.D.; Agapiou, A.; Tzouvaras, M.; Themistocleous, K.; Neocleous, K.; Michaelides, S.; Hadjimitsis, D. Integrated use of GIS and remote sensing for monitoring landslides in transportation pavements: The case study of Paphos area in Cyprus. Nat. Hazards 2013, 72, 119–141. [Google Scholar] [CrossRef]
  15. Yatoo, S.A.; Sahu, P.; Kalubarme, M.H.; Kansara, B.B. Monitoring land use changes and its future prospects using cellular automata simulation and artificial neural network for Ahmedabad city, India. GeoJournal 2020, 85, 1–22. [Google Scholar] [CrossRef]
  16. Vibhute, A.D.; Gawali, B.W. Analysis and modeling of agricultural land use using remote sensing and geographic information system: A review. Int. J. Eng. Res. Appl. 2013, 3, 81–91. [Google Scholar]
  17. Rahman, A.; Kumar, S.; Fazal, S.; Siddiqui, M.A. Assessment of land use/land cover change in the north-west district of Delhi using remote sensing and GIS techniques. J. Indian Soc. Remote. Sens. 2012, 40, 689–697. [Google Scholar] [CrossRef]
  18. Attri, P.; Chaudhry, S.; Sharma, S. Remote sensing & GIS based approaches for LULC change detection—A review. Int. J. Curr. Eng. Technol. 2015, 5, 3126–3137. [Google Scholar]
  19. Jiang, X.; Lu, D.; Moran, E.; Calvi, M.F.; Dutra, L.V.; Li, G. Examining impacts of the Belo Monte hydroelectric dam construction on land-cover changes using multitemporal Landsat imagery. Appl. Geogr. 2018, 97, 35–47. [Google Scholar] [CrossRef]
  20. Hazarika, N.; Das, A.K.; Borah, S.B. Assessing land-use changes driven by river dynamics in chronically flood affected Upper Brahmaputra plains, India, using RS-GIS techniques. Egypt. J. Remote Sens. Space Sci. 2015, 18, 107–118. [Google Scholar] [CrossRef] [Green Version]
  21. Lopez-Granados, E.; Mendoza, M.E.; Gonzalez, D.I. Linking geomorphologic knowledge, RS and GIS techniques for analyzing land cover and land use change: A multitemporal study in the Cointzio watershed, Mexico. Rev. Ambiente Agua 2013, 8, 18–37. [Google Scholar]
  22. Lambin, E.F.; Geist, H.; Lepers, E. Dynamics of land-use and land-cover change in tropical regions. Annu. Rev. Environ. Resour. 2003, 28, 205–241. [Google Scholar] [CrossRef] [Green Version]
  23. Serra, P.; Pons, X.; Sauri, D. Land-cover and land-use change in a Mediterranean landscape: A spatial analysis of driving forces integrating biophysical and human factors. Appl. Geogr. 2008, 28, 189–209. [Google Scholar] [CrossRef]
  24. Chowdhury, M.; Hasan, M.E.; Al Mamun, M.M.A. Land use/land cover change assessment of Halda watershed using remote sensing and GIS. Egypt. J. Remote. Sens. Space Sci. 2018, 23, 63–75. [Google Scholar] [CrossRef]
  25. Mohamed, M.A. Monitoring of temporal and spatial changes of land use and land cover in metropolitan regions through remote sensing and GIS. Nat. Resour. 2017, 8, 353–369. [Google Scholar] [CrossRef] [Green Version]
  26. El Gammal, E.A.; Salem, S.M.; El Gammal, A.E.A. Change detection studies on the world’s biggest artificial lake (Lake Nasser, Egypt). Egypt. J. Remote Sens. Space Sci. 2010, 13, 89–99. [Google Scholar] [CrossRef] [Green Version]
  27. Akinyemi, F.O. Land change in the central Albertine rift: Insights from analysis and mapping of land use-land cover change in north-western Rwanda. Appl. Geogr. 2017, 87, 127–138. [Google Scholar] [CrossRef]
  28. Cheruto, M.C.; Kauti, M.K.; Kisangau, P.D.; Kariuki, P. Assessment of land use and land cover change using gis and remote sensing techniques: A case study of Makueni County, Kenya. J. Remote. Sens. GIS 2016, 05, 1000175. [Google Scholar] [CrossRef]
  29. Yusof, F.M.; Jamil, N.R.; Laew, N.I.C.; Aini, N.; Manaf, L.A. Land use change and soil loss risk assessment by using geographical information system (GIS): A case study of lower part of Perak River. IOP Conf. Ser. Earth Environ. Sci. 2016, 37, 12065. [Google Scholar] [CrossRef] [Green Version]
  30. Noh, N.S.M.; Sidek, L.M.; Wayayok, A.; Abdullah, A.F.; Basri, H.; Farhan, S.A.; Sulaiman, T.; Ariffin, A.B. Erosion and sediment control best management practices in agricultural farms for effective reservoir sedimentation management at Cameron Highlands. Int. J. Recent Technol. Eng. 2019, 8, 6198–6205. [Google Scholar] [CrossRef]
  31. Hanif, M.F.; ul Mustafa, M.R.; Hashim, A.M.; Yusof, K.W. Spatio-temporal change analysis of perak river basin using remote sensing and GIS. In Proceeding of the 2015 International Conference on Space Science and Communication, Lakawi, Malaysia, 10–12 August 2015. [Google Scholar]
  32. Wang, C.; Lei, S.; Elmore, A.J.; Jia, D.; Mu, S. Integrating temporal evolution with cellular automata for simulating land cover change. Remote. Sens. 2019, 11, 301. [Google Scholar] [CrossRef] [Green Version]
  33. Arsanjani, J.J. Characterizing, monitoring, and simulating land cover dynamics using GlobeLand30: A case study from 2000 to 2030. J. Environ. Manag. 2018, 214, 66–75. [Google Scholar] [CrossRef]
  34. Brown, D.; Band, L.E.; Green, K.O.; Irwin, E.G.; Jain, A.; Lambin, E.F.; Pontius, R.G.; Seto, K.C.; Turner, B.L., II; Verburg, P.H. Advancing Land Change Modeling: Opportunities and Research Requirements; The National Research Council Press: Amsterdam, The Netherlands, 2013. [Google Scholar]
  35. Jiang, X.; Lin, M.; Zhao, J. Woodland cover change assessment using decision trees, support vector machines and artificial neural networks classification algorithms. In Proceedings of the 2011 Fourth International Conference on Intelligent Computation Technology and Automation, Shenzhen, China, 28–29 March 2011. [Google Scholar]
  36. Zhang, P.; Gong, M.; Su, L.; Liu, J.; Li, Z. Change detection based on deep feature representation and mapping transformation for multi-spatial-resolution remote sensing images. ISPRS J. Photogramm. Remote Sens. 2016, 116, 24–41. [Google Scholar] [CrossRef]
  37. Sivakumar, V. Urban mapping and growth prediction using remote sensing and GIS techniques, Pune, India. ISPRS Int. Arch. Photogramm. Remote. Sens. Spat. Inf. Sci. 2014, XL-8, 967–970. [Google Scholar] [CrossRef] [Green Version]
  38. Sinha, S.; Sharma, L.K.; Nathawat, M.S. Improved land-use/land-cover classification of semi-arid deciduous forest land-scape using thermal remote sensing. Egypt. J. Remote Sens. Space Sci. 2015, 18, 217–233. [Google Scholar]
  39. Subedi, P.; Subedi, K.; Thapa, B. Application of a hybrid cellular automaton—Markov (CA-Markov) model in land-use change prediction: A case study of Saddle Creek Drainage Basin, Florida. Appl. Ecol. Environ. Sci. 2013, 1, 126–132. [Google Scholar] [CrossRef] [Green Version]
  40. Zhan, X.; Sohlberg, R.; Townshend, J.; DiMiceli, C.; Carroll, M.; Eastman, J.; Hansen, M.; DeFries, R. Detection of land cover changes using MODIS 250 m data. Remote. Sens. Environ. 2002, 83, 336–350. [Google Scholar] [CrossRef]
  41. Zeng, Y.; Wu, G.; Zhan, F. Modeling spatial land use pattern using autologistic regression. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2008, 37, 115–119. [Google Scholar]
  42. Schneider, L.C.; Pontius, R.G., Jr. Modeling land-use change in the Ipswich watershed, Massachusetts, USA. Agric. Ecosyst. Environ. 2001, 85, 83–94. [Google Scholar] [CrossRef]
  43. Harun, N.; Yaacob, W.Z.W.; Simon, N. Potential areas for the near surface disposal of radioactive waste in Pahang. AIP Conf. Proc. 2016, 1784, 060021. [Google Scholar] [CrossRef] [Green Version]
  44. Omar, M.N.; Rahaman, Z.A.; Hashim, M. the development of a soil erosion risk map for Perak, Malaysia. Int. J. Acad. Res. Bus. Soc. Sci. 2018, 8, 1108–1123. [Google Scholar] [CrossRef]
  45. USGS. United States Geological Survey. Available online: https://www.usgs.gov/faqs/what-are-band-designations-landsat-satellites?qt-news_science_products=0#qt-news_science_products (accessed on 27 February 2021).
  46. Sun, Y.; Ren, H.; Zhang, T.; Zhang, C.; Qin, Q. Crop leaf area index retrieval based on inverted difference vegetation index and NDVI. IEEE Geosci. Remote. Sens. Lett. 2018, 15, 1662–1666. [Google Scholar] [CrossRef]
  47. McFeeters, S.K. Using the normalized difference water index (NDWI) within a geographic information system to detect swimming pools for mosquito abatement: A practical approach. Remote. Sens. 2013, 5, 3544–3561. [Google Scholar] [CrossRef] [Green Version]
  48. Li, W.; Du, Z.; Ling, F.; Zhou, D.; Wang, H.; Gui, Y.; Sun, B.; Zhang, X. A comparison of land surface water mapping using the normalized difference water index from TM, ETM+ and ALI. Remote. Sens. 2013, 5, 5530–5549. [Google Scholar] [CrossRef] [Green Version]
  49. Delbart, N.; Kergoat, L.; Le Toan, T.; Lhermitte, J.; Picard, G. Determination of phenological dates in boreal regions using normalized difference water index. Remote. Sens. Environ. 2005, 97, 26–38. [Google Scholar] [CrossRef] [Green Version]
  50. Mas, J.F. Monitoring land-cover changes: A comparison of change detection techniques. Int. J. Remote. Sens. 1999, 20, 139–152. [Google Scholar] [CrossRef]
  51. Karakuş, C.B. The Impact of Land Use/Land Cover (LULC) changes on land surface temperature in sivas city center and its surroundings and assessment of urban heat island. Asia-Pac. J. Atmos. Sci. 2019, 55, 669–684. [Google Scholar] [CrossRef]
  52. Ibrahim, F.U.; Rasul, G. Urban land use land cover changes and their effect on land surface temperature: Case study using Dohuk City in the Kurdistan Region of Iraq. Climate 2017, 5, 13. [Google Scholar] [CrossRef] [Green Version]
  53. Gao, J.; Liu, Y. Determination of land degradation causes in Tongyu County, Northeast China via land cover change detection. Int. J. Appl. Earth Obs. Geoinf. 2010, 12, 9–16. [Google Scholar] [CrossRef]
  54. Bolstad, P.T.M. Rapid maximum likelihood classification. Photogramm. Eng. Remote Sens. 1991, 57, 67–74. [Google Scholar]
  55. El-Tantawi, A.M.; Bao, A.; Chang, C.; Liu, Y. Monitoring and predicting land use/cover changes in the Aksu-Tarim River Basin, Xinjiang-China (1990–2030). Environ. Monit. Assess. 2019, 191, 480. [Google Scholar] [CrossRef]
  56. Tallón-Ballesteros, A.J.; Riquelme, J.C. Data mining methods applied to a digital forensics task for supervised machine learning. In Computational Intelligence in Digital Forensics: Forensic Investigation and Applications; Springer: Berlin/Heidelberg, Germany, 2014; Chapter 17; pp. 413–428. [Google Scholar]
  57. Zhang, J.; Goodchild, M.F. Uncertainty in Geographical Information; CRC Press: Boca Raton, FL, USA, 2002. [Google Scholar]
  58. Kumar, M.; Mondal, I.; Pham, Q.B. Monitoring forest landcover changes in the Eastern Sundarban of Bangladesh from 1989 to 2019. Acta Geophys. 2021, 69, 561–577. [Google Scholar] [CrossRef]
  59. Han, H.; Yang, C.; Song, J. Scenario simulation and the prediction of land use and land cover change in Beijing, China. Sustainability 2015, 7, 4260–4279. [Google Scholar] [CrossRef] [Green Version]
  60. MOLUSCE. Modules for Land Use Change Evaluation. Available online: https://wiki.gis-lab.info/w/Landscape_change_analysis_with_MOLUSCE_-_methods_and_algorithms (accessed on 21 March 2021).
  61. Al-Rubkhi, M.N.A.G. Land Use Change Analysis and Modeling Using Open Source (QGis)-Case Study: Boasher Willayat. Ph.D. Thesis, College of Arts and Social Science, Department of Geography, Sultan Qaboos University, Muscat, Oman, 2017. [Google Scholar]
  62. Gasarovic, M.; Jogun, T. The effect of fusing Sentinel-2 bands on land-cover classification. Int. J. Remote Sens. 2018, 39, 822–841. [Google Scholar] [CrossRef]
  63. Pijanowski, B.C.; Brown, D.; A Shellito, B.; A Manik, G. Using neural networks and GIS to forecast land use changes: A Land Transformation Model. Comput. Environ. Urban Syst. 2002, 26, 553–575. [Google Scholar] [CrossRef]
  64. Rumelhart, D.; Hinton, G.; Williams, R. Learning Internal Representations by Error Propagation; California University of San Diego, La Jolla Institute for Cognitive Science: San Diego, CA, USA, 1985. [Google Scholar]
  65. Lau, K.H.; Kam, B.H. A cellular automata model for urban land-use simulation. Environ. Plan. B Plan. Des. 2005, 32, 247–263. [Google Scholar] [CrossRef] [Green Version]
  66. Lantman, J.v.S.; Verburg, P.H.; Bregt, A.; Geertman, S. Core principles and concepts in land-use modelling: A literature review. In Land-Use Modelling in Planning Practice; Koomen, E., Borsboom-van, B.J., Eds.; Springer: Dordrecht, The Netherlands, 2011; pp. 35–57. [Google Scholar] [CrossRef] [Green Version]
  67. Jogun, T.; Lukić, A.; Gašparović, M. Simulation model of land cover changes in a post-socialist peripheral rural area: Pozga-Slavonia County, Croatia. Croat. Geogr. Bull. 2019, 81, 31–59. [Google Scholar] [CrossRef] [Green Version]
  68. GIS-Lab. Landscape Change Analysis with Methods of Land Use Change Evaluation (MOLUSCE) Methods and Algorithms 2018. Available online: https://wiki.gislab.info/w/Landscape_change_analysis_with_MOLUSCE_methods_and_algorithms (accessed on 8 March 2021).
  69. Pontius, R.G.; Huffaker, D.; Denman, K. Useful techniques of validation for spatially explicit land-change models. Ecol. Model. 2004, 179, 445–461. [Google Scholar] [CrossRef]
  70. Wright, J.; Lillesand, T.M.; Kiefer, R.W. Remote sensing and image interpretation. Geogr. J. 1980, 146, 448. [Google Scholar] [CrossRef]
  71. Jakovljevic, G.; Govedarica, M.; Álvarez-Taboada, F. Waterbody mapping: A comparison of remotely sensed and GIS open data sources. Int. J. Remote. Sens. 2018, 40, 2936–2964. [Google Scholar] [CrossRef]
  72. Ali, T.; Shahbaz, B.; Suleri, A. Analysis of myths and realities of deforestation in Northwest Pakistan: Implications for forestry extension. Int. J. Agric. Biol. 2006, 8, 107–110. [Google Scholar]
  73. Butt, A.; Shabbir, R.; Ahmad, S.S.; Aziz, N. Land use change mapping and analysis using remote sensing and GIS: A case study of Simly watershed, Islamabad, Pakistan. Egypt. J. Remote Sens. Space Sci. 2015, 18, 251–259. [Google Scholar]
  74. Pakistan, I.U.C.N. Rapid Environmental Appraisal of Developments in and Around Murree Hills; Technical Report; IUCN: Karachi, Pakistan, 2005. [Google Scholar]
  75. Gong, P.; Li, X.; Zhang, W. 40-Year (1978–2017) human settlement changes in China reflected by impervious surfaces from satellite remote sensing. Sci. Bull. 2019, 64, 756–763. [Google Scholar] [CrossRef] [Green Version]
  76. Al-Faraj, F.A.; Scholz, M. Impact of upstream anthropogenic river regulation on downstream water availability in transboundary river watersheds. Int. J. Water Resour. Dev. 2014, 31, 28–49. [Google Scholar] [CrossRef]
  77. Drieschova, A.; Giordano, M.; Fischhendler, I. Governance mechanisms to address flow variability in water treaties. Glob. Environ. Chang. 2008, 18, 285–295. [Google Scholar] [CrossRef]
  78. Veldkamp, T.; Wada, Y.; Aerts, J.; Döll, P.; Gosling, S.N.; Liu, J.; Masaki, Y.; Oki, T.; Ostberg, S.; Pokhrel, Y.; et al. Water scarcity hotspots travel downstream due to human interventions in the 20th and 21st century. Nat. Commun. 2017, 8, 15697. [Google Scholar] [CrossRef]
  79. Villarini, G.; Smith, J.A.; Baeck, M.L.; Sturdevant-Rees, P.; Krajewski, W.F. Radar analyses of extreme rainfall and flooding in urban drainage basins. J. Hydrol. 2010, 381, 266–286. [Google Scholar] [CrossRef]
  80. Bronstert, A. Floods and climate change: Interactions and impacts. Risk Anal. 2003, 23, 545–557. [Google Scholar] [CrossRef]
  81. Saputra, M.H.; Lee, H.S. Prediction of land use and land cover changes for North Sumatra, Indonesia, using an artificial-neural-network-based cellular automaton. Sustainability 2019, 11, 3024. [Google Scholar] [CrossRef] [Green Version]
  82. Memarian, H.; Balasundram, S.K.; Bin Talib, J.; Sung, C.T.B.; Sood, A.M.; Abbaspour, K. Validation of CA-Markov for Simulation of land use and cover change in the Langat Basin, Malaysia. J. Geogr. Inf. Syst. 2012, 04, 542–554. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Description of the study area.
Figure 1. Description of the study area.
Water 13 02286 g001
Figure 2. The elevation map of the study area.
Figure 2. The elevation map of the study area.
Water 13 02286 g002
Figure 3. Slope map of the study area.
Figure 3. Slope map of the study area.
Water 13 02286 g003
Figure 4. Detailed flow chart of methodology.
Figure 4. Detailed flow chart of methodology.
Water 13 02286 g004
Figure 5. Flow chart of future LULC prediction mapping.
Figure 5. Flow chart of future LULC prediction mapping.
Water 13 02286 g005
Figure 6. Land use and land cover classifications of Perak River basin for: (a) year 2000; (b) year 2010; and (c) year 2020.
Figure 6. Land use and land cover classifications of Perak River basin for: (a) year 2000; (b) year 2010; and (c) year 2020.
Water 13 02286 g006
Figure 7. LULC change graph for 2000, 2010, and 2020 for Perak River basin.
Figure 7. LULC change graph for 2000, 2010, and 2020 for Perak River basin.
Water 13 02286 g007
Figure 8. Validation graph of predicted and referenced data points.
Figure 8. Validation graph of predicted and referenced data points.
Water 13 02286 g008
Figure 9. Predicted LULC of Perak River basin for: (a) year 2030; (b) year 2040; and (c) year 2050.
Figure 9. Predicted LULC of Perak River basin for: (a) year 2030; (b) year 2040; and (c) year 2050.
Water 13 02286 g009
Table 1. Detailed data for the Landsat images used for the study area.
Table 1. Detailed data for the Landsat images used for the study area.
YearLandsat Scene IDPathDate AcquiredResolution (m)RowEarth-Sun Distance
2000LT51270562000237BKT0012724 August 200030561.0109543
LT51270572000237BKT0012724 August 200030571.0109543
LT51280562000020DKI0012820 January 200030560.9839503
LT51280572000148BKT0112827 May 200030571.0132773
2010LT51270562010056BKT0012725 February 201030560.9898357
LT51270572010360BKT0012726 December 201030570.9834906
LT51280562010047BKT0012816 February 201030560.9878930
LT51280572010047BKT0012816 February 201030570.9878930
2020LC81270562020228LGN0012715 August 202030561.0128081
LC81270572020228LGN0012715 August 202030571.0128080
LC81280562020315LGN0012810 November 202030560.9902608
LC81280572020267LGN0012823 September 202030571.0034227
Table 2. The classification scheme for land use and land cover (LULC).
Table 2. The classification scheme for land use and land cover (LULC).
ClassDescription
WaterbodiesRivers, open water, lakes, ponds, and reservoirs.
Barren and urban landsLand areas of exposed soil and barren areas influenced by humans.
Dense forestsContinuous stands of trees, many of which may attain a height of 50 m, including natural forest, mangrove, and plantation forest.
Agricultural landsMainly composed of grass, vegetation, crop plants, cultivated lands, and shrub lands.
Table 3. Accuracy assessment of the LULC classifications in the Perak River basin.
Table 3. Accuracy assessment of the LULC classifications in the Perak River basin.
2000ALWBLDFTotalUA (%)
AL360134090
W047114996
BL124945687
DF302394489
Total40495447189
PA (%)90959183
OA (%)90
K 0.86
2010ALWBLDFTotalUA (%)
AL550135993
W042114495
BL425716489
DF300323591
Total62445937202
PA (%)89949786
OA (%)92
K 0.88
2020ALWBLDFTotalUA (%)
AL610026397
W0450045100
BL302703090
DF201363992
Total66452838177
PA (%)921009695
OA (%)95
K 0.91
Note: AL = Agricultural lands; W = Waterbodies; BL = Barren and urban lands; DF = Dense forests; UA = User’s accuracy; PA = Producer’s accuracy; OA = Overall accuracy; K = Kappa coefficient.
Table 4. Outcomes of the classification of LULC for images from the years 2000, 2010, and 2020, demon-strating the area of each class and its percentage of the Perak River basin.
Table 4. Outcomes of the classification of LULC for images from the years 2000, 2010, and 2020, demon-strating the area of each class and its percentage of the Perak River basin.
Years200020102020
LULC (Area)Area (km2)%Area (km2)%Area (km2)%
Dense Forests12,512.3056.029739.5443.597925.0435.47
Agricultural Lands6195.9827.738714.0839.049603.3942.46
Barren and Urban Lands3270.0114.633510.3915.713817.0720.29
Waterbodies363.181.62377.471.68391.791.75
Table 5. Changes in area and percentage cover of LULC classes between the years 2000–2010, and 2010–2020 of the Perak River basin.
Table 5. Changes in area and percentage cover of LULC classes between the years 2000–2010, and 2010–2020 of the Perak River basin.
Years2000–20102010–2020
LULC (Area)Area (km2)%Area (km2)%
Dense Forests−2772.77−12.4108−1812.85−8.11
Agricultural Lands2518.09 11.27889.90 3.46
Barren and Urban Lands240.38 1.07307.364.59
Waterbodies14.290.0615.310.068
Table 6. Projected land use and land cover changes for the years 2030, 2040, 2050.
Table 6. Projected land use and land cover changes for the years 2030, 2040, 2050.
Class NameArea (km2) 2030Area (km2) 2040Area (km2) 2050Change from 2030 to 2040 (km2)Change from 2040 to 2050 (km2)
Waterbodies388.72387.16384.633−1.56−2.53
Barren and urban lands3830.353896.33948.2465.9551.24
Dense forests7461.27125.397030.62−335.91−94.77
Agricultural lands9788.39847.89894.760.646.9
Table 7. Land use transition matrix of the year 2030, 2040, and 2050.
Table 7. Land use transition matrix of the year 2030, 2040, and 2050.
2020–2030
WaterbodiesBarren and Urban LandsDense ForestsAgricultural Lands
Waterbodies0.900.060.030.01
Barren and urban lands0.000.940.010.05
Dense forests0.000.060.750.19
Agricultural lands0.000.120.050.83
2030–2040
WaterbodiesBarren and Urban LandsDense ForestsAgricultural Lands
Waterbodies0.960.000.000.04
Barren and urban lands0.000.980.010.01
Dense forests0.000.020.870.11
Agricultural lands0.000.030.060.91
2040–2050
WaterbodiesBarren and Urban LandsDense ForestsAgricultural Lands
Waterbodies0.970.000.010.02
Barren and urban lands0.000.950.010.04
Dense forests0.000.020.850.13
Agricultural lands0.000.070.050.88
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zeshan, M.T.; Mustafa, M.R.U.; Baig, M.F. Monitoring Land Use Changes and Their Future Prospects Using GIS and ANN-CA for Perak River Basin, Malaysia. Water 2021, 13, 2286. https://doi.org/10.3390/w13162286

AMA Style

Zeshan MT, Mustafa MRU, Baig MF. Monitoring Land Use Changes and Their Future Prospects Using GIS and ANN-CA for Perak River Basin, Malaysia. Water. 2021; 13(16):2286. https://doi.org/10.3390/w13162286

Chicago/Turabian Style

Zeshan, Muhammad Talha, Muhammad Raza Ul Mustafa, and Mohammed Feras Baig. 2021. "Monitoring Land Use Changes and Their Future Prospects Using GIS and ANN-CA for Perak River Basin, Malaysia" Water 13, no. 16: 2286. https://doi.org/10.3390/w13162286

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