CES Technical Report 131,   August 2013
¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
Ecological Modelling and Energy DSS
Energy & Wetlands Research Group, Centre for Ecological Sciences, Indian Institute of Science, Bangalore - 560012, INDIA
Land Surface Temperature Analysis of Himachal Pradesh through Multi-Resolution, Spatio – Temporal data
                 BACK  «  HOME  »  NEXT

Abstract

Rapid changes in the land use and land cover of a region have become a major environmental concern in recent times.  This has lead to unsustainable  development with the reduction of green spaces and also changes in local climate and formation of urban heat islands (UHI). Monitoring and management of land use dynamics would help in land use planning and mitigation of environmental impacts. The main goal of this paper is to quantify the changes in the land cover and consequent changes in land surface temperature. Land use and land cover dynamics were assessed using temporal remote sensing data (Landsat Thematic Mapper and Enhanced Thematic Mapper data) of Himachal Pradesh, India.  The thermal  infrared bands of  the Landsat data  were used  to retrieve Land Surface Temperature. The results revealed that there was a huge increase in Urban Area (including Barren land), which is the causal factor for the changes in Land Surface Temperature. Overall, Remote  sensing and Geographic Information System technologies were effective approaches for monitoring and analyzing urban growth patterns and evaluating their impacts on Land Surface Temperature.

Keywords : Land Surface Temperature, Landsat, Land Use,  Land cover, NDVI.

Introduction

Land cover changes induced by human and natural processes play a major role in global as well as at regional scale patterns of the climate and biogeochemistry of the Earth system. Studies have revealed of changes in water cycling process between land and atmosphere due to the large scale land cover changes, affecting the local to regional climate. Geoninformatics technologies such as remote  sensing and geographic  information systems are very  effective in measuring monitoring and predicating the land use/cover  changes (Qihao Weng, 2009; Enner  et al., 2010). Timely information with higher accuracy of landuse (LU)  Land cover (LC) changes is crucial for long-term planning, economic development, and sustainable management(Zhang et al., 2010) of natural resources. The analysis of temporal remote sensing data helps in understanding the land cover changes and its impact on the environment.  The thermal  infrared bands of  remote sesning data of space borne sensors help to retrieve Land Surface Temperature (LST). Land surface temperature is the measure of the heat emision from land surface due to various activities associated with the land surface. Increase in paved land cover  is an indication of concentrated human activities, which often leads to increased LST’s (Ramachandra  and Kumar, 2009). Increased LST in certain urban pockets in comparison  to its surroundings consequent to the increase in paved surfaces is known as urban heat island (UHI)  phenomenon (Landsbeg, 1981; Gallo and Owen, 1998, Li et al., 2004).

Detection of the thermal charecteristics of land surface using remote sensing data of space borne sensors and the analysis of land surface temperature (LST) has been reported earlier (Rao, 1972).  Spatio-temporal data were used to develop models of land surface atmosphere exchange, and to analyze the relationship between temperature and land use and land cover (LULC) in urban areas (Voogt and  Oke, 2003) highlighting the relationship between LST and surface characteristics such as vegetation indices (Carlson et al., 1994; Carlson 2007; Owen et al., 1998). Also, the studies reveal the effect of biophysical factors on LST by using vegetation fraction instead of qualitative LULC classes (Gallo and Tarpley, 1996; Owen et al., 1998; Dousset and Gourmelon, 2003). The vegetation index–LST relationship has also been used to retrieve surface biophysical parameters (Carlson et al., 1994), to extract sub-pixel thermal variations (Kustas et al., 2003), and to analyze land cover dynamics (Lambin and Ehrlich, 1996). Many investigators have observed a negative relationship between vegetation index and LST, leading to further research into two major pathways, namely, statistical analysis of the relationship and the temperature/vegetation index (TVX) approach. TVX is a multi-spectral method of combining LST and a vegetation index (VI) to monitor their associations.

