A new multiscale approach for monitoring vegetation using remote sensing-based indicators in laboratory, field, and landscape

Remote sensing is an important tool for studying patterns in surface processes on different spatiotemporal scales. However, differences in the spatiospectral and temporal resolution of remote sensing data as well as sensor-specific surveying characteristics very often hinder comparative analyses and effective up- and downscaling analyses. This paper presents a new methodical framework for combining hyperspectral remote sensing data on different spatial and temporal scales. We demonstrate the potential of using the “One Sensor at Different Scales” (OSADIS) approach for the laboratory (plot), field (local), and landscape (regional) scales. By implementing the OSADIS approach, we are able (1) to develop suitable stress-controlled vegetation indices for selected variables such as the Leaf Area Index (LAI), chlorophyll, photosynthesis, water content, nutrient content, etc. over a whole vegetation period. Focused laboratory monitoring can help to document additive and counteractive factors and processes of the vegetation and to correctly interpret their spectral response; (2) to transfer the models obtained to the landscape level; (3) to record imaging hyperspectral information on different spatial scales, achieving a true comparison of the structure and process results; (4) to minimize existing errors from geometrical, spectral, and temporal effects due to sensor- and time-specific differences; and (5) to carry out a realistic top- and downscaling by determining scale-dependent correction factors and transfer functions. The first results of OSADIS experiments are provided by controlled whole vegetation experiments on barley under water stress on the plot scale to model LAI using the vegetation indices Normalized Difference Vegetation Index (NDVI) and green NDVI (GNDVI). The regression model ascertained from imaging hyperspectral AISA-EAGLE/HAWK (DUAL) data was used to model LAI. This was done by using the vegetation index GNDVI with an R ² of 0.83, which was transferred to airborne hyperspectral data on the local and regional scales. For this purpose, hyperspectral imagery was collected at three altitudes over a land cover gradient of 25 km within a timeframe of a few minutes, yielding a spatial resolution from 1 to 3 m. For all recorded spatial scales, both the LAI and the NDVI were determined. The spatial properties of LAI and NDVI of all recorded hyperspectral images were compared using semivariance metrics derived from the variogram. The first results show spatial differences in the heterogeneity of LAI and NDVI from 1 to 3 m with the recorded hyperspectral data. That means that differently recorded data on different scales might not sufficiently maintain the spatial properties of high spatial resolution hyperspectral images.     

Change vector analysis to categorise land cover change processes using the tasselled cap as biophysical indicator

The continuous extraction of wood and the conversion of forest to small- and large-scale agricultural parcels is rapidly changing the land cover of the mount Cameroon region. The changes occur at varying spatial scales most often not more than 2ha for the small-scale subsistence farms and above 10ha for the extensive agricultural plantations of cocoa and palm. Given the importance of land use and land cover data in conservation planning, accurate and efficient techniques to provide up-to-date change information are required. A number of techniques for realising the detection of land cover dynamics using remotely sensed imagery have been formulated, tested and assessed with the results varying with respect to the change scenario under investigation, the information required and the imagery applied. In this study the Change Vector Analysis (CVA) technique was implemented on multitemporal multispectral Landsat data from the Thematic Mapper (TM) and Enhanced Thematic Mapper (ETM) sensors to monitor the dynamics of forest change in the mount Cameroon region. CVA was applied to multi-temporal data to compare the differences in the time-trajectory of the tasseled cap greenness and brightness for two successive time periods - 1987 and 2002. The tasseled cap was selected as biophysical indicator because it optimises the data viewing capabilities of vegetation, representing the basic types of land cover - vegetation, soil and water. Classes were created arbitrarily to predict the technique's potential in monitoring forest cover changes in the mount Cameroon region. The efficiency of the technique could not be fully assessed due to the inavailability of sufficient ground truth data. Assessment was based on the establishment of an error matrix of change versus no-change. The overall accuracy was 70%. The technique nevertheless demonstrated immense potentials in monitoring forest cover change dynamics especially when complemented with field studies.     

Vegetation Mapping for Change Detection on an Arid-Zone River

