1. Introduction {#sec1}
===============

Land surface temperature (LST) is an important biophysical parameter in the processes of surface energy and water balance at regional and global scales ([@bib3]; [@bib18]; [@bib29]; [@bib47], [@bib48]). LST is used in a large scale to determine soil moisture content ([@bib45]; [@bib26]; [@bib57]), to analyze the effect of urban heat island ([@bib55]; [@bib42]; [@bib25], [@bib23]; [@bib24]; [@bib59],[@bib61]), to evaluate diurnal temperature variation ([@bib50]; [@bib15]; [@bib1], [@bib2]; [@bib43]) to calculate surface longwave radiation ([@bib52]; [@bib34]), to compute different types of evapotranspiration ([@bib40]; [@bib35]; [@bib38]; [@bib22]), and to estimate surface albedo and thermal inertia ([@bib32]).

Current satellite imageries, such as the Landsat, Moderate Resolution Imaging Spectroradiometer (MODIS) and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), can provide LST at different spatial, spectral, radiometric, and temporal resolutions ([@bib31]; [@bib27]; [@bib20]). Due to technical limitations, these current satellite thermal sensors reflect a balance between spatial and temporal resolutions; i.e., the high spatial resolution sensors generally have a low temporal resolution, and vice versa ([@bib49]). In order to obtain LST at high spatial and high temporal resolutions downscaling technique is considered as an effective method ([@bib17]).

Various types of downscaling techniques have been developed to get high-resolution LST from a coarse resolution thermal infrared band ([@bib46]; [@bib12]; [@bib37]; [@bib56]; [@bib51]). Extended reviews on LST downscaling methods have been assessed so far ([@bib57]; [@bib14]). The most frequently used algorithm is the DisTrad algorithm ([@bib27]) which was modified as the TsHARP algorithm ([@bib2]) and the scale invariant relationship between the normalized difference vegetation index (NDVI) and LST is the foundation of these algorithms. The performance of these different algorithms was assessed over different land use/land cover (LULC) categories, including cultivated land ([@bib26]) or mixed urban land ([@bib20]). The main restriction of the TsHARP algorithm is that this relationship is not unique, and thus, provides an extensive range of LST for a single NDVI value. [@bib30] further enhanced the TsHARP algorithm by considering the effect of photosynthetically and non-photosynthetically active cells of the green plants within the spatial variability of LST. Furthermore, [@bib5] developed a nonlinear DisTrad (NL-DisTrad) model where NDVI-LST polynomial relationship generates from the hot-edge pixels at a coarse resolution which is also valid at high-resolution pixels.

The NDVI is not so appropriate for LST downscaling procedure performed in a heterogeneous urban land ([@bib43]; [@bib34]; [@bib16]; [@bib55]; [@bib20]). Some recent research works reflected modern image classification techniques for detail LULC classification ([@bib10], [@bib11]; [@bib59],[@bib61]). Different LULC types, found in a complex urban land, produce a unique emission rate. Emissivity ([@bib43]; [@bib34]) has been used as one of the most significant environmental factors for highly heterogeneous urban landscapes. [@bib41] observed a close relationship between LST and surface albedo in urban areas. [@bib16] developed high-resolution urban thermal sharpener algorithm for downscaling LST in an urban area by integrating NDVI and surface albedo. Impervious surface percentage ([@bib20], [@bib19]) and pure pixel index ([@bib53]) have also been applied as basic parameters in urban areas. Impervious surface provides a better result than NDVI in LST downscaling in the mixed urban land by comparing 15 different parameters ([@bib20]). [@bib20] and [@bib54] also established a strong linear relationship between LST and impervious surface irrespective of the seasonal influence. Contrary to that, the LST-NDVI relationship changes with the seasonal variability ([@bib27]). For a complex urban area with mixed LULC types, multiple environmental parameters must be integrated to achieve a high downscaling precision. Although various types of LST downscaling methods have been suggested, they have also some limitations in the available remotely sensed data and LULC categories. Besides, selection of appropriate environmental predictors for LST downscaling in complex LULC surface is still considered as an important task.

Basically, the statistical downscaling of LST is based on the correlation between LST and other environmental factors or remote sensing based various LULC indices. Generally, most of the popular LST downscaling methods apply NDVI observed at a fine resolution ([@bib27]), but the NDVI alone cannot be able to explain all the variations in LST for a complex urban surface. [@bib33] evaluated seasonal variation in downscaled LST in DisTrad, TsHARP, and local model using Landsat TM 5 data over a heterogeneous agricultural land in India. [@bib6] attempted to retrieve LST using Landsat TM data for Florence city in Italy and proposed a traditional downscaling framework analyzing its performances using high-resolution LST airborne image. A combination of built-up and vegetation spectral indices was adopted for the LST downscaling methods using MODIS and Landsat TM data performed in Milan city, Italy ([@bib7]). A multiple regression-based LST downscaling technique using spectral mixture analysis over the heterogeneous urban area of Aprilia, Italy was also evaluated to examine the effectiveness of multiple environmental parameters ([@bib8]). Another analytical research was performed in the Zhangye oasis and Beijing city to assess the scale effect in LST downscaling from medium to high-resolution satellite data ([@bib60]).

