Can satellite AI show where a town stays cooler?

The study asks a practical question in an under-studied setting

Heat-risk research is dominated by large cities with dense sensor networks. Oued Zenati is different: a medium-sized agricultural municipality in Guelma province, northeastern Algeria, with a Mediterranean semi-arid climate, varied elevation and a seasonal wadi. The paper asks whether globally available satellite and geographic data can indicate where vegetation and landscape conditions are most compatible with cooler surfaces.

That is a useful research question. Smaller North African towns may have limited monitoring budgets even as hotter summers, land conversion and urban growth increase heat exposure. A repeatable remote-sensing workflow could help local teams decide where to investigate, preserve vegetation or place sensors. It would be especially useful if it transferred to places that cannot afford dense field campaigns.

The important word, however, is potential. The study did not observe people, indoor conditions, street-level shade or measured air temperature in candidate locations. Its output is a susceptibility map: a statistical estimate of places resembling a rule-defined ‘cool’ class. It is not a verified map of safe refuges during a heatwave.[1]

Seven satellite scenes describe one summer, not a changing climate

The authors used seven Landsat 9 Collection 2 Level-2 scenes acquired between 11 June and 15 September 2023. Each had less than 10% reported cloud cover. They applied pixel-level cloud masking and combined the scenes into a median composite. Reflective bands were available at 30-metre resolution; the 100-metre thermal product was resampled to 30 metres.

The model predictors were land-surface temperature, albedo, elevation, slope, distance to roads, distance to the hydrographic network and population density. Elevation came from the Shuttle Radar Topography Mission, built-up land from ESA WorldCover, and population density from a WorldPop product dated 2019. The authors also used ERA5 monthly aggregates to describe the 2023 climatic context.

This mixture creates a spatial snapshot from sources with different dates and resolutions. A 30-metre cell can contain roofs, trees, roads and soil at once; resampling thermal data does not create new 30-metre temperature detail. The seven scenes cover one warm season. They cannot establish whether a cool-looking area stays cooler across years, seasons, droughts or extreme heat events.[1]

The answer was defined by vegetation because no field labels existed

The study had no direct field observations of urban cool islands. Instead, the researchers generated reference labels from satellite indicators. A pixel could enter the cool class if its normalised difference vegetation index was at least 0.20 or its water index was at least zero. Non-cool candidates came from ESA’s built-up class or a combined bare-soil rule. Pixels meeting conflicting rules were removed.

During the selected summer, the water-index threshold identified no pixels: values stayed at or below minus 0.06 and the wadi was dry. Every positive reference point therefore came from the vegetation rule. Stratified random sampling produced 50 potential-cool points and 50 potential-non-cool points.

The authors excluded the exact label-generating vegetation, water and built-up indices from the final predictor set to reduce direct information leakage. That is a sensible safeguard, but it does not make the labels independent observations. The model is still being asked to recover a vegetation-based class using variables such as land-surface temperature, albedo and terrain that are environmentally related to that class.

In plain language, the analysis shows which model most successfully reproduces a remotely defined vegetation-cooling proxy. It does not show that the selected locations produced lower air temperatures for residents or reduced heat illness. Calling the 30 test points ‘independent’ means they were held out from model fitting, not that they came from an independent sensor campaign or another city.[1][2]

High scores rest on 30 test points

The researchers split the 100 proxy-labelled points 70/30 and compared Classification and Regression Trees, Random Forest, Gradient Tree Boosting and MaxEnt. They used the same seven predictors for each model and report that they did not conduct an exhaustive hyperparameter search.

Gradient Tree Boosting achieved the highest area under the receiver-operating-characteristic curve, 0.916, followed by Random Forest at 0.909, MaxEnt at 0.853 and the single decision tree at 0.667. After choosing model-specific probability thresholds, Random Forest recorded 90% accuracy, 0.80 Cohen’s kappa, 0.929 precision and a 0.897 F1 score.

Those percentages need their denominator. The test set contained 15 proxy-cool and 15 proxy-non-cool points. Random Forest labelled 13 of the 15 positive points and 14 of the 15 negative points correctly: 27 correct classifications in total. Gradient Tree Boosting made four errors, one more than Random Forest. With such a small test set, changing one prediction moves accuracy by 3.3 percentage points and can change the model ranking.

The split also appears to be random points within one municipality and one composite image. The paper does not report spatially blocked cross-validation that would force a model to predict a geographically separate area. Nearby pixels can share land cover, elevation and temperature, so a random split may overstate how well the method would work in an unseen neighbourhood or another town.[1]

The maps find plausible landscape patterns, not causal cooling effects

Across the tree-based models, SHAP attribution ranked elevation as the most influential predictor, followed by slope, albedo and land-surface temperature. Distance to roads, distance to the river network and population density contributed less. The susceptibility maps concentrated higher cooling potential in vegetated and agricultural areas and lower potential in built-up or bare areas.

These patterns are environmentally plausible, and agreement across several algorithms adds some internal reassurance. It does not prove that elevation causes a local cooling benefit or that planting vegetation at a low-scoring site would deliver the predicted effect. Terrain, land use and vegetation are correlated, while shade, irrigation, canopy structure, building height, street geometry, anthropogenic heat and wind were not directly represented.

The paper’s 2019 population layer is also not a measure of who is exposed at a particular hour. A cooler agricultural fringe may score well while being inaccessible to residents, and a small shaded street could disappear inside a coarse satellite pixel. Planning for human heat safety requires air temperature, humidity, radiant heat, walking access and information about vulnerable residents—not surface susceptibility alone.[1]

The disclosure text needs clarification

The accepted manuscript reports no competing interests and, in its formal funding statement, says the research received no funding. An acknowledgement immediately above it nevertheless thanks King Khalid University’s deanship of research and graduate studies for funding under grant RGP.2/4/47. That internal inconsistency should be corrected or explained in the final version; readers should not have to choose between two conflicting disclosure statements.

The authors provide a Google Earth Engine code link and say data are available from the corresponding author on reasonable request. Public code is valuable, particularly for a workflow intended to transfer across lower-resource settings. Reproducibility would be stronger with a versioned archive, exact study boundary, exported sample coordinates or a privacy-safe equivalent, and complete parameter records.

The study is an accepted, peer-reviewed article in press, not a preprint. The publisher says the final edited version will replace it. Editorial production may clarify wording and disclosures, but it will not by itself supply missing field validation.[1]

What would make the map actionable

A stronger next study would prespecify the model, then compare predictions with fixed and mobile temperature sensors across streets, parks, agricultural edges and built-up areas. Measurements should include daytime and night-time air temperature, humidity and mean radiant temperature over several seasons and at least one extreme heat event. Spatially blocked validation would test whether performance survives geographic separation.

Transfer tests in other Algerian and Mediterranean towns should keep the model and thresholds fixed before evaluation. Reporting by urban form, elevation and land-cover type would show where errors concentrate. Higher-resolution imagery and three-dimensional building and canopy information could help with street-scale conditions, but they should be judged against observations rather than another proxy map.

For policy, the useful role today is hypothesis generation. The map could help decide where to inspect, install sensors or investigate green-corridor protection. It should not determine where residents are directed during dangerous heat, where public money is committed, or which neighbourhood is declared adequately cooled without ground measurements and community knowledge.

The broader lesson is positive but restrained: accessible satellite data and machine learning can organise evidence about landscape conditions in under-studied places. In Oued Zenati, the workflow produces a plausible picture of vegetation-associated cooling potential. Its small, proxy-labelled test does not yet demonstrate reliable heat protection for people.[1][2]