A State-of-the-Art Vegetation Map for Jordan: A New Tool for Conservation in a Biodiverse Country

In many countries, including Jordan, the updating of vegetation maps is required to aid in formulating development and management plans for agriculture, forest, and rangeland sectors. Remote sensing data contributes widely to vegetation mapping at different scales by providing multispectral information that can separate and identify different vegetation groups at reasonable accuracy and low cost. Here, we implemented state-of-the-art approaches to develop a vegetation map for Jordan, as an example of how such maps can be produced in regions of high vegetation complexity. Specifically, we used a reciprocal illumination technique that combines extensive ground data (640 vegetation inventory plots) and Sentinel-2 satellite images to produce a categorical vegetation map (scale 1:50,000). Supervised classification was used to translate the spectral characteristics into vegetation types, which were first delimited by the clustering analyses of species composition data from the plots. From the satellite image interpretation, two maps were created: an unsupervised land cover/land use map and a supervised map of present-day vegetation types, both consisting of 18 categories. These new maps should inform ecosystem management and conservation planning decisions in Jordan over the coming years.


Introduction
If natural resources are to be managed sustainably, it is crucial to assess land cover/land use and the spatial distribution of vegetation types. This is also helpful for protecting habitats and for providing data for vegetation cover modelling. Our goal here is to derive accurate maps and cover estimates for vegetation in Jordan through the use of both field study and satellite imagery. Jordan is a region of high vegetation complexity and provides an ideal testing ground for methods in land cover and vegetation type mapping.
Since Alexander von Humboldt first began categorising vegetation zones in 1807, biogeographers have classified the planet into chorological units. Prior to his work, it was typical to focus on individual plants; Humboldt shifted this towards a more collective focus on species-their distribution, growth forms, and how they relate to their surroundings [1]. This new approach meant that communities and associations among plants were able to be identified [2]. August Grisebach utilised this new approach to develop what would become known as phytogeographic regions, including the Mediterranean and Irano-Turanian regions, according to the species composition and physiognomy of individual vegetation groups [2,3].
Expanding on the work of Long [9], Al-Eisawi [10] determined that there are nine bioclimatic areas in Jordan, which are categorised into four overarching clusters: Saharan Mediterranean (cool, warm, very warm); arid Mediterranean (cool, warm, very warm); semi-arid Mediterranean (cool, warm); and sub-humid Mediterranean. Later, Al-Eisawi [11] also created a map based on the delineations made by Zohary [3] of the regions of vegetation in the Middle East.
The flora in the region under study is very species rich due to the area being situated at the meeting point of three phytogeographic regions [3,[11][12][13]. The land changes from Mediterranean maquis and forests to grassland, steppes, and extremely arid desert biomes within a relatively short space. While edaphic influences do have an effect on vegetation, it is the climate that has a greater impact [13]. As well as vegetation being influenced by environmental pressures, human activities also affect the way plants are distributed. Some research has accounted for forestry, agriculture, and grazing potential along with human activity when classifying vegetation [9,[14][15][16].
The vegetation types of the region have also been mapped. While earlier work (e.g., [3,5,6]) was conducted prior to modern development and experienced fewer restrictions on political boundaries, other modern works were either restricted by political boundaries or the study was conducted across limited borders based on extensive fieldwork [11,14,16,17]. The most recent and updated map of vegetation for Jordan was produced by Albert et al. [14]. The map included 19 vegetation types and was primarily based on topography and coarse-resolution satellite images to describe vegetation cover of the country at a scale of 1:1 M. However, no map to date has combined state-of-the-art earth observation approaches with extensive field data to ground-truth remote sensing data.
Remote sensing data have contributed to land cover and vegetation type mapping by applying advanced digital classification techniques to derive both land cover and vegetation types at different spatial resolutions [18]. For vegetation mapping, the outputs from digital classification of remotely sensed data are intensively used to map land cover dynamics [19], floristic diversity, phytogeography [20], and vegetation [21] over large geographic areas. In order to identify and map land cover and vegetation in a continuous manner, remote sensing is vital. Data from sensors such as Landsat, MODIS, Sentinel, SPOT, and ASTER can be obtained cheaply or free of charge, with some providing wide swath-width images (over 290 km).
Vegetation mapping of Jordan is now urgently required so that appropriate and effective planning can be made for managing agriculture, afforestation, and rangeland, and for the degradation of habitats to be reduced. Prior research has been less than ideal for a multitude of reasons; for instance, studies did not utilise satellite imaging, were not comprehensive enough, or did not include field observations to support mapping and definition of the areas under investigation. Further, much of the existing research fails to include data from maps based on land use, geological or meteorological information, vegetation sampling, or GIS methodologies, and often does not utilise herbarium specimens or pictures to identify the species found in each region.
Taking all of the abovementioned points into consideration, it is clear that no study currently exist that describes Jordan's vegetation in a spatial context for use in mapping or restoring particular kinds of vegetation. This would be extremely useful in the management and conservation of the country's biodiversity [22]. Hence, the aim of the current study is to map and classify the vegetation cover of Jordan using a combination of the latest satellite data, extensive ground data, and state-of-the-art statistical approaches.

Processing of Remote Sensing Data
For the main floral groups to be displayed and a map of the vegetation to be produced (1:50,000 scale), it was necessary to define the types of vegetation present by utilising a reciprocal illumination approach between on-the-ground sampling and satellite imagery, as illustrated in Figure 2. The first step was to examine Sentinel-2 satellite images of Jordan to assess how well they correlate with known land cover classes, using Al-Bakri et al.'s land cover map [23]. Then, a more up-to-date land cover/land use map was created. The third step was more specialised and involved assessing how compatible the satellite images were with existing maps of the vegetation. Finally, decisions were made regarding where field sampling and ground-based verification were needed. The study of sufficient samples was confirmed within the known vegetation types in satellite images, and a minimum of 30 sample areas were collected for each hypothesised class of the produced land cover/land use map. A digital elevation model (DEM) was used to represent the bare ground topographic surface of the earth excluding trees, buildings, and any other surface objects. The information from the geological map was useful in this area, as the map indicated shallow soil depth in this area as well as the presence of a marble layer close to the surface. This shallow soil has a high salt content as these areas receive marginal rainfall. The climatic data (i.e., isohyets of rainfall and temperature) were very useful to classify some types of vegetation (i.e., thermophilous vegetation in Jordan Valley).
Several data sources can be used for obtaining remotely sensed data. The current study utilised the United States Geological Service (USGS) EarthExplorer website (https: //earthexplorer.usgs.gov/ (accessed on 30 October 2017)) to obtain Sentinel-2A images and prepare them for ArcMap usage. Eighteen cloud-free Sentinel-2A images were selected for use in this study. Individual bands of Sentinel-2 images were examined and used to create a natural colour composite, which is an important step to understand how to combine these bands for use in visualisation and analysis. Other spectral bands contain information that is very useful for visualisation and is used in later stages (e.g., false colour combinations were used in this study to identify urban area boundaries). The Sentinel-2 Multi-Spectral Instrument (MSI) has a span of 13 spectral bands from the visible and the near-infrared to the shortwave infrared at different spatial resolutions ranging from 10 to 60 metres on the ground. Of all those bands, only the four bands at 10 m resolution that gave information about vegetation cover were composited and used: bands 2, 3, 4, and 8. Images taken from Sentinel-2 were pre-processed, mosaicked, and clipped so that the images were calibrated. The images were pre-processed so that land cover/land use could be interpreted visually. Two techniques can be used to classify vegetation from images: supervised and unsupervised, with each method selected according to whether or not ground data are included. In supervised classification, training data collected from ground observations are used to initiate spectral signatures of the land cover types. The method has the advantage of the final map being comprised of known classes. Supervised classification approaches tend to give more accurate results in homogenous landscapes and when enough training data are available. Unsupervised classification methods have the advantage of separating spectral classes without prior knowledge of the area and without training data. Therefore, when mapping thematically from imagery, unsupervised methods are commonly employed [24].
In this study, unsupervised classification was used for mapping the land cover/land use of the country, while supervised classification was carried out to derive a vegetation map based on vegetation type clusters obtained from ground surveys of the species composition of vegetation. In this hybrid approach, classification accuracy was improved using both auxiliary data and a wide range of expert knowledge [25,26]. Hybrid supervised and unsupervised classification methods were adopted, along with the field experience of the first author, so that a greater degree of accuracy could be obtained within a shorter timeframe. Further, ArcGIS 10.5.1 software and the Sentinel Toolbox Sen2cor plugin were utilised for classification of images taken from the satellite and for simplifying the spectral data.

Using Unsupervised Classification to Produce Land Use/Land Cover Map
An unsupervised technique was employed to produce a baseline map from which vegetation sampling sites could be located. This was also used to update the land cover/land use map for Jordan. Ancillary datasets were used to develop an initial overview of the area under investigation, and then land cover/land use maps from the Jordan Land Cover Atlas [27], Al-Bakri et al. [23], and geological and soil maps from the Ministry of Agriculture of Jordan were also consulted [28]. Unsupervised classifications are determined by automatic software calculation of pixel similarity. For instance, forest pixels can be differentiated from water pixels by automatically grouping these pixels based on their spectral characteristics.
First, the visible spectrum (RGB) was utilised to visualise the landscape in natural colour. Figure 3 shows the resulting images. However, it was necessary to assess the full spectral capabilities of Sentinel-2 to acquire a better understanding of the nature of the area under study through the different colour composites of the multispectral data. Therefore, the images of Sentinel-2 were displayed in both true colour (TCC) and false colour (FCC) composites. Both combinations showed different characteristic of vegetation and land cover. The TCC reflected the variations in soils and rocks (Figure 3), while the FCC reflected the variations in natural vegetation and agricultural areas ( Figure 4). In addition, spectral indices were produced for the study. Such indices are used for highlighting the landscape's geological features, bodies of water, and vegetation. Here, the data were cut down and changed into more meaningful information using ratios between bands to transform the spectral data: NDVI = (NIR-Red)/(NIR+Red).
Typically, many classes are identified during analysis, and these are then honed and grouped more precisely to obtain more meaningful classes of features. We did this using the unsupervised method known as the Iterative Self-Organising Data Analysis Technique (ISODATA). The ISODATA method was preferred over the K-means method as splitting and merging of clusters was possible, so that the final number of classes would not be limited to the number identified by the trainer [18]. In this study, unsupervised classification was initiated with 30 spectral clusters, and the image pixels were then clustered into groups using ISODATA. Ground reference data and researchers' expertise were then used to judge some of the clusters before a thematic urban land cover/land use map was produced. The classes were interpreted using two methods. The first involved interpreting images visually using field experience and feature attributes from high-resolution images taken from Google Earth. The second method was an automatic process using the vegetation index thresholds of feature attributes in each polygon. A visual interpretation was used for most of the mapping and the final results were checked using global positioning systems (GPS), topographic maps, and a number of field visits.

Field Data
We sampled the vegetation in order to validate vegetation types in the field. GPS data were gathered from the centre of homogeneous areas with respect to vegetation. Field surveys were conducted based on the previously produced land cover/land use map. A plan was made to sample 640 (50 × 50 m) plots during autumn 2017, spring and autumn 2018, and spring 2019. For each land cover class, approximately 30 sample areas were gathered [29], sometimes including multiple vegetation types. The percentage canopy cover for each layer (tree, shrub, and herbaceous layers) found in the sampling units was recorded. For all trees, shrubs, and perennial herbaceous plant species, the percentage of cover across the entire plot was estimated.
All of the vascular plants that could be identified to the sub-species level were listed for each of the vegetation classification plots. To accomplish this task, the team relied upon the experience of the first author, gained from his work as a botanist at the Royal Botanic Garden of Jordan for nearly 15 years. The Angiosperm Phylogeny Group (APG IV 2016) classification was followed along with a checklist of Jordan's vascular plants [30].

Statistical Analysis
Hierarchical clustering of field plots was executed based on the percent cover of each plant species. Annual species were excluded, because presence and percent cover could not be reliably recorded, as some plots were visited in different seasons. A dissimilarity matrix [31] was produced using the Bray-Curtis floristic distance approach, an index that helps to determine the differences in species' relative abundance profiles among sites [32]. Ward's hierarchical clustering algorithms were used to maximise the coherence of clusters, as per Murtagh and Legendre [33]. A number of approaches have been suggested to determine the optimal number of clusters including Silhouette, Elbow, and k-means [34,35]. The problem is that most apply arbitrary thresholds that do not account for knowledge of the vegetation. Thus, here, the appropriate number of clusters was estimated based on biological and field knowledge of the sites and vegetation. Analyses were carried out with the 'vegan' and 'recluster' packages in the Rstatistical environment [36][37][38].
For each identified habitat, an indicator species list was developed based on indicator values (IndVals) [39]. Two components make up IndVals, namely fidelity and specificity (percentage of sites containing the target species, and the probability of a species being abundant in a particular habitat, respectively) [39,40]. The species considered significant had IndVals indicators over 0.2, and a p-value lower than 0.05 for each particular habitat [39,41]. The 'indicspecies' R package was used for analysis [42].

Using Supervised Classification to Produce the Vegetation Map
The supervised classification method requires input training data for each class created by the researcher to calculate the classified image. Relying mainly on land cover classes, vegetation cover was classified using a combination system between physiognomic and phytosociological criteria [43]. All of the polygons/map units were given a vegetation type category according to the clustering analyses and training data, and this was revised and, if necessary, reassigned once the field data were analysed. In the interpretation phase, mosaics were reported, and the percentage of each type of vegetation found within the map unit boundary was provided.
In addition to mapping vegetation types, the boundaries of phytogeographical regions were determined using the approach described by White and Léonard [7]. A map of Jordan's phytogeographical regions was developed using a combination of remotely sensed images and field data to delineate the biogeographical boundaries.

Accuracy Assessment of Vegetation Maps
There were four key steps taken to assess the accuracy of our maps: (i) visually inspecting the maps, (ii) comparing the classes of the thematic maps, (iii) using accuracy metrics based on comparing the class labels from ground data and thematic maps, and (iv) using an error or confusion matrix [44]. The technique most often utilised for assessing accuracy is the confusion matrix. Here, the interpreted classification is compared to the referenced classification that is known to be correct [45,46].
As per [47], a Kappa index was worked out for all typology levels by crossing the classification raster layers and validation points. Kappa = 1 when classification is perfect; Kappa > 0 when the observed correct proportion is greater than the proportion expected by chance; Kappa < 0 when the observed correct proportion is less than the proportion expected by chance [47]. If a classification has a Kappa index over 0.8, this means the classification accuracy is very good. An index of 0.6-0.8 means the accuracy is good; anything lower than 0.6 indicates a low classification accuracy [48].
The accuracy assessment was undertaken using random stratified ground-truth 'square samples' validation. Squared samples of 2500 m 2 (the size of our field plots) were generated. Overall accuracy was calculated using the percentage of land cover that was classified correctly (sum of the correct classifications (diagonal elements) divided by the number of samples) [48].
ArcGIS 10.5.1 software was used to generate the reference maps from ground observations. The use of this approach enabled the agreement between the ground truth and the current land cover raster to be reflected through the use of four measures of accuracy: the user's accuracy, the producer's accuracy, the overall accuracy, and the Kappa coefficient.

Land Cover/Land Use Map
The land cover map is the result of the integration of remotely sensed data with thematic features from land cover models. According to this interpretation of satellite images, the map shows 18 classes of land cover ( Figure 5, Table 1).

Vegetation Map
Sixteen groups, or vegetation type classifications, were produced from hierarchical cluster analysis ( Figure 6). Sampling sites were grouped together based on the composition of the perennial vegetation present in each site. Two overarching vegetation groups emerged from the cluster analysis (Table 2).  From the species composition analysis, 54 indicator species were identified that could then be used to classify species of vegetation. Annual plant species were not included in the cluster analyses, but one vegetation type, specifically Gravel Hammada vegetation, lacks perennials. For this type, annual indicator species were added based on researcher experience and field records.

Using Remote Sensing to Classify Vegetation Types
Based on the land cover/land use map and as a result of using the supervised classification method in interpreting the satellite images, a current vegetation map was produced at a scale of 1:50,000 that contains 18 vegetation types (Figure 7, Table 3).
Two vegetation types-anthropogenic pine forest and gravel hammada vegetationwere added to the map, although they were not included in the cluster analysis. The reason for this is that there were no perennial indicator species for the gravel hammada vegetation, and there were no unique indicator species that distinguished both the natural and anthropogenic pine forest. According to the interpretations of the satellite images, there are 68 areas in which riparian habitats occur, mostly along tributaries of the Jordan River and Yarmouk River or wadis flowing towards the Dead Sea.
The Mediterranean region is restricted to the Highlands extending from Irbid in the north to Ras An-Naqab in the south, in addition to some isolated representation in the southern mountains. The Irano-Turanian zone separates the Mediterranean from the Saharo-Sindian-Arabian regional subzone on one side, and the Mediterranean from the Saharo-Sindian-Nubo-Sindian subzone on the other side [7]. The Saharo-Sindian-Arabian regional subzone is located at the eastern desert and accounts for the greatest proportion of land cover in Jordan-about 72%. The Saharo-Sindian-Nubo-Sindian subzone penetrates into Jordan, following Red Sea coastal areas, from tropical regions of the Arabian Peninsula and northeast Africa. A summary of each bio-geographic region and their associated vegetation types is presented in Table 4.

Accuracy Assessment of Vegetation Maps
In this study, four accuracy measures of each typology of the maps were calculated: the user's accuracy, the producer's accuracy, the overall accuracy, and the Kappa coefficient. These measurements were extracted using the specific confusion matrices listed in Appendix A. The overall accuracy and Kappa coefficient statistics were used to examine the quality of land cover and vegetation maps as a whole. The Kappa coefficient was used to determine the agreement between the map and reality, where 0 = complete disagreement and 1 = complete agreement.
All vegetation maps previously produced by experts and reviewed in this study were of poor to moderate accuracy, where the confusion matrices showed a diagonal line that illustrates a weak to medium agreement between the layers (land cover layer and groundtruth points). The classification accuracy of these maps indicated by the Kappa index also showed a low classification accuracy. On the other hand, the results of the analysis showed high overall accuracy for both the unsupervised land cover/land use map and the supervised vegetation map. The Kappa index also showed a high agreement for these two produced maps, as the result was more than 80%, indicating very good classification accuracy. Table 5 below summarises the accuracy and Kappa coefficient for all maps.

Land Cover/Land Use Map
The land cover/land use map of Jordan produced here includes detailed information pertaining to agricultural activity and urbanisation, which can be used to address the issues related to the country's increasing population. Overall, the accuracy of the land cover/land use map was relatively high and enabled the derivation of a vegetation map, after merging ground data sources. The land cover map was in agreement with other maps produced for Jordan [23,27] and with global land cover maps [49], particularly for the classes of dry lands, arable land, forests, and water bodies. For other classes, the differences in classification could be attributed to the type of remote sensing data and the classification scheme used in this study, which was intended for derivation of a vegetation map. The output map was produced as a series of digital layers, which were then overlain and analysed using a GIS in order to obtain land cover and land use percentage information. The map shows that 5.4% of the land is arable land, and dry lands account for 77.9% of the country. Highly productive rangelands are represented by the steppes, desert scrubs, and wadis located in the eastern desert, and constitute 14.1% of the total area. These findings can be used to update the percentages given in the Atlas of Jordan [27].

Vegetation Map
The clustering analyses delimited 16 vegetation types, all of which were included in the four older vegetation maps [11,14,16,17], with the exception of three vegetation types that were added during this study, namely runoff hammada vegetation, sandy gravel hammada vegetation with Haloxylon scoparia, and sandy gravel hammada vegetation with Vachellia gerrardii and Artemisia judaica. Al-Eisawi [11] indicated that the hammada is divided into four types, but he did not classify them on his map. The 16 vegetation types can be largely separated into two higher-level groups: mesic vegetation types growing in the highlands, and xeric vegetation growing in the eastern desert and the Jordan Valley [13,50]. Juniper forest and steppe vegetation form their own group distinct group, however, which represents a transition between the mesic and xeric areas [51]. While juniper is considered a Mediterranean species, it is also found in habitats alongside steppes shrubs.
Within the higher-level 'dry group', sand dune and acacia woodland, which grow in Wadi Araba and Wadi Rum, cluster together as their ground flora is compositionally similar; the shrub stratum dominates sand dune vegetation, and the tree stratum dominates acacia types [11,52]. Granite and sandstone scrubland form a group, likely because they are similar in having rock crevices that allow for water accumulation and can therefore support the germination of Mediterranean species typically growing in areas with higher rainfall [17]. Sandy gravel hammada vegetation with Vachellia gerrardii and Artemisia judaica, sandy gravel hammada vegetation with Hammada scoparia, and mudflat vegetation are clustered together because there is considerable overlap between these types of vegetation [11]. They also occur together in the eastern desert. The last sub-group within the higher-level dry group contains vegetation that is edaphically influenced such as riparian vegetation, saline vegetation, and thermophilous vegetation. All of these vegetation types occur in the Rift Valley, and frequently exhibit similar bush and herbaceous cover [12].
We mapped the distribution of these 16 vegetation types, plus two additional types (gravel hammada and cultivated pine forest), which we included based on remotely sensed information. These additional types of vegetation were only able to be determined through the use of high-quality satellite images, in which these latter two vegetation types were clearly visible along contour lines.
In our final map, the locations of Quercus ithaburensis-dominated deciduous oak forests were identical to those posited by both Al-Eisawi and Kasapligil, but not to locations outlined in Albert et al. or Danin [11,14,16,17]. In Albert et al.'s study, coverage of the Wadi Essir region is indicated to be wholly deciduous oak forest [14], while Danin showed no forested areas at all in either Wadi Essir or Al-Alouk [17]. Neither of these studies reflects the true coverage, which in Wadi Essir involves two oak species: deciduous oak Quercus ithaburensis and evergreen oak Quercus coccifera. In addition, the deciduous oak Quercus ithaburensis forest covers the Al-Alouk region. On the other hand, the findings of the current study support those of Kasapligil that suggested trees of Pistacia atlantica and Ceratonia siliqua are present [16]. As would be expected, the supervised current vegetation map gives a more accurate result than that obtained from the unsupervised map. Jordan's forested regions account for 0.71% of the total land area, with acacia woodland and forests combined totalling 2.3%. These figures can be considered an update to the estimates put forth by the Ministry of Agriculture [53].
All of the prior vegetation maps represent the Quercus coccifera-dominated evergreen oak forests fairly well, although there is substantive variation in their accuracy [11,14,16,17]. Evergreen oak constitutes the densest forest in the country, with some areas achieving 95% canopy cover. This forest is also extremely important in terms of providing shelter for biodiversity; indeed, Zohary noted the richness of the habitat for many different species of flora and asserted it to be characteristic of the east Mediterranean maquis [54]. The current study found Quercus coccifera trees in areas neighbouring Quercus ithaburensis deciduous oak forests, with areas where the two species coexist in Wadi Essir and to the west of Ajloun-a finding that correlates to those of both Kasapligil and Al-Eisawi [11,16].
Pine forests account for 18.29 km 2 of Jordan's land cover. The Ministry of Agriculture report an area greater than this, but this is because they have a tendency to also include cultivated pine forests in their estimations. Pine trees are mainly found in Zay and Dibeen, but this study found some small pine forest areas in Jarash and Ajloun, lending further support to prior research [3,55]. These forests are found anywhere between 500 masl and 1000 masl, contradicting Al-Eisawi's assertion that they only occur at elevations greater than 700 masl [11]. The Dibeen Forest Reserve is home to approximately 46% of the country's total cover of natural pine forest.
The study's findings did not correlate fully with existing perspectives regarding what constitutes the climax in Mediterranean forests, as natural regeneration takes place in both evergreen oak and pine forests. One argument put forth by Atkinson and Beaumont is that Pinus halepensis forests represent the climax [55], while others such as Liphschitz and Biger believe that it is either Quercus coccifera or a combination of Quercus coccifera and Pistacia palaestina [56]. However, this study found that, sometimes, Quercus coccifera replaces Pinus halepensis in places where the primary vegetation has been degraded, and this is in agreement with prior research by Al-Eisawi [11].
In support of prior research (i.e., [11]), we found that herbaceous coverage in natural pine forests is far greater than in anthropogenic pine forests. Coverage can reach up to 60% in natural pine environments, yet in cultivated forests only around 20% is recorded. Al-Eisawi reported that the natural vegetation that grows under Pinus halepensis trees on yellow rendzina soil is not present under the same trees cultivated in terra rosa [11]. The reason for this is that rendzina soil is alkaline, and pine needles are acidic; thus, when the pine needles fall onto the soil they neutralise it, improving its quality. Conversely, red soil is acidic, so when the needles fall here the condition of the soil worsens, making it less favourable for many types of vegetation.
All of the previous vegetation maps of Jordan illustrate well how juniper forests are distributed [11,14,16,17]. Since the regions of Petra and Dana experience a higher annual mean temperature range than other regions in the mountains, juniper forests occur in habitats that are more steppe-like, with shrubs such as Noaea mucronata and Artemisia sieberi growing over an Irano-Turanian understory [17].
Garrigue and Batha regions are Mediterranean areas, apart from the cultivated and forest areas [11], and their distribution matches that observed in previous studies.
The steppe region does not appear as one large, uniform area, but rather as a meeting place where the Arabian regional and Mediterranean phytogeographical regions come together. All of the remaining steppe forests are located in the south of Jordan at Beer Ad-Dabghat and Beer Khdad, the areas neighbouring the eastern border of the juniper forest in the Irano-Turanian region. The trees found in these forests are typically historical Pistacia atlantica, and it is thought that these trees have existed here for at least a thousand years. Among the pistachios are Rhamnus disperma and Prunus korshinskyi, although in much smaller numbers. It is believed these species are a relict species from the Mediterranean [17]. Trees of Pistacia atlantica are also found in the eastern desert at Wadi Al-Botom, accompanied by Retama raetam and Prunus arabica shrubs [11]. The growth of these plant species in the desert may be due to the presence of some types of soil or rocks that may provide the necessary moisture requirements for species that are usually found in more humid areas [57].
Granite-sandstone regions are full of fissures, meaning that water infiltration after rainfall is prevalent. The water is able to permeate the rock fissures and reach any soil there. As the water is protected between and under the rocks, it is less likely to evaporate in the heat. This environment is suitable for endemic and rare species to flourish, even in the hotter desert regions. One example is a shrub endemic to Jordan, Daphne mucronata subsp. linearifolia.
Coverage of riparian vegetation in the canyons, rivers, and wadis differs depending on the edaphic and climatic conditions. In the north, trees of Platanus orientalis, Populus euphratica, Tamarix palaestina, and Salix spp. are associated with Phragmites australis, Nerium oleander, and Arundo donax shrubs. In the south, Populus euphratica, Dalbergia sissoo, and Phoenix dactylifera trees are associated with the same shrubs plus Myrtus communis. On the banks of the Jordan River, trees of Populus euphratica populate the area closest to the water [3,6]. Many riparian vegetation sites were able to be identified from the satellite images, and a few of them appeared in the map produced by Al-Eisawi [11].
There are concentrated enclaves of thermophilous vegetation surrounding the Dead Sea. These enclaves were referred to by Zohary as 'Sudano-Decanian enclaves' [3,13], while Al-Eisawi named them 'Sudanian tropical vegetation' [11]. Degradation and fragmentation of this vegetation type is widespread due to the expansion of farming; this was first reported by Al-Eisawi and the current study reinforces this assertion [11]. Alongside the thermophilous vegetation is saline vegetation close to the Dead Sea and along the Jordan Valley, as well as in some parts of the eastern desert. Tamarix trees are able to tolerate high salt levels, while for shrubs, Atriplex spp. can withstand saline soils [11]. Herbaceous vegetation is particularly salt intolerant, however, and may not grow here at all.
Sand dune vegetation is very sparse, and generally trees cannot grow in pure sand. Some botanists delineated this vegetation growing in independent areas surrounded by different vegetation types in Wadi Rum and Wadi Araba in the south [11,14,16], whereas Danin considered it connected to desert savannoid vegetation in the south and extreme desert in the east [17]. Large sand dunes are frequently much better vegetated than the surrounding plains or sand flats. The huge sand sheet areas are populated with Hammada salicornica, Haloxylon persicum, Caroxylon tetrandrum, and Anabasis articulata. The only restricted plant of all the species in this environment is Haloxylon persicum, which can form quite dense shrubland on dunes [17].
All of the previously produced vegetation maps of Jordan have included acacia woodland. Kasapligil refers to it as "scattered acacia grasslands", Al-Eisawi names it "acacia and Sudanian rocky vegetation", Danin calls it "desert savannoid vegetation", and Albert et al. refer to it as "acacia woodland" [11,14,16,17]. This type of vegetation is represented by sporadic Vachellia species as well as other desert semi shrubs like Anabasis articulata, Hammada salicornica, and Retama raetam.
There are five types of vegetation that make up the eastern desert, namely runoff hammada, gravel hammada, mudflat, and two types of sandy gravel hammada-one dominated by Hammada scoparia and another dominated by Vachellia gerrardii and Artemisia judaica. While mudflats frequently occur in the Badia, only one is found in Wadi Araba. There is typically no vegetation growing on the mudflats, but some species can be found around its periphery [11].
Most of the eastern desert is covered with gravel hammada. In this type of vegetation, only desert annuals generally grow. Desert shrubs belonging to the families Asteraceae and Chenopodiaceae occur sporadically next to and within the wadi beds, forming Runoff Hammada Vegetation. There are two sections of the Badia that have been identified as having soil that consists of both sand and gravel: one in the south, along the border with Saudi Arabia; the other in the northeast, close to the Iraqi-Saudi border. In the south, Vachellia gerrardii trees and Artemisia judaica shrubs are dominant, while in the northeast, shrubs of Hammada scoparia are the most commonly occurring species.

Phytogeographical Regions
Jordan is in a unique location geographically, where four separate phytogeographical zones meet. The Mediterranean region appears to be the least difficult to delineate out of all the zones. The Mediterranean region is home to the vast majority of the country's population, and subsequently is the area that has experienced the most severe effects from human activities [53]. Here, disparate vegetation types are found intersecting with each other, such as garrigue, batha, maquis, and forests, which is related to vegetation degradation. The maquis vegetation type is characterised by Quercus coccifera-Pistacia palaestina, likely preceding climax forest in terms of vegetation succession. Formations of forest and maquis are visible at Wadi Essir and along the Ajloun mountains. Formations of garrigue and batha occur in degraded areas of the Mediterranean region that do not support forests under current physical and climatic conditions, such as many parts of Amman and Irbid. These regions are represented by Sarcopoterium spinosum, among other small shrubs [7,16].
In the Irano-Turanian region, the climax vegetation lacks substantive tree cover, but some tree species such as Pistacia atlantica and Ziziphus lotus do occur [11]. Tafila's oak and juniper forests and Shoubak's open forests of Pistacia atlantica are surrounded by steppes. Trees are uncommon in this habitat though, and Danin posits that these two areas with trees are likely to be relics from a climate that once had more rain [17].
In the Saharo-Sindian-Arabian regional subzone, the availability of water is the main determinant of the type of vegetation found here-often from wadis. The natural vegetation of the Arabian regional subzone has been destroyed in many areas due to ploughing and other agricultural practices for cultivating barley, particularly in runoff areas and wadi beds [11].
Finally, the vegetation of the Saharo-Sindian-Nubo-Sindian local centre of endemism is characterised by species such as Vachellia gerrardii and Ziziphus nummularia [7]. These species are at their northernmost limits of their distribution and represent the northernmost limits of tropical African vegetation in southwest Arabia.

Conclusions
Sentinel-2 imaging and extensive field work were utilised to produce two maps of Jordan. The first map is an unsupervised land cover/land use map including 18 distinct classes, and the second is a supervised present-day vegetation map, also including 18 classes of vegetation, which partially overlap with the 18 landcover classes. Remote sensing technology enabled the identification of two types of vegetation that were not represented in the hierarchical cluster analysis based on species composition, namely anthropogenic pine forest and gravel hammada vegetation. The hierarchical cluster analysis suggested that the vegetation was categorised into two major groups: mesic and xeric. The diversity assessment showed that the most species-rich plant communities are found in the evergreen oak forests, with mudflats representing the least species-rich habitat. Analysis of the species composition of the sampling sites enabled the production of a list containing 54 indicator species, which could then be used to identify the type of vegetation at currently unsampled sites. The country is uniquely located at the point where three continents converge, resulting in the presence of four phytogeographical regions: Mediterranean, Arabian regional subzone, Irano-Turanian, and Nubo-Sindian subzone. The vegetation maps produced by previous researchers were tested with regards to their accuracy; the findings revealed the overall accuracy ranged from weak, at 47% accuracy, to medium, at 61% accuracy. For the unsupervised land cover/land use map developed here, the overall accuracy was calculated at 91%, while the supervised current vegetation map had an overall accuracy of 95%. Jordan's forested regions account for 0.71% of the total land area, with acacia woodland and forests combined totalling 2.3%.  Acknowledgments: Appreciation goes to The Royal Botanic Garden, Edinburgh, and The Royal Botanic Garden, Jordan, for their support over the study period. We thank Edward Mitchard and Antje Ahrends for their valuable advice and discussions during the execution of this research.