In the present study, linear regression based TsHARP downscaling technique has been examined with some LULC indices along with NDVI. Basically, the TsHARP technique used fractional vegetation cover and it was analyzed in agricultural field ([@bib1],[@bib2]). Here, the technique has been evaluated on Landsat 8 OLI and TIRS data for four different seasons (pre-monsoon, monsoon, post-monsoon, and winter) in the year of 2014. The main objectives of the research were (1) to estimate the accuracy of the downscaled LST in a heterogeneous urban landscape (Raipur City of India) of humid tropical region using TsHARP technique for separate LULC indices (normalized difference vegetation index (NDVI), normalized difference water index (NDWI), normalized difference built-up index (NDBI), and normalized multi-band drought index (NMDI)) and (2) to assess the seasonal variation of the result in 2014.

2. Materials & methods {#sec2}
======================

The present study has been conducted in Raipur, a tropical city of India. Raipur is considered as a heterogeneous urban area characterized by dry tropical monsoon climate. Raipur has an average elevation of 298 m, is located between 21°11′22″N to 21°20′02″N and 81°32′ 20″E to 81°41′50″E ([Fig. 1](#fig1){ref-type="fig"}). Average annual temperature varies between 34 °C (in summer), and 20 °C (in winter). The area receives 90% of annual rainfall of 130 cm during June to September. This urban area is characterized by four main LULC types, namely, vegetation, water bodies, bare land, and built-up area. A high rate of land conversion (other lands to the built-up area) has been taken place in the recent time periods.Fig. 1Location of the study area with false colour composite image.Fig. 1

Four separate cloud-free Landsat 8 OLI and TIRS images of 5 June 2014, 25 September 2014, 12 November 2014, and 30 December 2014 were used for the evaluation of the TsHARP downscaling technique. Green, red, near infrared, shortwave infrared 1, and shortwave infrared 2 bands of 30 m spatial resolution in OLI sensor were used in determining the LULC indices and thermal infrared band of 100 m spatial resolution in TIRS sensor was used for estimating the LST.

For the convenience of LST downscaling calculation, Landsat 8 TIRS data (100-m resolution) were resampled into 30 m by the nearest neighbor method by simple raster data aggregation to coincide with Landsat 8 OLI data (30-m resolution). These aggregated OLI and TIRS data were further resampled by aggregation into 960 m. The 960-m aggregated coarse-resolution data were used to build the relationship model, whereas the 30-m reflective data were used to downscale the LST at 30 m. The 30 m retrieved LST was used to validate the 30-m downscaled LST. Subsequently, the four indices (NDVI, NDWI, NDBI, and NMDI) with 30-m resolution were then estimated from the 30-m OLI images. The final available output is the downscaled LST of 30-m resolution. The overall accuracy of the entire downscaled LST in various LULC was compared and evaluated.

2.1. Retrieval of LST using Landsat 8 OLI and TIRS data {#sec2.1}
-------------------------------------------------------

Retrieving LST from Landsat 8 OLI and TIRS data includes the brightness temperature calculation and the emissivity correction. The at-sensor brightness temperature assumes that the Earth\'s surface is a black body and includes atmospheric effects (along path absorption and emissions) ([@bib9]). The conversion formula from the at-sensor spectral radiance $L$ to the at-sensor brightness temperature $T_{B}$ is shown in [Eq. (1)](#fd1){ref-type="disp-formula"}. LST is retrieved by [Eq. (2)](#fd2){ref-type="disp-formula"}:$$T_{B} = \frac{K_{2}}{ln\left( {\left( {K_{1}/L_{\lambda}} \right) + 1} \right)}\ $$$$LST = \frac{K_{2}}{1 + \left( {{\lambda\sigma T_{B}}/\left( hc \right)} \right)\mathit{\ln}\text{ε}}\ \ $$Where, $T_{B}$is the brightness temperature in Kelvin (K), $L_{\lambda}$ is the spectral radiance in Wm^−2^sr^−1^mm^−1^; $K_{2}$ and $K_{1}\ $are calibration constants. For Landsat 8 TIRS data, $K_{1}\ $is 774.89, $K_{2}$ is 1321.08. $\lambda$ is the effective wavelength (11.335 μm for band 10 in Landsat 8 OLI data), $\sigma$ is Boltzmann constant (1.38 × 10^−23^ J/K), $h$ is Plank\'s constant (6.626 × 10^−34^ J), $c$ is the velocity of light at a vacuum (2.998 × 10^−8^ m/s), $\text{ε}$ is emissivity. Based on earlier studies, $\text{ε}$ is determined as follows: if NDVI is $\geq$0.157 and $\leq$0.727, $\text{ε}$ for that pixel is calculated by [Eq. (3)](#fd3){ref-type="disp-formula"} ([@bib44]):$$\text{ε} = 1.0094 + 0.0047\ \ln\left( {NDVI} \right)$$

The land cover type should be vegetation if NDVI values $\geq$0.727, and then a constant value of 0.99 is assumed ([@bib42]) and $\text{ε}$ is set to 0.995 for water bodies and 0.92 for the other land use and land cover types.

2.2. Determination of NDVI, NDWI, NDBI, and NMDI as LULC indices {#sec2.2}
----------------------------------------------------------------

NDVI ([@bib39]) is considered as 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. NDWI ([@bib21]) is generally used for water body extraction. NDBI ([@bib58]) is another spectral index which was applied in this study for built-up area extraction. NMDI ([@bib54]) was also used to extract the dry soil or bare land. The formulation of these four indices was presented by [Table 1](#tbl1){ref-type="table"}.Table 1LULC indices used for LST downscaling.Table 1AcronymDescriptionFormulationReferenceNDVINormalized difference vegetation index$\frac{NIR - Red}{NIR + Red}$[@bib39]NDWINormalized difference water index$\frac{Green - NIR}{Green + NIR}$[@bib21]NDBINormalized difference built-up index$\frac{SWIR1 - NIR}{SWIR2 + NIR}$[@bib58]NMDINormalized multi-band difference index$\frac{\left\lbrack NIR - \left( {SWIR1 + SWIR2} \right) \right\rbrack}{\left\lbrack NIR - \left( {SWIR1 + SWIR2} \right) \right\rbrack}$[@bib54]

2.3. TsHARP downscaling technique based on the regression of LST and LULC indices {#sec2.3}
---------------------------------------------------------------------------------

LST can be derived using thermal infrared images with coarse spatial resolutions. Regression models between ancillary environmental predictors and LST have been widely established to enhance LST resolution. If the relationships between LST and the predictors do not change with the variation in the spatial resolution, a detailed LST with a high resolution can be calculated by the predictors using such relationships. [Table 1](#tbl1){ref-type="table"} shows the used LULC indices in the present study for estimating downscaled LST using TsHARP technique.

The original TsHARP algorithm is primarily based on the regression model of LST and fractional vegetation cover. This technique was proposed by [@bib1],[@bib2]). Basically, it was a modification of the DisTrad algorithm ([@bib27]) and was evaluated on agricultural land of Central Iowa, USA. [@bib26] showed that NDVI should be used as a covariate in TsHARP technique. The present study uses NDVI, NDWI, NDBI, and NMDI as the covariates with the LST. The TsHARP algorithm was outlined in Eq. [(4](#fd4a){ref-type="disp-formula"}a--[4](#fd4d){ref-type="disp-formula"}d) where $a_{0}$ is the intercept and $a_{1}$ is the slope of the regression equations. The Eq. [4](#fd4a){ref-type="disp-formula"}a was followed by Eq. [(4](#fd4b){ref-type="disp-formula"}b--[4](#fd4d){ref-type="disp-formula"}d) by replacing the NDVI with NDWI, NDBI, and NMDI, respectively. The fine-resolution (30 m) LST, ($LST_{fine}$) could be determined by the Eq. [(4](#fd4a){ref-type="disp-formula"}a--[4](#fd4d){ref-type="disp-formula"}d):$$LST_{fine} = a_{0} + a_{1}.NDVI_{fine}$$$$LST_{fine} = a_{0} + a_{1}.NDWI_{fine}$$$$LST_{fine} = a_{0} + a_{1}.NDBI_{fine}\ \ $$$$LST_{fine} = a_{0} + a_{1}.NMDI_{fine}\ $$

Therefore, the coarse-resolution LST ($LST_{coarse}$) can be determined by the Eq. [(5](#fd5a){ref-type="disp-formula"}a--[5](#fd5d){ref-type="disp-formula"}d):$$LST_{coarse} = a_{0} + a_{1}.NDVI_{coarse} + \Delta T_{coarse}\ $$$$LST_{coarse} = a_{0} + a_{1}.NDWI_{coarse} + \Delta T_{coarse}\ $$$$LST_{coarse} = a_{0} + a_{1}.NDBI_{coarse} + \Delta T_{coarse}$$$$LST_{coarse} = a_{0} + a_{1}.NMDI_{coarse} + \Delta T_{coarse}$$

Then, a residual of LST ($\Delta T_{coarse}$) was computed as the difference between the retrieved LST ($LST_{coarse}$) and the corresponding observed LST ($LST_{ref}$) by [Eq. (6)](#fd6){ref-type="disp-formula"} ([@bib27]):$$\Delta T_{coarse} = LST_{ref} - LST_{coarse}$$

The residual $\Delta T_{coarse}$ was introduced in the algorithm to take into account part of LST spatial variability that depends on the environmental factors other than the applied predictors, such as soil moisture, emissivity or other LULC indices. The $a_{0}$, $a_{1}$ and residual are different when different indices are used. Finally, downscaled fine-resolution (30 m) LST ($LST_{down}$) was estimated by the Eq. [(7](#fd7a){ref-type="disp-formula"}a--[7](#fd7d){ref-type="disp-formula"}d):$$LST_{down} = a_{0} + a_{1}.NDVI_{fine} + \Delta T_{coarse}\ $$$$LST_{down} = a_{0} + a_{1}.NDWI_{fine} + \Delta T_{coarse}\ $$$$LST_{down} = a_{0} + a_{1}.NDBI_{fine} + \Delta T_{coarse}\ $$$$LST_{down} = a_{0} + a_{1}.NMDI_{fine} + \Delta T_{coarse}\ $$Where, the coarse-resolution regression coefficients were applied to fine-spatial resolution spectral indices, adding the residual error of the corresponding coarse-resolution image to increase the accuracy.

The NDVI-LST relationship always tends to be varied over the mixed urban landscape. In order to overcome the problem, some remote sensing indices may be tested along with NDVI to obtain higher accuracy in the downscaled LST. The multiple least-squares linear regression downscaling method with a number of predictors was applied in some recent studies ([@bib7]; [@bib6]). A number of remote sensing indices have also been tested individually to obtain a better downscaled LST ([@bib20]). In the present study, NDVI, NDWI, NDBI, and NMDI have been examined separately in four different dates for the downscaling method. Finally, the overall procedure of downscaling has been applied at 30 m, 120 m, 240 m, and 480 m spatial resolutions.

2.4. Accuracy assessment and validation {#sec2.4}
---------------------------------------

The LST of 100-m resolution is aggregated to 960, 480, 240, and 120-m resolution and these new aggregated data are known as reference data ($LST_{ref}$). The aggregated 960-m resolution LST was downscaled ($LST_{down}$) to 480, 240, 120, and 30-m resolution. Root mean square error (RMSE) and mean error (ME) statistics were applied to estimate the error in downscaled LST with respect to the reference LST. RMSE and ME have been calculated by the following [Eq. (8)](#fd8){ref-type="disp-formula"} and [Eq. (9)](#fd9){ref-type="disp-formula"}, respectively.$$RMSE = \sqrt{\left\lfloor {n^{- 1}\sum\limits_{i = 1}^{n}\left( {LST_{down} - LST_{ref}} \right)^{2}} \right\rfloor}\ $$$$ME = \sum\limits_{i = 1}^{n}\left( {LST_{down} - LST_{ref}} \right)\ $$

A complete flowchart of the entire study has been shown in [Fig. 2](#fig2){ref-type="fig"}. The computation cost of the proposed method was small as it required less time and storage space for the entire computation process.Fig. 2Flowchart of methodology.Fig. 2

3. Results & discussion {#sec3}
=======================

[Fig. 3](#fig3){ref-type="fig"} represented the spatial distribution of reference LST at 30 m resolution and [Table 2](#tbl2){ref-type="table"} showed the statistical information of LST at 30-m resolution in four multi-date images. Seasonal variation in the LST distribution shows a specific thermal pattern. The mean LST values in pre-monsoon, monsoon, post-monsoon, and winter season are 33.59 °C, 29.56 °C, 23.40 °C, and 22.26 °C respectively. The range of temperature is found as 13.50 °C in pre-monsoon, 12.28 °C in monsoon, 8.83 °C in post-monsoon, and 10.12 °C in the winter season, respectively. Basically, this type of heterogeneity in LST was observed due to the changes in vegetation abundance and soil moisture content. Monsoon and post-monsoon seasons are characterized by healthy vegetation and wet soil. Winter remains comparatively dry and having least standard deviation value in LST ([Table 2](#tbl2){ref-type="table"}).Fig. 3Spatial distribution of reference LST at 30-m spatial resolution: (a) pre-monsoon; (b) monsoon; (c) post-monsoon; (d) winter.Fig. 3Table 2Spatial distribution of LST (°C) for four multi-date reference images at 30 m resolution.Table 2Date of AcquisitionTimePath/RowLST (Minimum)LST (Maximum)LST (Mean)LST (Standard deviation)05-JUN-201404:55:45142/04425.7739.2733.591.6125-SEP-201404:56:11142/04424.5136.7929.561.7312-NOV-201404:56:21142/04419.3928.2323.401.1230-DEC-201404:56:09142/04417.6427.7622.261.11

Figs. [4](#fig4){ref-type="fig"}, [5](#fig5){ref-type="fig"}, [6](#fig6){ref-type="fig"}, and [7](#fig7){ref-type="fig"} and [Table 3](#tbl3){ref-type="table"} presents the downscaled LST using various LULC indices based TsHARP algorithm along with the retrieved reference LST in four different seasons (pre-monsoon, monsoon, post-monsoon, and winter). It is observed that the first panel (LST reference) has a smoothed pattern since the 30-m LST resolution is not an actual measurement but a resampling, and that the downscaling improve the detail at this spatial scale. It is clearly revealed that NDBI and NMDI based TsHARP algorithm generates almost similar nature in the downscaled LST while NDVI and NDWI based TsHARP algorithm has almost identical spatial distribution of downscaled LST. This particular type of spatial pattern is reflected in each and every season. The spatial pattern of downscaled LST generates from NDBI and NMDI based TsHARP technique has a similarity with retrieved reference LST where the study area achieves a higher temperature. Monsoon season indicates most identical scenario when reference LST compares with the downscaled LST, irrespective of any LULC indices based TsHARP technique.Fig. 4Spatial distribution of retrieved LST and downscaled LST at 30-m resolution for 5-JUN-14: (a) LST~ref~; (b) NDBI-based LST~down~; (c) NDVI-based LST~down~; (d) NDWI-based LST~down~; (e) NMDI-based LST~down~.Fig. 4Fig. 5Spatial distribution of retrieved LST and downscaled LST at 30-m resolution for 25-SEP-14: (a) LST~ref~; (b) NDBI-based LST~down~; (c) NDVI-based LST~down~; (d) NDWI-based LST~down~; (e) NMDI-based LST~down~.Fig. 5Fig. 6Spatial distribution of retrieved LST and downscaled LST at 30-m resolution for 12-NOV-14: (a) LST~ref~; (b) NDBI-based LST~down~; (c) NDVI-based LST~down~; (d) NDWI-based LST~down~; (e) NMDI-based LST~down~.Fig. 6Fig. 7Spatial distribution of retrieved LST and downscaled LST at 30-m resolution for 30-DEC-14: (a) LST~ref~; (b) NDBI-based LST~down~; (c) NDVI-based LST~down~; (d) NDWI-based LST~down~; (e) NMDI-based LST~down~.Fig. 7Table 3Seasonal variation in different LULC indices based downscaled LST (°C) at various resolutions.Table 3LULC Indices- based downscaled LST (°C)at 30 mat 120 mat 240 mat 480 mJune-05-2014MinMaxμσMinMaxμσMinMaxμσMinMaxμσNDBI-based25.7739.2733.591.6126.3338.1433.621.5827.0138.0633.641.5227.9837.0933.691.42NDVI-based25.1245.7933.641.3425.6938.4133.661.1228.8736.5233.670.9830.2636.2933.700.88NDWI-based28.6137.1433.740.9129.0036.6133.740.7730.5535.8733.740.6931.4135.5733.740.62NMDI-based29.5437.6533.740.8329.8237.2933.740.7231.3136.5133.740.6431.9036.0833.740.58NDBI-based26.2838.7733.640.9326.2838.7833.640.9330.9635.9833.670.7231.3935.6433.690.66September-25-2014MinMaxμσMinMaxμσMinMaxμσMinMaxμσNDBI-based24.5136.7929.561.7324.5935.8929.541.7124.6734.7629.521.6525.2634.4929.461.57NDVI-based21.2038.2829.502.0223.3834.7729.491.7323.6833.0629.471.5124.9632.2929.441.40NDWI-based23.3337.2229.492.0024.0635.3829.471.7024.3834.8729.461.5125.4933.1929.431.34NMDI-based23.4031.0229.471.9024.3736.5729.461.6124.7236.5129.451.4225.5934.0729.421.25NDBI-based18.8939.4629.471.8722.3835.5629.451.5822.7533.4429.441.4124.5433.2229.421.26November-12-2014MinMaxμσMinMaxμσMinMaxμσMinMaxμσNDBI-based19.3928.2323.401.1219.5327.9223.391.1119.6627.3723.381.0720.2626.9623.361.02NDVI-based15.2730.1523.341.2815.5927.5823.331.0319.8126.1623.330.8819.8125.2623.320.74NDWI-based18.3727.4723.381.0519.3826.4123.370.9320.5226.0823.360.8621.4825.5923.340.79NMDI-based18.1527.5823.381.0219.6426.8423.370.9120.7726.8023.360.8421.5925.9523.340.78NDBI-based20.0025.1023.300.2821.3124.3923.310.2222.4524.0723.310.1922.4823.9723.310.15December-30-2014MinMaxμσMinMaxμσMinMaxμσMinMaxμσNDBI-based17.6327.7622.261.1117.6826.8622.271.0817.7825.7922.281.0319.0825.4122.290.94NDVI-based13.4828.5622.291.0615.0726.2622.290.8418.7125.2922.300.7019.7125.2922.300.60NDWI-based20.7823.2922.310.2521.0323.1422.310.2221.5623.0722.310.2021.8322.8022.310.19NMDI-based21.4722.9322.310.1421.6822.9022.310.1321.9522.8622.310.1222.0222.6422.300.11NDBI-based20.3823.3322.290.1721.3022.9322.290.1421.7922.7922.290.1321.8922.7522.300.11[^1]

[Table 4](#tbl4){ref-type="table"} presents the seasonal variation of the estimated errors (RMSE and ME) produced in downscaled LST by various types of LULC indices based TsHARP algorithm. Downscaled LST at lower spatial resolution generates a greater error than downscaled LST at higher spatial resolution. It is a common phenomenon observed in each and every season and it remains constant for any LULC indices based TsHARP downscaling technique. NDBI-based TsHARP algorithm showed the best result among all the LULC indices-based TsHARP algorithm for all the season and at every resolution level as the values of RMSE of downscaled LST lie significantly below the standard deviation for the corresponding actual native LST reference images at 100 m spatial resolution in each and every occasion. NDBI proves as the most successful parameter for LST downscaling in almost every season. Since it is an urban area having so many mixed landscapes built-up area has the most dominant land use in controlling LST. Previously, NDVI-LST relationship based TsHARP algorithm proved as one of the most effective downscaling methods for single LULC dominated region ([@bib1],[@bib2]; [@bib26]). In the present study, NDVI also proves as an important indicator in determining the downscaled LST values in a heterogeneous urban landscape. Along with NDBI and NDVI, two more remote sensing indices (NDWI and NMDI) were also examined to compare the accuracy level in downscaled LST. Post-monsoon and winter season shows a better accuracy in terms RMSE and ME for any of the indices based downscaled LST. NDBI and NDMI reflect a smaller ME for pre-monsoon and winter seasons while NDVI and NDWI indicates a smaller ME in monsoon and post-monsoon seasons as an indicator used to develop TsHARP downscaling technique. Figs. [8](#fig8){ref-type="fig"}, [9](#fig9){ref-type="fig"}, [10](#fig10){ref-type="fig"}, and [11](#fig11){ref-type="fig"} and Figs. [12](#fig12){ref-type="fig"}, [13](#fig13){ref-type="fig"}, [14](#fig14){ref-type="fig"}, and [15](#fig15){ref-type="fig"} show the spatial distribution of RMSE and ME, respectively.Table 4Seasonal variation in mean RMSE and mean ME for downscaled LST (^o^C) at various resolution.Table 4LULC indices based LST~down~at 30 mat 120 mat 240 mat 480 mJune-05-2014 (the standard deviation of the native LST reference image at 30 m spatial resolution is 1.61)RMSEMERMSEMERMSEMERMSEMENDBI-based1.090.050.960.040.910.030.870.01NDVI-based1.220.141.170.121.130.091.050.05NDWI-based1.230.151.190.131.150.101.070.05NMDI-based1.090.051.030.040.990.020.960.01September-25-2014 (the standard deviation of the native LST reference image at 30 m spatial resolution is 1.73)RMSEMERMSEMERMSEMERMSEMENDBI-based1.11-0.050.81-0.050.69-0.050.61-0.03NDVI-based1.21-0.070.97-0.060.86-0.060.75-0.04NDWI-based1.25-0.081.08-0.080.98-0.070.86-0.04NMDI-based1.18-0.090.92-0.080.82-0.070.75-0.05November-12-2014 (the standard deviation of the native LST reference image at 30 m spatial resolution is 1.12)RMSEMERMSEMERMSEMERMSEMENDBI-based0.91-0.060.72-0.050.65-0.050.60-0.05NDVI-based0.75-0.010.66-0.020.60-0.020.55-0.03NDWI-based0.74-0.010.66-0.020.60-0.020.55-0.02NMDI-based0.81-0.090.79-0.080.78-0.070.76-0.06December-30-2014 (the standard deviation of the native LST reference image at 30 m spatial resolution is 1.11)RMSEMERMSEMERMSEMERMSEMENDBI-based0.850.030.710.030.650.020.590.01NDVI-based0.860.050.840.050.790.040.730.01NDWI-based0.860.050.840.040.800.030.720.01NMDI-based0.810.030.790.030.760.020.690.01Fig. 8Spatial distribution of ME at 30-m resolution for 5-JUN-14: (a) NDBI-based LST~down~; (b) NDVI-based LST~down~; (c) NDWI-based LST~down~; (d) NMDI-based LST~down~.Fig. 8Fig. 9Spatial distribution of ME at 30-m resolution for 25-SEP-14: (a) NDBI-based LST~down~; (b) NDVI-based LST~down~; (c) NDWI-based LST~down~; (d) NMDI-based LST~down~.Fig. 9Fig. 10Spatial distribution of ME at 30-m resolution for 12-NOV-14: (a) NDBI-based LST~down~; (b) NDVI-based LST~down~; (c) NDWI-based LST~down~; (d) NMDI-based LST~down~.Fig. 10Fig. 11Spatial distribution of ME at 30-m resolution for 30-DEC-14: (a) NDBI-based LST~down~; (b) NDVI-based LST~down~; (c) NDWI-based LST~down~; (d) NMDI-based LST~down~.Fig. 11Fig. 12Spatial distribution of RMSE at 30-m resolution for 5-JUN-14: (a) NDBI-based LST~down~; (b) NDVI-based LST~down~; (c) NDWI-based LST~down~; (d) NMDI-based LST~down~.Fig. 12Fig. 13Spatial distribution of RMSE at 30-m resolution for 25-SEP-14: (a) NDBI-based LST~down~; (b) NDVI-based LST~down~; (c) NDWI-based LST~down~; (d) NMDI-based LST~down~.Fig. 13Fig. 14Spatial distribution of RMSE at 30-m resolution for 12-NOV-14: (a) NDBI-based LST~down~; (b) NDVI-based LST~down~; (c) NDWI-based LST~down~; (d) NMDI-based LST~down~.Fig. 14Fig. 15Spatial distribution of RMSE at 30-m resolution for 30-DEC-14: (a) NDBI-based LST~down~; (b) NDVI-based LST~down~; (c) NDWI-based LST~down~; (d) NMDI-based LST~down~.Fig. 15

The high values of RMSE (positive) and ME (positive or negative) have mostly been found in the areas where LST is very high or abnormally low. Due to the presence of mineral-based industrial agglomeration or thermal power plants in the northwest portions of the study area, LST generally remains high. This portion has come under one of the most erroneous zones regarding the downscaled LST. It happens mainly due to the anthropogenic activities and no such hypothetical relationship is built to support the particular behaviour.

[Table 5](#tbl5){ref-type="table"} presents the percentage of RMSE in downscaled LST at 30 m spatial resolution for different seasons. In pre-monsoon season, NMDI and NDBI based TsHARP downscaling technique provide almost 73% of the total pixels in the study area having less than 1.5 °C RMSE in downscaled LST. NDBI again proves as the best parameter for TsHARP downscaling algorithm in monsoon season as 72.7% pixels having less than 1.5 °C RMSE. Post-monsoon season reflects the best result regarding RMSE in downscaled LST. In post-monsoon season, NDVI and NDWI based TsHARP algorithm have proved as the best output at 30 m resolution irrespective of any season as they generates downscaled LST with more than 90% pixels having less than 1.5 °C RMSE. In winter season, NDVI and NMDI based TsHARP techniques give the best result (\>86% pixels have less than 1.5 °C RMSE).Table 5Percentage of root mean square error (RMSE) in downscaled LST (^o^C) at 30 m resolution.Table 5RMSE in downscaled LST (°C)0.0--0.50.5--1.01.0--1.51.5--2.02.0--2.5\>2.5June-05-2014NDBI-based28.70%24.43%19.52%13.66%7.34%6.35%NDVI-based24.00%23.13%19.00%16.33%9.43%6.46%NDWI-based23.94%23.10%20.60%16.46%9.54%6.36%NMDI-based29.46%25.80%18.61%11.86%6.86%7.41%September-25-2014NDBI-based28.82%25.27%18.61%11.91%7.28%8.13%NDVI-based29.49%24.26%17.51%11.49%7.03%10.22%NDWI-based28.66%24.09%17.78%11.87%7.06%10.54%NMDI-based28.29%23.80%17.92%12.29%7.72%9.98%November-12-2014NDBI-based30.07%27.52%17.58%9.91%5.19%3.74%NDVI-based43.58%31.52%14.98%5.50%1.81%2.61%NDWI-based46.43%30.42%13.77%4.91%1.60%2.87%NMDI-based41.28%30.01%14.26%7.18%3.47%3.80%December-30-2014NDBI-based36.98%28.71%18.75%9.13%3.74%2.70%NDVI-based35.47%31.16%19.47%7.30%3.29%3.35%NDWI-based36.28%30.67%18.31%7.54%3.62%3.58%NMDI-based39.58%30.55%16.15%7.14%3.64%2.94%

Seasonal variation in the percentage of mean error (ME) in downscaled LST at 30 m resolution has been shown in [Table 6](#tbl6){ref-type="table"}. It is very clear from [Table 6](#tbl6){ref-type="table"} that all the LULC indices are almost equally suitable for TsHARP downscaling technique. Post-monsoon season simply generates the best result for downscaled LST due to low pollution level, high vegetation coverage and high moisture content. NDWI and NDVI based TsHARP downscaling algorithm generates lowest level of mean error \[\>75% pixels have a ME from 0 °C to 1 °C (positive/negative)\].Table 6Percentage of mean error (ME) in downscaled LST (°C) at 30 m resolution.Table 6ME in downscaled LST (°C)\< -4-4---3-3---2-2---1-1--00--11--22--33--4\>4June-05-2014NDBI-based0.40%1.61%5.90%13.79%23.35%29.78%19.39%4.80%0.61%0.36%NDVI-based0.01%0.29%5.00%19.48%25.89%21.25%17.53%7.89%1.06%1.63%NDWI-based0.01%0.32%5.47%19.18%25.36%21.67%17.88%7.18%1.01%1.92%NMDI-based0.57%2.26%6.93%10.70%21.39%33.86%19.77%3.59%0.58%0.35%September-25-2014NDBI-based0.61%2.15%6.65%14.35%24.53%29.56%16.16%4.77%1.01%0.20%NDVI-based0.62%2.29%6.99%16.35%27.62%26.13%12.65%4.06%1.09%2.20%NDWI-based0.47%2.06%7.22%18.24%29.85%22.90%11.41%4.00%1.23%2.62%NMDI-based1.11%2.68%6.79%14.37%25.42%26.68%15.84%5.47%1.36%0.29%November-12-2014NDBI-based0.14%0.85%4.57%13.97%30.98%32.6%13.52%2.88%0.41%0.06%NDVI-based0.01%0.04%1.03%12.45%41.24%33.86%8.03%1.52%0.78%1.04%NDWI-based0.01%0.04%0.95%11.73%43.41%33.49%6.95%1.33%0.74%1.41%NMDI-based0.18%1.55%3.80%10.12%34.98%36.31%11.32%1.73%0.01%0.01%December-30-2014NDBI-based0.19%0.80%3.68%11.13%29.75%35.93%16.75%1.62%0.12%0.03%NDVI-based0.03%0.37%3.11%12.55%30.79%35.84%14.17%1.86%0.74%0.54%NDWI-based0.04%0.48%3.43%12.30%30.19%36.76%13.55%2.02%0.75%0.46%NMDI-based0.07%0.65%3.62%11.19%28.49%41.64%12.10%1.68%0.51%0.05%

The study was evaluated the downscaling LST in a humid tropical city of India and reflected a good result. The similar types of studies were also performed successfully in the cities of other climatic zones, e.g., humid subtropical zone ([@bib36]), temperate zone ([@bib4]; [@bib6]).

The proposed work also have some limitations. The 30 m LST resolution is not an actual measurement but a resampling, and the downscaling improve the detail at this spatial scale. Only one date of a particular season may not be appropriate for seasonal assessment of downscaling technique. The other relevant LULC indices may also be incorporated to improve the downscaling results.

4. Conclusion {#sec4}
=============

The aim of this research is to estimate remote sensing indices (NDBI, NDVI, NDWI, and NMDI) based downscaled LST at 30, 120, 240, and 480 m resolution and compare the results in different seasons (pre-monsoon, monsoon, post-monsoon, and winter) over a heterogeneous urban area like Raipur city. The NDVI-LST linear relationship is the primary basis of TsHARP downscaling technique. NDBI, NDWI, and NMDI are tested separately in place of NDVI to form new linear relationships, i.e., NDBI-LST, NDWI-LST, and NMDI-LST. The downscaled LST based on NDBI-LST, NDVI-LST, NDWI-LST, and NMDI-LST was finally evaluated through error estimation. The seasonal variation of the results was also determined. The results were examined at various spatial resolutions. It is clear from the various results that NDBI based TsHARP technique provides the lowest RMSE in downscaled LST at 30 m spatial resolution, irrespective of all seasons. In pre-monsoon, monsoon, and winter seasons, NDBI and NMDI based TsHARP technique have a smaller range of ME, but, in post-monsoon season, NDVI and NDWI indicate small ME due to the presence of a higher percentage of chlorophyll and moisture content. The RMSE and ME became gradually smaller with the increase of spatial resolution in downscaled LST. The LULC based TsHARP models may also be examined in determining the downscaled LST for daytime and nighttime thermal data. Different satellite sensors may also be used to assess the capability of this technique in a separate study area. Any further modification in the present TsHARP algorithm will be expected by the future researchers.

Declarations {#sec6}
============

Author contribution statement {#sec6.1}
-----------------------------

Guha, S. Govil, H. Dey, A. Gill, N.: Conceived and designed the analysis; Performed the experiments; Analyzed and interpreted the data; Contributed analysis tools or data; Wrote the paper.

Funding statement {#sec6.2}
-----------------

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Competing interest statement {#sec6.3}
----------------------------

The authors declare no conflict of interest.

Additional information {#sec6.4}
----------------------

No additional information is available for this paper.

The authors thank to United States Geological Survey and National Institute of Technology Raipur, India. The authors also thank the anonymous reviewers for their relevant comments.

[^1]: μ = mean, and σ = standard deviation.
