Estimate Land Surface Temperature in Relation to Land Use Types and Geological Formations Using Spectral Remote Sensing Data in Northeast Jordan

Land Surface Temperature (LST) is one of the important indicators to understand the spatial changes and surface processes on the earth surface that leads to actual assessment of environmental quality from local to global scales. The relation between spatial analysis of the land surface temperature and existing land use/land cover changes is important to evaluate the climate processes. Monitoring of this relation in the arid and semi-arid regions is necessary to make an appropriate decision about Land surface temperature and environmental status. In this paper, generally the split-window algorithm is used to estimate LST from thermal bands of the Landsat Operational Land Imager (OLI) and the Thermal Infrared Sensor (TIRS) using remote sensing and Geographic Information System (GIS) techniques as well as meteorological data through Moderate Resolution Imaging Spectroradiometer (MODIS). The results show the relationships between land use types and land surface temperature. MODIS data were analyzed. The relationship between MODIS and Landsat data temperature is moderate relation and the (R2 = 0.5109) according on 200 random points were selected. This research concludes that the maximum temperatures of the land use types in MODIS and Landsat data for the rock formation are 59 ̊ and 45 ̊ respectively, whereas the maximum temperatures of the geological formation in MODIS and Landsat data for the basalt are 59 ̊ and 45 ̊ respectively. In conclusion, the MODIS and Landsat OLI and TIRS Data have high ability to distinguish the land use types. The correlation coefficient of the relation between the surface temperature with rock was (R2 = 0.6197). Therefore, it is found that there is an ability to monitor the How to cite this paper: Ibrahim, M. and Abu-Mallouh, H. (2018) Estimate Land Surface Temperature in Relation to Land Use Types and Geological Formations Using Spectral Remote Sensing Data in Northeast Jordan. Open Journal of Geology, 8, 174-185. https://doi.org/10.4236/ojg.2018.82011 Received: January 15, 2018 Accepted: February 23, 2018 Published: February 26, 2018 Copyright © 2018 by authors and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access DOI: 10.4236/ojg.2018.82011 Feb. 26, 2018 174 Open Journal of Geology M. Ibrahim, H. Abu-Mallouh spatial and temporal changes for land surface temperature and thus it can be useful to environmental studies.


Introduction
Global environmental processes to understand climate changes has been studied for over a decade using Land Surface Temperature considered as a good indicator to monitor physical properties of surface processes and climate changes for this reason also LST is an important model of the energy balance at the surface, which means that the LST has a significant role to the effect on environmental processes [1] [2] [3] [4] and this effect occurs in the land surface by different natural and human activities.
LST has a significant impact on analyzing the environmental issues such as urban heat islands, soil moisture, urban heat islands and vegetation which play role in exchanging bioprocesses of water and energy between land surface and air [1] [5]- [10], which are related to the surface temperature (surface heat) [5].
This leads to having different land-use/cover types or surface materials which have several effects on LST from local to global scales.
Moreover, estimate LST is useful to study areas that witness changes and activities in land types and surface materials, which located in the arid and semi-arid areas [5] [11], because high temperature decreases energy activity used in cooling building in hot months [5] [12].On the other hand, high temperature in these arid and semi-arid regions evaporates water, decreases probability of vegetation growing and increases soil erosion activity.
Furthermore, it's found that the relation between Spatial distribution of the LST and surface materials and land use/land cover types is good indicator to evaluate general ecosystem and environmental status.LST is estimated using different remote sensing data in many fields such as; forestry, vegetation, climate change hydrology, oceanography, urban heat island and climate [13] [14] [15], where the remote sensing has high techniques which evaluate and monitor environmental changes [1].In this case, Landsat is used as one of the widely remote sensors used, Landsat 8 follow the Landsat series which provides thermal data [1] [16].LST is investigated by meteorological data, Moderate Resolution Imaging Spectroradiometer (MODIS) one of the most remote sensing data range which is considered one of the important meteorological data from space, where MODIS used to evaluate estimated LST from TIRS of Landsat.

M. Ibrahim, H. Abu-Mallouh
with provided data of TIRS Landsat 8.This study aims at using spectral remote sensing data (TIRS and MODIS data) to estimate Land surface temperature in relation to Land use types and geological formation in northeast Jordan.