Land surface and atmospheric temperatures rises by various  antropogenic activities like increased land surface coverage by artificial  materials, energy consumption, which have a high heat capacity and conductivity, and is also associated with the decreases in vegetation and water  surfaces, which are the major factors that reduce surface temperature  through evapo-transpiration (Kato and Yamaguchi., 2005). Temperatures  can be monitored through space borne remote sensing(rs) sensors, which account for the top of the atmosphere (TOA) radiances in the thermal infrared (TIR) region. TOA radiance is the net radiance of the emitted radiance from the earth’s surface upwelling radiance from the atmosphere, and downwelling radiance from the sky. The brightness temperatures (also known as blackbody temperatures) can also be derived from the TOA radiance (Dash et al., 2002). These brightness temperatures account for the various properties of the land surface , the amount and nature of vegetation cover , the thermal properties and moisture content of the soil (Friedl, 2002). However, lack of knowledge of spectral  emissivity can introduce an error which ranges  from 0.2 to 1.2k for mid-latitude summers and  0.8 to 1.4k for the winter conditions for an emissivity of 0.98 and at the ground height of 10km (Dash et al., 2002). Two approaches have been developed to recover LST from multispectral TIR imagery (Schmugge et al., 1998) as on date. The first approach utilises the radiative transfer equation to correct the at sensor radiance to surface radiance, followed by an emissivity model to separate the surfaces radiance into temperature and emissivity (Friedl, 2002). The second approach applies the split window technique for sea surfaces to land surfaces.  Assuming that the emmisivity in the channels used for the split window is similar (Dash et al., 2002). TIR region corresponding  to 8-14 µm in the electromagnetic spectrum is being used quantifying the thermal urban environment.

Landsat satellites, are one of the most  widely used for environmental studies. Landsat thematic mapper is composed by seven bands, six of them in the visible and near infrared, and only one band located in the thermal infrared region (with an effective wavelength of 11.457 µm) is used for LST retrieval. Possesing only one thermal band  might be stated as an disadvantage/limitation in order to obtain LST as it does not allow the application of a split-window method (Sobrino et al., 1996) neither a temperature/emissivity separation method (Gillespie et al., 1998; Snyder et al., 1998) to obtain information about the emissivity spectrum of natural surfaces.

The objective  here is to investigate the Land Surface Temperature with Land use dynamics to understand the Urban Heat  Island phenomenon in Himachal Pradesh considering  Multi -senor, Multi-resolution and temporal RS data acquired through space borne sensors.  This  involved:

  1. Temporal  LU change  analysis (during  1989 to 2005)
  2. Computation  of LST and NDVI (Normalised  Difference  Vegetation Index) from Landsat Tm(1989) and Landsat ETM(2000) and Landsat ETM+(2005) Data;
  3. Investigation of the  role  of NDVI in LST;

Study area:  This analysis has been carried out for Himachal Pradesh, which  lies between the Latitudes :30o 22' 40"N  to 30o 12' 40"N, Longitude:75o 47' 55"E to 79o 04' 20"E (figure 1). LST was computed for a region lies between the latitudes: 30o 18' 30"N to 30o 10' 30"N and Longitude: 76o 19' 35"E to 78o 59' 10"E and covers Shimla.

Figure 1: Study area - Himachal Pradesh, India

Himachal comprises of 12 districts in total, covering the area of 55,673 sq km. Its total population is 6,077,248 as per 2001 census. Shimla is the state capital of Himachal Pradesh, with its population around  7,21,745 as per 2001 census with geographical area of 5131 sq. km. Density of population in Himachal Pradesh is 109 per sq. kms. Himachal Pradesh is one of the Major Indian state undergoing rapid urbanisation. Labour Force has about 49.3% of total population employed in Industrial sectors. Industries located in Himachal as on 31/03/2010 are 36845 (micro, small, medium and large Enterprises) of which 444 are in medium and  large scale industries.

