Karachi / 2001 delimitation / 184 units

Karachi's fire risk,
finally on a map.

B-MCRI is an independent, open fire-risk map for Karachi. It takes the things that make a fire worse (packed housing, blocked roads, a slow response, a history of past incidents) and turns them into a single score for every block.

0
Union councils
0
Cantonment boards
0
Towns covered
0
Open data sources
The Index / interactive

Poke around the map

Tap any unit to see its structural index R̂, the Bayesian-network posterior, and the season it's in. Boundaries follow the 2001 city-district delimitation: 178 union councils across 18 towns plus 6 cantonment boards. The newer 246-committee structure isn't mapped yet.

184 units · model estimate
01 / Method

Four models, stacked.

No single model sees the whole picture, so B-MCRI stacks four. Each one answers a question the others can't, and together they fold into a single score you can actually read.

01
MCRI · the base layer
A multi-criteria risk score for every block

Start with what you can see on a map: how tightly the buildings are packed, how a fire truck would reach the street, how far the nearest station is, what the land is used for, and whether it's an informal settlement. Those roll up into four dials (Hazard, Exposure, Vulnerability, and protective Capacity) and combine into a score from 0 to 1. Everything after this is just a correction on top.

02
Bayesian network · the reasoning
What the evidence actually implies

Risk factors aren't independent. Old wiring, occupancy, and season all lean on each other. The Bayesian network writes those links as a graph and works out P(Fire | evidence), updating its belief as conditions change and staying sensible even when half the inputs are missing.

03
Hidden Markov model · the season
Which regime the city is in right now

Fire risk breathes with the calendar. From the run of past detections, an HMM picks out three hidden moods (Calm, Transitional, Dry), each nudging the structural score up or down. It's how a flat map learns to tell January from a dry, load-shedding July.

04
Hawkes process · the aftershock
One fire makes the next more likely

Urban fires cluster. One blaze pulls in nearby crews, jumps across shared walls, and briefly raises the odds of another next door. The Hawkes process captures that as a conditional intensity λ(t): a background rate plus a decaying jolt after every incident. Have a go below.

λ(t) · live conditional intensity
λ(t) now: 0.20 background μ: 0.20 events in view: 0

Click the plot (or the button) to start a fire. Watch the intensity jump, then cool back down toward the background rate. That's the same self-excitation the model applies across the city.

02 / Provenance

Where the numbers come from.

Five feeds, each open or public-domain. Some are live; some are frozen at a census epoch. All of it is traceable.

SourceWhat it feedsVariables derived
01
OpenStreetMap
OSM Foundation · community-mapped
Continuous
Road network, building footprints, fire-station locations, and land use across Karachi's districts and towns.
Road accessibility · building density per cell · response-time estimates · land-use multipliers
02
FIRMS
NASA · fire information for resource management
Near real-time
Thermal anomalies from MODIS (1 km) and VIIRS (375 m), with time, confidence, and radiative power. A proxy for activity, not the real thing (see Limits).
Detection frequency · clustering parameters · HMM training sequences · Hawkes baseline λ₀
03
PBS Census 2023
Pakistan Bureau of Statistics · national census
2023 epoch
District population and household totals from the published tables. Informal-settlement extent is not a census field; it's a derived proxy from built-up density and footprint shape.
Exposure weights · baseline denominators for the Exposure and Vulnerability terms
04
WorldPop
Univ. of Southampton · gridded population
2020 · 100 m
Constrained top-down gridded population for Pakistan, spreading census totals onto settled pixels at sub-district scale.
Per-unit population disaggregation · grid exposure surface beneath the Exposure term
05
GHSL
EC JRC · global human settlement layer
R2023A · 2020 · 100 m
Built-up surface from GHS-BUILT-S, quantifying built density per cell across the metro area.
Built-up density in the Hazard term · informal-settlement proxy (with OSM morphology)
03 / Honesty

What it can't tell you.

Every model has limitations. These are the main ones, set out plainly rather than buried in a footnote.

01
Detection, not incidence
FIRMS sees thermal anomalies, not building fires. Many structural fires stay below MODIS and VIIRS thresholds, while flaring, industrial heat, and open burning trip detections that aren't fires in our sense. The Hawkes baseline and HMM sequences inherit that bias.
02
No ground truth
There's no public dispatch record from Rescue 1122 or the KMC Fire Brigade to validate against. So R̂ is uncalibrated: it ranks units by modelled risk, but its absolute values carry no frequency meaning.
03
Census uncertainty
The 2023 census is contested, with acknowledged undercounts across all seven Karachi districts. Population enters the Exposure and Vulnerability terms directly, so that uncertainty flows into R̂ without being shown in it.
04
Point estimates only
The index reports posterior means without credible intervals. Two units separated by 0.02 in R̂ should not be treated as meaningfully different.
05
Old boundaries
Coverage is the 2001 delimitation: 178 union councils across 18 towns plus 6 cantonments (184 units). The boundary file predates the 2022 redistricting, so today's 246 union committees aren't represented, and units are keyed to the older town-level names.
0
Probabilistic components fused into one composite score
0
Primary data sources, all open or public-domain
0.0M
Karachi Division residents living with recurring, preventable fire risk
04 / About

An independent tool,
not a product.

B-MCRI started from a gap. Karachi keeps losing people to preventable fires in dense, hard-to-reach settlements, and there is almost no spatial risk data anyone can act on.

It combines classic spatial risk scoring with probabilistic methods (Bayesian networks, Hidden Markov Models, and self-exciting Hawkes processes) into an index that updates as new incident data comes in. The goal is a working tool, not only a paper.

Independent research Open access Karachi, Pakistan
05 / Contact

Get in touch.

Collaboration, data access, or a gap you've spotted in the method: email is the best way to reach me.

02
Location
Karachi, Sindh, Pakistan
03
Model framework
Publication in progress

“The point of a risk map isn't to predict the next fire. It's to argue, with evidence, about where the next truck, hydrant, or inspection should go.”