Open Access
13 May 2019 Analytical study of seasonal variability in land surface temperature with normalized difference vegetation index, normalized difference water index, normalized difference built-up index, and normalized multiband drought index
Author Affiliations +
Abstract
Remote sensing technique often analyzes the thermal characteristics of any area. Our study focuses on estimating land surface temperature (LST) of Raipur City, emphasizing the urban heat island (UHI) and non-UHI inside the city boundary and the relationships of LST with four spectral indices (normalized difference vegetation index, normalized difference water index, normalized difference built-up index, and normalized multiband drought index). Mono-window algorithm is used as LST retrieval method on Landsat 8 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS) data, which needs spectral radiance and emissivity of TIRS bands. The entire study is performed on 11 multidate Landsat 8 OLI and TIRS images taken from four different seasons; premonsoon, monsoon, postmonsoon, and winter, in a single-year time period. The Landsat 8 data derived LST is validated significantly with Moderate Resolution Imaging Spectroradiometer (MOD11A1) data. The results show that the UHI zones are mainly developed along the northern and southern portions of the city. The common area of UHI for four different seasons is developed mainly in the northwestern parts of the city, and the value of LST in the common UHI area varies from 26.45°C to 36.51°C. Moreover, the strongest regression between LST and these spectral indices is observed in monsoon and postmonsoon seasons, whereas winter and premonsoon seasons revealed comparatively weak regression. The results also indicate that landscape heterogeneity reduces the reliability of the regression between LST with these spectral indices.

1.

Introduction

Land surface temperature (LST) estimation is considered as one of the most significant tasks in thermal remote sensing study, where thermal impacts of changing land use/land cover (LULC) may regulate the generation of urban heat island (UHI) in mixed urban landscapes.15 The effect of UHI on some major cities (e.g., Beijing, Columbus, Shanghai, Baltimore, and Chicago) in the world is growing at an alarming rate and is directly related to the LULC categories.611 Various LULC indices (vegetation, impervious surface area, etc.) were applied in a UHI related remote sensing study to determine the changed pattern of LST in different types of natural and man-made landscape systems.1216 Some recent articles discussed the statistical linear regression between LST and some selected LULC indices for separate study areas like Brisbane,17 Raipur,18 Shanghai,19 Addis Ababa,20 Mexico,21 Philadelphia,22 Florence, and Naples.23

A more recent direction includes the application of some sophisticated statistical techniques to estimate the seasonal variability of LST in the urban area. Cui and de Foy24 showed that vegetation cover, daytime insolation, and atmospheric stability are related to seasonal UHI variations in Mexico City, Mexico. Zhou et al.25 made an attempt in Baltimore, USA, to build relationships between land cover and UHI and their seasonal variability. Haashemi et al.26 presented a seasonal variability in the regression between LST and fractional vegetation cover, albedo, impervious surfaces, and elevation in Tehran, Iran.

In later studies, seasonal variation in the UHI effect was critically analyzed. Seasonal effects in urban thermal patterns were determined by spatial regression analyses in Ohio, USA.27 Multiple statistical methods were integrated to show the seasonal contrast of some LULC indices for LST distribution in Shenzhen, China.28 A seasonal variation in LST and selected LULC indices was investigated in Jaipur, India.29 A trend and seasonal decomposition model for LST change over Beijing, China, was investigated.30 A regression-based model was performed to determine the average annual and seasonal trend of LST in peninsular Spain.31,32 The relationship of urbanization and climate variability with urban and rural LST for 31 cities of China was critically investigated.33 An analytical research work performed in more than 86 major cities of China presented that UHI is seasonally different and is negatively related to the cloud cover percentage in transitional seasons.34 A landscape source-sink distance index was applied to relate the landscape connectivity and LST variation in Beijing, China.35 A variation of LST due to elevation change in summer, winter, and monsoon season for Jaipur, India, was assessed.36 In summary, recent researchers tried to make a true attempt in the seasonal variability of LST in UHI of the major cities using common statistical methods.

This paper was an analytical study based on the seasonal variability of LST distribution inside the UHI and the area outside the UHI zones in a typical Indian city. The area outside the UHI zones is described here as non-UHI.18,23 Basically, the term “non-UHI” is slightly different from the term “urban cool island (UCI).” Generally, the term UCI indicates that the low-temperature zones inside an urban area are characterized by urban vegetation, wetland, and water bodies, whereas the term “non-UHI” is used to present an area other than the UHI, which is actually denoted moderate-to low-temperature zones of an urban area. Hence, UCI can be considered as a part of non-UHI in an urban area. Another term “common-UHI” was used in this study to show the common area of UHI for each and every satellite image (total of 11 images) for four different seasons.18,23

Landsat TIR data are often used in the identification of UHI zones.3739 LST derivation using Landsat TIR data becomes popular through the introduction of some algorithms, such as mono-window algorithm,40,41 single-channel algorithm,4145 radiative transfer algorithm,41,45 split-window algorithm.41,4547 Landsat 8 TIRS dataset has two TIR bands (bands 10 and 11). Generally, TIR band 10 is used for LST derivation as TIR band 11 is associated with larger uncertainties and thus it is not recommended to use the split-window algorithm. In the radiative transfer equation, LST estimation requires atmospheric correction and emissivity correction, which are difficult to perform.41 Ground emissivity, atmospheric transmittance, and effective mean atmospheric temperature are the parameters needed to derive the LST using mono-window algorithm, whereas single-channel algorithm needs atmospheric transmissivity, upwelling, and downwelling atmospheric radiance. Single-channel algorithm and mono-window algorithm both present significant results in LST retrieval.18,23,43,48 In this study, LST derivation process was conducted using the mono-window algorithm.

The specific aims of the study were: (1) to describe the nature of LST in the UHI, non-UHI, common UHI, and the whole of Raipur City for Landsat 8 OLI and TIRS images selected from premonsoon, monsoon, postmonsoon, and winter seasons in a single-year span and (2) to examine the variation in nature and strength of regression between LST and normalized difference vegetation index (NDVI), normalized difference water index (NDWI), normalized difference built-up index (NDBI), and normalized multiband drought index (NMDI) for the UHI, non-UHI, and common UHI and for the whole of Raipur City in these four particular aforesaid images.

