Much research in recent years has focused on the impact of urbanization and land degradation on agricultural activities. These phenomena are important problems globally, regionally, and locally particularly in the rural landscape  . Many studies have reported that urban growth has the greatest impact on transforming agricultural lands into other uses    . However, it has been found that not only acceleration of urbanization has a significant impact on the amount of land used for food production, but also shifting sand, population increases and soil degradation   may have the same effect.
As a response to economic growth worldwide, several industrial countries have adopted new policies to meet that growth; one of these is industrialization  . For example, Saudi Arabia since the discovering of the petrol oil in 1923 has shifted from only agricultural dependence to industrial dependence. As a result, this shift has increased the urbanization process  . These increases had very consequences for the study area; accelerating urban growth (urbanization), led to high consumption of food and water, and increased internal migration.
According to  Eastern region is one of the highest concentrations of population in Saudi Arabia with approximately 3,799,773 populations, and 36.7% houses rented. A recent study by  reported that the transformation of agricultural lands can lead to a series of “social, cultural, environmental and economic problems” in recent times or even in the future. One of the most critical issues, when built-up-areas extend into arable lands, is decreasing of agricultural activities and potentially insecurity with respect food supply  . Brown  has mentioned that amount of valuable cultivation lands have been shifted to others’ land use including road, homes, industries as a result of the urban growth; in Kumasi Metropolis as an example arable lands were declined by 15.6% from 1986 and 2006. Likewise, in Malaysia, as Elhadary et al.  reported most of the croplands have been converted to urban uses by several reasons one of them is that most of the owners tend to sell or rent their lands rather than cultivation because of the low income and high-input cost of cultivation.
On the other hand, several studies have reported that urban growth is not the only factor leading to arable land reduction. Land degradation is another significant factor and the impact of this factor might increase in the future as the results of climate change. Reyonlds et al.  stated that “Land degradation is frequently associated with vegetation productivity decline over time”. The process of land degradation can lead to vegetation productivity decline by either gradual soil degradation  or by land transformation. Land degradation of different types  may degrade most of the croplands in many African’s countries. In Sudan for example, sand encroachment threatening croplands in many states (e.g., White Nile, Northern State, and River Nile),   for example Salih et al.  found that following sand encroachment the normal reaction is shrinking, reduction, and declining of cropland productivity. And therefore, there may be a food gap in the area.
In our study area, this phenomenon was demonstrated by few studies    , for instance, Allbed et al.  have documented that arable lands in the Al Hassa oasis have been suffering from shifting sand and soil salinity as a major threat. For example,  found out that there is a very positive strong correlation between salinity and degradation in vegetation cover, and they detected that change in concentration of soil salinity resulting in a change of vegetation cover density. Al Saud et al.  showed that the some of the Al Hassa oasis arable lands have been covered by sand dunes. Their study suggests that a possible cause for these dunes is that the location of Al Hassa oasis surrounded by many deserts (see Plate 1). It is clear that these deserts considered as the main source of these dunes. Holm  found out that (“the study area is located in the center of the Arabian peninsula with a total area estimated at 1,060,000 square miles, has a sand cover on at least 300,000 square miles surrounded by large deserts”). This meant that desertification and land degradation with its broad meaning is going to take a place in the study area.
In spite of the considerable amount of literature on the risks of urbanization and land degradation, there are still important gaps with respect information that documents these impacts. Few researchers have addressed and investigated the problem of agricultural land loss in the study area. Previous work has focused only on the causes, they don’t provide any statistical information about the quantity of the transformation of the land, therefore, the current work needed for future development and conservation of the lands. Until today, there are no interesting and significant findings to highlight the hotspot issues about the changing and persistence of agricultural lands in the study area.
To address the gap in data availability with respect to the impact of urbanization on the arable land conversion in Al Hassa’s oasis, Eastern region, Saudi Arabia, we employed well-known and widely used digital image processing techniques and GIS-based tools. We used different sources of data including remote sense satellite images as a primary source of information for land use/cover (LULC) classifications and change assessment  , along with other data sources for investigating, assessing and providing a wide preview for the different factors behind the urban expansion and arable lands conversion that have taken place in Al Hassa’s oasis. On the basis of the findings, this paper argues that urbanization and to great extent salinity is the major driving force for Al Hassa’s oasis arable land conversion. Next, we applied an accuracy assessment statistical model using ground truth to measure the degree of accuracy of the final land cover classification results.
Plate 1. The accumulation of the sand dunes alongside the study area with different size and shape.
This paper is organized as follows section one used to describe the data and methods used in this study, and section two is used to report our findings, while section three and four are used for the discussion of findings, and conclusions.
2.1. Study Area
The present work took place at Al Hassa oasis, Eastern Region in Saudi Arabia. According to  the term “oasis” has come to mean “a cluster of palm trees near a waterhole with perhaps a coupled of tents, all encircled by sandy desert”. Whereas Al Hassa name is “an area of vegetation surrounded by desert”. Al-Hassa Oasis is considered to be one of largest agricultural land. It is located in the eastern province of the Kingdom of Saudi Arabia. Al-Hassa is an L-shaped area  , lie approximately between 25˚21' to 25˚37' Latitude. N. and 49˚37' Longitude. E. (Figure 1), at an elevation of approximately 130 to 160 metres above sea level   ; covering an area of approximately 260 Km2. The selection of the area might be justified by that it is the only cultivated land in Al Hassa province. Moreover, it is reported that the area has long been affected by sand movement and human activities which led to decrease in the lands available for cultivation.
It has been reported in  that the area is the important and greatest agricultural area in the eastern region of Saudi Arabia; it is establishing an extension being due to the huge reserves of groundwater. The primary crops grown were dates, wheat, barley, millet, rice, various fruits and vegetables and alfalfa and the most important livestock raised were camels, sheep, goats, cattle and donkeys  . Currently, the total cultivated land area within the oasis is approximately 80 km2, of which 92% is planted by date palm  .
Furthermore, there are over 51 small towns scattered over the oasis that are either fully or partially surrounded by date palms. In relation to the climate, the area is mainly located between hot-dry and hot-humid conditions because of its closeness to the seashore. Around 12,000 hectare of vegetation cover is a major driving force for the improvement of the climate characteristics  . The temperature exceeds 45˚C during the summer season  with annual rainfall approximately 488 mm.
Figure 2 shows the population census (in thousands) of Al Hassa Oasis from 1935 to 2016. In 2010 1,067,691 persons were lived in the study area  . It increased to nearly 1,063,112 persons in 2014 while in 2016 the population increased to 1,220,000. Data from recent study by Abdelatti et al.  reported that the population of Al Hassa oasis rose during the period 1992-2016 from 445,000 in 1992 to 768,000 in 2016. Such an increase in population may require a significant amount of lands for urban development. Consequently, agricultural lands are under pressure and negatively may impact it. Accordingly, the present work is established and carried out by adapting remote sensing data as the main source of information.
Figure 1. Location of the study area indicated by the red square within the Eastern Region of Saudi Arabia (The background image is landsat standard false color composite (RGB: 5, 4, 3) of the year 2017).
Figure 2. Population numbers in Al Hassa oasis by population Census for the periods 1935-2016. Data obtained from   .
2.2. Data Used and Processing
The current investigation involves an in-depth analytical investigation of urbanization and its impact on agricultural lands. To perform this, a qualitative  and quantitative approaches to data collection and results interpretation were used. Integration of data from different sources (e.g., field word data, and remote sensing data, available reports and previous studies) to assess changes and trends of the agricultural lands at the study area has been achieved. The primary source of information was from satellite images which were obtained from the Landsat archive data center. Consequently, three satellite imageries were downloaded from the USGS archive website “https://earthexplorer.usgs.gov/”, these are 1) Landsat-5 (TM) for the 1984 and 2000 and 2) Landsat-8 (OLI) for the date 2017 as shown in (Figure 3). The characteristics of the satellite images are provided in (Table 1).
Figure 3. The three selected Landsat imagery: (TM) for 1988, 2000, and (OLI) for 2017, respectively, used to study the agricultural lands change in the study area.
Table 1. The characteristics of the Landsat satellite images used in the study.
In addition to all these the study referred back to different reports and previous studies to formulate and adjusting the objectives and hypothesis and to detect the knowledge gaps related to the research problem.
The initial step for properly accomplishing change detection comparison for unprocessed images and classified images is image preprocessing techniques  . Generally speaking, image processing is an essential processing for those photos with 2D elements using a computer capability  . Image preprocessing is an operation widely used prior to the main analysis of satellite images  . Accordingly, the necessary image preprocessing (i.e., geometric correction, and radiometric calibration) was accomplished. For the first one, we used the image to image registration method; as a result, the images can be overlapped with vector data in a geographic information system (GIS) and used for comparison of land use/cover classifications. While for the second one, we used the methods presented by  as a prerequisite for creating high-quality image data. Finally, the concept of tasseled cap transformation was applied to each of the four images using ERDAS 2014 for compressing spectral data into a few bands associated with physical scene characteristics and to make the visual interpretation and digital classification process for the satellite images as easier as possible.
It is a linear transformation of data that projects soil and vegetation information on a single plane in multispectral data space. The results of the TCT produce an image file that consists of three bands that are attributed to soil brightness, vegetation greenness, and soil or vegetation wetness. Detailed information about the idea behind this transformation can be found in the following literature    .
The output from the previous methods was used as input for the images classification. In this work a clustering system (i.e., Iterative Self-Organization Data Analysis Techniques (ISODATA)  known as unsupervised classification approaches were applied to classify the three preprocessed satellite images into three different land use/cover classes. Firstly this approach classifies the image into a number of unknown classes as an initial step. It is then the analyst’s responsibility  to labels each homogenous group of pixels with the suitable class name as described in (Table 2). Managing took place when two or three classes have the same pixel values which mean that they belong to the same feature class. The reason behind choosing this method over other methods it is simplicity comparing to others methods (i.e., sub-pixel; supervised classification methods). All the preprocessing and analysis of data have been carried out using ERDAS IMAGIN software version 2014.
Table 2. The definition and description of the mainland use/cover used in this study.
According to  remote sense data can’t be used without ground truth datasets. Ground truth can be used alongside remote sense data for several reasons as follows: Firstly, assessing the visual interpretation process, secondly, geometric correction of satellite images, and finally, to verify information extracted from remote sense data  . Accordingly, fieldwork data were collected in January 2018. GPS device was used for collecting, locating and positioning of ground truths and ground control points (GCPs). A set of about 100 training points was used for several purposes these are a geometric correction, accuracy assessment representing the corresponding feature on the ground, where overall accuracy, commission, and omission  statistical information were calculated for each classified images, and for signature evaluation. In the present work to identifying differences in the state of agricultural lands and to get a clear picture about the impact of urban expansion on agricultural activities in the study area a comparison of land cover classifications  in between the classified images was carried out
The Normalized Difference Vegetation Index (NDVI) was used in this study as described by  for assessment of vegetation cover change. The NDVI is a widely used method for detecting the change in vegetation greenness rate. It is a simple measurement of the differences between red channel reflectance and infrared channel reflectance, by which the area under study can be assessed for agricultural drought  . The values of the NDVI generally ranged from −1 to +1. It is clear that the area tend to be covered by vegetation when the value of NDVI is near (+1), if the area is covered by bare soil or bare rock, the NDVI values will be near zero, while if the area is covered by water or sand, the values will be near (−1). In this way, the study area was classified to the desired land cover classes (Table 2). Differencing  on two dates of NDVI were subtracted one from another to determine change ∆NDVI image with the positive increase and negative decrease.
People are interested in change detection for several reasons. These are for example; determine the changes that have occurred and extrapolate them to the future, decide on a course of action, or to assess the results of an action or policy  . By mapping where and how things move over a period of time, you can gain insight into how they behave in the future. Another reason for mapping change is to predict future needs.
By mapping earth features conditions before and after an action or event, you can predict the impact. You can map change by showing the location and condition of a feature within a different date. Knowing the type of change and the type of features you’re dealing with, how you’re measuring time, and the type of information you need from the analysis will help you decide how to map  . Accordingly, there are two types of change:
1) Change in location: mapping change in location helps you see how features behave in the future.
2) Change in character or magnitude: Mapping change in character or magnitude shows you how conditions in a given place have changed  . The change can be in the type of feature in a place, for example, the different categories of land cover/use in the study area now, compared to 15 years ago. Or, the change can be in a quantity associated with each feature, for example, the amount of feature class has increased or decreased in each area over the past 15 years.
According to the above description change detection based on dividing the classified images into different size and shape of watersheds was accomplished using QGIS 2.8.2 software and topographic data which was obtained by the Shuttle Radar Topographic Mission (SRTM). These data are available at a 3''-(90 meter) horizontal resolution and 16 meter vertical accuracy  .
3.1. Land Use/Cover Conversions in Al Hassa Oasis
In earlier studies attempts were made to classify and evaluate the land use/cover in the study area using remote sensing techniques (e.g. supervised classification approaches). These studies indicated that the Al Hassa oasis’s land cover consists of different land cover classes.   classified the land use/cover of the study area into 5 to 7 cover classes including buildings, vegetation cover, sabakh, sand dunes, sand sheets, barren land. They documented that the Al Hassa oasis is dominated by barren and shifting sand feature.
In the present study, as the main purpose was to use satellite images to assess to what extent the arable lands in Al Hassa oasis were affected by urban development, three land use/cover categories were analyzed and mapped which include: agricultural lands, urban lands, and barren and sand dunes using unsupervised classification method. The analysis results reveal that the dominant land cover class in the study area is urban besides the barren and sand dunes as shown and presented in (Table 3) and (Figures 4(a) to (c)). This class has increased in area over the years. Comparison of classified images revealed that different situations were found in relation to the conversion of agricultural lands in the study area. In 1988, urban land-use covered around 185, 1 km2 of the total area of Al Hassa oasis, and then increased to 254.52 km2 in 2017, an increase of 69 km2 during thirty-three years of time. The high rate of urban growth (12%) that exceeds the world average rate of between 4 and 4.5 per annum  occurred as a response to the increase in the population of Al Hassa oasis due to the oil discovery. The area covered by arable lands, however, decreased from 160 km2 in 1988 to 156 km2 in 2000. This decrease can be justified by the development of the first oilfields to the north and east of the oasis of 1950 to1989 which resulted in the social and economic change impact on arable lands in the Al Hassa oasis  and in the same time led to the increase in urban use that took place between the period 2000 to 2017. By different words, the developments that took place in the Al Hassa oasis since 1950 have encouraged the transformation from the agricultural sector to other sectors (e.g., oil utilization, industry, trade, construction, and civil service). These non-agricultural sectors have opened the door to new employment opportunities and accelerated the immigration process to the urban area  which has the greatest negative effect on agricultural activities. In 1974 it is estimated that about 695,000 persons were practice agricultural as the main source of income. This number of persons was declined by about 96,000 persons in 1979  . In addition to all these factors, increasing in soil salinity due to inefficient use of water and insufficient drainage has also led to a decline in the area covered by crops  .
It is important, however, to note that science 2000 there has been an improvement in a cropland area, it slightly increased by 22 km2 The increase in agricultural area could be justified by many factors, including (According to local respondents: 1) establishment of new irrigation projects to the north and south of the study area  , 2) the stabilization of sand dunes project by planting different type of species in the north of the study area, 3) the construction and improvement of irrigation and drainage systems increased the cropland area to the south of the study area.
In the period 2000-2017, the largest decrease in the areas covered by barren lands and sand dunes can be detected. It decreased by almost −85.41 Km2 in 2017 compared to the period 1988-2000. However, the obtained result is not consistent with the findings of other earlier studies e.g.,   and  . It is immediately obvious that there are some quite significant differences, which are areas covered by sand dunes are decreased in 2017 which explain the high NDVI values obtained for 2017 compared to other years. It is important to note that according to this result the study area has been affected by a sand movement which may lead to loss of some arable lands in addition to the losses due to urban expansion as the main driving force about arable land conversion between 2000 and 2017.
It should, however, be noted that in the land cover classification comparison, the accuracy of the classification itself may affect the results. Nevertheless, this effect is only local. Future work should, therefore, include other methods based on pixel comparison rather than on classes.
Data from (Table 4 and Table 5) respectively, show the land cover change matrix with more detailed change information  . For the period 1988-2000, it revealed that around 80.31% of the area did not exhibit a change in the three land covers. While around 39.7% of the area did not exhibit a change in the same three classes. It can be noted that the arable land showed a slight decrease in the area for the period 1988-2000 and at the same time exhibit greatest area with no change (83.81%); this might be explained by the large area occupied by this cover type. Whereas the same cover type is shown a slight increase for the period 2000-2017 by (69%). The area covered by urban increased by 14% from 1988 to 2000, and by 11% from 2000 to 2017.
To verify our classification results and confirm what we did and found, an accuracy assessment was accomplished using ground truth data. The results of this verification can be found in section 3.2.
Table 3. Change in the Arable land, urban lands, and barren and sand dunes cover in the study area, Al Hassa 1988-2017.
Table 4. Land-Cover change matrix and transitions in study area, Al Hassa 1988-2000.
Table 5. Land-Cover change matrix and transitions in study area, Al Hassa 2000-2017.
Bold numbers indicate not changed pixels between two dates.
Figure 4. Land-cover classifications for three time periods; 1988 (a/upper), 2000 (b/upper right) and 2017 (c/bottom).
3.2. Accuracy Assessment Results
A quantitative analysis to determine analysis uncertainties was applied, based on an accuracy assessment model described by  . Based on these approaches, the obtained accuracies for all classified images slightly exceed the acceptable limit defined by the expertise in the area of remote sensing application. Nevertheless, these results suggest that data obtained from satellite images supported by data obtained from fieldwork may provide more accurate and up to date information for assessing the impact of urban development on the agricultural future. Tables 6-8 provide statistical information for the overall map accuracy assessment results. We found that the map accuracy of all classified images is 74%, 87%, and 89% for the classified images of 1988, 2000, and 2017. These tables also provide information about the producers and users accuracy of the land use/cover classifications.
The producer’s accuracy of the agricultural category is about 100 percent, for all classified images, while in the urban category is in the range of 60 percent for the image 1988 to less than 80 percent for the image 2000, which is quite acceptable. For the barren lands and sand dunes category, the producer’s accuracy seems to be imponderable. By reference to the user’s accuracy of all classified image, we revealed that about 86% of the agricultural lands of image 1988, 89% of image 2000, and 100% of image 2017 has been correctly identified as agricultural lands on the ground.
The user’s accuracy of 1988 for the urban and barren land reveals that there is significant confusion in discriminating urban lands from barren lands. This can be justified as follows: 1) the study area is dominated by a sand dune which has an effect on the reflectance of the others land cover and land use materials  , 2) in the study area the materials that buildings have been made from are sand particles, which have the same spectral reflectance of the sand dunes category. All these factors reduced the accuracy of the classification process and provided misleading results. Generally speaking, the obtained accuracy assessment results revealed that the dataset and methods that we used and the results that were obtained from the analysis can help us in assessing and making a judgment on the current environmental problem related to the conversion of the agricultural lands to other uses. This could lead to the fact that the moderate resolution data such as Landsat-5 and 7 can provide useful information to assess human impacts on agricultural lands as reported by  .
Table 6. Statistical report of the classification results of 1988 classified image.
Table 7. Statistical report of the classification results of 2000 classified image.
Table 8. Statistical report of the classification results of 2017 classified image.
3.3. Normalized Difference Vegetation Index Results
The main aim of this section is to provide information about the general status of vegetation cover in Al Hassa oasis based on satellite images. Data obtained in previous studies   using remote sensing data indicated that different indices can be used as an effective tool for change detection in land use/cover studies. For example, according to Allbed et al.  , the NDVI index can give optimal results when used for change detection studies compared to post-classification comparison methods. In the present work, the well-known and wieldy used NDVI index was used to detect a change in vegetation cover along the study area. NDVI results have positively increased the accuracy of the final products obtained from the unsupervised classification approach. Figure 5 shows the results obtained using the NDVI index. Statistical information for areas covered by health vegetation is summarized in (Table 9 and Table 10). As can be seen, the maximum mean value of the NDVI is 0.933 of the image 1988 and −0.029 for the image 2000. The maximum and minimum NDVI values were 2.300 of the image 1988 and 0.493 of the image 2000. Also, it can be noticed that the percent of very health vegetation was significantly lower in 1988 than in 2000 and 2017, which isn’t consistent with results obtained in previous studies e.g.,  . However, the obtained results are in line with the recent situation of the vegetation cover in the study area which is that it experiencing remarkable change over the study area. In other words, there has been very strong evidence that the vegetation in the study area is more affected by drought and other factors than ever before especially in the northern and southern parts. Although, the previous efforts undertaken by  towards this effect, it is still the area suffer from many aspects of drought, for example, sand movement (Plate 1).
It can be observed in (Figure 5) that in 1988 compared to 2000 and 2017 there was a considerable decline in vegetation cover. In spite of the considerable efforts made by  toward vegetation cover in the study area, stress in vegetation cover was clear in 1988 compared to years 2000 and 2017. Furthermore, between 1988 and 2000 a decrease in vegetation cover was observed especially in the north of the study area. Whereas the increase in vegetation cover was observed between the years 2000 to 2017 mostly towards the southern part of the study area. From (Figure 5) it can be seen that the vegetation greenness was improved where the high NDVI value was 0.30, representing around 7% as health vegetation in 2017.
Table 9. The area covered by very health vegetation along the study area based on NDVI derived data, using Landsat images of three years.
Table 10. Statistical Information of the NDVI values of vegetation covers along the study area for different study periods of time derived from Landsat satellite images.
Figure 5. The area covered by very health vegetation along the study area (1988, 2000 and 2017).
Based on this approach, it was difficult to judge that if the change and decline in vegetation cover occur by urban expansion or by other human activities or by land degradation (i.e. soil salinity). Nevertheless, these results suggest that data obtained using NDVI index to evaluate to what extent vegetation cover degraded may provide more sensitive and accurate information for assessing the impact of drought on vegetation cover in general than other traditional methods.
Data obtained in recent study by  using NDVI image differencing technique as efficient method, which is a process of pixel-by-pixel subtraction of different NDVI images for land cover change detection indicated that a decrease in vegetation cover was detected and concentrated in the northern part of the study area, whereas the increase of vegetation cover was observed throughout the study area for the 1985-2000 period, and was only concentrated in the southern part between 2000-2013 period.
In our study, the same technique was used to obtain information in relation to vegetation cover change. (Figure 6) demonstrates the results obtained using the subtraction method. Data in (Figure 6) indicate that there were a significant decrease and increase in vegetation along the northern part of the study area. These trends are in line with the previously discussed results obtained by  . The results seem to indicate that there was no change (0 values) in vegetation cover in the center of the study area. The overall volume changes are presented in (Figure 6 (upper)). Comparison of (Figure 6 (upper and bottom)) shows that the vegetation cover was increased positively between the periods of 1988-2017 and negatively between the periods of 1988-2000. Overall, however, data in Figure 7 suggest that the study area has suffered a severe change in the vegetation cover, particularly between the periods of 1988 to 2000. Overall, it is believed that the general trend in the study area refers to urban expansion in the regions of ALuyun, ALjafor, Hofuf, and Joatha. It should, however, be noted that in the subtraction method, as pointed out by  , mixed features within one pixel may cause some problem in the final result and despite the change in pixel value the changed pixel might still represent the same land cover class.
Figure 6. Vegetation cover change detection using Landsat data (TM, ETM+, and OLI) between the years of (1988 to 2000 upper), and (2000 to 2017 bottom). (-) indicates for negative change, (0) indicates for No change, and (+) indicates for positive change.
Figure 7. Croplands, urban areas, and barren sand dunes cover change. Notice: Red = Urban Lands, Green = Agricultural lands, and Gray = Barren land and sand dunes.
Data obtained in SRTM using GIS-based tools indicated that change detection is more obvious and easy to detect when the study areas were divided into different size and shape of watersheds. The main necessity here is to highlight the areas where urban and cultivated lands have decreased or increased or remained unchanged during the study period. Figure 7 show the results obtained from change detection procedures based on GIS tools. As can be seen in Figure 7 there are great differences in terms of the size and shape between the three watersheds in all years. As can be seen in each watershed the differences between the three lands cover classes for each date. In all watersheds, it is clear that the urban use was prevailing in 2017 compared to 1988 and 2000. There are some similarities, however. For example, in both 1988 and 2000, the proportion of cultivated land was roughly the same area, particularly in Watersheds no 1, 4 and 5. Moreover, from these figures, a significant expansion in the urban area was detected in all directions. This meant that there is no control and a clear plan to restrict this increase in urban development along the study area.
Prior works have documented that the effectiveness of sand movement can lead to significant losses of agricultural lands especially in arid and semi-arid environments where the water shortage is severe, and there is drought;   for example, reported that if legislation protecting arable areas is not enacted then they may be transformed to others uses. Using visual comparison it can be observed that the sand dunes and barren lands are the dominant cover features in watersheds numbers 2, 5 and 6. Although the agricultural activities in the study area have been receiving much attention from the government during the last five years, it is still suffering from sand encroachment. The decline in area covered by crops are very clear in 2017 specially in watersheds number 2, 3, and 4, whereas the increase is more obvious in other watersheds for instance 6 and 7.
Prior work has documented the impact of urbanization on the future of agriculture particularly in developing countries;  for example, report that rapid urbanization might generate great pressures on arable lands, also they documented that urban expansion has a significant effect on the countryside by changing their landscape and lifestyle of rural communities. Other study carried by  focused on the impact of sand movement on the agricultural lands in Al Hassa oasis without considering the impact of urban expansion. However, these studies have either focused on one factor affecting the future of the agriculture or not based on satellite images and GIS techniques as these methodologies facilitate to a greater extent the achievement of the present study.
This study has attempted by employing remote sense and GIS techniques to evaluate and assess the impact of urban development on Al Hassa oasis’s agricultural land using different data sources; where significant results were found. The urban areas increased during the period of 1988-2017 while the arable lands were decreased. In 1988, urban land-use covered around 185, 1 km2 of the total area of Al Hassa oasis, and then increased to 254.52 km2 in 2017, an increase of 69 km2 during thirty-three years of time. This increase in urban use can be supported by recent study findings carried out by  has addressed the nature and trend of urban in the Al Hassa oasis they revealed that as the socio-economic situation has improved, the land use pattern has changed and urban area has increased rapidly for the periods 1957-2011. They also pointed out that the huge increase fraction was occurred in the built-up area by (366%) for the period 1985-1994 with concentration mainly in Hofuf and Mubarraz cities. In contrast, the area covered by arable lands however decreased from 160 km2 in 1988 to 156 km2 in 2000. This decrease can be justified by the development of the first oilfields to the north and east of the oasis of 1950 to1989 which resulted in the social and economic change impact on arable lands in the Al Hassa oasis  and in the same time led to the increase in urban use that took place in the period 2000 to 2017. In the same context it is important, however, to note that since 2000 there has been an improvement in a cropland area, it slightly increased by 22 km2. This increase have been justified by many facts, which includes: 1) establishment of new irrigation projects to the north and south of the study area  , 2) the stabilization of sand dunes project by planting different type of species in the north of the study area, 3) the construction and improvement of irrigation and drainage systems increased the cropland area to the south of the study area [According to the respondents]. The increase in crop area has been supported by a recent study of  ; they pointed out that the change in NDVI values increased by 7% for the period 2000-2013. They also revealed that approximately 7% of the total the Al Hassa area experienced positive change (i.e., increase in general vegetation cover) between 2000-2013, whereas nearly 3% of showed negative change, that is a decline in vegetation cover for the same period of time. The increase in vegetation cover in general may reflect the increase in crop area.
This increase and decrease in urban and agricultural areas can be justified for many reasons. Firstly, lands owners seem to have largely transformed their arable lands to other investment because of the low income of the agricultural productivity compared to the expenses on the land cultivation. That is to say, in terms of recent development in the area, it has been shown that most of the owners of the land tend to convert their land to urban uses because it saves their many and in the same time can lead to more outcomes. Secondly, it is clear that the study area falls in the desert, accordingly, it has long been threatening by shifting sand. This shifting sand changed most of the soil characteristics and has become a big challenge for the farmers to keep their lands with fertile soil. Moreover, a recent study carried out by  has revealed that the most of the arable lands in Al Hassa oasis have been affected or affecting by soil salinity, and they justify that by different factors one of them is poor quality irrigation of groundwater. The same study and other several studies   have found that the soil in Al Hassa oasis’s arable lands contents high EC ranging between 1.4 and 16 ds/m and high salinity. From all these studies it can conclude that there is a strong relationship between the changes in vegetation cover and the salinity of the soil, increase in soil salinity might lead to a decrease in vegetation cover.
Furthermore, based on the survey undertaken in the study area it was found that the farmers preferred to work in other sectors rather than cultivating their land. Moreover, in some cases, they can be employed in nonagricultural jobs. It is important to mention that agricultural lands in the study area are facing great challenging in relation to the current economic situation. That because of the oil operation in the Eastern Province, in general, has led to the establishment and rapid growth of new settlements  .
The findings of this study are consistent with the study by Elhadary et al.  carried out in Seberang Perai region, Malaysia (Despite the fact that both areas are located in the Asian countries and almost have the same population number. It is important to note that the pattern of land use system and the way of managing it are slightly different. Malaysia depends on growing rice and oil palm while Saudi Arabia is based on growing date palm. Both countries have witnessed rapid expansion of urban building at the area of agricultural land and this phenomenon found in most of the developing countries) which showed that many arable lands are transformed to urban use for many reasons one of them is the urban development in the region, while the second one is that most of landowners tend to sell or rent their land to others who might keep it without cultivation, and finally they prefer to transform their land to others uses as a response to the economic investment.
Overall, these studies suggest that there will be no future for the agricultural practices in the study area if all these challenges continue. As consequences, the owners of the land probably preferred to choose others activities rather than keep their land without any cultivation.
This article has assessed and evaluated the issue of urbanization and its impact on the arable lands conversion in the Al Hassa oasis of Eastern Region in Saudi Arabia using remote sensing and Geographical Information System (GIS). The analysis results from this study can be considered as an alarm for the arable land conversion in Al Hassa oasis. It might lead to social and environmental impacts in the near future. In addition, the findings presented in this study expand prior work. The study provides a framework for future studies to assess current problem solutions. However, some limitations are worth noting. Although our findings were supported by accuracy assessments, it still needs more validation. Future work should, therefore, be achieved to evaluate the impact of urbanization on agricultural lands using sub-pixel classification comparisons and data with high resolution.
We wish to thanks the anonymous reviewers for their valuable and useful suggestions and comments on the manuscript and their valuable editing and enhancement of the English language. Also, we would like to thanks the USGS organization for providing the remote sensing data for free.
Data processing and writing of the methods section have been achieved by the third author. The interpretation of the results, the discussion, and conclusions sections were written by the second author. Finally, the introduction and the abstract sections were written by the first author.