Modeling Past, Present, and Future Urban Growth Impacts on Primary Agricultural Land in Greater Irbid Municipality, Jordan Using SLEUTH (1972–2050)

: Urban expansion and loss of primarily agricultural land are two of the challenges facing Jordan. Located in the most productive agricultural area of Jordan, Greater Irbid Municipality (GIM) uncontrolled urban growth has posed a grand challenge in both sustaining its prime croplands and developing comprehensive planning strategies. This study investigated the loss of agricultural land for urban growth in GIM from 1972–2050 and denoted the negative consequences of the amalgamation process of 2001 on farmland loss. The aim is to unfold and track historical land use/cover changes and forecast these changes to the future using a modiﬁed SLEUTH-3r urban growth model. The accuracy of prediction results was assessed in three different sites between 2015 and 2020. In 43 years the built-up area increased from 29.2 km 2 in 1972 to 71 km 2 in 2015. By 2050, the built-up urban area would increase to 107 km 2 . The overall rate of increase, however, showed a decline across the study period, with the periods of 1990–2000 and 2000–2015 having the highest rate of built-up areas expansion at 68.6 and 41.4%, respectively. While the agricultural area increased from 178 km 2 in 1972 to 207 km 2 in 2000, it decreased to 195 km 2 in 2015 and would continue to decrease to 188 km 2 by 2050. The district-level analysis shows that from 2000–2015, the majority of districts exhibited an urban increase at twice the rate of 1990–2000. The results of the net change analysis of agriculture show that between 1990 and 2000, 9 districts exhibited a positive gain in agricultural land while the rest of the districts showed a negative loss of agricultural land. From 2000 to 2015, the four districts of Naser, Nozha, Rawdah, and Hashmyah completely lost their agricultural areas for urbanization. By 2050, Idoon and Boshra districts will likely lose more than half of their high-quality agricultural land. This study seeks to utilize a spatially explicit urban growth model to support sustainable planning policies for urban land use through forecasting. The implications from this study conﬁrm the worldwide urbanization impacts on losing the most productive agricultural land in the outskirts and consequences on food production and food security. The study calls for urgent actions to adopt a compact growth policy with no new land added for development as what is available now exceeds what is needed by 2050 to accommodate urban growth in GIM.