2.

Study Area and Data

Raipur, the capital city and the largest city of Chhattisgarh State in east-central India, located along the west of Mahanadi River, was selected for the present research work. Raipur is characterized by comparatively dry climate having no influence of coastal or maritime climatic condition. Moreover, Raipur is a rapidly growing city in India where a lot of new conversions of land have been seen in the recent time period. It has a wide variety in the population. The city has also seen exponential growth in industrial sectors and has emerged as a major business hub in central India. It is also ranked seventh in Ease of Living Index 2018 by Union Ministry of Housing and Urban Affairs.49 Raipur is considered to be one of the modern and smart technological cities in India due to the addition of new technology to the older system, power plant retrofit, home energy retrofit, seismic retrofit, eco-friendly, and sustainable nature; and thus the city was selected as the study area of the research work. Due to the presence of various types of surface materials, the nature of LST distribution is too dynamic. The latitudinal and longitudinal extents of Raipur City are 21°11’22”N-21°20’02”N and 81°32’20”E-81°41’50”E, respectively (Fig. 1). The mean elevation of the city is around 298 m. Raipur City is characterized by a tropical wet and dry climate. There are four seasons observed in Raipur according to the India Meteorological Department, i.e., premonsoon or summer (March to May), monsoon (June to September), postmonsoon (October to November), and winter (December to February). Generally, the summer months (March to May) remain hot and dry. The mean annual temperature of Raipur City lies between 21°C to 34°C. Amount of rainfall is moderate (average annual rainfall is 120 to 150 cm.). Rain occurs mainly in the monsoon season. At that time, vegetation was looking green and healthy, containing more chlorophyll component. In the rainy monsoon season, temperature falls at a significant rate from the summer or premonsoon season. A pleasant weather condition persists throughout the postmonsoon season. Winter season remains cool and dry. December is considered the coldest month (average temperature 12°C) of the city.

Fig. 1

The study area: (a) Chhattisgarh in India, (b) Raipur City in Chhattisgarh, and (c) Raipur City.

JARS_13_2_024518_f001.png

A total of eleven (three from the premonsoon, two from the monsoon, three from the postmonsoon, and three from winter season) Landsat 8 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS) data (<10% cloud coverage) were selected to determine the UHI zones over all of Raipur City (Table 1). The Landsat 8 datasets were freely downloaded from the website of the United States Geological Survey.50 Only Landsat TIR band 10 data were used for LST retrieval process because TIR band 11 data faces some calibration uncertainty. Optical bands 3, 4, 5, 6, and 7 datasets were used in developing NDVI, NDWI, NDBI, and NMDI. The entire study was performed in the environment of some sophisticated image processing and GIS software packages.

Table 1

Specification of multidate Landsat 8 OLI and TIRS satellite images.

Date of acquisitionSeasonTimePath/rowSun elevation (deg)Sun azimuth (deg)Cloud cover (%)Earth–Sun distance (astronomical unit)
09-Nov-2013Postmonsoon14:27:51142/04447.1003152.51470.030.9905
25-Nov-2013Postmonsoon14:27:43142/04443.2895153.97500.000.9871
11-Dec-2013Winter14:27:42142/04440.6309153.52450.000.9847
27-Dec-2013Winter14:27:33142/04439.5066151.60831.690.9835
28-Jan-2014Winter14:27:13142/04442.3669144.86230.490.9849
17-Mar-2014Premonsoon14:26:36142/04455.9542129.38720.000.9949
02-Apr-2014Premonsoon14:26:19142/04460.9189121.72160.000.9995
20-May-2014Premonsoon14:25:38142/04468.563890.40855.461.0118
05-Jun-2014Monsoon14:25:45142/04468.382183.30980.021.0146
25-Sep-2014Monsoon14:26:11142/04459.2100134.18040.811.0030
12-Nov-2014Postmonsoon14:26:21142/04446.2266152.46647.590.9899

3.

Methodology

3.1.

Determination of NDVI, NDWI, NDBI, and NMDI

NDVI51 is considered one of the most frequently used vegetation indices in remote sensing study. It is also applied in deriving LST and normally shows a negative regression with LST. NDWI52 is generally used for water body extraction. NDBI53 is another spectral index, which was applied in this study for built-up area extraction. NMDI54 was also used to extract the dry soil or bare land. The formulation of these four spectral indices was presented in Table 2. Integration of suitable threshold values of these spectral indices is useful in identifying several LULC types.55

Table 2

Description of NDVI, NDWI, NDBI, and NMDI.

AcronymDescriptionFormulationReference
NDVINormalized difference vegetation indexNIRRedNIR+Red51
NDWINormalized difference water indexGreenNIRGreen+NIR52
NDBINormalized difference built-up indexSWIR1NIRSWIR1+NIR53
NMDINormalized multiband difference index[NIR(SWIR1+SWIR2)][NIR+(SWIR1+SWIR2)]54

3.2.

LST Derivation Using Landsat 8 OLI and TIRS Band

For retrieving LST using mono-window algorithm, the original Landsat 8 TIR band (100-m resolution) was resampled into 30 m. The entire procedure included the following equations:56

Eq. (1)

Lλ=0.0003342×DN+0.1,
where Lλ is spectral radiance (Wm2sr1mm1)

Eq. (2)

Tb=K2ln(K1Lλ+1),
where Tb is the at-sensor brightness temperature (K),57 K2 and K1 are the calibration constants. Here K1 is 774.89 (Wm2sr1mm1) and K2 is 1321.08 (Wm2sr1mm1)

Eq. (3)

Fv=(NDVINDVIminNDVImaxNDVImin)2,
where NDVImin is the minimum value (0.2) of NDVI for bare soil pixel and NDVImax is the maximum value (0.5) of NDVI for healthy vegetation pixel.5860

