Anaplian.com · urban remote sensing · synthesis

City ranking of the IVC (surface climate-vulnerability index)

April–May 2026

v 0.1 · 15 September 2026

Mexico City · heat island
Mexico City · heat island

What this ranking compares

This report compares surface heat in 18 cities. It does not name the hottest place on Earth. It does not predict heat deaths. The cities have more than two million inhabitants and sit between 18° and 42° north. All of them are measured in April and May 2026.

The IVC combines four measures, each with equal weight:

  1. LST (land surface temperature), peak: urban fraction whose per-pixel maximum exceeds 43 °C.
  2. LST mean: urban fraction whose composite mean exceeds 43 °C.
  3. NDVI (normalized difference vegetation index), inverted: less greenness raises the index.
  4. SUHI (surface urban heat island): urban mean LST minus that city’s rural temperature. A negative SUHI adds nothing.

The urban mask and the rural reference come from GHSL (Global Human Settlement Layer) and its SMOD (settlement model).

We chose them for hazard, exposure and vulnerability — not for temperature records. That is why an index of 60 in Ahmedabad can be read next to a 60 in Iztapalapa.

One list, one window

City reports already exist one by one. Here they share a single list, using each city’s official polygon, so the rank can be read next to the district or borough reports.

Headline result. Ahmedabad leads (61.1/100). Rome is last (4.1/100). City median 37.2. 6 cities have urban SUHI below 0 °C.

City ranking

IVC 0–100. Equal weights. Not renormalised. Source: Landsat 8/9, April–May 2026.
Each bar is 25 × z. Dark red: peak LST. Gold: mean LST. Green: greenness deficit. Blue: SUHI.
Contribution of peak LST to the IVC. Each layer can add at most 25 points (equal weights).
Contribution of mean LST to the IVC. Ahmedabad stays high on both layers; Phoenix and Mexico City diverge.
Urban mean LST − rural T. Negative values mean the GHSL rural reference is hotter than the city.
Rural GHSL mean. High values (Ahmedabad, Delhi, Karachi, Riyadh) flatten the SUHI component.
#CityPIVCPeak > 43 °CMean > 43 °CSUHI °CT ruralNDVIUnits
1Ahmedabad261.199%89%-1.748.10.277
2Phoenix255.499%39%+1.340.30.1715
3Mexico City154.164%13%+13.025.60.2816
4Riyadh350.9100%9%-0.039.90.0814
5Karachi150.191%35%-3.244.10.1824
6Baghdad146.095%1%+0.837.60.144
7Delhi144.989%42%-2.545.10.3111
8Casablanca342.689%0%+2.627.10.2316
9Cairo139.569%3%-0.639.00.1241
10Tehran235.047%0%+2.630.20.1722
11Lahore234.979%2%+1.736.60.312
12Barcelona331.152%0%+3.026.30.2910
13Monterrey228.766%0%-0.527.80.3113
14Dhaka224.433%0%+4.726.80.3841
15Madrid116.922%0%+0.322.80.3321
16Houston314.114%0%+2.923.80.4411
17Wuhan311.86%0%+2.924.40.4413
18Rome34.12%0%+0.926.40.5115

Patterns

Ahmedabad leads (61.1/100) because almost all of its urban fabric exceeds 43 °C at both the pixel maximum (~99%) and the composite mean (~89%). The GHSL rural reference is itself a furnace (48 °C), so SUHI is negative (-1.7 °C) and adds zero to the index. Ahmedabad ranks high on three layers, not on a classic city-versus-countryside island. Inside the city the ranking is almost flat (about 58–63): every taluka (subdistrict) is already a furnace. The 2013 Heat Action Plan is a public-health pioneer; the IVC measures remaining surface exposure, not lives the plan saves.

Mexico City is not first on the IVC (54.1/100, third), but it has the strongest heat island in the sample (SUHI +13.0 °C). Southern conservation land sits near 26 °C while the paved east exceeds 43 °C. That is the harshest urban heat wave in this set: not the hottest desert, but the city that warms the most relative to its own surroundings. 6 cities have SUHI below 0 °C; that component is clipped at zero.

Rome (4.1) and Madrid (16.9) sit at the bottom because April–May is spring, not their annual peak. Houston and Wuhan enter on the same fixed window. The scale is kept so that a 60 in a Houston district remains comparable to Iztapalapa.

Casablanca sits between Cairo and Barcelona: same latitude as Phoenix, Atlantic softening, less adaptive capacity than Europe. Lahore is close to Tehran at city scale; two tehsils (subdistricts) hide the intra-city range that 41 qisms (district sections) in Cairo or 16 alcaldías (boroughs) in Mexico City show. Grain is not population-weighted.

Ahmedabad — first on the index