Introduction
Since the late 20th century, urbanization has become one of the leading anthropogenic activities transforming the Earth's surface and consuming agricultural land; by 2000 about 88% of these areas were destroyed by urban conversion [1][2][3]. Cities are home to more than 54% of the world's population, and by 2050 this percentage is expected to reach 66% [4]. To provide a favorable-living environment and housing units for the growing population, cities have rapidly grown beyond their peripheries, mostly on the primary agricultural lands [5]. A low-intensity urban sprawl pattern has evolved as the reliance on the automobile and the desire for suburban living increase. The spatial extent of urban areas worldwide is expanding twice as fast as their population [6,7]. development and sustaining limited agricultural resources, especially farmlands. From 1970 to 2010, the Arab urban population grew four times, and the number is expected to double by 2050 [42]. Here, city planning is often challenged by the disconnection between national strategic planning cycles and spatial planning for urban development. This becomes more problematic when the Arab countries neglect the importance of proper management through balancing between available agricultural lands and conversion for urban use [44]. Nevertheless, and due to increasing conversion rates of farmlands to urbanization, countless efforts have been paid towards modeling long-term land-use/cover dynamics, in general, and urbanization, in particular, in these countries [44][45][46][47]. Jordan is categorized as a food deficit Arab country, where 0.5% of all households suffer from food insecurity. Jordan's arable area has undergone considerable decline since 1975 mainly due to rapid urban growth and land degradation [48][49][50]. Although the contribution of the agricultural sector to the GDP is small, farming remains a major source of food for nearly 25% of Jordan's poor [51,52]. While numerous studies were conducted to establish a linkage between urbanization and agricultural land consumption, the focus was yet on the capital city of Amman [53][54][55][56].
Located to the north of Amman in the most fertile and productive agricultural zone of Jordan, Irbid city was neglected for years. In recent years and as land-use/cover began to change dramatically, built-up urban area rapidly has grown on the productive peripheral farmlands, forming a major metropolitan area in the north. To facilitate the problem of governance of such a vast area, the number of municipalities and village councils around the city of Irbid were amalgamated to form what is known today as the Greater Irbid Municipality [57]. Accordingly, the municipal boundary of Irbid city was adjusted resulting in integrating agriculturally dominant rural districts to the city limits. Since 2001, the GIM area has undergone significant land cover changes to support the growing urban population resulting from (1) economic development (2) domestic rural-urban migration, (3) changing demographic structure, (4) changing lifestyle, and (5) the influx of Syrian refugees during the Arab Spring of 2011 [58].
While numerous case-studies investigated urban growth in Irbid using remote sensing and GIS [59][60][61][62][63], no study has yet used the SLEUTH model to predict urban dynamics and quantify the past, present, and future urban expansion on agricultural land of Jordan. The SLEUTH model is applied for the first time in Jordan, providing a robust approach for making long-term sustainable urban planning strategies that integrate spatially explicit models and geospatial technologies. In this study, the impact of past, present, and future urban development patterns and agricultural land consumption across a 324 km 2 GIM area are examined. These transformations are analyzed at the district level for the first time to highlight the most affected districts and their agricultural values. The GIM was selected as a study area because (1) nearly 40% of Jordan's primary croplands located in Irbid, (2) it is the second-fastest-growing urban area in the country, (3) it has been experiencing rapid socioeconomic and demographic transformations, and (4) its zoning and urban planning policies are not well designed, not fully enforced, and their sustainability perspectives are not pronounced.  (47 cm). Its relatively level topography made it suitable for settlement and farming activities. The major land covers are Urban and Agriculture. Rain-fed farming mainly dominates the eastern and southeastern parts of the study area where the Vertisols are the predominant soil type (Figure 1). Vertisols cover more than 60% of the study area and are found in the undulating eastern and southeastern plains. Vertisols are deep fertile soils containing high clay [64]. Inceptisols are found in the less fertile parts of GIM along the northern and western rolling hills of GIM.

History of GIM Evolution
Greater Irbid Municipality (GIM) is the second-largest populated metropolitan municipality in Jordan, after Amman, with a population estimate of 907,675 which represents nearly half of the Irbid governorate population [65]. It is located 80 km north of Amman on the Highland Plateau close to the Syrian boundaries, forming an important transportation and business hub not only in Jordan but in the entire region. Irbid city was established in 1884 with a total population of 1300. It is considered one of the oldest cities, established in 1884, and municipalities in Jordan that grew from a small village with a total population of 1300 surrounded by vast areas of high-quality farmland into a metropolitan agglomeration center in the 21st century. In its early development stages, urban areas expanded smoothly from their center around Al-Tal toward the surrounding wheat fields. That was until 1948 when the city, along with other Jordanian cities, received the first wave of Palestinian refugees. A new camp was established in 1950 on the northern edge of the city center [66]. The second wave of Palestinian refugees happened in 1967 and a new camp was built outside the city in Hoson town. Afterward, Irbid's urban area expanded substantially, and its population reached 112, 864 inhabitants. Since then, Irbid has represented the only modern city in northern Jordan attracting domestic migration from rural communities. Besides its strategic location and agricultural value, Irbid has become an educational elite center with the establishment of Yarmouk University and Jordan University for Science and Technology in 1976 and 1986, respectively. Nowadays, the metropolitan area of Greater Irbid Municipality covers an area of 324 km 2 spreading over 23 municipal districts.