dε is the geometric distribution effect for the natural surface and internal reflection. The value of dε may be 2% for mixed and elevated land surfaces58

Eq. (4)

dε=(1εs)(1Fv)Fεv,
where εs is the soil emissivity, Fv is the fractional vegetation, F is a shape factor (0.55),58 and εv is the vegetation emissivity

Eq. (5)

ε=εvFv+εs(1Fv)+dε,
where ε is land surface emissivity. The value of ε is determined by the following equation:58

Eq. (6)

ε=0.004×Fv+0.986.

Water vapour content is determined by the following equation:61

Eq. (7)

w=0.0981×{10×0.6108×exp[17.27×(T0273.15)237.3+(T0273.15)]×RH}+0.1697,
where w is the water vapor content (g/cm2), T0 is the near-surface air temperature (K), and RH is the relative humidity (%). These parameters of the atmospheric profile were the average values of 14 stations around Raipur, which were obtained from the Meteorological Centre, Raipur and the Regional Meteorological Centre, Nagpur.

Eq. (8)

τ=1.0314120.11536w,
where τ is the total atmospheric transmittance.62

Raipur City is located in a tropical region. Thus the following equations were applied to compute the effective mean atmospheric transmittance of Raipur:

Eq. (9)

Ta=17.9769+0.91715T0,

Eq. (10)

Ts={a(1CD)+[b(1CD)+C+D]TbDTa}C,

Eq. (11)

C=ετ,

Eq. (12)

D=(1τ)[1+(1ε)τ],
where Ta is the mean atmospheric temperature, Ts is the LST, a=67.355351, and b=0.458606.

3.3.

Mapping UHI and Non-UHI Inside a City Area

UHI and non-UHI zones were demarcated using the following methods:18,23

Eq. (13)

LST>μ+0.5*σ,

Eq. (14)

0<LSTμ+0.5*σ,
where μ and σ are the mean and standard deviation values of LST for the whole area of the city.

4.

Results and Discussion

4.1.

Variation in the Distribution of NDVI, NDWI, NDBI, and NMDI

The descriptive statistics (Table 3) were presented the reliable nature of NDVI, NDWI, NDBI, and NMDI (Fig. 2) for the whole of Raipur City. The highest mean NDVI value (0.1468) was found in the monsoon season. Premonsoon season reflected the minimum value (0.1466) of mean NDWI. The lowest mean value (0.0569) of NDBI was observed in the postmonsoon season. The highest (0.2743) mean NMDI was found in premonsoon season and the lowest (0.2866) mean NMDI value was found in winter. In the winter season, all the four indices had the least value of standard deviation.

Table 3

Descriptive statistics of NDVI, NDWI, NDBI, and NMDI for the whole of Raipur City in four seasons.

SeasonMinMaxMeanStandard deviation
NDVI values
Premonsoon (mean)0.10820.47020.14280.0645
Monsoon (mean)0.10890.46640.14680.0635
Postmonsoon (mean)0.10250.41830.13920.0659
Winter (mean)0.12650.42840.09540.0561
NDWI values
Premonsoon (mean)0.41120.15140.14260.0576
Monsoon (mean)0.41020.14340.14660.0564
Postmonsoon (mean)0.35960.13710.11920.0602
Winter (mean)0.40060.15820.08770.0398
NDBI values
Premonsoon (mean)0.31950.16680.01980.0474
Monsoon (mean)0.32190.23730.01580.0477
Postmonsoon (mean)0.31530.17540.05690.0426
Winter (mean)0.35080.18570.01990.0374
NMDI values
Premonsoon (mean)0.46350.07730.27430.0518
Monsoon (mean)0.54830.07140.27610.0511
Postmonsoon (mean)0.53610.06020.23850.0499
Winter (mean)0.59460.09580.28660.0398

Fig. 2

Seasonal variations of NDVI, NDWI, NDBI, and NMDI values in the whole of Raipur City: (a) premonsoon, (b) monsoon, (c) postmonsoon, and (d) winter.

JARS_13_2_024518_f002.png

4.2.

Variation of LST Distribution

The LST maps retrieved from satellite image were shown in Fig. 3. Seasonal variation in the LST distribution showed a specific thermal pattern. The mean LST values in the premonsoon, monsoon, postmonsoon, and winter seasons were 34.40°C, 31.70°C, 25.22°C, and 24.71°C, respectively. The ranges of temperature (maximum temperature–minimum temperature) were found as 12.71°C in premonsoon image, 10.20°C in monsoon image, 8.81°C in the postmonsoon image, and 10.03°C in the winter image, respectively. For UHI zones of the study area, the highest threshold LST (35.23°C) was observed in the premonsoon season. Postmonsoon season reflected the least standard deviation value (1.06°C) of LST (Table 4). Basically, this type of LST variation was noticed due to the change in vegetation abundance and soil moisture content. Monsoon and postmonsoon seasons were characterized by healthy vegetation and wet soil.

Fig. 3

Seasonal variations in LST maps of Raipur City: (a) premonsoon, (b) monsoon, (c) postmonsoon, and (d) winter.

JARS_13_2_024518_f003.png

Table 4

Seasonal variations in spatial distribution of LST (°C) for the whole of Raipur City.

SeasonLST (minimum)LST (maximum)LST (mean)LST (standard deviation)Threshold LST for UHI
Premonsoon (mean)26.9739.6834.401.6535.23
Monsoon (mean)26.4336.6331.701.1632.28
Postmonsoon (mean)21.5330.3425.221.0625.75
Winter (mean)20.1230.1524.711.1525.29

4.3.

LST Validation with Respect to MODIS Satellite Sensor

A validation of retrieved LST using Landsat 8 data with in situ measurement or with any other satellite data is mandatory to conduct any kind of further analysis. Due to the unavailability of in situ measurement, Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data were applied as a reference image for the validation of LST retrieved from Landsat 8 data. The MODIS LST daily product is derived from bands 31 and 32 by the emissivity of land cover types, atmospheric column water vapor, and lower boundary air surface temperature.63,64

