Urban Heat Island Analysis of Prayagraj

Remote Sensing & GIS-Based Assessment of Surface Thermal Patterns During the 2026 North Indian Heatwave

Remote Sensing · GIS · Land Surface Temperature · NDVI · Python · QGIS · 2026

Abstract

This study investigates the spatial distribution of Urban Heat Islands (UHI) across Prayagraj, Uttar Pradesh, India, using multi-sensor satellite remote sensing and geospatial analysis. Land Surface Temperature (LST) was derived from Landsat 9 Collection 2 Level-2 thermal data, and vegetation density was assessed using the Normalized Difference Vegetation Index (NDVI) computed from Sentinel-2 Level-2A imagery. Both datasets were acquired during the peak of the 2026 North Indian heatwave, one of the most intense heat events recorded in the region in over a century , providing a temporally critical and scientifically relevant observational window. The analysis employed 1 km² grid-based zonal statistics to characterize spatial variation in surface temperature and vegetation cover across the study area. Results demonstrate a statistically significant negative correlation between NDVI and LST (Pearson r = −0.386, p < 0.001), confirming that vegetation cover exerts a measurable cooling influence on surface temperatures. Mean LST across the study area reached 45.7°C, with extreme hotspots recording values above 51°C in the dense northwestern urban core, while river channels of the Ganga and Yamuna registered temperatures as low as 38–39°C. The study provides spatial evidence in support of evidence-based urban greening interventions and contributes to the growing body of literature on UHI dynamics in rapidly urbanizing South Asian cities.

Methodology

NDVI was computed from Sentinel-2 Level-2A Band 4 (Red) and Band 8 (NIR) imagery, and LST was retrieved from Landsat 9 Collection 2 Level-2 ST_B10 using USGS scaling parameters. Both datasets were reprojected to UTM Zone 44N, clipped to the Prayagraj study boundary, and cloud-masked prior to analysis. A 1 km² fishnet grid was overlaid on the study area and zonal statistics (mean NDVI and mean LST) were extracted per cell. Pearson correlation and ordinary least-squares regression were applied across approximately 280 valid grid cells to quantify the NDVI–LST relationship.

NDVI Map

Classified NDVI map of Prayagraj derived from Sentinel-2 Level-2A imagery, peak summer 2026. Blue: water bodies; orange: bare soil/sand/built-up; yellow-green: sparse vegetation; medium green: moderate vegetation; dark green: dense vegetation. The Prayagraj Cantonment is identifiable as the concentrated dark green cluster in the west-central portion of the study area. Ganga (north) and Yamuna (south/west) rivers are clearly delineated in blue.

Classified NDVI Map of Prayagraj

Land Surface Temperature Map

Land Surface Temperature map of Prayagraj derived from Landsat 9 Collection 2 Level-2 ST_B10, peak summer 2026. Colour scale from blue (cool, ~34°C) through orange to deep red (hot, ~53°C). River channels of the Ganga and Yamuna are clearly visible as blue cooling corridors. The northwest and northeast bare riverbanks show persistent deep red temperatures. The Cantonment zone (centre-left) is distinguishable as a lighter thermal signature relative to adjacent dense urban areas.

Land Surface Temperature Map of Prayagraj

LST Grid Map

Mean Land Surface Temperature per 1 km² grid cell across Prayagraj, peak summer 2026. Each cell is labelled with its mean LST value in degrees Celsius. The northwest hotspot cluster (50–51°C cells) and the river cooling corridors (38–40°C cells) are directly readable from the labelled grid. The Cantonment zone is identifiable as the cooler cell cluster (41–43°C) in the west-central area.

Mean LST per 1 km² Grid Cell, Prayagraj

LST Frequency Distribution

Frequency distribution of Land Surface Temperature pixel values across the Prayagraj study area, peak summer 2026. Mean (45.7°C) and median (45.3°C) are annotated. The pronounced high-temperature shoulder at 49–50°C reflects the thermal signature of the bare riverbanks.

LST Frequency Distribution Histogram

NDVI Frequency Distribution

Frequency distribution of NDVI pixel values across the Prayagraj study area, peak summer 2026. The mean NDVI of approximately 0.15 falls at the boundary between sparse and moderate vegetation classes, confirming the overwhelmingly low-vegetation character of the urban surface during peak summer.

NDVI Frequency Distribution Histogram

NDVI vs LST Scatter Plot

Scatter plot of mean NDVI versus mean Land Surface Temperature per 1 km² grid cell across Prayagraj, peak summer 2026. Each point represents one grid cell. The linear regression line confirms a negative relationship between vegetation density and surface temperature. Pearson r = −0.386, p < 0.001 (n ≈ 280 grid cells).

NDVI vs LST Scatter Plot with Regression Line

Key Findings

  1. Mean LST of 45.7°C was recorded across the study area, with extreme hotspots exceeding 51°C in the northwestern and northeastern bare riverbanks of the Ganga river, and minimum values of 34–35°C in river channels, an overall thermal range of approximately 18–19°C driven by land cover variation.
  2. A statistically significant negative correlation (Pearson r = −0.386, p < 0.001) was confirmed between mean NDVI and mean LST across 1 km² grid cells. Each 0.10 increase in mean NDVI is associated with approximately 1.5–2.0°C reduction in mean LST.
  3. River cooling is real but spatially confined, the Ganga and Yamuna channels provide 10–13°C of cooling relative to adjacent built-up surfaces, but this effect is restricted to the water channel itself. The sandy reh land flanking both rivers records LST values of 46–49°C, comparable to the densest urban surfaces.
  4. The Prayagraj Cantonment demonstrates a persistent 6–9°C thermal advantage over the dense urban core, a direct, measurable consequence of its historically planned, vegetated, low-density urban fabric, providing empirical evidence that deliberate green urban planning delivers durable climate resilience over decadal timescales.

Tools & Data

Dataset / Software Provider / Version Role
Sentinel-2 Level-2A ESA Copernicus Data Space NDVI computation (B04, B08 at 10 m)
Landsat 9 Collection 2 Level-2 USGS EarthExplorer Land Surface Temperature retrieval (ST_B10 at 30 m)
QGIS 3.34 LTS Pre-processing, spatial analysis, map production
Python · rasterio · numpy · scipy · matplotlib · geopandas 3.10+ LST conversion, NDVI calculation, statistical analysis, chart generation

Skills Demonstrated

  • Satellite remote sensing: Sentinel-2 and Landsat 9 data acquisition and processing
  • Land Surface Temperature retrieval using USGS Collection 2 scaling methodology
  • NDVI computation, classification, and interpretation
  • Raster pre-processing: reprojection, clipping, cloud masking
  • Grid-based zonal statistics and spatial aggregation
  • Pearson correlation and ordinary least-squares regression analysis
  • Python scripting for geospatial raster processing and scientific visualization
  • Cartographic design and thematic map production in QGIS
  • Urban climate analysis and scientific interpretation