The Estimation of Soil Organic Matter Variation in Arid and Semi-Arid Lands Using Remote Sensing Data ()
1. Introduction
Soil organic matter (SOM) has a significant effect on the soil processes which in turn affects the ability and productivity of the soil. Moreover, it has an impact on water capacity of soil, soil structure and nutrients which increases soil quality [1] .
Soil organic matter is a mix component which leads to change soil characteristics and affected by land use, climate, vegetation and soil type. That means the productive ability of agriculture will be compromised by deterioration in soil physical properties if SOM concentration decreased in soil [2] . SOM is carbon-rich material that includes plant, animal and microbial residue in various decomposition stages [3] . SOM is the most important indicator of soil quality and productivity. It also consists of a complex and varied mixture of organic substances [4] . SOM enhances soil function and environmental quality because it binds soil particles together into stable aggregates, thus 1) improving porosity, infiltration and root penetration; 2) reducing runoff and erosion so that the soil in arid, semi-arid, and hot humid regions commonly contains less organic matter than the one in other environments [3] . Therefore, SOM contents affect the soil color where, the dark-colored soil has a higher content of organic matter (OM), and vice versa is true [5] .
The quantity of soil organic matter determines the input of organic matter and its rate of oxidation, soil texture, climate, mineralization rate and organic matter decomposition. The clay helps to settle down and limit SOM decomposition. Therefore, soil texture is considered an important factor to agricultural crops as well as rainfall and temperature is a significant factor to the agricultural crops because of its high impact on the OM decomposition rate [6] .
SOM includes all OM in the soil except in decayed material. Thus, SOM supports plant because OM acts as an important source of soil nutrients to plants. The OM is reduced by agricultural and economic development [7] . SOM has a key function in exploring the soil and ecological quality [8] .
The prime interest for soil experts and environmental directors is studying soil characteristics effectively and opportunely. With the growing request of locative and temporal decision in soil characteristics in different agriculture sectors, traditional laboratory techniques were proved to be unsuitable [5] [9] . Furthermore, the costs of soil test with agriculture systems are costly compared to traditional techniques [5] [10] . In this regard, the range of remote sensing can be explored as substitutional technique for studying soil characteristic [5] [11] . Remote sensing (RS) carries the possibility for identifying soil attributes onto field [5] [12] . This technology assisted to soil exploration and decreased the fieldworks of soil to a large level [5] .
Moreover, traditional methods are not always available because of lacking both time and cost connected with sample aggregation. So, it was focused on the employ of remotely sensed data to quantify variation in soil physical characteristic [13] . Soil characteristics description is one of the oldest usage of remotely sensed data [13] .
Soil scientists are interested in using remote sensing and its applications in soil surveys and soil mapping especially after the evolution of optical sensor in conjunction with field measurements [14] [15] [16] . Remote sensing and Geographic Information System (GIS) techniques have modern and useful techniques to assess SOM at high scales due to the quantities of data necessary and the big sites coverage. For this purpose, these mechanisms have been openly used, adopted and displayed in SOM applications [17] .
This study aims to use simple method of results from satellite data and determine the relationship between measured SOM in field and obtained data using remote sensing techniques, find correlation coefficients which applied to evaluate the strength and direction of the linear relationships between variable couples [18] . This paper is characterized by takes into account the soil degradation parameter, where is this parameter effects to the OM particularly in the arid and semi-arid areas.
2. Materials and Method
2.1. Study Area
Study area is located northeast and northwest Jordan and distributed in two governorates: Mafraq and Irbid governorates Figure 1. The Mediterranean climate is dominant in the region (arid and semi-arid climate). Winter season can be described as cold with low temperature through the night and warm arid summer. Mean early rainfall 350 mm [19] .
The average monthly overflow volumes are in December, January, February
and March. Evaporation rate constitutes about 90% of the aggregate rainfall [20] . Average predestined infiltration rate is about 4% - 10% [20] [21] . The climate characteristic of eastern highlands in Jordan is hot in summer and moderate in winters. Low quantities of precipitation happen through comparatively wet months from October to May and they are normally related with the inland steam transport from the Mediterranean Sea. The annual moderate temperature is 17.3˚C. The daily moderate temperature amplitude is from about 8˚C in the winter to 25˚C in the summer. Monthly evaporation average varies from about 85 up to 90 millimeters in December and January to 210 millimeters in April, and up to 300 millimeters in May [22] , and the evaporation extent 567 millimeters in August [23] . The atmospheric dust, arid climate and low precipitation impact the precipitation water quality thus growing salt content [24] [25] .
2.2. Geological Setting
The study area consists fundamentally of Paleogene to quaternary continental tuffs and basalts, outcrops back to the Azraq Formation, existing within a rift zone, limited by the NW–SE trending Fuluq and Sirhan faults [26] . Azraq formation base is unconformable on below (older strata) and the top (younger strata) is known by the present-day depositional and erosional surface [27] . Study area contains two groundwater basins Amman-Zerqa basin in Mafraq area and Yarmouk basin in Irbid area [20] . Most of the geologic formations consist of marly layers [24] [25] .
2.3. Methodology
The methods implemented in this study can be broadly divided into three subsections; fieldwork, image processing and statistics analysis.
2.3.1. Fieldwork and Lab Work
Soil Samples from 22 sites located in two zones (semi-arid and arid zones) were carefully collected in plastic bags and appropriately labeled. The soil samples were collected from depth 30 cm. The fieldwork visit includes three types of lands; available agricultural in the study period, abandoned land, and protected land. The purpose of this fieldwork was to collect soil samples from these locations in validation of the land degradation index by assessing soil quality and structure in addition to chemical analyses of these samples. Samples were then classified to different soil texture type by sieve analysis method [28] . In a similar way, the SOM content was measured using the TPDM (Titration of Potassium Dichromate Method) [5] [29] . TPDM actually was used to estimate carbon percentage, Organic Carbon (OC) and Organic Matter (OM) percentages which were calculated based on the following Equation (1):
(1)
B is the volume of Fe2+ solution used to titrate blank (mL). S is the volume of Fe2+ solution used to titrate sample (mL), and 12/4000 is milliequivalent weight of C in grams to convert easily oxidizable organic C to total C, divided by 0.77 (or multiply by 1.30) or other experimentally determined correction factor. Therefore, the percentage of organic matter can be estimated by the following Equation (2):
(2)
SOM was classified for five classes and used to implement this study after a modified soil-organic-matter test [30] [31] that shown in Table 1.
2.3.2. Data Acquisition
The satellite used Landsat 8 OLI (Operational Land Imager) and TIRS (Thermal Infrared Sensor) used in this study in May 2017 path 169 and row 45. The image was obtained from the United States Geological Survey (USGS) Global Visualization (GloVis) site and geometrically corrected and rectified to UTM zone 36. In order to prepare data, managing and analyzing software were conducted as follows: ERDAS IMAGINE 2014 in this study ERDAS was applied in importing and enhancement, as well as ArcGIS 10.5 used to digitizing, indexing and image analysis, geo-referencing, creation of database. Arc map was used for the composition and generation of maps.
2.3.3. Image Analysis
Three indices were applied in this research to achieve the aims, where these indices are widely used and these were selected to suit the acquisition data that used in this paper. The Normalized Difference Vegetation Index (NDVI)—which is a normalized ratio of red and near-infrared reflectance. Rouse et al., 1973 have been used in many phonological studies. NDVI is calculated as Equation (3) [32] :
(3)
where the (λNIR) represent reflection in near infrared, IR (λ ≈ 0.8 µm) regions of the spectrum. (λR) represent surface reflectance in the visible light, R (λ ≈ 0.6 µm, “red”). The large variations in NDVI can be due to variations in soil brightness [33] . The wide regions of bare soil in arid and semi-arid regions may be the
Table 1. Percentages of organic matter in the soil.
reason of considerable differences in NDVI values and that it’s not represent of the actual vegetation cover [34] [35] [36] .
The NDVI is measured by satellite sensors. Rouse et al. (1973) [32] were the first to write a report of NDVI [37] . Many factors affect NDVI values other than the vegetation itself such as total plant cover, plant photosynthetic activity, biomass, and plant stress. Soil and atmosphere have been shown to have a great effect on the NDVI [38] [39] [40] .
The Soil Adjusted Vegetation Index (SAVI) (Huete, 1988) [41] was also used in this study. SAVI is an adjustment on the (NDVI) that decreases soil effect, it was developed by (Huete 1988; Huete et al., 1987; Kaufman et al., 1902; Epiphanio et al., 1995) [38] [39] [40] [41] , and is expressed by the following Equation (4):
(4)
where the NIR represents reflection in near infrared (730 - 1000 nm), RED represents surface reflectance in the visible light (550 - 700 nm) and L is soil adjustment factor (Ray et al., 2001). The value of L varies by the amount or cover of green vegetation and varying from 0 (close canopy cover) to 1 (open canopy cover): in very high vegetation regions L = 0; and in areas with no green vegetation, L = 1. Generally, an L = 0.5 works well in most situations and is the default value used.
However, removal of plant covers jeopardizes the carbon cycle and decreases the percentage of organic matter in the soil, as well as resulting in a degradation of physical soil properties and structure [42] . The removal of vegetation processes increases erosion and land degradation, particularly in the semi-arid and dry sub-humid areas that have a high sensitivity to desertification and degradation. Table 2 shows three classes of vegetation cover that were used in this study [42] [43] [44] .
Bare Soil Index (BSI) value was estimated using combines blue, red, green and near-infrared bands to detect the soil changes of the Landsat 8 OLI and TM data to determine the condition of the uncovered soil by vegetation which is helpful to determine the state of the soil organic matter on the region. BSI is described by the following Equation (5) according to Jamalabad & Abkar (2004) [45] :
(5)
Table 2. vegetation removal percentages.
3. Results and Discussions
Results of laboratory measurements of soil indicate that organic matter percentage in the soil (Organic Carbon) in study area varies from low to very low based on Table 1. SOM was ranging between (0.38% to 0.9%) in the areas with relatively arid climate condition with mean value of 0.68 ± 0.156, while in areas with relatively semi-arid (humid) climate condition SOM ranged from (0.59% to 1.68%) with mean value of 0.89 ± 0.32. Statistical analysis of this study relies on the five variables which were used; SOM, water content, NDVI, SAVI and BSI in the study area shown in Table 3.
3.1. Relation between Vegetation Index and SOM
It was found that the relationship between Vegetation Index (NDVI, SAVI) and SOM was statistically significant and the (R2 = 0.19, p < 0.05), as shown in Table 4 and Figure 2, the form of the model that appears this relationship is as following Equation (6):
SOM = 2.96 × NDVI + 0.36 (6)
This model describes the relation between (NDVI, SAVI) and SOM in both arid and semi-arid regions. However, the relationship is relatively different in arid areas from that in semi-arid areas (R2 = 0.01, p < 0.05), (R2 = 0.13, p < 0.05), respectively. It was also found that the relationship between (NDVI, SAVI) and SOM increases in semi-arid areas that have diversity in vegetation cover and humidity more than arid areas that have dry climate change but do not have variety in vegetation cover. The means of laboratory measurement of soil in both arid and semi-arid or humidity regions were (0.68, 0.89), respectively.
3.2. Relation between Vegetation Index and BSI
The relationship between vegetation index such as NDV and BSI is represented
Table 3. The statistical analysis for four variables used; SOM, NDVI, SAVI and BSI in the study area.
Table 4. The statistically significant relationship for four variables used.
Figure 2. Relation between vegetation index and SOM.
in Figure 3. The results shown a statistically significant relation between NDVI and BSI (R2 = 0.81, p < 0.05). That increasing degradation conditions mean decrease in rainfall and increase in dry climate that lead to decrease in NDVI.
3.3. Relation between BSI and SOM
It was found that there are differences in the relationship between measured SOM and Bare soil Index (BSI). This relation is shown in Figure 4. The results showed a statistically significant relation between measured SOM and BSI (R2 = 0.114, p < 0.05). The form of the model that appears this relationship is as following Equation (7):
SOM = 0.09 × BSI − 9.99 (7)
The outcomes also indicated that SOC increases with increasing NDVI and decreasing BSI. Finally, the multiple regression models were fitted into ArcGIS to generate the digital SOC map using raster calculator Figure 5. The statistical correlation of the relationship between the estimated SOM and measured SOM was calculated (R2 = 0.23) as shown in Figure 6.
It was also found that the regions, that have both a high rainfall (250 - 400 mm annually) and a good vegetation percentage, have the high concentrations relatively of SOM. These regions were distributed in eastern part boarders and some regions in the central part of study area which is semi-arid. The region covered by urban, rock and sand has extremely low SOC values about zero, while the central part and southern-east boarders have low of SOC values and have a low rainfall relatively (less than 200 mm annually) and a high temperature.
4. Conclusion
The SOM measured values showed a good correlation with used indices, where
Figure 3. Relation between vegetation index and BSI.
Figure 6. The relationship between the estimated SOM and measured SOM.
this study suggested that these vegetation indices (NDVI and SAVI) and bare soil index (BSI) are helpful to detect the spatial distribution of SOM concentrations. In addition the statistical analysis results were useful to analyse measured values in the field and to understand the relationships that helped to estimated SOM concentrations by remote sensing data, although the pixel size (30 m × 30 m) used in the study area was relatively large with this area, which reflected the variation in SOM concentrations. The results of this study support the approach that the remote sensing data are useful to detect and investigate the environmental changes. This approach was used through mapping and determines the spatial distribution of SOM concentrations using spectral indices processed remote sensing data; techniques could be effective to monitor and manage soil.
Acknowledgements
Authors are thankful to anonymous reviewers for their constructive comments and suggestions to improve the manuscript. Authors also are thankful to Department of Earth and Environmental Sciences members in Al Al-Bayt University for providing all necessary support.