Landsat and MODIS sensors could not provide images of the same study area for any particular date. Thus the acquisition dates of MODIS data were either one day before or one day after the acquisition dates of Landsat 8 data. A total of eleven (three from the premonsoon season, two from the monsoon season, three from the postmonsoon season, and three from winter season) MOD11A1 data were taken for the validation of estimated LST using Landsat 8 data. No such atmospheric disturbances were noticed in the acquisition dates of the selected satellite images for both the sensors. For MOD11A1 and Landsat 8 data, the spatial resolution of retrieved LST was 1000 and 100 m, respectively. Therefore, before integrating the data, 1000-m pixel size was resampled into 100-m pixel size (Fig. 4). There was found a small difference between Landsat derived mean LST and the corresponding MODIS derived mean LST values due to (a) 30 min intervals between the Landsat 8 and MODIS sensors (b) water vapor content, and (c) scale effect in resampling method.63 Figure 4 compared four Landsat 8 data derived LST images with the four corresponding MODIS LST images for four different seasons.

Fig. 4

Retrieved LST from Landsat 8 data and MODIS data: (a) 02-Apr-2014 (Landsat 8), (b) 03-Apr-2014 (MODIS), (c) 25-Sep-2014 (Landsat 8), (d) 24-Sep-2014 (MODIS), (e) 09-Nov-2013 (Landsat 8), (f) 08-Nov-2013 (MODIS), (g) 27-Dec-2013 (Landsat 8), and (h) 28-Dec-2013 (MODIS).

JARS_13_2_024518_f004.png

For the validation of the results, the pixel-wise regression analysis method was applied between Landsat 8 retrieved LST values and the corresponding MODIS retrieved LST values. Significant positive correlation coefficients (0.70 for the premonsoon, 0.75 for the monsoon, 0.69 for the postmonsoon, and 0.59 for winter) were found between Landsat 8 data derived LST and MODIS data derived LST (Table 5). Hence, it can be said that the results were reliable and consistent in spite of lacking some LST retrieval parameters.

Table 5

Validation of LST retrieved from Landsat 8 data with MODIS data in different dates of acquisition.

Landsat 8 dataMODIS dataCorrelation coefficient
09-Nov-201308-Nov-20130.69
25-Nov-201326-Nov-20130.66
11-Dec-201312-Dec-20130.61
27-Dec-201328-Dec-20130.68
28-Jan-201427-Jan-20140.45
17-Mar-201418-Mar-20140.74
02-Apr-201403-Apr-20140.68
20-May-201421-May-20140.68
05-Jun-201404-Jun-20140.73
25-Sep-201424-Sep-20140.77
12-Nov-201413-Nov-20140.71
Premonsoon (mean)0.70
Monsoon (mean)0.75
Postmonsoon (mean)0.69
Winter (mean)0.59

4.4.

UHI, Non-UHI, and Common UHI of the Study Area

The UHI intensity (difference between the mean LST of UHI and the mean LST of non-UHI) of the study area was determined in Table 6. In premonsoon and winter seasons, the UHI zones were mainly generated in the north, west, and south-east periphery (Fig. 5). But in monsoon and postmonsoon seasons, the northern and central parts (the main built-up areas and bare lands of the city) were considered as the UHI zones. The UHI intensity for the whole of Raipur City was 2.54°C, 1.87°C, 1.76°C, and 1.96°C in the premonsoon, monsoon, postmonsoon, and winter seasons, respectively. The mean LST values of the common UHI of the city in all four seasons were ranged between 26.45°C (winter image) and 36.51°C (premonsoon or summer image). Regardless of any particular season, the common UHI zones were developed mainly in the northwest portion (bare lands and built-up areas) of the city (Fig. 6).

Table 6

Seasonal variations in LST (°C) of UHI, non-UHI, and common UHI in Raipur City.

SeasonLST (minimum)LST (maximum)LST (mean)LST (standard deviation)
UHINon-UHICommon UHIUHINon-UHICommon UHIUHINon-UHICommon UHIUHINon-UHICommon UHI
Premonsoon (mean)35.2326.9735.2339.6835.2339.6836.2333.6936.510.781.320.78
Monsoon (mean)32.2826.4332.2836.6332.2836.6333.0731.2033.570.760.830.81
Postmonsoon (mean)25.7521.5325.7530.3425.7530.3426.4924.7326.990.710.700.79
Winter (mean)25.2920.1225.2930.1525.2930.1526.1724.2126.450.760.780.78

Fig. 5

Seasonal variations in UHI and non-UHI of Raipur City: (a) premonsoon, (b) monsoon, (c) postmonsoon, and (d) winter.

JARS_13_2_024518_f005.png

Fig. 6

Seasonal variations in common UHI of Raipur City: (a) premonsoon, (b) monsoon, (c) postmonsoon, and (d) winter.

JARS_13_2_024518_f006.png

4.5.

Relationship of LST with NDVI, NDWI, NDBI, and NMDI for the Whole Area, UHI, non-UHI, and Common UHI of Raipur City

Generally, LST built negative relationships with NDVI and NMDI, whereas it built positive relationships with NDWI and NDBI. This particular pattern was seen in the whole of Raipur City, irrespective of dates (Table 7). NDBI and NMDI equally built strong regression with LST for four seasons. Their relationships became stronger in the premonsoon, monsoon, and postmonsoon seasons compared to winter. NDVI and NDWI showed weak relationships with LST compared to NDBI and NMDI. NDVI presented a slightly better relationship than NDWI. Winter season was the least reliable season for the relationships between LST and the four LULC indices. LST-NDWI relationship remained almost neutral in winter. This is partially due to the presence of high amount of dust particles in the air in winter.

Table 7

Correlation coefficients for LST-NDVI, LST-NDWI, LST-NDBI, and LST-NMDI relationships.