Study Area and Data Used
The study area is located in north-east of Jordan as shown in Figure 1, where the study area has a semi-arid climate of the Mediterranean region with a limited amount of rainfall and high temperatures.Its weather is cold and rainy from October to April, dry and hot from May to September.The annual rainfall is 152.1 mm.The mean minimum temperature is 12˚C and the mean maximum temperature is 26.9˚C.The annual relative humidity is 49.7% according WMO (2017) [18] [19].
The main source of data in this research are satellite images from Landsat 8 OLI (Operational Land Imager), TIRs (Thermal Infrared Sensor) and MODIS imagery (Moderate Resolution Imaging Spectroradiometer) which downloaded free from the United States Geological Survey (USGS) Global Visualization (GloVis) site.These data were acquired in the same seasonal period (April 2015).Moreover, geological map was used to find relationship between LST and geological formation that covered study area which was obtained from Natural Resources Authority (NRA).Many software were used in order to pre-processing these data software such as: ERDAS IMAGINE 2014 to pre-processing and processing images, GIS techniques to get land use/land cover using one of the classification approaches, accuracy assessment and to apply LST equation and find correlation R studio to statistics analysis.

Methods
The methodology applied in this research is illustrated in Figure 2. First, the raw data of remote sensing are used to carry out some processes.Second, the spectral indices are applied to carry out equations of LST, the land use type is extracted from Landsat 8 OLI, MODIS is used as reference data, all these were used to estimate land surface temperature.

Image Processing
There are a several ways of calculating and estimating land surface temperature (LST), different algorithms to do that like Split-Window (SW), Dual-Angle (DA), Single-Channel (SC).In this study we will use the Split-Window, the thermal bands and NDVI are acquired for the study area [20] [21] [22].The (Equation ( 1)) of LST is: ) where: LST: Land Surface Temperature.BT: Brightness Temperature.W: Wavelength.P: 1438.e/LSE: Land Surface Emissivity.

Brightness Temperature (BT)
Brightness temperature is the temperature of a blackbody which is used to produce the radiance perceived by the sensor, according NASA 2012 [23].Moreover,  it is the temperature that been received by the satellite at the time that the image was taken.Therefore, this is not the real temperature on the ground; it's the temperature at satellite [22].TIRS band data can be converted from spectral radiance to brightness temperature using the thermal constants provided in the metadata file, the (Equation ( 2)) used to convert from spectral radiance to brightness temperature. where: K 1 and K 2 : Thermal constant.Lλ: Top of Atmospheric spectral radiance (TOA).
Landsat 8 provides some constant to estimate LST such as thermal constant and rescaling factor; it found it in metadata file of Landsat satellite images (Table 1, Table 2).Now, to find the Brightness Temperature (BT), the Top Of Atmospheric spectral radiance (TOA) is acquired.The Equation (3) states: where: ML: Band specific multiplicative rescaling factor (radiance_mult_band_10/11). AL: Band specific additive rescaling factor (radiance_add_band_10/11). Qcal: band 10/11 image.

Land Surface Emissivity
Land surface emissivity is a relative factor that estimates blackbody radiance that is measured from land surface temperature [24].There are several ways to estimate land surface emissivity (Equation ( 4)).One of them is Normalized Difference Vegetation Index (NDVI) method, taking into account the proportion vegetation (Pv), then LST in Celsius is determined.The formula of land surface emissivity is "Equation ( 5)".

Image Classification
Satellite data was downloaded from USGS Earth explorer.Supervised classification was used to land use map, which is considered most common technique of classification approaches [25].These techniques are required to define the training samples for each class and creating the spectral signature where 200 points were samples.

Accuracy Assessments
The accuracy assessments provide more information on where the errors of classification happened.To know how much a classification is accurate, a set of random points must be created to evaluate the data, at the location of each random point.The result would be finding the type of land use of that spot using Google earth (truth points) and comparing it to land use of the classified raster [26].Three standard criteria were used to assess the accuracy of the classifications (overall accuracy, producer's accuracy and user's accuracy) [27].the spatial resolution between sources data and surface roughness effects on the surface temperature [28].Points out that atmospheric impurities can be effect to the surface temperature.