Data:  The  data  used were from Landsat series thematic mapper (28.5m) and enhanced thematic mapper (28.5m) data  acquired over years 1989 to 2010.  Collateral data includes google earth imagery (http://earth.google.com) and Survey of India top sheets. 

Methodology: The RS data used to study the temporal changes in landscape pattern were Landsat Thematic Mapper (TM), Landsat Enhance TM Plus (ETM+) of 1989 to 2006. The data were georeferenced, rectified and cropped pertaining to the study area. Landsat ETM+ bands of 2010 were corrected for the SLC-off by using image enhancement techniques, followed by nearest-neighbour interpolation.

  1. Land use and land cover analysis: This was carried out using data of Landsat satellite using Supervised Pattern classification using Maximum likelihood classifier. This method has already been proved as a superior method as it uses various classification decisions using Probability and Cost functions (Duda et al., 2000). Mean and covariance matrix are computed using estimate of Maximum likelihood estimator. Application of this method resulted in accuracy of about 75% in all the datasets. For the purpose of accuracy assessment, a confusion matrix was calculated. Land Use was computed using the temporal data through open source program GRASS: Geographic Resource Analysis Support System (http://grass.fbk.eu/)
  2. Calculation of Land surface Temperature from Landsat data:  LST was computed (Weng etal.,2004) from TIR bands (Landsat  TM and ETM). Emissivity corrections for specified  LC is carried out using surface Emissivity  as per Synder et al.,(1998); Stathopoplou et al., (2007) and land surface temperature is calculated as per Artis and Carnahan (1982) and Carnahan and Larson, 1990.

                                                       ……. ………. 1
Where, λ the wavelength of emitted radiance for which the peak response and average of the limiting wavelength (λ = 11.5µm) (Markham and Barkar, 1985) were used, ρ=1.439 x 10-2 mk and E=Spectral Emissivity.

LC was determined through the computation of Normalised Difference Vegetation Index (NDVI) using   Landsat visible Red (0.63 – 0.69 µm) and near-infrared band (0.76 – 0.9 µm) bands of Landsat TM/ETM. NDVI was computed in order to calculate emissivity for computing LST.

NDVI is given by

                                      ……………. (2)

Results and Discussion

The Classified images from 1989 – 2006 showed an overall accuracy of 76%. Land Use changes were more prominent in the area during the last 2 decades consequent to increase in Barren land as indicated in Table 1 and Represented in Figure 2. The LU analysis shows that there has been 55% increase in urban  and open areas during 1989-2000 and 39% increase during 2000-2005, and 18.92% increase during 2005-06. This land cover changes has also influenced the local climate. The minimum (min) and maximum (max) temperature was found to be -2º C and 31º C from Landsat Data as tabulated in the Table 2. The analysis showed that there has been an increase in Temperature from 1989 to 2006 as evident from temporal analysis.

Figure 2. Land use classification in 1989, 2000, 2005 and 2006

Table 1: Land use changes in Himachal Pradesh

Class Area in Hectares
(1989)
Area in Hectares (2000) Area in Hectares (2005) Area in Hectares
(2006)
Water 16485.679 16039.310 10669.536 9454.235
Vegetation 1774602.387 1537166.846 1250625.535 1050875.265
Snow 142757.716 107648.808  99234.989 98756.325
Rock and Urban 491566.388 764532.924 1064898.297 1266326.257
Total 2425412.172 2425412.172 2425412.172 2425412.172

Table 2: Land surface Temperature changes in Himachal Pradesh

Year Temperature (Maximum) in Degree Celsius Temperature (Minimum) in Degree Celsius
1989 27  3
2000 28 -1
2005 30 -1
2006 31 -2

The temperature drawn by the study are,  in the year 1989 minimum  30 C and  maximum  270 C, where as  in the year 2006 the minimum -20 C and  maximum 310 C ( given in Table 2), which highlights that LU characteristics play a significant role in maintaining the ambient temperature and also in the regional heat island phenomenon. Correlation analysis shows that the vegetation and water bodies are negatively correlated with temperature suggesting that these LU aid as heat sinks and hence maintains the regional climate.   The population of this region exceeds its carrying capacity and is exerting pressure on the local natural resources such as land, water, etc. (MOEF,1992; Graymorea et al., 2009).

Conclusion:

The LU analysis shows that there has been 55% increase in urban  and open areas during 1989-2000 and 39% increase during 2000-2005, and 18.92% increase during 2005-06. This land cover changes has also influenced the local climate. The minimum (min) and maximum (max) temperature was found to be -2º C and 31º  C. This study clearly shows that the rate of increase in the urbanisation leads  to change in the Land Surface Temperature. The increase in the temperature is in the range of   30 C to 40 C  during 1989 to 2010.  Change detection techniques from multi-resolution images integrating spectral, structural and textural features to generate changed patches and change attribute is also desirable and a challenging area of research.

Acknowledgement

We are grateful to NRDMS Division, The Ministry of Science and Technology, Government of India for the financial support.

References

  1. Artis, D.A., and Carnahan, W.H., 1982, Survey of emissivity variability in thermography of urban areas. Remote Sensing of Environment, 12, 13-329.
  2. Carnahan, W.H., and Larson, R.C., 1990, An analysis of an urban heat sink. Remote Sensing of Environment, 33, 65-71.
  3. Carlson, T. N., Gillies, R. R., and Perry, E. M., 1994., A method to make use of thermal infrared temperature and NDVI measurements to infer surface water content and fractional vegetation cover. Remote Sensing Reviews, 9, 161−173.
  4. Carlson, T. N., 2007. An overview of the so-called “Triangle Method” for estimating surface evapotranspiration and soil moisture from satellite imagery. Sensors, 7,1612−1629.
  5. Dousset, B., and Gourmelon, F., 2003., Satellite multi-sensor data analysis of urban surface temperature and land cover. ISPRS Journal of Photogrammetry and Remote Sensing, 58,43−54.
  6. Duda, R.O., Hart, P.E., and Stork, D.G., 2000, Pattern Classification, A Wiley-Interscience Publication, Second Edition, ISBN 9814-12-602-0.
  7. Friedl, M.A., 2002, Forward and inverse modeling of land surface energy balance using surface temperature measurements. Remote Sensing of Environment, 79, 344-354.
  8. Gallo, K. P., and Tarpley, J. D., 1996., The comparison of vegetation index and surface temperature composites for urban heat island analysis. International Journal of Remote sensing, 17, 3071−3076.
  9. Gallo, K.P., and Owen, T.W., 1998, Assessment of urban heat island: A multi-sensor perspective for the Dallas-Ft. Worth, USA region. Geocarto International, 13, 35– 41.
  10. Kato, S., and Yamaguchi, Y., 2005, Analysis of urban heat-island effect using ASTER and ETM+ Data: Separation of anthropogenic heat discharge and natural heat radiation from sensible heat flux. Remote Sensing of Environment, 99, 44–54.
  11. Kustas,W. P., Norman, J.M., Anderson,M. C., and  French, A. N., 2003., Estimating subpixel surface temperatures and energy fluxes from the vegetation index-radiometric  temperature  relationship. Remote sensing of Environment, 85, 429−440.
  12. Lambin, F. F., and Ehrlich, D., 1996., The surface temperature–vegetation index space for land use and land cover change analysis. International Journal of Remote Sensing, 17, 463−487.
  13. Landsat 7 science data user’s handbook, Landsat Project Science Office, Goddard Space Flight Center, 2002, URL: http://ltwww.gsfc.nasa.gov/IAS/handbook/handbook_toc.html (last accessed: 25 March, 2008).
  14. Li, F., Jackson, T.J., Kustas, W., Schmugge, T.J., French, A.N., Cosh, M.L., Bindlish, R., 2004, Deriving land surface temperature from Landsat 5 and 7 during SMEX02/SMACEX. Remote Sensing of Environment, 92, 521 – 534.
  15. Markham, B.L., and Barker, J.K., 1985, Spectral characteristics of the LANDSAT Thematic Mapper Sensors. International Journal of Remote Sensing, 6, 697-716.
  16. Owen, T. W., Carlson, T. N., and Gillies, R. R., 1998., Assessment of satellite remotely sensed land cover parameters in quantitatively describing the climatic effect of urbanization. International Journal of Remote sensing, 19, 1663−1681.
  17. Ramachandra, T.V., and Kumar, U., Land Surface Temperature with land cover dynamics: Multi resolution, spatio – temporal data analysis of greater Bangalore, International journal of geoinfomatics, vol5, No 3. September 2009.
  18. Rao, P. K., 1972, Remote sensing of urban heat islands from an environmental satellite, Bulletin of the American Meteorological Society, 53, 647−648.
  19. Schmugge, T., Hook, S.J., and Coll, C., 1998, Recovering surface temperature and emissivity from thermal infrared multispectral data. Remote Sensing of Environment, 65, 121– 131.
  20. Snyder, W.C., Wan, Z., Zhang, Y., and Feng, Y.Z., 1998, Classification based emissivity for land surface temperature measurement from space. International Journal of Remote Sensing, 19, 2753– 2774.
  21. Sobrino J.A et al., 2004,Land surface temperature retrieval from LANDSAT TM 5. Remote Sensing of Environment, 90, 434– 440.
  22. Stathopoplou, M., Cartalis, C., and Petrakis, M., 2007, Integrating CORINE Land Cover data and landsat TM for surface emissivity definitions: an application for the urban areas of Athens, Greece. International Journal of Remote Sensing, 28, 3291-3304.
  23. Voogt, J. A., and Oke, T. R., 2003., Thermal remote sensing of urban climate. Remote Sensing of Environment, 86, 370−384.
  24. Weng, Q., Lu, D., and Schubring, J., 2004, Estimation of land surface temperature – vegetation abundances relationship for urban heat island studies. Remote Sensing of Environment, 89, 467-483.
  25. WUP (World Urbanization Prospects): The 2005 Revision, Population Division, Department of Economic and Social Affairs, UN.
  26. http://himachal.nic.in/industry/getinfon.htm
  27. Zhang, J., and Zhang, Y., 2007, Remote Sensing research issues of the National Land Use Change Program of China. Photogrammetry & Remote Sensing, 62(6), 461-472.

BACK  «  TOP  »  NEXT

E-mail     |     Sahyadri     |     ENVIS     |     GRASS     |     Energy     |     CES     |     CST     |     CiSTUP     |     IISc     |     E-mail