SeasonLST-NDVILST-NDWILST-NDBILST-NMDI
(Whole of Raipur City)
Premonsoon (mean)0.42350.29370.69840.6988
Monsoon (mean)0.49780.41230.45280.5291
Postmonsoon (mean)0.51240.42720.62390.6647
Winter (mean)0.17480.04140.56760.5509
(UHI of Raipur City)
Premonsoon (mean)0.32980.26090.28830.1622
Monsoon (mean)0.49980.5755−0.08500.0876
Postmonsoon (mean)0.48660.51440.09050.2367
Winter (mean)0.31040.25720.22130.3657
(Non-UHI of Raipur City)
Premonsoon (mean)0.28340.12340.68770.6553
Monsoon (mean)0.26080.08590.54370.5509
Postmonsoon (mean)0.20320.08220.58270.5637
Winter (mean)0.0141−0.16330.54970.4782
(Common UHI of Raipur City)
Premonsoon (mean)0.35910.32670.17590.2477
Monsoon (mean)0.51700.59080.19920.0183
Postmonsoon (mean)0.56070.54520.15790.3367
Winter (mean)0.33140.27260.06960.1892

But these relationships tend to be changed within the UHI of the city where NDVI and NDWI reflected a stronger regression with LST for all the four seasons (Table 7). Again, the scenario became different within the non-UHI portions of the city where NDBI and NMDI presented a stronger regression with LST compare to the other two spectral indices. In the common UHI region, NDVI and NDWI had a much better regression with LST, but these relationships became gradually weak with the increase of the heterogeneous surface features. Monsoon and postmonsoon images revealed a better regression of LST with NDVI and NDWI due to the higher percentage of water and vegetation.

In the UHI and common UHI, high LST values were found in such areas where NDVI values were low but NDWI values were high. NDWI showed the strongest positive regression with LST, whereas NDVI showed the strongest negative regression. The monsoon and postmonsoon season revealed the most consistent relationship. In the non-UHI, the high LST values were corresponding to high NDBI and low NMDI values. Due to the presence of more moisture content in soil and air, the relationships between LST and these spectral indices became more consistent in the monsoon and postmonsoon images. NDBI showed the strongest positive regression with LST, whereas NMDI showed the strongest negative regression. The premonsoon season reflected the most significant relationship.

The best regression between LST and the four spectral indices was found for the whole of Raipur City, irrespective of any date. This relationship became weaker with the increase in heterogeneity in an urban landscape. Common UHI of all seasons simply indicated the built-up area and semibare lands, which were more heterogeneous. Thus the least regression was found in the common UHI of the city area. UHI and non-UHI of the city reflected the moderate range of the correlation coefficient values.

5.

Conclusion

In this paper, eleven Landsat 8 data of four different seasons in a single-year span were selected for analyzing the effect of UHI intensity in Raipur City of India and to interpret the seasonal variations for the relationships of NDVI, NDWI, NDBI, and NMDI with LST. The above relationships were examined for the whole area, UHI, non-UHI, and common UHI of the Raipur City. UHIs were identified through the spatial distribution of LST, which mainly existed in the northern and southern parts of the city. Generally, high LST values were generated in the bare land and built-up area. LST level was reduced significantly in the vegetal covered areas and water bodies.

Moreover, the relationships of LST to NDVI, NDWI, NDBI, and NMDI were analyzed using linear regression per-pixel level. In whole Raipur City, LST showed a strong positive regression with NDBI and a moderate to strong negative regression with NMDI, irrespective of any season. Inside the UHI, NDVI and NDWI showed a stronger regression (NDVI-negative and NDWI-positive) with LST in comparison with the other indices. Conversely, inside the non-UHI zone, NDBI and NMDI presented a stronger regression (NDBI-positive and NMDI-negative). Again, NDVI-LST and NDWI-LST regression became stronger within the common UHI in all seasons.

Another objective of the research was to estimate the variations in the regression analysis for the satellite images of four different seasons. Monsoon and postmonsoon seasons were more prominent in showing the regression between LST and spectral indices due to the higher percentage of healthy green vegetation and soil moisture. The premonsoon season was less dominant compared to the monsoon and postmonsoon seasons while this regression became weakest in the winter image.

Acknowledgments

The authors were indebted to the United States Geological Survey.36 The authors also thank the anonymous reviewers whose feedback improved this manuscript significantly.

References

1. 

J. Arnfield, “Two decades of urban climate research: a review of turbulence, exchanges of energy and water, and the urban heat island,” Int. J. Climatol., 23 (1), 1 –26 (2003). https://doi.org/10.1002/joc.859 IJCLEU 0899-8418 Google Scholar

2. 

P. A. Mirzaei, “Recent challenges in modeling of urban heat island,” Sustain. Cities Soc., 19 200 –206 (2015). https://doi.org/10.1016/j.scs.2015.04.001 Google Scholar

3. 

A. M. Rizwan et al., “A review on the generation, determination and mitigation of the urban heat island,” J. Environ. Sci., 20 (1), 120 –128 (2008). https://doi.org/10.1016/S1001-0742(08)60019-4 Google Scholar

4. 

C. Rinner and M. Hussain, “Toronto’s urban heat island exploring the relationship between land use and surface temperature,” Remote Sens., 3 1251 –1265 (2011). https://doi.org/10.3390/rs3061251 Google Scholar

5. 

M. Zhao et al., “Influence of urban expansion on the urban heat island effect in Shanghai,” Int. J. Geogr. Inf. Sci., 30 (12), 2421 –2441 (2016). https://doi.org/10.1080/13658816.2016.1178389 Google Scholar

6. 

B. Chun and J.-M. Guldmann, “Spatial statistical analysis and simulation of the urban heat island in high-density central cities,” Landsc. Urban Plann., 125 76 –88 (2014). https://doi.org/10.1016/j.landurbplan.2014.01.016 LUPLEZ 0169-2046 Google Scholar

7. 

Z. Dai et al., “Spatial regression models of park and land-use impacts on the urban heat island in central Beijing,” Sci. Total Environ., 626 1136 –1147 (2018). https://doi.org/10.1016/j.scitotenv.2018.01.165 Google Scholar