A vegetation mapping system for change detection was tested at the Havasu National Wildlife Refuge (HNWR) on the Lower Colorado River. A low-cost, aerial photomosaic of the 4200 ha, study area was constructed utilizing an automated digital camera system, supplemented with oblique photographs to aid in determining species composition and plant heights. Ground-truth plots showed high accuracy in distinguishing native cottonwood (Populus fremontii) and willow (Salix gooddingii) trees from other vegetation on aerial photos. Marsh vegetation (mainly cattails, Typha domengensis) was also easily identified. However, shrubby terrestrial vegetation, consisting of saltcedar (Tamarix ramosissima), arrowweed (Pluchea sericea), and mesquite trees (Prosopis spp.), could not be accurately distinguished from each other and were combined into a single shrub layer on the final vegetation map. The final map took the form of a base, shrub and marsh layer, which was displayed as a Normalized Difference Vegetation Index map from a Landsat Enhanced Thematic Mapper (ETM+) image to show vegetation intensity. Native willow and cottonwood trees were digitized manually on the photomosaic and overlain on the shrub layer in a GIS. By contrast to present, qualitative mapping systems used on the Lower Colorado River, this mapping system provides quantitative information that can be used for accurate change detection. However, better methods to distinguish between saltcedar, mesquite, and arrowweed are needed to map the shrub layer.     

Analysis of agricultural drought using vegetation temperature condition index (VTCI) from Terra/MODIS satellite data

The most commonly used normalized difference vegetation index (NDVI) from remote sensing often fall short in real-time drought monitoring due to a lagged vegetation response to drought. Therefore, research recently emphasized on the use of combination of surface temperature and NDVI which provides vegetation and moisture conditions simultaneously. Since drought stress effects on agriculture are closely linked to actual evapotranspiration, we used a vegetation temperature condition index (VTCI) which is more closely related to crop water status and holds a key place in real-time drought monitoring and assessment. In this study, NDVI and land surface temperature (T _s) from MODIS 8-day composite data during cloud-free period (September–October) were adopted to construct an NDVI–T _s space, from which the VTCI was computed. The crop moisture index (based on estimates of potential evapotranspiration and soil moisture depletion) was calculated to represent soil moisture stress on weekly basis for 20 weather monitoring stations. Correlation and regression analysis were attempted to relate VTCI with crop moisture status and crop performance. VTCI was found to accurately access the degree and spatial extent of drought stress in all years (2000, 2002, and 2004). The temporal variation of VTCI also provides drought pattern changes over space and time. Results showed significant and positive relations between CMI (crop moisture index) and VTCI observed particularly during prominent drought periods which proved VTCI as an ideal index to monitor terminal drought at regional scale. VTCI had significant positive relationship with yield but weakly related to crop anomalies. Duration of terminal drought stress derived from VTCI has a significant negative relationship with yields of major grain and oilseeds crops, particularly, groundnut.     

Water quality monitoring and assessment of an urban Mediterranean lake facilitated by remote sensing applications

Degradation of water quality is a major problem worldwide and often leads to serious environmental impacts and concerns about public health. In this study, the water quality monitoring and assessment of the Koumoundourou Lake, a brackish urban shallow lake located in the northeastern part of Elefsis Bay (Greece), were evaluated. A number of water quality parameters (pH, temperature, dissolved oxygen concentration, electrical conductivity, turbidity, nutrients, and chlorophyll-a concentration) were analyzed in water samples collected bimonthly over a 1-year period from five stations throughout the lake. Moreover, biological quality elements were analysed seasonally over the 1-year period (benthic fauna). Statistical analysis was performed in order to evaluate the water quality of the lake and distinguish sources of variation measured in the samples. Furthermore, the chemical and trophic status of the lake was evaluated according to the most widely applicable classification schemes. Satellite images of Landsat 5 Thematic Mapper were used in order for algorithms to be developed and calculate the concentration of chlorophyll-a (Chl-a). The trophic status of the lake was characterized as oligotrophic based on phosphorus and as mesotrophic–eutrophic based on Chl-a concentrations. The results of the remote sensing application indicated a relatively high coefficient of determination (R ²) among point sampling results and the remotely sensed data, which implies that the selected algorithm is reliable and could be used for the monitoring of Chl-a concentration in the particular water body when no field data are available.     

Quantitative analysis of shoreline changes at the Mediterranean Coast in Turkey

