Detecting and Analyzing Land Use and Land Cover Changes in the Region of Al-Jabal Al-Akhdar, Libya Using Time-Series Landsat Data from 1985 to 2017

The region of Al-Jabal Al-Akhdar in northeastern Libya has undergone rapid, wide-ranging changes in the land use and land cover (LULC) intensified by the conversion of natural resources for food purpose, urbanization, and other socioeconomic benefits. This study examined the use of geographic information system (GIS) and remote sensing techniques to gain a quantitative understanding of the spatiotemporal dynamics of LULC. In addition, the major factors behind LULC changes and decline of natural vegetation in the region were analyzed. A post-classification comparison approach was used to detect LULC changes in the study area between 1985 and 2017 using four Landsat images from 1985, 2000, 2010, and 2017. The observed changes were indicative of a decrease in the expanse of the natural Mediterranean forest which lost 9018 ha over 32 years, 39% of its total area, with the highest deforestation rate registered between 2010 and 2017 estimated at 513 ha. year−1. Orchards and rain-fed agriculture lands gained 4095 ha, which matches 55% of initial area, whereas the land under irrigated crops increased by 2266 ha, about 85% of the original area. The area of urban and built-up land in 2017 was more than double in 1985 and achieved the highest urbanization rate between 2010 and 2017 at 203 ha.year−1. Results indicate an unstable trend of bare and low vegetation lands which generally increased by about 50%. From the outcomes of this research, it is strongly recommended that urgent measures be taken to conserve the natural forest and to achieve a rational use of agricultural land in the region of Al-Jabal Al-Akhdar.


Introduction
Many countries in the world are faced with rapid, wide-ranging changes in land use (e.g., human activities and various utilization of the land) and land cover (e.g., physical characteristics of the land surface) [1]. Numerous studies have shown that only a few landscapes on earth, in peripheral sites and remote locations, are still remaining in their natural state [2]. The rapid increase in the world population, which necessitates an increase in anthropogenic activities, has resulted in speedy alteration in the land use and land cover (LULC), leading to forest deterioration and transformation of fertile land to urban construction with significant impact on the ecosystem [3]. LULC change has generated considerable interest as it focuses on the central issue in the field of global change study. Management of the earth's surface transformation due to LULC changes represents one of the pressing

Study Area
Libya lies along the southern side of the Mediterranean coast, approximately between the latitudes of 19 • 30 N and 33 • 10 N and the longitudes of 09 • 30 E and 25 • 00 E. Libya is mostly a desert country with a dry climate and a rare vegetation cover. However, the region of Al-Jabal Al-Akhdar in the northeastern part is an exception with a relatively high annual rainfall (250-650 mm) combining intense farming areas, natural forests, and shrubs. This region is the most favored part of the country from the climatic aspect with human settlements and land use dating back to the classical period and where agricultural activities are concentrated and expanding until today. As shown in Figure 1, the region is located in the northeastern part of Libya immediately south of the coastal belt extending along with it for about 300 km at the latitudes of 32 • 00 N to 32 • 56 N and the longitudes of 20 • 19 E to 23 • 08 E, encompassing approximately 1.149 million ha. It is a high plateau, largely rocky and stony, rising to about 882 m above sea level with an undulant surface gradually inclining to the south and dissected frequently by many deep ravines (rivers) streaming toward the east, west, and north through a steep escarpment toward the north to either the sea or a narrow coastal plain. This work was carried out in an area covering the middle northern part the of Al-Jabal Al-Akhdar region between the latitudes of 32°42'N and 32°56'N and the longitudes of 21°30′ and 21°50′ E ( Figure 1) and with a total area of about 36,860 ha. The study area included Al-Bayda, the capital city of the Jabal al Akhdar region, and three villages, Massah, Al-Hamamah, and Al-Haniya, in addition to the Lussaitah area that is very intensively used for agriculture. The rationale behind choosing this area as a research site emanates from the fact that urbanization and agricultural expansion have adversely affected the natural ecosystem, particularly after the Libyan revolution in 2011.