8. 

J.-P. Kim and J.-M. Guldmann, “Land-use planning and the urban heat island,” Environ. Plann. B, 41 1077 –1099 (2014). https://doi.org/10.1068/b130091p Google Scholar

9. 

H. Zhang et al., “Analysis of land use/ land cover change, population shift, and their effects on spatiotemporal patterns of urban heat islands in metropolitan Shanghai, China,” Appl. Geogr., 44 121 –133 (2013). https://doi.org/10.1016/j.apgeog.2013.07.021 Google Scholar

10. 

W. Q. Zhou et al., “Does spatial configuration matter? Understanding the effects of land cover pattern on land surface temperature in urban landscapes,” Landsc. Urban Plan., 102 (1), 54 –63 (2011). https://doi.org/10.1016/j.landurbplan.2011.03.009 LUPLEZ 0169-2046 Google Scholar

11. 

P. Coseo and L. Larsen, “How factors of land use/land cover, building configuration, and adjacent heat sources and sinks explain Urban Heat Islands in Chicago,” Landsc. Urban Plann., 125 117 –129 (2014). https://doi.org/10.1016/j.landurbplan.2014.02.019 LUPLEZ 0169-2046 Google Scholar

12. 

J. Li et al., “Impacts of landscape structure on surface urban heat islands: a case study of Shanghai, China,” Remote Sens. Environ., 115 3249 –3263 (2011). https://doi.org/10.1016/j.rse.2011.07.008 Google Scholar

13. 

J. Peng et al., “Urban thermal environment dynamics and associated landscape pattern factors: a case study in the Beijing metropolitan region,” Remote Sens. Environ., 173 145 –155 (2016). https://doi.org/10.1016/j.rse.2015.11.027 Google Scholar

14. 

R. Amiri et al., “Spatial-temporal dynamics of land surface temperature in relation to fractional vegetation cover and land use/ cover in the Tabriz urban area, Iran,” Remote Sens. Environ., 113 2606 –2617 (2009). https://doi.org/10.1016/j.rse.2009.07.021 Google Scholar

15. 

J. Song et al., “The relationships between landscape compositions and land surface temperature: quantifying their resolution sensitivity with spatial regression models,” Landsc. Urban Plann., 123 145 –157 (2014). https://doi.org/10.1016/j.landurbplan.2013.11.014 LUPLEZ 0169-2046 Google Scholar

16. 

W. Kuang et al., “What are hot and what are not in an urban landscape: quantifying and explaining the land surface temperature pattern in Beijing, China,” Landsc. Ecol., 30 357 –373 (2015). https://doi.org/10.1007/s10980-014-0128-6 LAECEH 0921-2973 Google Scholar

17. 

K. Deilami and M. Kamruzzaman, “Modelling the urban heat island effect of smart growth policy scenarios in Brisbane,” Land Use Policy, 64 38 –55 (2017). https://doi.org/10.1016/j.landusepol.2017.02.027 Google Scholar

18. 

S. Guha et al., “Dynamic analysis and ecological evaluation of urban heat islands in Raipur City, India,” J. Appl. Remote Sens., 11 (3), 036020 (2017). https://doi.org/10.1117/1.JRS.11.036020 Google Scholar

19. 

Q. Nie et al., “Spatiotemporal impact of urban impervious surface on land surface temperature in Shanghai, China,” Can. J. Remote Sens., 42 (6), 680 –689 (2016). https://doi.org/10.1080/07038992.2016.1217484 CJRSDP 0703-8992 Google Scholar

20. 

G. L. Feyisa et al., “Locally optimized separability enhancement indices for urban land cover mapping: exploring thermal environmental consequences of rapid urbanization in Addis Ababa, Ethiopia,” Remote Sens. Environ., 175 14 –31 (2016). https://doi.org/10.1016/j.rse.2015.12.026 Google Scholar

21. 

J. M. R. Lopez et al., “Frontiers of urbanization: identifying and explaining urbanization hot spots in the south of Mexico City using human and remote sensing,” Appl. Geogr., 79 1 –10 (2017). https://doi.org/10.1016/j.apgeog.2016.12.001 Google Scholar

22. 

H. Pearsall, “Staying cool in the compact city: vacant land and urban heating in Philadelphia, Pennsylvania,” Appl. Geogr., 79 84 –92 (2017). https://doi.org/10.1016/j.apgeog.2016.12.010 Google Scholar

23. 

S. Guha et al., “Analytical study of land surface temperature with NDVI and NDBI using Landsat 8 OLI and TIRS data in Florence and Naples city, Italy,” Eur. J. Remote Sens., 51 (1), 667 –678 (2018). https://doi.org/10.1080/22797254.2018.1474494 Google Scholar

24. 

Y. Y. Cui and B. de Foy, “Seasonal variations of the urban heat island at the surface and the near-surface and reductions due to urban vegetation in Mexico City,” J. Appl. Meteorol. Clim., 51 855 –868 (2012). https://doi.org/10.1175/JAMC-D-11-0104.1 Google Scholar

25. 

W. Zhou et al., “Relationships between land cover and the surface urban heat island: seasonal variability and effects of spatial and thematic resolution of land cover data on predicting land surface temperatures,” Landsc. Ecol., 29 153 –167 (2014). https://doi.org/10.1007/s10980-013-9950-5 LAECEH 0921-2973 Google Scholar

26. 

S. Haashemi et al., “Seasonal variations of the surface urban heat island in a semi-arid city,” Remote Sens., 8 352 (2016). https://doi.org/10.3390/rs8040352 Google Scholar

27. 

B. Chun and J. M. Guldmann, “Impact of greening on the urban heat island: seasonal variations and mitigation strategies,” Comput. Environ. Urban Syst., 71 165 –176 (2018). https://doi.org/10.1016/j.compenvurbsys.2018.05.006 CEUSD5 Google Scholar

28. 

J. Peng et al., “Seasonal contrast of the dominant factors for spatial distribution of land surface temperature in urban areas,” Remote Sens. Environ., 215 255 –267 (2018). https://doi.org/10.1016/j.rse.2018.06.010 Google Scholar