This research is focused on the coastline evolution monitoring and its potential change estimation by remote sensing techniques using multi-temporal Landsat images at the southeast coasts of the Mediterranean Sea in Turkey. The study area includes the coastal zone located in the Cukurova Delta coasts. The Cukurova Delta has accreted toward the Mediterranean Sea as a result of sediment discharge and transport from Seyhan and Ceyhan rivers. These processes have caused the morphological changes (accretion or erosion) of coastline along some parts of the southeast coasts of the Mediterranean Sea. In this study, coastline changes were researched by using radiometrically and geometrically corrected multi-temporal and multi-spectral data from Landsat Multispectral Scanner dated 1972, Thematic Mapper dated 1987, and Enhanced Thematic Mapper dated 2002. In the image processing steps, mosaicing, subset, Iterative Self-Organizing Data Analysis Technique classification, band ratioing (B5/B2), edge detection, and overlay techniques were used to carry out coastline extraction and the Digital Shoreline Analysis System was used to calculate rate of coastline changes. As a result of the analysis, in some parts of the research area, remarkable coastline changes (more than 2,900 m withdrawal and ??24.50 m/year erosion) were observed for a 30–year period.     

Evaluating the difference between the normalized difference vegetation index and net primary productivity as the indicators of vegetation vigor assessment at landscape scale

Both the net primary productivity (NPP) and the normalized difference vegetation index (NDVI) are commonly used as indicators to characterize vegetation vigor, and NDVI has been used as a surrogate estimator of NPP in some cases. To evaluate the reliability of such surrogation, here we examined the quantitative difference between NPP and NDVI in their outcomes of vegetation vigor assessment at a landscape scale. Using Landsat ETM+ data and a process model, the Boreal Ecosystem Productivity Simulator, NPP distribution was mapped at a resolution of 90 m, and total NDVI during the growing season was calculated in Heihe River Basin, Northwest China in 2002. The results from a comparison between the NPP and NDVI classification maps show that there existed a substantial difference in terms of both area and spatial distribution between the assessment outcomes of these two indicators, despite that they are strongly correlated. The degree of difference can be influenced by assessment schemes, as well as the type of vegetation and ecozone. Overall, NDVI is not a good surrogate of NPP as the indicators of vegetation vigor assessment in the study area. Nonetheless, NDVI could serve as a fairish surrogate indicator under the condition that the target region has low vegetation cover and the assessment has relatively coarse classification schemes (i.e., the class number is small). It is suggested that the use of NPP and NDVI should be carefully selected in landscape assessment. Their differences need to be further evaluated across geographic areas and biomes.     

Seagrass resource assessment using remote sensing methods in St. Joseph Sound and Clearwater Harbor,Florida, USA

In the event of a natural or anthropogenic disturbance, environmental resource managers require a reliable tool to quickly assess the spatial extent of potential damage to the seagrass resource. The temporal availability of the Landsat 5 Thematic Mapper (TM) imagery provided a suitable option to detect and assess damage of the submerged aquatic vegetation (SAV). This study examined Landsat TM imagery classification techniques to create two-class (SAV presence/absence) and three-class (SAV estimated coverage) SAV maps of the seagrass resource. The Mahalanobis Distance method achieved the highest overall accuracy (86%) and validation accuracy (68%) for delineating the seagrass resource (two-class SAV map). The Maximum Likelihood method achieved the highest overall accuracy (74%) and validation accuracy (70%) for delineating the seagrass resource three-class SAV map. The Landsat 5 TM imagery classification provided a seagrass resource map product with similar accuracy to the aerial photointerpretation maps (validation accuracy 71%). The results support the application of remote sensing methods to analyze the spatial extent of the seagrass resource.     

Vegetation NDVI Linked to Temperature and Precipitation in the Upper Catchments of Yellow River

Vegetation in the upper catchment of Yellow River is critical for the ecological stability of the whole watershed. The dominant vegetation cover types in this region are grassland and forest, which can strongly influence the eco-environmental status of the whole watershed. The normalized difference vegetation index (NDVI) for grassland and forest has been calculated and its daily correlation models were deduced by Moderate Resolution Imaging Spectroradiometer products on 12 dates in 2000, 2003, and 2006. The responses of the NDVI values with the inter-annual grassland and forest to three climatic indices (i.e., yearly precipitation and highest and lowest temperature) were analyzed showing that, except for the lowest temperature, the yearly precipitation and highest temperature had close correlations with the NDVI values of the two vegetation communities. The value of correlation coefficients ranged from 0.815 to 0.951 (p?<?0.01). Furthermore, the interactions of NDVI values of vegetation with the climatic indicators at monthly interval were analyzed. The NDVI of vegetation and three climatic indices had strong positive correlations (larger than 0.733, p?<?0.01). The monthly correlations also provided the threshold values for the three climatic indictors, to be used for simulating vegetation growth grassland under different climate features, which is essential for the assessment of the vegetation growth and for regional environmental management.     