Data Collection
Satellite imagery and ancillary data were used for investigating the historical changes of LULC in the study area over 32 years from 1985 to 2017. Four Landsat satellite scenes for Path/Row 183/37 from three types of sensors were downloaded from the United States Geological Survey (USGS) website (https://earthexplorer.usgs.gov/). Table 1 shows the characteristics of the Landsat images used in this study. To assure the best comparability, we chose images that were captured between May and September (dry season in Libya) when the sky is usually clear, which enabled us to obtain cloud-free imagery with the best visibility.
The dates of imagery were defined based on events that had considerable impacts in expanding deforestation activities and LULC changes in Libya. The image of Landsat 5 Thematic Mapper (TM) This work was carried out in an area covering the middle northern part the of Al-Jabal Al-Akhdar region between the latitudes of 32 • 42'N and 32 • 56'N and the longitudes of 21 • 30 and 21 • 50 E ( Figure 1) and with a total area of about 36,860 ha. The study area included Al-Bayda, the capital city of the Jabal al Akhdar region, and three villages, Massah, Al-Hamamah, and Al-Haniya, in addition to the Lussaitah area that is very intensively used for agriculture. The rationale behind choosing this area as a research site emanates from the fact that urbanization and agricultural expansion have adversely affected the natural ecosystem, particularly after the Libyan revolution in 2011.

Data Collection
Satellite imagery and ancillary data were used for investigating the historical changes of LULC in the study area over 32 years from 1985 to 2017. Four Landsat satellite scenes for Path/Row 183/37 from three types of sensors were downloaded from the United States Geological Survey (USGS) website (https://earthexplorer.usgs.gov/). Table 1 shows the characteristics of the Landsat images used in this study. To assure the best comparability, we chose images that were captured between May and September (dry season in Libya) when the sky is usually clear, which enabled us to obtain cloud-free imagery with the best visibility.
The dates of imagery were defined based on events that had considerable impacts in expanding deforestation activities and LULC changes in Libya. The image of Landsat 5 Thematic Mapper (TM) Sustainability 2020, 12, 4490 5 of 24 captured in 1985 was used to represent the status before the first abolishment of the Ministry of Agriculture in 1986. Unfortunately, we could not get a usable image from the same year. The Landsat 7 Enhanced Thermal Mapper Plus (ETM+) image of 2000 was captured in May, about two months after the second abolishment of the Ministry of Agriculture, and thus was selected to picture the LULC status at the beginning of the second period. The popular uprising in Libya which started on 17 February 2011, was an important event, after which the whole country experienced dramatic changes in many political, economic, and social aspects that significantly contributed to LULC change; therefore, the Landsat 7 ETM+ image acquired on 5 September 2010, just five months before this event, was used as a reference for the beginning of the third period.
Ancillary data comprised a collection of Ground Control Points (GCP) and a combination of digital and paper maps covering the study area including a series of topographic maps at scale (1:50,000) and land use and land cover maps of 2005. In addition, aerial photographs covering a part of the study area taken on 30 May 1979, with a scale of 1:6000, were obtained and used for identification of areas subjected to LULC change. Ground truthing points were collected using a hand-held Geographical Positioning System (GPS) through field visits organized in 2016 during the same season in which satellite images were acquired. Points were separated into two groups, one for training samples used in the classification stage and another for test samples for classification accuracy assessment.

Image Processing and Analysis
Image data were treated using ENVI 5.3 and Arc-GIS 10.3 software packages. The data processing steps are given in Figure 2 and explained in the following section. , was an important event, after which the whole country experienced dramatic changes in many political, economic, and social aspects that significantly contributed to LULC change; therefore, the Landsat 7 ETM+ image acquired on 5 September 2010, just five months before this event, was used as a reference for the beginning of the third period. Ancillary data comprised a collection of Ground Control Points (GCP) and a combination of digital and paper maps covering the study area including a series of topographic maps at scale (1:50,000) and land use and land cover maps of 2005. In addition, aerial photographs covering a part of the study area taken on 30 May 1979, with a scale of 1:6000, were obtained and used for identification of areas subjected to LULC change. Ground truthing points were collected using a hand-held Geographical Positioning System (GPS) through field visits organized in 2016 during the same season in which satellite images were acquired. Points were separated into two groups, one for training samples used in the classification stage and another for test samples for classification accuracy assessment.

Image Processing and Analysis
Image data were treated using ENVI 5.3 and Arc-GIS 10.3 software packages. The data processing steps are given in Figure 2 and explained in the following section.   Comparative analysis using different Landsat data (i.e., TM, ETM+, and OLI (Operational Land Imager)) can be challenging due to the instrumental errors related to ache sensor, noise from several sources, and uncertainty in scale and geometric conditions [32]. Preprocessing satellite imagery before conducting image classification and change detection is very vital to minimize those errors and to build a more thorough association between the obtained data and biophysical features on the ground [33]. In this study, a number of sequential steps of data preprocessing was performed using Arc GIS 10.3 and ENVI 5.3 software packages. These steps comprised radiometric, atmospheric, and geometric corrections in addition to image gap-filling, sub-setting, and enhancement. Preprocessing also involved the selection of the appropriate band combinations to be used in image classification.
To reduce radiometric errors, images were calibrated using the radiometric correction tool in ENVI 5.3, where raw data from the sensors (DNs) were converted to top-of-atmosphere reflectance. For atmospheric correction, a dark object subtraction model was applied in the current study as it is relatively the simplest and most widely used empirical method for classification and change detection applications [34]. This technique operates on the assumption that radiances received at the satellite for dark objects (zero or small surface reflectance) within the image is attributed to atmospheric scattering (path radiance), and thus, the minimum radiance value from the entire image is subtracted from all pixel values [35].
The imageries were projected to the Universal Transverse Mercator (UTM) coordinate system, Datum LGD2006, zone 34 North and corrected for geometric errors from the sources. Image-to-map rectification method was applied to geometrically correct the Landsat 8 OLI/TIRS (Operational Land Imager and Thermal Infrared Sensor) satellite image dated 11 May 2017 using well-distributed ground control points (GCPs) taken from topographic maps at a scale of 1:50000. The other three Landsat images (1985 TM, 2000 TM, and 2010 ETM+) were then referenced by performing the image-to-image registration method using the corrected Landsat 8 OLI/TIRS 2017 image as a reference image. The first-degree polynomial method and nearest neighbor resampling technique were applied when performing the geometrical rectification and image registration.
Because of scan line corrector (SLC) failure (known as the SLC-Off problem) that occurred on May 31, 2003, for the Enhanced Thermal Mapper Plus (ETM+) sensor of Landsat 7, all Landsat 7 images collected from that date till 2013 have had line gaps appear as long lines on both sides of each image, resulting in a loss about 22% of data. However, it was found that the effect of those missing data was slight or passable in some applications, such as monitoring land use change and crop quality [36]. For this study, only the Landsat 7 image dated 05 September 2010 was affected by the SLC-Off error where line gaps covered a portion of the study area in the northwestern side. Therefore, the image was corrected in such a way that the missing data caused by the failure of the SLC were filled. A correction was made to the gap-filling based on a single image by applying the Landsat Gap Fill extension in ENVI 5.3 using the single file gap fill (triangulation) method. This procedure performs correlations between the pixels of the image based on geostatistical interpolation techniques, which use data from the same image where empty pixels are replaced by the mean of their neighbors in the nearest two lines. Figure 3 shows the scene from LSC-Off-affected image before and after correction.
Landsat 8 OLI and Landsat 7 ETM+ imagery were sharpened to increase the resolution from 30 m to 15 m, where the spatial information in the high-resolution panchromatic bands and the color information in the multispectral bands were merged using the nearest neighbor diffusion (NNDiffuse) pan sharpening technique in order to create a more valuable color image with high resolution that improves mapping and classification accuracy [37,38]. However, it was not possible to perform this on the Landsat 5 TM image captured in 1985 since the panchromatic band is not available. Therefore, data were resampled from 30 to 15 m using the nearest neighbor technique to ensure consistency with OLI and ETM+ data used in this study. Furthermore, a collection of enhancement techniques including histogram stretching and spatial filtering and creating band compositions were performed to improve Sustainability 2020, 12, 4490 7 of 24 the clarity and visualization of the Landsat images used in this study. All images were clipped to the study area using Arc GIS 10.3 to speed up the data processing. Landsat 8 OLI and Landsat 7 ETM+ imagery were sharpened to increase the resolution from 30 m to 15 m, where the spatial information in the high-resolution panchromatic bands and the color information in the multispectral bands were merged using the nearest neighbor diffusion (NNDiffuse) pan sharpening technique in order to create a more valuable color image with high resolution that improves mapping and classification accuracy [37,38]. However, it was not possible to perform this on the Landsat 5 TM image captured in 1985 since the panchromatic band is not available. Therefore, data were resampled from 30 to 15 m using the nearest neighbor technique to ensure consistency with OLI and ETM+ data used in this study. Furthermore, a collection of enhancement techniques including histogram stretching and spatial filtering and creating band compositions were performed to improve the clarity and visualization of the Landsat images used in this study. All images were clipped to the study area using Arc GIS 10.3 to speed up the data processing.

Image Classification
Image classification is an important remote sensing technique used to catalog all pixels in an image or raw remotely sensed data into a finite number of individual LULC classes to produce beneficial thematic maps and information [39]. The image classification process is normally conducted to assign different spectral signature from the dataset to several classes based on reflectance attributes of the diverse types of LULC.
For this study, a classification scheme was developed after Anderson et al. [40] based on prior knowledge about the study area, reconnaissance field survey, and information from previous studies. The field observations and available map of land cover/vegetation in the northeastern zone of Libya, prepared in 2005 by FAO and The United Nations Development Programme (UNDP), were the basic references for identification of the prevalent LULC classes. Accordingly, five LULC classes, including the natural Mediterranean forest (NMF), orchards and rain-fed agriculture (OR), irrigated crops (IC), urban and built-up lands (UB), and bare and low vegetation lands (BLV), were defined to form the adopted general LULC classification scheme for this study as described in Table 2. Table 2. Definition of land use and land cover (LULC) classes of classification scheme adopted in this study.

LULC Code Description
Natural Mediterranean Forest NMF Areas covered by a dense evergreen natural woody vegetation known as Maquis which characterizes the Mediterranean regions, with cover reaching 70%-80%, generally composed of shrubs, bushes, and trees (usually less than 6 m)

Image Classification
Image classification is an important remote sensing technique used to catalog all pixels in an image or raw remotely sensed data into a finite number of individual LULC classes to produce beneficial thematic maps and information [39]. The image classification process is normally conducted to assign different spectral signature from the dataset to several classes based on reflectance attributes of the diverse types of LULC.
For this study, a classification scheme was developed after Anderson et al. [40] based on prior knowledge about the study area, reconnaissance field survey, and information from previous studies. The field observations and available map of land cover/vegetation in the northeastern zone of Libya, prepared in 2005 by FAO and The United Nations Development Programme (UNDP), were the basic references for identification of the prevalent LULC classes. Accordingly, five LULC classes, including the natural Mediterranean forest (NMF), orchards and rain-fed agriculture (OR), irrigated crops (IC), urban and built-up lands (UB), and bare and low vegetation lands (BLV), were defined to form the adopted general LULC classification scheme for this study as described in Table 2.
There are several different methods of imagery data classification; however, unsupervised and supervised techniques are the two main approaches commonly used for mapping land cover [41]. Most researchers prefer using the supervised technique because it produces more accurate classification than the unsupervised technique [42][43][44]. There are several statistics-based supervised classifiers. However, the maximum likelihood, minimum distance, parallopiped, and Mahalanobis are the most popular classifiers [45]. While reviewing the literature, it transpired that the maximum likelihood classifier is widely used in remote sensing research as it is quick, is simple to implement, allows a clear interpretation of the outcomes, and mostly delivers a satisfactory accuracy [46][47][48]. The technique presumes that the statistics for each LULC category in each band are subject to the normal distribution and computes the likelihood for each pixel to belong to a particular category [4]. In this work, supervised classification was applied for LULC classification. The process began with the most important step, which is selecting more than 200 training samples in each Landsat image by delineating polygons around representative areas of each LULC class. True and false composites were used to improve the visualization of different features to easily delineate LULC classes in the imagery. Training sites for the predefined LULC classes ( Table 2) were developed by the visual interpretation of the Landsat images (1985,2000,2010, and 2017) based on ground truthing information, Google Earth imagery, familiarity with the study area, field observations, and historical information derived from interviews conducted with local landholders. Once the training sites were developed, the maximum likelihood classifier was applied to create spectral signatures and then used to classify all pixels in the image. After the initial classification, a majority filter was used to smoothen classification results by removing "salt and pepper" noise from the classified maps, and thus, final LULC maps were produced and used for further analyses. Table 2. Definition of land use and land cover (LULC) classes of classification scheme adopted in this study.

LULC Code Description
Natural Mediterranean Forest NMF Areas covered by a dense evergreen natural woody vegetation known as Maquis which characterizes the Mediterranean regions, with cover reaching 70-80%, generally composed of shrubs, bushes, and trees (usually less than 6 m) Orchards and Rain-fed Agriculture OR Agricultural lands that totally depend on rainfall for growing cereal crops and some fruit crops Irrigated Crops IC Agricultural lands that depend on groundwater for irrigation, often cultivated as small fields for production of different vegetable crops Urban and Built-up Lands UB Areas of intensive use where most land is covered with structures including residential, commercial, and industrial areas; buildings and farmhouses spread in the countryside are involved in this class Barren and Low Vegetation BLV Areas with low or no vegetation including rocky areas with thin soil layer on steep slopes or in areas of severe erosion, exposed soil resulting from overgrazing or poor cultivation practices, or construction activities and areas that have been stripped of vegetation due to deliberate or accidental fires

Classification Accuracy Assessment
The accuracy assessment in this study was done using historical high-resolution google imagery obtained from Google Earth Pro 7.3.1.4507 (Google Inc., Mountain View, CA, USA). Two google images captured in October 2010 and July 2017 were used to carry out the accuracy validation for LULC maps of 2010 and 2017, respectively. The random sampling technique was used to identify verification points based on the density and distribution of each LULC class which was performed in ArcGIS 10.3 using the "create random points" tool. A total number of 280 verification points were randomly selected with a minimum sample number of at least 50 for each class [49]. However, no accuracy assessment was conducted for the land use map of 1985 and 2000 due to the dearth of adequate reference data, and it was assumed that accuracy validation for classification maps of 2010 and 2017 should be adequate to give information on the general accuracy of the adopted procedures in this research for mapping LULC.
Classification accuracy was assessed using measures derived from the error matrices, and these included overall accuracy and kappa coefficient. In addition, user's and producer's accuracies were also obtained for each class, which respectively measure the commission and omission errors. The value of overall accuracy greater than 70% is regarded as an acceptable classification accuracy [50]. For most applications, values of Kappa statistic more than 0.75 can be accepted as an excellent or very good agreement while values in the range of 0.40 to 0.75 are indicative of fair or good correspondence and values less than 0.4 represent poor correspondence [51].

Change Detection
Change detection is defined as "the process of identifying differences in the state of an object or phenomenon by observing it at different times" [52]. Change detection analysis makes LULC study a subject of interest because it not only explore changes that have taken place but also identifies their nature and determines their spatial extent and pattern [53,54]. There are several methods that have been developed for change detection analysis of remote sensing data, and they have been classified into Sustainability 2020, 12, 4490 9 of 24 two major categories: pre-classification and post-classification comparison techniques [55]. However, the post-classification comparison approach has been proven to be the most commonly used method for LULC change detection [56].
The post-classification approach involves the classification of multiple date images separately to generate thematic maps, after which a pixel-based comparison of the corresponding classes is used to produce tables and maps of changes that have occurred [57]. There are several advantages to using post-classification comparison approach, and this includes the following: 1) it minimizes the possible effects of sensor and atmospheric and environmental differences between images as they are independently classified, 2) it provides more detailed information obtained from a complete matrix of LULC change, and 3) it quantifies the change rates and magnitude of change [58]. However, the degree of change detection success depends on the accuracy of individual classifications [59]. Many studies have employed post-classification techniques based on maximum likelihood classification for detecting LULC changes in different areas around the world using Landsat imagery data and have attained high accuracy [47,[60][61][62][63]. The capability of the post-classification comparison techniques is not limited to just identifying the size and distributions of areas that have changed, but it goes beyond that to determine the individual participation rate of each LULC class in changing another particular class. Therefore, the post-classification comparison technique was chosen in this study to be employed for change detection. Pairs of the produced classified thematic maps for the four dates (1985,2000,2010, and 2017) were compared by applying the detection algorithm using ENVI 5.3. Change maps and matrices were produced to understand the change status and its magnitude rate for the periods of 1985 to 2000, 2000 to 2010, and 2010 to 2017.
The magnitude of change (MC), the percentage of change (PC), and the annual rate of change (ARC) for each LULC class during each period of time were computed based on the following equations: where A i is the class area (ha) at the initial time, A f is the class area (ha) at the final time, and n is the number of years of the time period.

Accuracy Assessment
The use of maximum likelihood classifier to categorize the Landsat images (TM, ETM+, and OLI/TIRS) has produced maps showing the distribution of the five prevalent LULC classes in the study area for the years 1985, 2000, 2011, and 2017. Accuracy assessment was carried out for the classified maps of 2010 and 2017 to examine the agreement between the produced maps and what actually exists on the ground. Table 3 illustrates the summarized results of classification accuracy. The classified maps of 2010 and 2017 achieved overall accuracy values of 81% and 83%, respectively, indicating that both LULC maps met the criterion. Furthermore, values of the kappa coefficient for both the two maps were more than 0.75, indicating good agreements. Results also showed that means of commission and omission errors were less than 20%, inferring an acceptable classification accuracy.
The LULC classification of Landsat ETM+ 2010 and OLI 2017 recorded the largest commission and omission error within the BLV, IC, and OR classes (Figure 4). Commission errors varied among these classes between 22% and 34%, while omission errors varied from 20% to 27%. These classes were confused with each other and with other classes. These errors might be due to the spatial pattern of vegetation growth which might lead to spectral confusion. For example, the fields of IC class in the study area often are small areas adjacent to lands of other land uses, especially OR. Also, leaving the agricultural lands (IC and OR) exposed without any vegetation cover after harvesting and before the next planting season might be a source of spectral confusion with BLV class. However, there are some other sources of errors, such as mixed pixels and atmospheric reflection effects, that may have resulted in confusion between LULC classes. Generally, the accuracy achieved by the maximum likelihood classification for the different LULC in this study is in agreement with previous studies that have worked on LULC classification in semiarid areas and is hence adequate for further investigation.

Year
Overall Accuracy (%) The classified maps of 2010 and 2017 achieved overall accuracy values of 81% and 83%, respectively, indicating that both LULC maps met the criterion. Furthermore, values of the kappa coefficient for both the two maps were more than 0.75, indicating good agreements. Results also showed that means of commission and omission errors were less than 20%, inferring an acceptable classification accuracy.
The LULC classification of Landsat ETM+ 2010 and OLI 2017 recorded the largest commission and omission error within the BLV, IC, and OR classes (Figure 4). Commission errors varied among these classes between 22% and 34%, while omission errors varied from 20% to 27%. These classes were confused with each other and with other classes. These errors might be due to the spatial pattern of vegetation growth which might lead to spectral confusion. For example, the fields of IC class in the study area often are small areas adjacent to lands of other land uses, especially OR. Also, leaving the agricultural lands (IC and OR) exposed without any vegetation cover after harvesting and before the next planting season might be a source of spectral confusion with BLV class. However, there are some other sources of errors, such as mixed pixels and atmospheric reflection effects, that may have resulted in confusion between LULC classes. Generally, the accuracy achieved by the maximum likelihood classification for the different LULC in this study is in agreement with previous studies that have worked on LULC classification in semiarid areas and is hence adequate for further investigation.

Supervised Classification
The processing and supervised classification of Landsat imagery (TM, ETM+, and OLI/TIRS) acquired in 1985, 2000, 2010, and 2017 produced four classification maps for the five LULC types prevalent in the study area ( Figure 5). Classification results showed that the areas of LULC classes in the four produced maps (1985, 2000, 2010, and 2017) are in descending order as follows: NMF > OR > IC > UB > BLV ( Figure 6). NMF is the dominating class, especially in the northern latitudes of the study area, covering more than half of the total landscape in 1985 and 2000 but decreased to less than 40% in 2017. The OR class occupied an area ranging from 20% to about 32% of the total area over the study time period, while the maximum area of each of the other classes represents less than 15%. With regard to agricultural uses, lands of the OR class are distributed over different parts within the study area, whereas lands under the IC class are concentrated in the middle latitudes of the study area (Lussaitah area) and in the narrow coastal strip in the north. These areas are considered as the most important agricultural areas in the region of Al-Jabal Al-Akhdar due to the sufficient amount of precipitation and soil fertility in addition to the availability of groundwater. On the other hand, areas of the UB class are concentrated in the southeastern part, where the city of Al-Bayda lies.

Supervised Classification
The processing and supervised classification of Landsat imagery (TM, ETM+, and OLI/TIRS) acquired in 1985, 2000, 2010, and 2017 produced four classification maps for the five LULC types prevalent in the study area ( Figure 5). Classification results showed that the areas of LULC classes in the four produced maps (1985, 2000, 2010, and 2017) are in descending order as follows: NMF > OR > IC > UB > BLV ( Figure 6). NMF is the dominating class, especially in the northern latitudes of the study area, covering more than half of the total landscape in 1985 and 2000 but decreased to less than 40% in 2017. The OR class occupied an area ranging from 20% to about 32% of the total area over the study time period, while the maximum area of each of the other classes represents less than 15%. With regard to agricultural uses, lands of the OR class are distributed over different parts within the study area, whereas lands under the IC class are concentrated in the middle latitudes of the study area (Lussaitah area) and in the narrow coastal strip in the north. These areas are considered as the most important agricultural areas in the region of Al-Jabal Al-Akhdar due to the sufficient amount of precipitation and soil fertility in addition to the availability of groundwater. On the other hand, areas of the UB class are concentrated in the southeastern part, where the city of Al-Bayda lies.

Supervised Classification
The processing and supervised classification of Landsat imagery (TM, ETM+, and OLI/TIRS) acquired in 1985, 2000, 2010, and 2017 produced four classification maps for the five LULC types prevalent in the study area ( Figure 5). Classification results showed that the areas of LULC classes in the four produced maps (1985, 2000, 2010, and 2017) are in descending order as follows: NMF > OR > IC > UB > BLV ( Figure 6). NMF is the dominating class, especially in the northern latitudes of the study area, covering more than half of the total landscape in 1985 and 2000 but decreased to less than 40% in 2017. The OR class occupied an area ranging from 20% to about 32% of the total area over the study time period, while the maximum area of each of the other classes represents less than 15%. With regard to agricultural uses, lands of the OR class are distributed over different parts within the study area, whereas lands under the IC class are concentrated in the middle latitudes of the study area (Lussaitah area) and in the narrow coastal strip in the north. These areas are considered as the most important agricultural areas in the region of Al-Jabal Al-Akhdar due to the sufficient amount of precipitation and soil fertility in addition to the availability of groundwater. On the other hand, areas of the UB class are concentrated in the southeastern part, where the city of Al-Bayda lies.

LULC Change Detection
The outcomes of the post-classification comparison give an account of the amount and type of change that has taken place for each LULC class. Areal change of LULC types in all imageries was obviously visible and occurred as either an expansion or a decline between sequential dates. Attempts to further study these outcomes can give insight to the land use behavior, pattern, and rate at which the changes occur in the study area in particular and in the region of Al-Jabal Al-Akhdar in general. This essential information would be vital for development planners and decision makers.
Change In 1985, the NMF class constituted the major type of LULC in the study area at 64.17% of the total area, followed by OR (20.62%), IR (7.44%), and UB (4.79%), whereas BLV only covered 3% and occupied the smallest area. During the 15-year time frame from 1985 to 2000, the area covered by NMF declined by 2724 ha, viz. 11.79% of the initial area, while other LULC types increased at different rates. Table 4 and Figure 7 summarize the major LULC conversion that occurred within the study area between 1985 and 2000.  It is important to remember that, in 1986, under the previous political regime, the Ministry of Agriculture in Libya was abolished with all its departments and various control and enforcement agencies at the national level including the forest guard; thus, the agricultural sectors in all regions regrouped under the management of the sector directors, and this encouraged forests encroachment and failure of protection efforts. The weakness of institutions and absence of accountability among other factors encouraged the local people to engage in forest encroachment by either clearance of vegetation cover for agricultural expansion or by cutting trees for gathering firewood, especially for making the popular Libyan bread, in addition to producing charcoal for heating, barbecue, and cooking. In addition, some landowners have tried exploiting forest land in an attempt to regain their lands that were confiscated and turned into state property by law (No. 142 of 1970). Although the Ministry of Agriculture was reinstated in 1990 and a new forest law was enacted in 1992 to correct the accesses committed, the disaster of forest loss had already taken place and its negative effects As shown in Table 4, based on 36,010 ha, the total size of the study area, NMF covered 23,109 ha in 1985 (64.17% of the total area), but this changed to 20,385 ha (56.61%) in the year 2000, with a negative annual rate of change of 182 ha.year −1 (0.79%). As shown in Figure 7, out of the 23 It is important to remember that, in 1986, under the previous political regime, the Ministry of Agriculture in Libya was abolished with all its departments and various control and enforcement agencies at the national level including the forest guard; thus, the agricultural sectors in all regions regrouped under the management of the sector directors, and this encouraged forests encroachment and failure of protection efforts. The weakness of institutions and absence of accountability among other factors encouraged the local people to engage in forest encroachment by either clearance of vegetation cover for agricultural expansion or by cutting trees for gathering firewood, especially for making the popular Libyan bread, in addition to producing charcoal for heating, barbecue, and cooking. In addition, some landowners have tried exploiting forest land in an attempt to regain their lands that were confiscated and turned into state property by law (No. 142 of 1970). Although the Ministry of Agriculture was reinstated in 1990 and a new forest law was enacted in 1992 to correct the accesses committed, the disaster of forest loss had already taken place and its negative effects were visible and tangible.
On the other hand, forests in Libya, like others in the Mediterranean region, which is one of the most affected regions by forest fires, have suffered excessive degradation by fires which are mainly due to manmade activities [64]. Destructive forest fires in Al-Jabal Al-Akhdar have always been present since the last few decades, often attributed to human activities and sometimes to natural factors. People in the region use fires to remove the natural vegetation, either to exploit the lands for agricultural use, or to induce annual plants to grow intensively to improve site characteristics as pasturelands for grazing their livestock. According to official reports from the Ministry of Justice and Public Security and the Ministry of Agriculture cited in OMU [19], forest fires devoured as much as 33,230 ha of the natural Mediterranean forest in the region of Al-Jabal Al-Akhdar in the period between 1985 and 2003, of which 4000 hectares were located in the focus area of this study (in Massah, Lussaitah and Al-Haniya).
The area of agricultural land in 2000, in comparison to 1985, increased by 2050 ha, where the area of OR increased from 7426 ha (20.62%) to 8659 ha (24.05%) with an annual rate of 82 ha.year −1 , while the area of IC increased from 2679 ha (7.44%) to 3496 (9.71%) with an annual rate of 55 ha.year −1 ( Table 4). The results shown in Figure 7 illustrate that about 5382 ha (72.48%) of OR land in 1985 did not change in 2000 whereas about 2044 ha (27.52%) were converted to other land uses, especially for IC, which acquired an area of 1202 ha from OR (58.81% of OR loss). However, OR gained new areas in 2000, estimated at 3277 ha. Much of OR gains were obtained from the clearance of NMF and conversion of IC land, which account for about 55.36% and 37.56% of OR gains, respectively. Clearance of NMF also contributed 900 ha to IC (42.08% of IC gain).
Libya has given the agricultural sector the largest share of the allocation of its transformation agenda through developmental projects during the 1970s and 1980s such as vertical and horizontal expansion in crop and fodder production. The region of Al-Jabal Al-Akhdar received a large share of these allocations that were used to reclaim large tracts of land estimated at about 160,000 ha. New farms were established in the Al-Jabal Al-Akhdar region to produce fruit and grain crops to attain self-sufficiency in food production; however, substantial areas of natural cover have been left without change for conservation purposes or to be used for natural grazing, particularly on hilly sites. The local people later exploited these marginal unused areas for establishing more private fields to produce cereal crops as well as fruits and irrigated vegetables in order to capitalize on the increasingly lucrative market.
The results shown in Table 4 indicate that the urban expansion in this period was not very large in the study area. The areas of UB increased by only 295 ha, with an annual rate of 20 ha.year −1 , and most of the UB gains were from the conversion of OR land, which represented 73.71% of the UB gains (Figure 7). This slight rate of increase in urban expansion may be attributed to several factors, one of them is the decrease in the annual rate of population growth in the Al-Jabal Al-Akhdar region and the whole country in general. As shown in Figure 8, the annual rate of population growth in the Al-Jabal Al-Akhdar region was 4.53% during the period 1964-1973 and dropped to 3.92% during the period 1973-1984 and then to 2.35% during the period 1984-1995. Another factor that strongly limited the urban expansion is the suspension of development programs and housing projects in the whole country due to many reasons such as the global economic depression of the mid-80s, the decline in oil prices and sales, and the trade sanctions against Libya during the 90s.
produce cereal crops as well as fruits and irrigated vegetables in order to capitalize on the increasingly lucrative market.
The results shown in Table 4 indicate that the urban expansion in this period was not very large in the study area. The areas of UB increased by only 295 ha, with an annual rate of 20 ha.year −1 , and most of the UB gains were from the conversion of OR land, which represented 73.71% of the UB gains (Figure 7). This slight rate of increase in urban expansion may be attributed to several factors, one of them is the decrease in the annual rate of population growth in the Al-Jabal Al-Akhdar region and the whole country in general. As shown in Figure 8, the annual rate of population growth in the Al-Jabal Al-Akhdar region was 4.53% during the period 1964-1973 and dropped to 3.92% during the period 1973-1984 and then to 2.35% during the period 1984-1995. Another factor that strongly limited the urban expansion is the suspension of development programs and housing projects in the whole country due to many reasons such as the global economic depression of the mid-80s, the decline in oil prices and sales, and the trade sanctions against Libya during the 90s.  (Table 4). As shown in Figure 7, the clearance of NMF land and abandonment of agricultural lands, particularly OR land, were the main factors that led to the increase in BLV area which represented 48.58% and 43.41% of BLV gains, respectively, by a combined area of about 873 ha. However, forest restoration occurred on about 303 ha (58.83% of the BLV loss) and about 208 ha (40.39% of the BLV loss) were converted to OR.

LULC Change between 2000 and 2010
During the 10 years between 2000 and 2010, more forest land was cleared or changed to other land uses. NMF areas, which amounted to 20,385ha (56.61%) in 2000, decreased to 17,685 ha (49.11%) in 2010, with an annual rate of 270 ha.year −1 (1.32%) ( Table 5). From the result in Figure 9, it is apparent that, out of 20,385 ha of NMF in 2000, about 17,035 ha (83.57%) of NMF did not change in the subsequent ten years whereas 3350 ha (16.43%) were cleared or converted to other land uses. Most of the lost area of NMF changed to agricultural land uses, where about 1547 ha (46.18% of NMF loss)  (Table 4). As shown in Figure 7, the clearance of NMF land and abandonment of agricultural lands, particularly OR land, were the main factors that led to the increase in BLV area which represented 48.58% and 43.41% of BLV gains, respectively, by a combined area of about 873 ha. However, forest restoration occurred on about 303 ha (58.83% of the BLV loss) and about 208 ha (40.39% of the BLV loss) were converted to OR.

LULC Change between 2000 and 2010
During the 10 years between 2000 and 2010, more forest land was cleared or changed to other land uses. NMF areas, which amounted to 20,385ha (56.61%) in 2000, decreased to 17,685 ha (49.11%) in 2010, with an annual rate of 270 ha.year −1 (1.32%) ( Table 5). From the result in Figure 9, it is apparent that, out of 20    It was argued that the administrative instability and the permanent state of change in the Ministry of Agriculture, which was reabolished in March 2000, and legal inconsistencies in the way of allocating the liabilities have contributed significantly to rising encroachments and violations on NMF land within the Al-Jabal Al-Akhdar region. In 2006, the ministerial resolution (No. 13; 2006) was enacted through which an agricultural policing body was created for the purpose of It was argued that the administrative instability and the permanent state of change in the Ministry of Agriculture, which was reabolished in March 2000, and legal inconsistencies in the way of allocating the liabilities have contributed significantly to rising encroachments and violations on NMF land within the Al-Jabal Al-Akhdar region. In 2006, the ministerial resolution (No. 13; 2006) was enacted through which an agricultural policing body was created for the purpose of monitoring and limiting encroachment on forest land. However, the results shown in Figures 7 and 9 revealed that the area of deforested land, which has been converted into other land uses, increased from 3238 ha (216.87 ha.year −1 ) through the period 1985-2000 to 3350 ha (335 ha.year −1 ) through the period 2000-2010, which suggests that those employed for the purpose of forest protection were incapable of executing their work properly. Such incapability could be due to the lack of financial and technical support. Moreover, according to many local residents, deforestation was often carried out within private property in the cover of the night where agricultural police control is absent.
During the second period, each of the OR and IC classes increased by 1214 ha, with an annual rate of 121 ha.year −1 ( Table 5). The results in Figure 9  Population growth and the consequent increase in food demand encouraged farmers in the study area to extend their cultivated areas through the clearance of more forests in order to benefit from the increased profits by marketing the agricultural produce in neighboring towns and cities as well as in the remote areas in the west of the country. Drilling more underground wells by the local population in the study area, either at their own expense or through state support in form of agricultural and pastoral loans, has greatly augmented the rate of deforestation and land conversion for irrigated agriculture.
Classification results shown in Table 5 reveal that the annual rate of urban expansion has doubled during the period 2000-2010 compared to the period 1985-2000, where it increased from 20 ha.year −1 in the first period to 41 ha.year −1 in the second period. The area of UB increased from 2020 ha (5.61%) in 2000 to 2430 ha (6.75%) in 2010. Figure 9 shows that most of the UB gains were obtained from converting OR land which represented about 48.43% of the UB gains, followed by clearance of NMF land, which represented about 32.49% of the UB gains.
The increase in urban expansion might be attributed to population growth within the region of Al-Jabal Al-Akhdar through the period 1995-2006 (Figure 8), which required constructing more housing units to accommodate the growing numbers of families. On the other hand, the period 2000-2010 experienced a recovery in the country's economic situation. The housing crisis was beginning to ease, especially after lifting economic sanctions imposed against Libya in 2003 after which oil incomes dramatically increased. As a result of this positive event, significant economic and political changes were made in the country. Thus, many residential projects were constructed and more housing loans were granted in an attempt to solve the housing crisis. This progress was accompanied by the expansion of construction of networks of paved and unpaved roads across forests over mountains and valleys of the Al-Jabal Al-Akhdar region as in all areas of the country. However, although urban expansion and road construction are extremely crucial for national development with the growing population, they are critical factors that significantly affect the rate of deforestation and decline of agricultural land as they facilitate the access to more land for forests clearance, construction, and exploitation for different purposes [66].
The area of BLV in the period between 2000 and 2010, unlike in the previous period (1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000), decreased by 138 ha with an annual rate of 9.52% (Table 5). Although the BLV in the second period gained about 845 ha, mostly from the clearance of NMF land and abandonment of OR land, which represent 50.94% and 38.76% of the BLV gains, respectively, there is about 992 ha of BLV that was converted in 2010 to other LULC types. As seen in Figure 9, during the period 2000-2010, an area of 504 ha of BLV land was exploited for OR land use which represents 50.81% of the BLV loss. NMF cover was restored in some areas (about 372 ha) which were classified as BLV in 2000, representing 37.50% of the BLV loss. The data in Figures 7 and 9 indicate that the rate of exploitation of BLV land by conversion into OR increased in the period 2000-2010 compared to the past (1985-2000), where the first period witnessed a conversion of 208 ha (about 14 ha. year −1 ) of BLV land to be used for OR land use while in the second period more BLV land, estimated at 504 ha (about 51 ha. year −1 ), was converted into OR. Moreover, the rate of change from NMF to BLV increased from about 34.16 ha. year −1 through the first period to about 50 ha. year −1 through the second period, indicating the inability of the agricultural police to effectively perform its function to protect forests and to limit clearance activities. The anthropogenic violations of any kind against these forests, at whatever size, definitely represent a real threat to the natural forests in the Al-Jabal Al-Akhdar region.

LULC Change between 2010 and 2017
The results shown in Table 6 illustrate that, during this period, more changes occurred in the areas of LULC types within the study area. After the Libyan uprising in 2011, the natural vegetation cover in the study area was subjected to massive deforestation in the forms of forest cutting, frequent fires, and overgrazing. In general, the natural forest in Al-Jabal Al-Akhdar continues to be deteriorated in terms of quality and density in the absence of proper management and silvicultural activities [67]. According to the supervised classification results shown in Table 6 (Tables 4 and 5). The results in Figure 10   Comparing the results of this study to those obtained from Achour et al. [68] in the northern part of Tunisia to evaluate forest cover 4 years before and after the 2011 Tunisian popular uprising, it was found that the current deforestation rate in our study area is much higher than that reported in Tunisia. As shown in Table 6, the annual deforestation rate in our study area that was covered by natural forest on about 17,685 ha of its area before the Libyan uprising in 2011 was estimated at 513 ha.year −1 between 2010 and 2017, representing a loss of 2.9% of forest area per year. The annual rate of forest loss in the north of Tunisia that was covered by forest in 72,000 ha before the revolution in 2011 was estimated at about 2245 ha.year −1 during the period 2010-2014, representing a loss of 0.31% of forest area per year [68]. However, the deforestation rate in our study area appears to be lower compared to that documented in the Algarabulli district in the northwest part of Libya [16], where Comparing the results of this study to those obtained from Achour et al. [68] in the northern part of Tunisia to evaluate forest cover 4 years before and after the 2011 Tunisian popular uprising, it was found that the current deforestation rate in our study area is much higher than that reported in Tunisia. As shown in Table 6, the annual deforestation rate in our study area that was covered by natural forest on about 17,685 ha of its area before the Libyan uprising in 2011 was estimated at 513 ha.year −1 between 2010 and 2017, representing a loss of 2.9% of forest area per year. The annual rate of forest loss in the north of Tunisia that was covered by forest in 72,000 ha before the revolution in 2011 was estimated at about 2245 ha.year −1 during the period 2010-2014, representing a loss of 0.31% of forest area per year [68]. However, the deforestation rate in our study area appears to be lower compared to that documented in the Algarabulli district in the northwest part of Libya [16], where 3883 ha of forest was lost through the period 2010-2013, with an annual rate of decline estimated at 1295 ha.year −1 (16.2% per year). This difference can be attributed to the fact that most of the forests cover in the Algarabulli district are planted forests located in accessible sites near residential communities with a high population unlike the natural forest of Al-Jabal Al-Akhdar that is largely situated in hard-to-reach areas such as steep slopes and valley bottoms. Being located away from main roads and residential communities gives the forest greater protection against human encroachments. However, deforestation rate in the study area still remains at a significant rate, given the fact that the country is very lightly forested, with forest only covering about 217,000 ha [7] approximately representing 0.1% of the country's total area and mostly situated in the region of Al-Jabal Al-Akhdar, which comprises 90% of the Libyan natural forests. Indeed, if the present trend of deforestation continues, more than 50% of the natural forest cover in the study area will be lost over the next 15 years. Therefore, deforestation in Al-Jabal Al-Akhdar and Libya in general presents a huge environmental challenge since the governments after the revolution were unable to act effectively to protect country's natural resources including natural forests.
It should be recalled that the deforestation during the third period was induced by some factors entirely unlike those during the previous periods from 1985 to 2010. Presently, increased deforestation might be attributed to a range of factors including the total breakdown of the Al-Qaddafi regime in 2011 and the accompanying chaos and spread of weapons, the inability of subsequent governments to prevent forest clearance, the lack of housing projects, and the significant increase in land price. Some citizens exploited the fall of the government regime, the cessation of regulatory agencies, and the complete absence of agricultural policing to clear more forest land for agricultural expansion or to sell cleared land building houses or construction of holiday homes and resorts.
The supervised classification analysis indicated an increase in the area of agricultural land in the study area during the third period (2010-2017), which was estimated at 1883 ha, where the area of agricultural land use increased from 14,583 ha (40.50% of total area) in 2010 to 15,466 ha (44.72% of total area) in 2017 ( Table 6). The expansion of agriculture due to increased demand for food is one of the major factors that contributed to the reduction of the natural vegetation cover area, where more than 70% of deforested areas in the study area between 2010 and 2017 was used for cultivating irrigated and rainfed crops. This is consistent with the report by Kissinger et al. [69] which identified agriculture as being responsible for roughly 80% of worldwide deforestation. The report rated commercial agriculture and subsistence agriculture as the main direct drivers of deforestation in Africa and both account for one-third of all deforestation, while in Latin America, two-thirds of forest degradation is a result of the expansion of commercial agriculture. Although agriculture is important in the development of any nation, it has some undesirable impacts on natural environments. As the area of agricultural land grows and increases, the area of natural vegetation declines and the soil is exposed to erosion and deterioration of its physical and chemical properties.
The results shown in Table 6 illustrate that the area of OR land increased from 9873 ha (27.42%) in 2010 to 11,521 ha (about 32%) in 2017 by an annual rate of 235 ha.year −1 (2.38%), whereas the area of IC increased in the same period from 4710 ha to 4954 ha, registering a lower annual rate of change estimated at 34 ha.year −1 . Data presented in Figure 10 show that about 7254 ha (73.47%) of OR land in 2010 did not change in the ensuing 7 years, whereas about 2619 ha (26.53%) was changed to other land uses. About 51.78% of OR loss was converted to IC land use, while 28.56% was altered to UB and 14.20% was altered to BLV. However, during the 7 years between 2010 and 2017, OR land use gained about 4267 ha, representing 37.04% of its area in 2017. Most gains of OR were obtained from converting IC and NMF land, which represented 48% and 40.83% of OR gains, respectively. The change matrix in Figure 10 also shows that about 51.50% of IC land in 2010 did not change in 2017 whereas the remaining was converted into other agricultural and even nonagricultural land uses. Most the lost area of IC land was converted to OR, which represents 89.67% of IC net loss. IC land in the third period gained 2519 ha mostly from the conversion of OR land, which contributed about 1356 ha (53.83% of IC gains), in addition to clearance of about 1122 ha NMF land (44.54% of IC gains). The gained area of IC in this period was only 10% greater than the area was lost (2284 ha), which is different from the first and second periods when IC gains were 62% and 70% higher than losses, respectively. From Figure 10, it can be also seen that the conversion from IC to OR during the recent period (2010-2017) was greater than conversion from OR to IC by about 50%, unlike the previous periods when both conversions were similar to some extent. The annual rate of increase in the area of IC land use, during the period from 2010 to 2017, which was estimated at 34 ha.year −1 , considers the lowest rate of this land use when compared to previous periods, where the annual rate of increase reached 55 ha.year −1 during the period 1985-2000 and increased to 121 ha.year −1 during the period 2000-2010 (Tables 4-6). Several factors have contributed to reducing the annual rate of increase in the total area of IC land in the study area during the period 2010-2017. The high increase in agricultural input cost, such as fertilizers, pesticides, seeds, and equipment in addition to increasing labor wages and shrinking of skilled workforce, which often is nonlocal, were the main factors that have driven farmers to abandon the cultivation of irrigated crops and to turn to other activities.
During the third period, UB extended from 6.75% to 10.70% of the total land use and had an annual rate of increase of 8.36%, which was much higher compared to the previous periods. The UB area increased from 2430 ha in 2010 to 3852 ha in 2017, achieving a total increase of 1422 ha and representing about 58.52% of the initial area. The annual rate of increase during this period which was estimated at 203 ha.year −1 ( Table 6) is about 10 times the annual rate recorded between 1985 and 2000 as 20 ha.year −1 ( Table 4) and about 5 times the rate achieved during the second period from 2000 to 2010 and estimated at 41 ha.year −1 ( Table 5).
Results in Table 6 shows that the urban expansion at the expense of agricultural land during the third period saw a significant increase comparing the first and second periods, particularly around the city of Al-Bayda. An area of about 845 ha of land which was previously classified as agricultural land (OR + IC) was converted during the 7 years between 2010 and 2017 into UB, which approximate twice the converted agricultural area in 25 years between 1985 and 2010. It also can be deduced that conversion from NMF into UB during the third period (450 ha) equals twice the converted area in the first and second periods combined (218 ha). Additionally, exploitation of BLV for urban expansion was doubled during the third period compared to the previous periods, where about 206 ha of BLV was converted into UB between 2010 and 2017 compared to the 107 ha converted between 1985 and 2000.
Several factors were responsible for the increase in the urban expansion in the Al-Jabal Al-Akhdar region. The most important factors are administrative causes that are related to the state policy, the weak application of legislation, and urban laws relating to the protection of agricultural and forest lands. Furthermore, population concentration in cities is not in tandem with their needs for residential schemes due to poor planning and organization which have contributed to the emergence of social and financial motivation for unplanned urban sprawl represented in the dominance of some tribes on large areas of land who capitalized on the housing deficit and the lack of new housing projects to sell land as residential plots without approval from the state. On the other hand, after the fall of the previous political regime in 2011, some people misappropriated state funds and bought lands in or around cities and then constructed residential or holiday houses as a means of money laundering. Additionally, some people from outside the district, particularly displaced citizens and migrants from some Arabic countries who left their homelands fleeing the menace of ongoing wars and armed conflicts, have shown interest in acquiring land and building houses to live in the study area. These people perceive the study area as comparatively safer than other regions in Libya, so land ownership is regarded as a means of attaining stability.
Additionally, because the region of Al-Jabal Al-Akhdar includes the finest places for tourism and recreation in Libya, the demand for accommodations for tourists and visitors has increased, which encouraged some citizens to establish more holiday homes and resorts in the absence of government control and stringent regulations.
The results in Table 6 illustrate that BLV land during the third period increased by 289 ha, achieving a positive annual rate of change of 41 ha.year −1 (3.15%), which is greater when compared to the previous periods. The results shown in Figure 10 shows that clearance of NMF land and abandonment of OR land were the main contributors to the increase in the area of BLV land, where their contributions were estimated at 763 ha (64.88% of BLV gains) and 372 ha (31.63% of BLV gains), respectively.
On the other hand, forest fires were a serious occurrence that led to the increase in the area of BLV land and to the decline of natural vegetation. Government reports indicated a significant increase in forest fires in Libya after 2011. In May 2013, after a dry season with severe shortage of rainfall and severe heat wave, a total of 5924 ha was destroyed in Libya due to more than seven fires of at least 40 ha, most of which occurred in the Al-Jabal Al-Akhdar region, which includes a single very large fire of more than 3300 ha [70]. Some parts of our study area, especially the Lussaitah area, were affected by these fires. In 2017, more land was affected by indiscriminate fires that occurred simultaneously and stripped considerable areas of their natural vegetation cover. To this day, forest fires are still a real threat to the natural forests in the Al-Jabal Al-Akhdar region, where fires occur every year at several locations due to various reasons.
The increased demand for firewood during the years following the Libyan uprising as a result of the wars which broke out in some Libyan cities was another factor that contributed to the decline in forest cover and to the increase in BLV. A severe shortage of cooking gas and kerosene in addition to electricity outage were common during the period after 2011. This has forced citizens to seek alternatives for heating and cooking; regrettably, they found cutting forest trees to be the easier option for gathering firewood and for producing charcoal for use or selling in local markets. Cutting forest trees is a negative phenomenon that can lead to stripping plenty of spaces of their natural vegetation cover. In fact, in such a semiarid climate, the loss of natural trees is difficult to recover due to the inability of the existing tree species in the region to regenerate naturally in addition to other limiting factors such as insufficient moisture due to frequent droughts and soil shallowness resulting from increased erosion rate after removal of protective vegetation.

Conclusion
Remote sensing is an effective tool for mapping natural resources at a regional level. The outcomes of the investigation enumerated in this research highlight the importance of satellite digital image processing with the aid of GIS technology in mapping and detecting LULC changes. Maximum likelihood supervised classification of Landsat TM5, ETM+7, and OLI8 imagery and post-classification comparison approach can be used to predict and obtain fairly accurate maps of LULC changes. The overall accuracy and kappa coefficient values of the LULC classifications achieved in this research is encouraging. It was generally observed that there was a significant decline in natural Mediterranean forest (NMF) and a corresponding increase in other land uses such as orchards and rain-fed agriculture (OR), irrigated crops (IC), urban and built-up lands (UB), and bare and low vegetation (BLV). The changes observed were indicative of a decrease in the expanse of NMF land, which lost about 39% of its total area over the 32 years, and the highest deforestation rate was registered in the third period (2010-2017). The high deforestation was generally a result of citizen exploitation of the weakness of the state's institutions post 2011 revolution in addition to the frequency of forest fires. Most of the changed area of NMF were converted for agricultural land use. Land under the OR class increased by about 55% and achieved a greater rate of increase during the third period, whereas the land under the IC class increased by about 85% with a higher rate of increase achieved between 2000 and 2010, and the lowest rate was recorded during the third period from 2010 to 2017. The urban expansion in the study area was mostly at the expense of forest and agricultural lands. The annual rate of urbanization during the third period was approximately 10 and 5 times more than those recorded in the first period (1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000) and second period (2000-2010), respectively. Our results indicate an unstable trend of BLV land which generally increased by about 50%.
Indeed, if the present trend of deforestation continues, the forest in the study area will lose more than 50% of its current natural forest cover over the next 15 years. Deforestation in the study area and Al-Jabal Al-Akhdar in general is an important environmental challenge.

Conflicts of Interest:
The authors declare no conflict of interest.