Relationship Estimated LST and Land Use Changes
Four land use types resulted from classification of the image for the study area, shown in the Figure 6.These classes were considered rock, soil, agriculture and urban.For example, by comparing the reference and classified data in the Table 3, it was observed that 82.35% (user accuracy percent) of the rock were being correctly identified as rock, while 66.67% (produces accuracy percent) of the areas called rock are actually rock.Kappa coefficient values of the classified image were 0.73 and overall accuracy was 80%.This accuracy assessment test made to know how the classification is accurate and these values were obtained from   analysis of the reference data (truth points) and the classification output Figure 6 and Table 3.
In the classified image, classes have minimum, maximum mean and standard deviation (SD) values of observed surface temperature and estimated LST Table 4.For the investigate area, results show that the observed surface temperature and estimated LST from single-channel algorithm using Landsat OLI/TIRS sensor data have high accuracy with minimum , maximum mean and standard

Relationship Estimated LST and Geological Formations
Four geological formations in the study area were geologically exposed according NRA (1997) [29], as shown in Figure 7 and these formations were considered: basalt, WadiAlshalla, Muwqqar and Rijam.These formations have minimum, maximum mean and standard deviation (SD) values of observed surface temperature and estimated LST (Table 5).The results show that the observed surface temperature and estimated LST values for these formation has high accuracy with minimum , maximum mean and standard deviation (SD), where the R 2 are 0.98, 0.99, 0.97 and 0.99 for rock, basalt, WadiAlshalla, Muwqqar and Rijam, respectively (Table 5).
It's generally observed that surface temperature has been increased in all types of land use and all type of geological formations, but the good increase in types of land use which was registered is rock type, and it should be noted that agriculture and urban had experienced a slight decrease over average temperatures.This is may be due to source cooling like water for the agriculture or cultivated area that mean in general the relationship between vegetation cover was weak, by increasing vegetation cover the surface temperature will decrease, and materials of building that help to reflect energy or cooling devices in the summer season where this area is considered one of the arid and semi-arid for high temperature in most months of the year.
While highest increase in the geological formation which was registered is basalt, and the other formation had slight decrease over average temperatures.This is probably because of the chemical components of the basalt.Other formation had experienced a slight decrease over average temperatures This may be because of the weak relationship between the minerals content of these formation and energy which lead to reflect energy or minerals have ability to store water which lead to decrease land temperature by cooling.DOI: 10.4236/ojg.2018.82011

Figure 1 .
Figure 1.Location of study area.

Figure 2 .
Figure 2. Methodology flow chart of the study.

3. 1 .
Relationship between the Observation Surface Temperature and Estimated LST In the present study, land surface temperature was estimated by the single-channel from TIRS (Thermal Infrared Sensor) data of the Landsat 8 OLI and the observation surface temperature from MODIS satellite data.The spatial distribution of estimated LST in the Figure 3, surface temperature (LST image) ranged from 24˚C to 51˚C (mean = 38˚C) in the Figure 3, high surface temperatures are shown by the dark red areas which mean the far northeast part of the study area have a higher temperature than the other part of the study area.While the spatial distribution of observation surface temperature in the Figure 4, surface temperature (LST image) ranged from 42˚C to 59˚C (mean = 46˚C).In the Figure 4, high surface temperatures are shown by the dark red areas.This result matches with results of the estimated LST. Statistical regression analysis was used to established relationship between observation surface temperature values and estimated LST values.The result was analyzed using 200 control points, and the correlation coefficient between two different types of surface temperature values (observation surface temperature values and estimated LST values) was R 2 = 0.51 at the same time for the satellites passed over the study area as shown in Figure 5.The difference between the observation temperature values and estimated LST values is due to the difference in

Figure 3 .
Figure 3. Final estimated land surface temperature (LST) of the study area.

Figure 4 .
Figure 4. Observation surface temperature of the study area (MODIS image).

Figure 5 .
Figure 5. Relationship between the observed surface temperature and estimated LST.

Figure 6 .
Figure 6.Land use map of study area (Landsat April 2015).

Table 3 .
Classification accuracy assessment report.

Table 4 .
Land surface temperature of different land use classes of study area.

Table 4 .
On the other hand the results show that these data have low accuracy with 50 control points for each class except rock class has good accuracy where the R2 are 0.62, 0.23, 0.02 and 0.26 for the rock and soil, agriculture and urban respectively.