A change vector analysis technique for monitoring land cover changes in Copsa Mica,Romania, in the period 1985–2011

During the communist regime, Romania’s planned economy focused exclusively on production neglecting the environment protection. The lack of less polluting production technologies and of environmental protection measures led to excessive pollution in certain industrialized areas. This is the case of the town of Copsa Mica in Sibiu County, which in 1987 was considered one of the most polluted towns in Europe. The present study assesses the change vector analysis (CVA) technique using a Landsat Thematic Mapper (TM) image time series to monitor land cover changes caused by carbon black and heavy metal pollution. CVA was applied to the tasseled cap greenness (TCG) and tasseled cap brightness (TCB) indices, as well as to the Normalized Difference Vegetation Index (NDVI) and bare soil index (BI). Various maps were generated for the periods 1985–1994, 1994–2003, 2003–2011, and 1985–2011, and threshold values were determined for the detection of land cover change/no change. The change direction and magnitude values were cross-tabulated and classified. The technique was assessed based on the change versus no-change error matrix. The results show that in the area of Copsa Mica, land cover changes occurred because of a considerable decrease in the area affected by carbon black and heavy metal pollution. The CVA technique proved efficient in monitoring the land cover changes caused by pollution and especially by carbon black pollution. Soil pollution by heavy metals is reflected in the bare soil surfaces present in the imagery.     

Pattern of NDVI-based vegetation greening along an altitudinal gradient in the eastern Himalayas and its response to global warming

The eastern Himalayas, especially the Yarlung Zangbo Grand Canyon Nature Reserve (YNR), is a global hotspot of biodiversity because of a wide variety of climatic conditions and elevations ranging from 500 to > 7000 m above sea level (a.s.l.). The mountain ecosystems at different elevations are vulnerable to climate change; however, there has been little research into the patterns of vegetation greening and their response to global warming. The objective of this paper is to examine the pattern of vegetation greening in different altitudinal zones in the YNR and its relationship with vegetation types and climatic factors. Specifically, the inter-annual change of the normalized difference vegetation index (NDVI) and its variation along altitudinal gradient between 1999 and 2013 was investigated using SPOT-VGT NDVI data and ASTER global digital elevation model (GDEM) data. We found that annual NDVI increased by 17.58 % in the YNR from 1999 to 2013, especially in regions dominated by broad-leaved and coniferous forests at lower elevations. The vegetation greening rate decreased significantly as elevation increased, with a threshold elevation of approximately 3000 m. Rising temperature played a dominant role in driving the increase in NDVI, while precipitation has no statistical relationship with changes in NDVI in this region. This study provides useful information to develop an integrated management and conservation plan for climate change adaptation and promote biodiversity conservation in the YNR.     

Weed mapping in cotton using ground-based sensors and GIS

Site-specific weed management presupposes the careful monitoring and mapping of weed infestation areas. Cut-edge sensor technologies coupled with geographical information systems (GIS) provide the means for reliable decision-making concerning weed management even in sub-field level. In present research, two different spectral sensing systems were engaged in order to digitally map weed patches as grown in four different cotton fields in Central Greece. The systems used were a set of two Crop Circle multispectral sensors ACS-430 and a digital camera Nikon D300S. The spaces between cotton rows were scanned and photographed with the two systems accordingly. Raw recorded data were stored and analyzed in GIS environment producing spatially interpolated maps of red-edge normalized difference vegetation index (NDVI) and weed cover percentage values. Both mapping approaches were satisfactorily related to weed distribution as occurred in the fields; however, the photographic method tended to underestimate weed populations. Correlation of red-edge NDVI and weed cover values, at the points where photographs were taken, as revealed by Pearson’s correlation coefficient was high (r?>?0.83) and statistically significant at the 0.01 level. A first-degree linear equation adequately modeled (R²?>?0.7) the between value pair relations, strengthening the validity of the two methodologies in spatially monitoring weed patches. The methodologies and the technologies used in the study can be used for yearly mapping weed flora in cotton cultivation and potentially constitute a means of rationalizing herbicide application in terms of doses and spatio-temporal decision-making.     

Vegetation response to hydrologic and geomorphic factors in an arid region of the Baja California Peninsula