Land Cover and Urban Data
For land cover information and urban data, a historical land cover timeline was developed and compiled from different sources for the years: 1972, 1980, 2000, and 2015. The 1972 land cover map was developed by visually interpreting land cover types and manually digitizing them from a Landsat MSS scene (Path 187, Row 37) that was downloaded from the Landsat data portal (https://www.usgs.gov/core-science-systems/nli/landsat) (accessed on 16 April 2016) following a consistent mapping technique. The other land cover maps were extracted from the Jordan National Land Cover Database (JNLCD) product at 30 m resolution [67]. The JNLCD is the first decadal land cover database for Jordan that was developed using a total of 44 scenes for 4 decades (11 tiles in each year). These satellite remote sensing data (Table 1) were mostly cloud-free and were obtained during the growing season. They were geo-registered to Jordan Transfer Mercator (JTM) projection and radiometrically (surface reflectance) corrected.  A set of ancillary and derived spatial data (i.e., slope, potential land use map, soil map, and derived NDVI values) were used to help to delineate land cover types. Land cover classification methods used in developing the JNLCD dataset combined unsupervised with rule-based classifications to resolve the problem of spectral clusters with low separability. Thematic accuracy for the JNLCD dataset was performed for the 1980s and 2015s maps due to the availability of high-resolution reference data and were 94 and 90%, respectively. For this study, the modified Anderson Level 1 JNLCD land cover legend was aggregated ( Table 2) to enable comparison with other land cover maps that were developed for the study area. The JNLCD detailed agricultural-related land cover types, specifically, were aggregated to agricultural land.

. SLEUTH Inputs and Calibration
The inputs used for SLEUTH-3r were : 1972, 1980, 2000, and 2015 urban extents; 1972, 1984, 1998, and 2015 road networks; land surface slope; hillshade image; and exclusion layer (Table 3). Urban coverages were extracted from the land cover maps that were developed and compiled in the previous section. For the road network, only primary roads were included and digitized from different sources (Landsat images, topographic maps (for the 1984 network), and Google Earth ® (for the 2015 network)). In the exclusion layer, we obtained the map of the zoned area from GIM which shows the areas where development is permitted. These areas were given the value of 50 indicating neutrality, neither attracting nor preventing development, all other areas were given the value of 100 indicating the highest level of protection from development. A 30-m resolution ASTER dataset was used to create inputs for slope and hillshade. All map inputs for slope, exclusion, urbanization, roads, and hillshade were converted to an 8-bit GIF format used by the SLEUTH-3r model ( Figure 2).  The SLEUTH-3r calibration process is typically undertaken using a "brute force" methodology, in which many combinations of parameter values are tested automatically. In this process, both growth rules and self-modification rules are refined to the locale. The model has a "self-modification" function to simulate different rates of growth more realistically over time. When the urban growth rate exceeds a specified critical threshold, the growth coefficients are multiplied by a factor greater than one, simulating a development "boom" cycle. Similarly, when the urban growth rate falls below a specified critical threshold, the growth coefficients are multiplied by a factor less than one, simulating a development "bust" cycle. While the calibration process is usually performed in three phases; coarse, fine, and final, researchers proved negligible gains in performance among these phases.
Therefore, a coarse calibration in which all growth coefficient values ranged from 1 to 100 was performed with only increments of 25 were tested. The coarse calibration process for each study area was initiated to first set the appropriate Diffusion Multiplier value (DM) that can capture dispersive growth within GIM. In this process, the SLEUTH-3r's diffusion growth coefficients were set to produce the maximum level of dispersive growth (i.e., diffusion coefficient was set to 100 and all other growth coefficients set to 0). The proper DM value was when the cluster fractional difference first exceeded zero. The coarse calibration was then performed to develop growth coefficients for each area, including the four control years (1972, 1980, 2000, and 2015) and using all five growth coefficients with 25 Monte Carlo trials for each study (Figure 3). To better evaluate the performance of the simulation phase and choose the appropriate growth coefficient values, representative goodness of fit measures was selected. Selecting proper fit statistics is crucial for an accurate future forecast; however, there is no agreement on a standard set. Here, the Pixel Fractional Difference PFD and Cluster Fractional Differ-ence CFD ratios were chosen because they proved to be the most relevant to the application of SLEUTH-3r. The metrics of Compare and Clusters were also chosen because they proved to provide the most robust results.
After the best-fit parameters were identified for the study areas, the model was initialized in 1972 and ran it in predict mode to 2015, with 25 Monte Carlo trials. This resulted in a predicted development probability surface for 2015, which I then compared to the observed patterns in 2015 to assess the accuracy of the calibration process.

SLEUTH-3r Forecasts to 2050
To create forecasts with SLEUTH, the model initially was run in prediction mode. The model was initialized with the latest urban extent map, the year 2015, and the growth coefficient values derived during the calibration phase. The target year was stopped in 2050 by performing 25 Monte Carlo trials using the same exclusion layer that was used for calibration. To evaluate the accuracy of the SLEUTH-3r forecast, a series of photos were taken in the years of 2015, 2017, and 2020 in three sites that showed a future infill development in SLEUTH predictions. These sites were in Idoon district which is one of the highly affected districts by urbanization in GIM.

Land Cover Changes in GIM (1972-2050)
For the 1972-2050 land cover change results, change trends and rates at the municipal and district levels were analyzed in two ways. (1) The areal coverages and change extents of land cover types were calculated for five periods of 1972-1980, 1980-1990, 1990-2000, 2000-2015

SLEUTH-3r Results
The study found the DM value of 0.005 sufficient to capture dispersive growth around Irbid. The high score of the road gravity parameter demonstrates that historic urban growth in GIM was driven by a transportation network. Also, the high spread parameter reflected a high probability of urban growth outward the existing urban clusters. The final calibration correlations for compare_r 2 and cluster_r 2 were high indicating that the model reflected the evolution of urbanization around Irbid. Comparison of the simulated urbanization against the actual urbanization of the control years reflected in the compared_score show that the model reflected the evolution of urbanization in the study area (a score of 0.98). The final calibration correlation was 0.84 in the case of the cluster_r 2 ( Table 4). The forecasted urban growth in GIM shows a continuation and intensification of development along transportation corridors (Amman-Irbid road and Petra road) to the south and east.
The accuracy assessment for the SLEUTH-3r predictions showed urban expansion (mostly infill development pattern) in three sites in Idoon district. The series of photos that were taken for each site in the years 2015, 2017, and 2020 proved the prediction results ( Figure 4). The photos timeline reflects the present and most likely the future dominant housing type (residential multi-story buildings) in GIM offering tens of apartments in each building to meet the needs of the increasing urban dwellers.    (1980, 1990, 2000, and 2015) and the SLEUTH-3r forecasts to 2050 allowed for the development of a long-term land cover timeline for the second-largest populated city in Jordan ( Figure 5). This timeline featured remarkable landscape transformation driven by significant urban growth across the 78 years. During the study period, Irbid City grew from a concentrated localized city surrounded by scattered low-intensity rural communities and vast primary cropland into a major metropolitan area in Jordan. Overall, built-up areas extent increased across the 78 years, with a total increase of 77.8 km 2 or 72.7%. (Table 5). According to the results during the mapping period, the total built-up areas coverage was 29.2 km 2 (9% of relative coverage) in 1972 and 100 km 2 (31% of relative coverage) in 2015, and will likely increase to 107 km 2 (33% of relative coverage) by 2050. Although persistently increasing, built-up areas change across the 78 years was not uniform in either space or time. Spatially, the southern and eastern parts of GIM area exhibited urban growth faster than the northern or western parts. Temporally, the rate of increase showed a decline across the study period, with the third and fourth intervals of 1990-2000 and 2000-2015 having the highest rate of built-up areas expansion at 68.6 and 41.4%, respectively. However, the rate of increase will likely decline by 2050 at 6.6% reflecting the impacts of the economic and housing market recessions. The total 2015 GIM coverage of agricultural land was 194.6 km 2 and will likely decline to 188 km 2 by 2050. In 1972, the total extent of agricultural land was 178.2 km 2 representing a slight 9.8 km 2 expansion over 78 years. However, the rate of agricultural land change was not consistent across the fife intervals. The first three intervals of 1972-1980, 1980-1990, and 1990-2000 exhibited an increase in agriculture coverage, with the second interval of 1980-1990 having the highest rate of agricultural land increase at 10.7% (Table 5). The increase rate declined to 3% for the following period of 1990-2000. Afterward, the total agriculture extent in GIM steadily declined, especially during the 2000-2015 period, with a total loss of 12.4 km 2 at a decrease rate of 6%. This loss of agricultural land will likely continue to the future at a slower rate of 3.4%. The sparsely vegetated/barren coverage in GIM dramatically declined with a total loss of about 87.6 km 2 ( Table 5). According to the study results, the total sparsely vegetated/barre extent was 116.

Urban Development Change at District Level (1972-2050)
The previous analysis of land cover change in GIM showed a persistent and permanent urban increase in the study area over 78 years. However, the most urban change occurred after 1990 and was not spatially even across the study area. Therefore, urban change was examined at the district-level to identify where exactly and by how much urban development has transformed the landscape of GIM's districts in the past decades and will likely do in the next coming decades. Figure 6 illustrates proportions of built-up areas extend to the district area in 1972 (a) 2000 (b), 2015 (c), and 2050 (d). In general, districts within the main city limits (i.e., Hashmyah, Nozha, and Rawdah) and peripheral districts (i.e., Boshra, Idoon, and Bait Ras) exhibited larger proportions of built-up areas than the rural-dominant districts in the farthest north and south (i.e., Al'al and Nu'aymah). For example, the built-up area proportion in Hashmyah and Nozha districts increased from 75% and 6% in 1972, respectively, to 100% in 2000. To the east, the peripheral districts of Idoon and Boshra exhibited a substantial built-up expansion with proportion increased respectively from 0.6% and 0.9% in 1972 to 15% and 14.4% in 2000, and to 53.1% and 48.9% in 2015. By 2050, the built-up proportion will likely increase to 56.9% in Idoon and to 55.6% in Boshra.  In the other 18 districts, however, urban increase rates are expected to decline. This declining rate has been denoted globally and can be attributed to the global recession, especially during and following the current COVID-19 pandemic. Rawdah. This decline in agriculture proportions will most likely continue to 2050, reflecting a slowdown in agricultural land consumption mainly by urbanization. Figure 9 illustrates the net change in total agricultural land in a district for the intervals of 1990-2000, 2000-2015, and 2015-2050 expressed as a percentage of the total agricultural land in the base year (1990, 2000, and 2015) of each interval. The results of the net change analysis of agriculture show that between 1990 and 2000, 9 districts exhibited a positive gain in agricultural land, with the districts of Kafer Jayez, Fo'ara, and Hawar having the highest gain rate at 54.5%, 51%, and 36%, respectively. The rest of the districts showed a negative loss of agricultural land, with the districts of Naser, Nozha, and Rawdah exhibiting the highest negative loss at 80%, 36%, and 35, respectively. From 2000 through 2015, the four districts of Naser, Nozha, Rawdah, and Hashmyah completely lost their agricultural areas for urbanization. The agricultural land consumption by urban growth expanded towards the peripheral districts, with Idoon having the highest negative loss of agricultural land at 51.5%. From 2015 to 2050, no change in agricultural land expected in the districts within the main city limits due to either a complete consumption of agricultural land as the case in Nozha, Naser, Hashmyah, and Rawdah or to the topographic constraints of urban development as the case in Barha, Manarah, and Rabiah districts. The negative loss of land in agriculture will likely continue in all districts, with Hakama having the highest negative loss at 21% followed by Boshra at 13.6% and Sareeh and Idoon both at 11%

Discussion
The previously presented results provide vital insights on the magnitude of land cover transformation and rate of agricultural land loss not only to Greater Irbid Municipality officials but also to Jordan Ministries of Environment, Agriculture, and local administration officials. These results question the efficiency and sustainability of decision-making in Jordan to ensure land resources and food security meanwhile introducing a new means of economic development crossing multiple aspects of an anthropogenic environment to help reach the 2030 UN Sustainable Development Goals (SDGs). Transitions towards adopting sustainable development models in decision-making can be challenging to a developing country such as Jordan. However, the movement should start immediately to demonstrate the commitment towards global eradicating poverty and achieving sustainable development by 2030. This study provides spatial evidence on past and future land transitions that can be of great value to local and national decision-makers to draw more sustainable planning policies based on realistic LULCCs.

Agriculture Implications on Food and Food Security
The LULC forecasts revealed that agricultural land will be the cover type that will suffer under urban expansion. The recession in agricultural land, as a response to urbanization, has been a worldwide problem hindering sustainable development strategies [68,69]. Generally, agricultural land that is adjacent to main urban areas mostly affected by irreversible urbanization. Such permanent replacement of agricultural land has negative consequences on arable lands and hence food production and food security. According to World Bank data, Jordan's arable land declined from 3.4% to 2.6% in 2015 [68]. Besides, arable land per person declined from 0.1 ha/cap in 1980 to 0.0 in 2015. The negative consequences of the paradoxical competition between land for settlement and agricultural land for food. Such juxtaposition becomes more problematic when occurring in a limited resources country with agricultural land only accounts for 12% of its total areas [69]. Despite this steady decline of land in agriculture, the production increased in absolute terms. The absolute growth in agricultural products is mainly for cash export crops such as olives and fruit. In the late 1980s, Jordan's main agricultural development strategies were directed towards maximizing the production of irrigated agricultural lands in the Jordan River Valley and the arid region. Although between 1975 and 1980s cereal production reached a striking peak with an increased rate of 150%, current and future production are expected to remain way far from meeting the domestic needs due to the highly fluctuated production from year-to-year resulting from frequent drought [70] and increasing agriculture abandonment [71][72][73].
Moreover, prices of food continue to rise globally; the FAO Food Price Index (FFPI) increased 2.1 points in August 2020 higher than its level in the same month in 2019 [72]. A reduction in cereal crop yields and wheat production is expected worldwide, with the current COVID-19 pandemic harshly affecting the top 50 food-exporting countries [72]. Although statistics showed that Jordan's agricultural output significantly rose from 2011 to 2016, with 2015 exhibiting the highest growth rate at 16% [73], Jordan continues to be a food-importing country, with 98% of primary food products (i.e., grains) are imported [74]. Besides, nearly 3% of Jordanian households are food insecure and 53% are vulnerable to food insecurity [74]. This alarming situation of Jordan's agricultural sector and food security crisis requires urgent actions to preserve the available agricultural land from uncontrolled urban growth and to enhance households' resilience and food security by (1) incorporating geospatial technologies to simulate past and forecast future LULCC trajectories and their potential impact on the agricultural sector, and (2) providing financial and governmental support for farmers in governorates with high-quality agricultural land (i.e., Irbid governorate).
The SLEUTH simulation results allow for a thorough investigation on the foreseen status of agriculture given a case study from the most fertile agricultural plains in Jordan. The recession of rain-fed cropland in GIM confirms the observed trend not only in Jordan but in the MENA region as well, with an anticipated decline in agricultural production by 21% at the end of the 21st century [75] The forecasts presented in this research were based on assumption that reflects the "business as usual (BAU)" scenario indicating no change in the exclusion layer (only included land open for urban development). Nevertheless, the results reflect reality offering a vital baseline to evaluate other alternative scenarios in GIM and other areas in Jordan. Incorporating multiple socioeconomic and environmental variables in the exclusion layer of SLEUTH-3r along with the use of self-modifications of growth rates proved to be promising in simulating future trends under different scenarios [28,67].

Implications for Urban Planning in GIM
Urban growth forecasts presented in this study assure the worldwide megatrend of urbanization and growth of metropolitan areas due to accelerating urban population growth [4]. The forecasted urbanization also aligns with the global urban sprawl pattern on the rural-urban fringe. The continued urban growth pattern in GIM is characterized by a low-density suburban growth (Figure 10a), with new urban cores. This urban growth can be attributed to multiple municipal, socioeconomic, demographic, and political variables that played a major role in increasing the urban population of GIM, especially after 2000 [59]. The number of populations, here, increased from 370,998 in 2000 to 535,199 in 2015 and expected to reach 1,219,055 by 2050 [65]. As a result, housing demand and supply greatly increased to meet the needs for housing for the new urban dwellers. Such housing demand was triggered by an increase in the average annual income of households which, according to the Household Expenditure and Income Surveys, nearly tripled from 1,703 JOD in 1992 to 3,995 JOD in 2013, respectively. Statistics showed that the percentage of apartment units in Irbid increased from 51% in 2006 to 64.7% in 2013, and 74.4% in 2018 [76]. This may explain the widespread growth in multi-story residential buildings and the construction of new roads in the suburbs (Figures 4 and 10). The nationwide increase in housing units proved the ability of the housing sector to respond to the accelerating urban growth, especially after the influx of Iraqi and Syrian refugees. It was estimated that between 2004 and 2015, 1.1 million new units were built, exceeding the projected units (300,000) by 800,000 units, contributing to an average of 7.2% of GDP [77]. After 2015, nevertheless, the housing sector and real estate market, in general, collapsed because of many global economic shocks (i.e., Iraqi War, 2008 financial crisis, and Syrian conflict) [77,78]. A substantial slowdown in Jordan's GDP growth occurred and consequently resulted in a decline in urban growth. The forecasted results from 2015 to 2050 confirm this slowdown as the urban growth rate drastically decreased from 41.1% between 2000 and 2015 to only 6.6% between 2015 and 2050. The projected built-up areas will reach 107 km 2 by 2050-or 12 km 2 below the officially zoned areas for urban development. This indicates that the available land for urban development is much more than what is needed and can accommodate future population growth of GIM beyond the year 2050 with no new land is consumed. With that said, however, the problem is that sufficient land might not be the right type for housing (i.e., agricultural land as the case in GIM), causing a structural shortage and mismatch between the types of residential land available on the market and the actual types needed. This structural shortage has led to reduce the land available for affordable housing and consequently increase the cost per housing unit [77]. This dilemma of shortage in land available for affordable housing might be driven by the inappropriate zoning categories of residential areas (A for luxurious housing, B for low-density higher-end housing, C and D for higher-density more affordable housing). Since the amalgamation process in 2001, new areas in GIM have been added for housing yet mostly zoned for the higher-end housing (Zones A and B) with large plot size and hence the least affordable type for housing. Besides the structural shortage of land available for development, land zoning in Irbid represents a mixture of urban and rural land uses with an absence of environmental, political, and social planning. City planning policies, generally, are drawn in favor of a group of competitors and upon expectations rather than demographic and socioeconomic statistics.

Conclusions
This study investigated past, present, and future land cover transformations at both the municipality and district levels of Greater Irbid Municipality featuring the application of SLEUTH for the first time in Jordan. The magnitude of past land cover changes and the foreseen changes were examined. The findings of this research suggest that GIM landscape composition exhibited a significant change in the past decades and will continue to experience more changes mainly associated with urbanization. The amount of agricultural land lost for urban development is of great value for sustainable development. Both fragmentations of agricultural lands and urbanization have been recognized as the main challenges in sustaining and developing the agricultural sector in Irbid. GIM is in urgent need of the adoption of a compact policy by which no more new lands are added for urban development. Infill development patterns should be enforced, and the striking land prices should be regulated.
While past urban growth in GIM was mainly driven by changing lifestyles and ruralurban migration, present growth is highly driven by the influx of refugees from neighboring countries, changing municipal policies, and socioeconomic transformations. The present and future low-intensity growth patterns are shaped by arbitrarily urban planning policies, absence of an environmental land value, land speculation, and striking land prices. GIM suffers from uncontrolled land uses as there are no hierarchy of zones and most areas are mixed-uses. Irbid city is unique in its location, physiographic settings, and agricultural value. Yet, urban planning strategies did not preserve this value and its environmental assets. Urban planning strategies must be drawn within a more sustainable development model reflecting an official commitment to attaining the SDGs goals and targets towards achieving "sustainable urbanization and capacity for participatory, integrated and sustainable human settlement planning and management in all countries" (SDG 11 (11.3)). City municipal boundaries should have been adjusted at several stages over time and based on reliable demographic and socioeconomic variables. For a balanced development, municipal authorities need tools to monitor how the land is currently used, assess future demand, and take actions to assure adequacy of future supply. This research promotes and supports sustainable development planning through producing spatially explicit forecasts proving the efficiency and capability of forecasting to not only investigating the implications for agriculture, food security, and urban planning but also to use the results as input to other models to measure consequences on, for instance, land degradation and water quality. Such information is of great importance for decision-making in Jordan where there is an urgent need for proper land management practices and sustainable land-use policies.