29. 

A. Mathew et al., “Investigating spatial and seasonal variations of urban heat island effect over Jaipur city and its relationship with vegetation, urbanization and elevation parameters,” Sustain. Cities Soc., 35 157 –177 (2017). https://doi.org/10.1016/j.scs.2017.07.013 Google Scholar

30. 

J. Quan et al., “Time series decomposition of remotely sensed land surface temperature and investigation of trends and seasonal variations in surface urban heat islands,” J. Geophys. Res-Atmos., 121 (6), 2638 –2657 (2016). https://doi.org/10.1002/2015JD024354 Google Scholar

31. 

M. Khorchani et al., “Trends in LST over the peninsular Spain as derived from the AVHRR imagery data,” Global. Planet. Change, 166 75 –93 (2018). https://doi.org/10.1016/j.gloplacha.2018.04.006 GPCHE4 0921-8181 Google Scholar

32. 

M. Khorchani et al., “Average annual and seasonal land surface temperature, Spanish Peninsular,” J. Maps., 14 (2), 465 –475 (2018). https://doi.org/10.1080/17445647.2018.1500316 Google Scholar

33. 

R. Yao et al., “Interannual variations in surface urban heat island intensity and associated drivers in China,” J. Environ. Manage., 222 86 –94 (2018). https://doi.org/10.1016/j.jenvman.2018.05.024 Google Scholar

34. 

J. Lai et al., “Does quality control matter? Surface urban heat island intensity variations estimated by satellite-derived land surface temperature products,” ISPRS J. Photogramm. Remote Sens., 139 212 –227 (2018). https://doi.org/10.1016/j.isprsjprs.2018.03.012 IRSEE9 0924-2716 Google Scholar

35. 

R. Sun et al., “A landscape connectivity model to quantify contributions of heat sources and sinks in urban regions,” Landsc. Urban Plann., 178 43 –50 (2018). https://doi.org/10.1016/j.landurbplan.2018.05.015 LUPLEZ 0169-2046 Google Scholar

36. 

S. Khandelwal et al., “Assessment of land surface temperature variation due to change in elevation of area surrounding Jaipur, India,” Egypt. J. Remote Sens. Space Sci., 21 (1), 87 –94 (2018). https://doi.org/10.1016/j.ejrs.2017.01.005 Google Scholar

37. 

Q. Weng, “A remote sensing-GIS evaluation of urban expansion and its impact on surface temperature in Zhujiang Delta, China,” Int. J. Remote Sens., 22 (10), 1999 –2014 (2001). https://doi.org/10.1080/713860788 IJSEDK 0143-1161 Google Scholar

38. 

Q. Weng and S. Yang, “Managing the adverse thermal effects of urban development in a densely populated Chinese city,” J. Environ. Manage., 70 (2), 145 –156 (2004). https://doi.org/10.1016/j.jenvman.2003.11.006 Google Scholar

39. 

Z. Zhang et al., “Towards an operational method for land surface temperature retrieval from Landsat 8 data,” Remote Sens. Lett., 7 (3), 279 –288 (2016). https://doi.org/10.1080/2150704X.2015.1130877 Google Scholar

40. 

Z. Qin et al., “A mono-window algorithm for retrieving land surface temperature from Landsat TM data and its application to the Israel–Egypt border region,” Int. J. Remote Sens., 22 (18), 3719 –3746 (2001). https://doi.org/10.1080/01431160010006971 IJSEDK 0143-1161 Google Scholar

41. 

Z.-L. Li et al., “Satellite-derived land surface temperature: current status and perspectives,” Remote Sens. Environ., 131 14 –37 (2013). https://doi.org/10.1016/j.rse.2012.12.008 Google Scholar

42. 

J. C. Munoz and J. A. Sobrino, “A generalized single channel method for retrieving land surface temperature from remote sensing data,” J. Geophys. Res., 108 (D22), 4688 (2003). Google Scholar

43. 

J. C. Munoz et al., “Revision of the single-channel algorithm for land surface temperature retrieval from Landsat thermal-infrared data,” IEEE Trans. Geosci. Remote Sens., 47 (1), 339 –349 (2009). https://doi.org/10.1109/TGRS.2008.2007125 IGRSD2 0196-2892 Google Scholar

44. 

J. C. Munoz and J. A. Sobrino, “A single-channel algorithm for land-surface temperature retrieval from ASTER data,” IEEE Geosci. Remote Sens. Lett., 7 (1), 176 –179 (2010). https://doi.org/10.1109/LGRS.2009.2029534 Google Scholar

45. 

X. Yu et al., “Land surface temperature retrieval from Landsat 8 TIRS—comparison between radiative transfer equation-based method, split window algorithm and single channel method,” Remote Sens., 6 9829 –9852 (2014). https://doi.org/10.3390/rs6109829 Google Scholar

46. 

O. Rozenstein et al., “Derivation of Land Surface Temperature for Landsat-8 TIRS Using a Split Window Algorithm,” Sensors, 14 (4), 5768 –5780 (2014). https://doi.org/10.3390/s140405768 SNSRES 0746-9462 Google Scholar

47. 

Z. Qin et al., “Derivation of split window algorithm and its sensitivity analysis for retrieving land surface temperature from NOAA-advanced very high resolution radiometer data,” J. Geophys. Res., 106 (D19), 22655 –22670 (2001). https://doi.org/10.1029/2000JD900452 JGREA2 0148-0227 Google Scholar

48. 

A. Bendib et al., “Contribution of Landsat 8 data for the estimation of land surface temperature in Batna city, Eastern Algeria,” Geocarto Int., 32 (5), 503 –513 (2016). https://doi.org/10.1080/10106049.2016.1156167 Google Scholar

49. 

, “Raipur beats Indore, Bhopal in ease of living index,” (2018) https://www.dailypioneer.com/2018/state-editions/raipur-beats-indore-bhopal-in-ease-of-living-index.html August ). 2018). Google Scholar