Remotely sensed imageries were used to analyze the response of desert vegetation to physiographic factors and accumulated precipitation in drier and wetter years within a region of >16,500 km² sampled with 5,000 random pixels of 30 m. Vegetation development was indexed by the annual maximum values for greenness (SAVI) and canopy water content (NDII). Precipitation was interpolated from the 0.25° grid of the Tropical Rainfall Measurement Mission satellite-based estimates, showing a regional average of ～55 mm in the wetter year. The vegetation indices were only weakly related to total precipitation, often in a negative sense. Terrain factors that most often affected the vegetation indices, in multiple regression models, were Topographic Wetness Index, elevation, and slope gradient; these often had different signs for SAVI and for NDII. Models for NDII on intrusive igneous rocks gave better results than on extrusive igneous rocks. The strongest patterns in vegetation development were the contrast among Pacific coast, Cordillera, and Gulf coast subregions and the generally stronger results for NDII than SAVI.     

Estimation of Vegetation Coverage in Semi-arid Sandy Land Based on Multivariate Statistical Modeling Using Remote Sensing Data

The estimation of vegetation coverage is essential in the monitoring and management of arid and semi-arid sandy lands. But how to estimate vegetation coverage and monitor the environmental change at global and regional scales still remains to be further studied. Here, combined with field vegetation survey, multispectral remote sensing data were used to estimate coverage based on theoretical statistical modeling. First, the remote sensing data were processed and several groups of spectral variables were selected/proposed and calculated, and then statistically correlated to measured vegetation coverage. Both the single- and multiple-variable-based models were established and further analyzed. Among all single-variable-based models, that is based on Normalized Difference Vegetation Index showed the highest R (0.900) and R ² (0.810) as well as lowest standard estimate error (0.128024). Since the multiple-variable-based model using multiple stepwise regression analysis behaved much better, it was determined as the optimal model for local coverage estimation. Finally, the estimation was conducted based on the optimal model and the result was cross-validated. The coefficient of determination used for validation was 0.867 with a root-mean-squared error (RMSE) of 0.101. The large-scale estimation of vegetation coverage using statistical modeling based on remote sensing data can be helpful for the monitoring and controlling of desertification in arid and semi-arid regions. It could serve for regional ecological management which is of great significance.     

Determining the impact of urban components on land surface temperature of Istanbul by using remote sensing indices

For the past 60 years, Istanbul has been experiencing an accelerated urban expansion. This urban expansion is leading to the replacement of natural surfaces by various artificial materials. This situation has a critical impact on the environment due to the alteration of heat energy balance. In this study, the effect upon the urban heat island (UHI) of Istanbul was analyzed using 2009 dated Landsat 5 Thematic Mapper (TM) data. An Index Based Built-up Index (IBI) was used to derive artificial surfaces in the study area. To produce the IBI index, Soil-Adjusted Vegetation Index, Normalized Difference Built-up Index, and Modified Normalized Difference Water Index were calculated. Land surface temperature (LST) distribution was derived from Landsat 5 TM images using a mono-window algorithm. In addition, 24 transects were selected, and different regression models were applied to explore the correlation between LST and IBI index. The results show that artificial surfaces have a positive exponential relationship with LST rather than a simple linear one. An ecological evaluation index of the region was calculated to explore the impact of both the vegetated land and the artificial surfaces on the UHI. Therefore, the quantitative relationship of urban components (artificial surfaces, vegetation, and water) and LST was examined using multivariate statistical analysis, and the correlation coefficient was obtained as 0.829. This suggested that the areas with a high rate of urbanization will accelerate the rise of LST and UHI in Istanbul.     

Growing stock and woody biomass assessment in Asola-Bhatti Wildlife Sanctuary,Delhi, India

Biomass is an important entity to understand the capacity of an ecosystem to sequester and accumulate carbon over time. The present study, done in collaboration with the Delhi Forest Department, focused on the estimation of growing stock and the woody biomass in the so-called lungs of Delhi—the Asola-Bhatti Wildlife Sanctuary in northern Aravalli hills. The satellite-derived vegetation strata were field-inventoried using stratified random sampling procedure. Growing stock was calculated for the individual sample plots using field data and species-specific volume equations. Biomass was estimated from the growing stock and the specific gravity of the wood. Among the four vegetation types, viz. Prosopis juliflora, Anogeissus pendula, forest plantation and the scrub, the P. juliflora was found to be the dominant vegetation in the area, covering 23.43 km² of the total area. The study revealed that P. juliflora forest with moderate density had the highest (10.7 m³/ha) while A. pendula forest with moderate density had the lowest (3.6 m³/ha) mean volume. The mean woody biomass was also found to be maximum in P. juliflora forest with moderate density (10.3 t/ha) and lowest in A. pendula forest with moderate density (3.48 t/ha). The total growing stock was estimated to be 20,772.95 m³ while total biomass worked out to be 19,366.83 t. A strong correlation was noticed between the normalized difference vegetation index (NDVI) and the growing stock (R ²?=?0.84)/biomass (R ²?=?0.88). The study demonstrated that growing stock and the biomass of the woody vegetation in Asola-Bhatti Wildlife Sanctuary could be estimated with high accuracy using optical remote sensing data.     