Ahmedabad leads at 61.1/100. About 99% of the urban area exceeds 43 °C at the pixel peak and 89% at the two-month mean. Peak LST therefore contributes 24.9 points to the index and mean LST 22.2 points (each layer caps at 25). Rural GHSL land is 48 °C, so SUHI is -1.7 °C and adds zero. The ranking inside the city is almost flat (58–63 across 7 talukas / subdistricts). Sabarmati is highest; Dholka is lowest. The 2013 Heat Action Plan is a public-health pioneer; the IVC measures remaining surface exposure, not lives the plan saves.

Full city report

LST max
Peak LST, overlay 43 °C.
LST mean
Mean LST of the observations (map colours from 35 °C).
NDVI
Mean NDVI, same window.
SUHI
Surface heat island (pixel minus rural T).
IVC units
IVC of official subdivisions in the city.

Phoenix — second

Phoenix is second (55.4/100). Peak LST is almost saturated (99% above 43 °C, 24.7 points), but the mean layer is much weaker (39%, 9.6 points): April–May is not Phoenix’s July peak. SUHI is only +1.3 °C against a 40 °C rural desert. Greenness is low (NDVI 0.17). Maryvale leads the urban villages at 71.2; Desert View is last at 43.6.

Full city report

LST max
Peak LST, overlay 43 °C.
LST mean
Mean LST of the observations (map colours from 35 °C).
NDVI
Mean NDVI, same window.
SUHI
Surface heat island (pixel minus rural T).
IVC units
IVC of official subdivisions in the city.

Mexico City — third, strongest heat island

Mexico City is third (54.1/100) but has the sample’s strongest SUHI (+13.0 °C). Peak LST covers 64% of the urban area (15.9 points); mean LST only 13% (3.3 points). The heat-island layer adds 21.7 points — almost the whole 25-point cap — because southern conservation land is near 26 °C while the paved east exceeds 43 °C. Intra-city range is the widest in the sample (11.3–82.6). Venustiano Carranza leads; Cuajimalpa de Morelos is last.

Full city report

LST max
Peak LST, overlay 43 °C.
LST mean
Mean LST of the observations (map colours from 35 °C).
NDVI
Mean NDVI, same window.
SUHI
Surface heat island (pixel minus rural T).
IVC units
IVC of official subdivisions in the city.

Rome — lowest score

Rome is last (4.1/100). April–May is spring, not the Mediterranean summer peak. Only 2% of the urban area exceeds 43 °C at the pixel maximum and 0% at the mean, so the two LST layers add 0.4 and 0.0 points. NDVI is high (0.51) and SUHI is modest (+0.9 °C). Municipio V is the warmest municipio (municipality) at 12.9; Municipio XIV is 0.0. The low city score is the fixed window working as designed: a 60 here would still be comparable to Iztapalapa.

Full city report

LST max
Peak LST, overlay 43 °C.
LST mean
Mean LST of the observations (map colours from 35 °C).
NDVI
Mean NDVI, same window.
SUHI
Surface heat island (pixel minus rural T).
IVC units
IVC of official subdivisions in the city.

Karachi — what a negative SUHI means

Karachi ranks fifth (50.1/100) with SUHI -3.2 °C. That does not mean the city is cool. Rural GHSL land around Karachi averages 44 °C — hotter than the urban mean — so the “island” is inverted: the city is a relative oasis next to an arid belt, while 91% of the urban fabric still exceeds 43 °C at peak (22.7 points) and 35% at the mean (8.8 points). The SUHI component is clipped at zero, so a negative value neither raises nor lowers the IVC: Karachi is not rewarded for being cooler than a furnace, and it is not treated as an island city. NDVI is low ({k['ndvi']:.2f}), so greenness deficit still adds 18.5 points. New Karachi Town leads at 59.3; Manora Cantonment (coastal) is last at 17.0. Six cities in this sample have SUHI below 0 °C for the same reason: the rural reference is itself extreme heat.

Full city report

LST max
Peak LST, overlay 43 °C.
LST mean
Mean LST of the observations (map colours from 35 °C).
NDVI
Mean NDVI, same window.
SUHI
Surface heat island (pixel minus rural T).
IVC units
IVC of official subdivisions in the city.

Methodology

1. Peak LST

Landsat 8 and 9, Collection 2 Level 2, thermal band 10, 1 April–31 May 2026. Fill, cloud and shadow pixels are dropped, as are temperatures outside 0–80 °C. The peak component is the urban fraction (GHSL SMOD ≥ 21) whose per-pixel maximum LST exceeds 43 °C, the ISO 13732-1 prolonged-contact limit.

2. Mean LST of the observations

Mean of the valid scenes in the same window. The mean component is the urban fraction whose composite-mean LST exceeds 43 °C. The map colour starts at 35 °C; the index still uses 43 °C.

3. NDVI