50. 

, “Home-United States Geological Survey,” http://earthexplorer.USGS.gov Google Scholar

51. 

T. S. Purevdorj et al., “Relationships between percent vegetation cover and vegetation indices,” Int. J. Remote Sens., 19 (18), 3519 –3535 (1998). https://doi.org/10.1080/014311698213795 IJSEDK 0143-1161 Google Scholar

52. 

C. Gao, “NDWI-a normalized difference water index for remote sensing of vegetation liquid water from space,” Remote Sens. Environ., 58 (3), 257 –266 (1996). https://doi.org/10.1016/S0034-4257(96)00067-3 Google Scholar

53. 

Y. Zha et al., “Use of normalized difference built-up index in automatically mapping urban areas from TM imagery,” Int. J. Remote Sens., 24 (3), 583 –594 (2003). https://doi.org/10.1080/01431160304987 IJSEDK 0143-1161 Google Scholar

54. 

L. Wang and J. J. Qu, “NMDI: a normalized multi-band drought index for monitoring soiland vegetation moisture with satellite remote sensing,” Geophys. Res. Lett., 34 (20), L20405 (2007). https://doi.org/10.1029/2007GL031021 GPRLAJ 0094-8276 Google Scholar

55. 

X. C. Chen et al., “Remote sensing image-based analysis of the relationship between urban heat island and land use/cover changes,” Remote Sens. Environ., 104 (2), 133 –146 (2006). https://doi.org/10.1016/j.rse.2005.11.016 Google Scholar

56. 

B. L. Markham and J. K. Barker, “Spectral characteristics of the LANDSAT thematic mapper sensors,” Int. J. Remote Sens., 6 (5), 697 –716 (1985). https://doi.org/10.1080/01431168508948492 IJSEDK 0143-1161 Google Scholar

57. 

G. E. Wukelic et al., “Radiometric calibration of Landsat thematic mapper thermal band,” Remote Sens. Environ., 28 339 –347 (1989). https://doi.org/10.1016/0034-4257(89)90125-9 Google Scholar

58. 

J. A. Sobrino et al., “Land surface temperature retrieval from Landsat TM 5,” Remote Sens. Environ., 90 434 –440 (2004). https://doi.org/10.1016/j.rse.2004.02.003 Google Scholar

59. 

J. A. Sobrino et al., “A comparative study of land surface emissivity retrieval from NOAA data,” Remote Sens. Environ., 75 (2), 256 –266 (2001). https://doi.org/10.1016/S0034-4257(00)00171-1 Google Scholar

60. 

T. N. Carlson and D. A. Ripley, “On the relation between NDVI, fractional vegetation cover, and leaf area index,” Remote Sens. Environ., 62 241 –252 (1997). https://doi.org/10.1016/S0034-4257(97)00104-1 Google Scholar

61. 

J. Li, “Estimating land surface temperature from Landsat-5 TM,” Remote Sens. Technol. Appl., 21 (4), 322 –326 (2006). https://doi.org/10.11873/j.issn.1004-0323.2006.4.322 Google Scholar

62. 

Q. Sun et al., “An ERDAS image processing method for retrieving LST and describing urban heat evolution: a case study in the Pearl River Delta Region in South China,” Environ. Earth Sci., 59 1047 –1055 (2010). https://doi.org/10.1007/s12665-009-0096-3 Google Scholar

63. 

Z. Wan, Z. L. Li, “MODIS land surface temperature and emissivity,” Land Remote Sensing and Global Environmental Change NASA’s Earth Observing System and the Science of ASTER and MODIS, 567 Springer-Verlag, Heidelberg, New York (2011). Google Scholar

64. 

L. Yang et al., “Land surface temperature retrieval for arid regions based on Landsat-8 TIRS data: a case study in Shihezi, Northwest China,” J. Arid Land., 6 704 –716 (2014). https://doi.org/10.1007/s40333-014-0071-z Google Scholar

Biography

Subhanil Guha received his MSc degree in geography from the University of Calcutta, India, and his postgraduate diploma degree in geoinformatics from the International Institute for Geo-Information Science and Earth Observation (ITC), Enschede, The Netherlands, and IIRS, Dehradun, India. He is pursuing his PhD in the Department of Applied Geology, National Institute of Technology (NIT), Raipur, India. His research interests include applications of thermal infrared remote sensing in geography, environment, and geology.

Himanshu Govil received his MSc and PhD degrees in geology from Aligarh Muslim University, Aligarh, India, and his postgraduate diploma degree in geoinformatics from the ITC, Enschede, The Netherlands, and IIRS, Dehradun, India. Currently, he is an assistant professor in the Department of Applied Geology, NIT, Raipur, India. His research interests include application of hyperspectral remote sensing in mineral exploration and geological mapping.

Prabhat Diwan received his MTech and PhD degrees in geology from Dr. Harisingh Gour University, Sagar, India. Currently, he is an associate professor and a head of the Department of Applied Geology, NIT, Raipur, India. His research interests include structural geology, economic geology, and remote sensing.

© 2019 Society of Photo-Optical Instrumentation Engineers (SPIE) 1931-3195/2019/$25.00 © 2019 SPIE
Subhanil Guha, Himanshu Govil, and Prabhat Diwan "Analytical study of seasonal variability in land surface temperature with normalized difference vegetation index, normalized difference water index, normalized difference built-up index, and normalized multiband drought index," Journal of Applied Remote Sensing 13(2), 024518 (13 May 2019). https://doi.org/10.1117/1.JRS.13.024518
Received: 17 December 2018; Accepted: 17 April 2019; Published: 13 May 2019
Lens.org Logo
CITATIONS
Cited by 30 scholarly publications.
Advertisement
Advertisement
RIGHTS & PERMISSIONS
Get copyright permission  Get copyright permission on Copyright Marketplace
KEYWORDS
Earth observing sensors

Landsat

Vegetation

MODIS

Algorithm development

Satellites

Satellite imaging

Back to Top