Monitoring and identification of spatiotemporal landscape changes in multiple remote sensing images by using a stratified conditional Latin hypercube sampling approach and geostatistical simulation

This study develops a stratified conditional Latin hypercube sampling (scLHS) approach for multiple, remotely sensed, normalized difference vegetation index (NDVI) images. The objective is to sample, monitor, and delineate spatiotemporal landscape changes, including spatial heterogeneity and variability, in a given area. The scLHS approach, which is based on the variance quadtree technique (VQT) and the conditional Latin hypercube sampling (cLHS) method, selects samples in order to delineate landscape changes from multiple NDVI images. The images are then mapped for calibration and validation by using sequential Gaussian simulation (SGS) with the scLHS selected samples. Spatial statistical results indicate that in terms of their statistical distribution, spatial distribution, and spatial variation, the statistics and variograms of the scLHS samples resemble those of multiple NDVI images more closely than those of cLHS and VQT samples. Moreover, the accuracy of simulated NDVI images based on SGS with scLHS samples is significantly better than that of simulated NDVI images based on SGS with cLHS samples and VQT samples, respectively. However, the proposed approach efficiently monitors the spatial characteristics of landscape changes, including the statistics, spatial variability, and heterogeneity of NDVI images. In addition, SGS with the scLHS samples effectively reproduces spatial patterns and landscape changes in multiple NDVI images.     

Estimating and mapping forest biomass using regression models and Spot-6 images (case study: Hyrcanian forests of north of Iran)

Hyrcanian forests of North of Iran are of great importance in terms of various economic and environmental aspects. In this study, Spot-6 satellite images and regression models were applied to estimate above-ground biomass in these forests. This research was carried out in six compartments in three climatic (semi-arid to humid) types and two altitude classes. In the first step, ground sampling methods at the compartment level were used to estimate aboveground biomass (Mg/ha). Then, by reviewing the results of other studies, the most appropriate vegetation indices were selected. In this study, three indices of NDVI, RVI, and TVI were calculated. We investigated the relationship between the vegetation indices and aboveground biomass measured at sample-plot level. Based on the results, the relationship between aboveground biomass values and vegetation indices was a linear regression with the highest level of significance for NDVI in all compartments. Since at the compartment level the correlation coefficient between NDVI and aboveground biomass was the highest, NDVI was used for mapping aboveground biomass. According to the results of this study, biomass values were highly different in various climatic and altitudinal classes with the highest biomass value observed in humid climate and high-altitude class.     

Spatio-temporal footprints of urbanisation in Surat, the Diamond City of India (1990–2009)

Urbanisation is a ubiquitous phenomenon with greater prominence in developing nations. Urban expansion involves land conversions from vegetated moisture-rich to impervious moisture-deficient land surfaces. The urban land transformations alter biophysical parameters in a mode that promotes development of heat islands and degrades environmental health. This study elaborates relationships among various environmental variables using remote sensing dataset to study spatio-temporal footprint of urbanisation in Surat city. Landsat Thematic Mapper satellite data were used in conjugation with geo-spatial techniques to study urbanisation and correlation among various satellite-derived biophysical parameters, [Normalised Difference Vegetation Index, Normalised Difference Built-up Index, Normalised Difference Water Index, Normalised Difference Bareness Index, Modified NDWI and land surface temperature (LST)]. Land use land cover was prepared using hierarchical decision tree classification with an accuracy of 90.4 % (kappa?=?0.88) for 1990 and 85 % (kappa?=?0.81) for 2009. It was found that the city has expanded over 42.75 km² within a decade, and these changes resulted in elevated surface temperatures. For example, transformation from vegetation to built-up has resulted in 5.5?±?2.6 °C increase in land surface temperature, vegetation to fallow 6.7?±?3 °C, fallow to built-up is 3.5?±?2.9 °C and built-up to dense built-up is 5.3?±?2.8 °C. Directional profiling for LST was done to study spatial patterns of LST in and around Surat city. Emergence of two new LST peaks for 2009 was observed in N–S and NE–SW profiles.