NDVI uses Landsat near-infrared and red reflectance. The greenness-deficit component inverts urban NDVI on a fixed 0.05–0.55 scale: less vegetation raises the index. The scale is not stretched to each city’s own min and max.

4. SUHI

SUHI is urban mean LST minus rural temperature. Rural land is GHSL SMOD 11–13, excluding water; if internal rural land is under 80 km², a 20 km GHSL belt is used. A negative SUHI adds nothing: a city is not rewarded for sitting in a hotter desert.

Land surface temperature. LST (°C) = 0.00341802 × DN + 149 − 273.15

DN is the Landsat Collection 2 Level-2 digital number of band 10 (ST_B10).

Peak component. zmax = urban fraction with LSTmax > 43 °C

Urban mask: GHSL 2020 SMOD ≥ 21. Threshold: ISO 13732-1 prolonged contact.

Mean component. zmean = urban fraction with LSTmean > 43 °C

Normalized difference vegetation index. NDVI = (NIR − red) / (NIR + red)

Greenness-deficit component. zNDVI inv = 1 − clip((NDVI − 0.05) / (0.55 − 0.05), 0, 1)

NIR and red are Landsat SR_B5 and SR_B4. clip(x, 0, 1) = min(max(x, 0), 1).

Surface urban heat island. SUHI (°C) = LSTmean, urban − Trural

Rural T is the GHSL SMOD 11–13 mean, excluding water. If internal rural land is under 80 km², a 20 km GHSL belt is used.

Heat-island component. zSUHI = clip(max(SUHI, 0) / 15, 0, 1)

Climate-vulnerability index. IVC = 100 × (zmax + zmean + zNDVI inv + zSUHI) / 4

Equal weights. Caps are fixed, so an IVC of 60 is comparable across cities.

Limits

The index does not include population, age, informal work or access to cooling. Unit grain differs (2 tehsils / subdistricts in Lahore, 41 qisms / district sections in Cairo). Rural T is GHSL, not a weather station. Landsat overpass is late morning. April–May 2026 is a common window, not each city’s worst month.

Annotated bibliography

  1. IPCC, Sixth Assessment Report, Working Group II (2022). — Frames climate risk as hazard × exposure × vulnerability. The city sample follows that product, not a hottest-place list.
  2. WHO / WMO, Heatwaves and Health (2015, updates). — There is no universal air-temperature threshold for mortality. Heat-health guidance is relative to the local climate; our 43 °C cut is a contact standard, not a death threshold.
  3. Kariyawasam et al., Sustainable Cities and Society / Oxford Smith School (2026). — Ranking of 205 cities over one million people. More than 95% of the worst-ranked lie in South/Southeast Asia and Africa; Basra, next to Baghdad, is first. Justifies including Karachi, Delhi, Cairo and Baghdad.
  4. Turner, Higgs, Sun et al., GUHVI, Urban Climate 64:102716 (2025). — Global urban heat-vulnerability index. Combines heat exposure, sensitivity and adaptive capability at neighbourhood scale; used here as a sampling frame, not as the Landsat score.
  5. Tuholske et al., PNAS (2021). — Urban person-days with wet-bulb globe temperature above 30 °C tripled from 1983 to 2016. Growth is concentrated in South Asian cities in this sample (Lahore, Dhaka, Ahmedabad).
  6. Mora et al., Nature Climate Change (2017). — Deadly heat conditions could become annual at +2 °C. Cited for Karachi’s humid-coastal regime after the 2015 heatwave.
  7. Gasparrini et al., The Lancet (2015). — In temperate cities, mortality rises faster per extra degree because the population is less acclimatised. Why Madrid, Barcelona and Rome stay in the sample even though April–May is spring.
  8. ISO 13732-1:2006. — Ergonomics of the thermal environment — skin contact. 43 °C is the prolonged-contact limit used as the Landsat overlay and as z_max / z_mean.
  9. Chakraborty et al. (2025), “Satellite-derived land surface temperatures strongly mischaracterise urban heat hazard.” — LST maps materials and shade; it does not replace air, humidity or wind for mortality. The IVC is a surface-exposure score.
  10. European Commission, GHSL SMOD V2-0, 2020 epoch. — Urban mask (SMOD ≥ 21) and rural reference (SMOD 11–13, no water). The 20 km belt is used when internal rural land is under 80 km².
  11. USGS / NASA, Landsat 8 and 9 Collection 2 Level 2. — ST_B10 surface temperature and SR_B4/B5 reflectance, 1 April–31 May 2026, fill/cloud/shadow mask.
  12. Knowlton et al. / NRDC–IITM, Ahmedabad Heat Action Plan (2013–). — Pioneer South Asian heat plan. The IVC still ranks Ahmedabad first because it measures remaining surface exposure, not the plan’s lives saved.