QGIS · Estonian 15-minute accessibility data
Finding the empty parcels where new biodiversity reaches the most children and families — translating the ..._W15 walk-accessibility attributes into a defensible priority score and a do-this-first matrix.
We are finding the empty plots of land where planting new, biodiverse green space would do the most good — above all for children and families.
For every location in the city we ask four simple questions: do people live here, are children nearby, is nature already missing, and is it a lively walkable place? The answers become a single 0–100 priority score and a clear map of where to act first.
The approach rests on published health research and a transparent, repeatable method that anyone is free to audit.
Adding forest floor, turf & planting to daycare yards measurably improved children's immune regulation within 28 days. Contact with biodiversity, not just "green", is the active ingredient.
Roslund et al. 2020, Science AdvancesGIS multi-criteria siting consistently finds population density the single strongest factor (weight ≈0.36) for where new green space delivers most benefit.
Frontiers Public Health 2024; Env. Systems Research 2024The 3-30-300 rule: everyone within 300 m of quality green space. <15% of people in 862 European cities comply — so prioritise the under-served, not the already-green.
Konijnendijk 2023; Nature Communications 2026A Philadelphia randomised trial: greening vacant lots cut residents' depression by up to 69% in poorer areas, and lowered stress/heart rate — cheap and scalable.
South/Branas et al. 2018, JAMA Network OpenWhy these weights — the academic argument. The weights are not arbitrary; they encode a hierarchy that the evidence supports, while the exact magnitudes remain an explicit, testable calibration.
Biodiversity deficit (0.35, the largest). The programme's purpose is to add biodiversity where people currently lack it, so deprivation carries the heaviest weight. Two lines of evidence justify this: the health case — direct contact with environmental biodiversity measurably improves children's immune regulation [1, 2] — and the equity case — under the 3-30-300 standard, fewer than 15% of residents across 862 European cities have adequate green access [4, 5]. Targeting the gap is the single strongest lever.
Household demand (0.25) and children reach (0.20). GIS multi-criteria siting studies consistently find population the dominant determinant of where new green space delivers benefit (the leading factor, with weights near 0.36 in AHP analyses) [6, 7]; demand is weighted second. Children are a distinct, high-value subset of that population — the health evidence is child-specific — but a complement to, not a replacement for, the population term. Together the two “people” terms (0.45) slightly exceed the deficit, so the model balances where the need is against how many people are reached: need takes the largest single weight, people the largest combined weight.
Walkable context (0.20). Green space is used and cared for best when embedded in active, walkable neighbourhoods; the 15-minute-city accessibility literature motivates rewarding integration into daily life, set equal to children as a supporting rather than leading factor.
Status of the weights. The literature robustly supports the ranking and direction of the four terms; the precise values (35/25/20/20) are an expert judgement by the research team, offered as a transparent hypothesis to be stress-tested with a sensitivity analysis or a formal AHP elicitation.
Every location gets a score out of 100, built from four simple questions. The share next to each one is how much it counts towards the final score.
A new green space helps more people when more people live next to it. We simply count the homes in each area, so that investment lands where it benefits the most residents.
The health benefits of contact with nature are strongest for children. We look at how many kindergartens, crèches, schools, playgrounds and children's activities sit within a 15-minute walk — a good signal of where families with young children actually are.
This is the heart of the programme, so it carries the most weight. If people can already walk to a park, trail, swimming spot or outdoor sports ground, the area is well served — so it scores lower for new greening. Where there is nothing green within a 15-minute walk, it scores higher. This is how we steer investment towards the people who are missing out today, rather than adding more green where there is plenty already.
We favour places that are already part of everyday life — where people walk to the shop, the pharmacy, the doctor, the post office, the library or a bus stop. Greening here completes a living neighbourhood instead of landing in the middle of nowhere.
Here the wheel is the score. Each wedge is one weighted term in the regreening formula, and the size of the wedge is its weight — so the biggest slice, the biodiversity deficit, is what moves the score most. The outer ring lists the attributes that feed each term. ▲ means more of it raises priority; ▼ means the term is inverted, so low outdoor-nature access is what scores high. Everything else in the dataset (health, culture, civic, retail) is mapped in the reference table below but does not enter the formula.
| Category | Dataset attributes (Estonian original in brackets) | Role in regreening score |
|---|---|---|
| Children & learningBenefit driver ▲ | Kindergarten (Lasteaed) · Childcare (Lastehoid) · Basic school (Põhikool) · Upper-secondary (Gümnaasium) · Higher ed (Kõrgkool) · Vocational (Kutseõppeasutus) · Continuing ed (Täienduskoolitus) · Playground (Mänguväljak) · Kids' hobby & playroom (Huvi, Mängutuba, Huvikeskus) · Children's shop (Pood-Lapsed) · underlying demand: residents DUs_Count / DUs_Sum | Children reach ↑ + demand (DUs) |
| Health & careBenefit driver ▲ | Pharmacy (Apteek) · GP (Perearst) · Specialist (Eriarst) · Dentist (Hambaarst) · Ambulance (Kiirabi) · Care home (Hooldekodu) · Day centre (KOV-Päevakeskus) · Health service (Teenus-Tervis) | Context ↑ vulnerable-group benefit |
| Nature & recreationBiodiversity gap ▼ | Walking routes (Liikumine-Jalutamine) · Dog walking (Koduloomad) · Nature destination (Sihtkoht) · Swimming (Suplus) · Theme park (Teemapark) · Outdoor sport (Välisport: staadion, tennis, ballcourt, pool) | Biodiversity deficit ↓ scored inverted — the core |
| Everyday & foodLived-in context ▲ | Grocery (Toidupood) · Market (Turg) · Daily goods (Tarbekaubad) · Eatery (Toidukoht) · Bar / café / pub (Toidukoht-Bar) · Post (Post) | Context ↑ |
| Movement & accessLived-in context ▲ | Transit terminal (ÜT-terminal) · Bike parking & rental (Rattaparkla, Rattarent) · Car rental (Autorent) · Fuel (Kütus) | Context ↑ (transit) |
| Community & cultureLived-in context ▲ | Local government (KOV) · County centre (Maakonnakeskus) · Church (Kirik) · Culture centre (Kultuurikeskus) · Public library (Raamatukogu) · Cinema (Kino) · Museum (Muuseum) · Theatre (Teater) · Gallery (Galerii) · Hall & stage (Saal, Lava) · Nightclub (Ööklubi) · Police (Politsei) · Fire/rescue (Päästekomando) · Civil shelter (Varjumiskoht) · Cemetery (Kalmistu) · Public toilet & waste (WC, Konteinerid) | Context ↑ (library) |
| Retail & servicesLow relevance | Clothing & textile (Riided, Tekstiil) · Electronics (Tehnika) · Furnishings (Sisustus) · Books/stationery/music (Raamatud, Kirjatarbed, Muusika) · Plants (Taimed) · Hardware (Ehitus) · Sport & bikes (Sport, Rattad) · Hobby/pets/vehicles/optics (Hobi, Koduloomad, Mootorsõidukid, Prillid) · Beauty & care (Pood-Ilu, Teenus-Ilu, Massaaž) · Gym & indoor sport (Sisesport, Jõusaal, Ujula) · Home services (Kodu, Laenutus, Print, Autopesu) · Hotel (Majutus) | — |
Note: the role labels above describe the whole category. Only the specific attributes in the scored baskets actually enter the formula — e.g. within Health & care only pharmacy and GP feed walkable-context, while ambulance, dentist, specialist and care home are not scored. Car rental, fuel and transit sit under movement/services and are not part of the biodiversity (nature) term.
The four colours match the index map. As the greening case weakens (deep green → orange), the binding constraint shifts from missing nature to missing everyday services. The top two classes say “add greenery”; the bottom two say “greenery helps, but service access is the real gap.”
Dense, child-rich areas that lack nearby nature but already have everyday services within reach. New biodiverse green here reaches the most people and completes the neighbourhood.
→ Plant first.A clear greening gap with plenty of people to benefit. The case for new nature is strong; act here alongside the top tier.
→ Plant / pipeline.These areas score lower mainly because everyday services are sparse — not because nature isn’t wanted. Greenery still adds value, but improving access to shops, care and transit matters as much.
→ Pair greening with services.The dominant shortfall is missing everyday services (and often few residents). Greening is a welcome co-benefit, but the primary intervention is basic amenities and access.
→ Address services; green opportunistically.The precondition. Only empty (or convertible) parcels can actually be regreened — gates everything.
Public land = act directly. Private = strategic target needing negotiation.
Already in your data. Dense surroundings = more people reached per intervention; very low density caps impact.
During the workshop, two teams proposed their own regreening priorities. Each is rebuilt below as a wheel using only the attributes in our dataset, with equal weight per category (20% each) and a presence-based reading (more of a category within a 15-minute walk → higher priority). Where the data cannot represent a category, it is flagged in grey.
| Team category | Mapped to our data | Status |
|---|---|---|
| T1 · Streets (home streets) | Dwellings (DUs) as a residential proxy | ~ Proxy — no street layer |
| T1 · Education | Kindergarten, childcare, schools, higher & vocational | ✓ Full |
| T1 · Shops — food (parking lots) | Grocery, market, eatery, bar | ~ Partial — parking lots ✗ no data |
| T1 · Services | Pharmacy, post, bank, fuel, beauty, home, WC… (Teenus*) | ✓ Full |
| T1 · Dense housing | Dwellings (DUs_Sum) | ✓ Full (same signal as “streets”) |
| T2 · Home (kodu) | Dwellings (DUs) | ✓ Full |
| T2 · Teenused – TOIT (food) | Grocery, market, eatery, bar | ✓ Full |
| T2 · Liikumine (movement) | Walking, dog-walking, nature destinations, swimming, sport | ✓ Full * |
| T2 · Huvid (interests) | Hobby, hobby centre, culture (cinema, museum, theatre…) | ✓ Full |
| T2 · Haridus (education) | Kindergarten, childcare, schools, higher & vocational | ✓ Full |
team1_score and team2_score.Each wheel is really a different story about what a good city is and where green belongs. None is more correct than another — they optimise for different values. Here is what each narrative sees, and what changes when you follow it.
We scored all 5,966 parcels on each of the three indices and compared them directly. Two questions decide the strategy: how strongly do the indices agree (correlation across every parcel), and how much do their shortlists overlap (the top 20% each index would actually fund). The answer is not a clean three-way split.
Across every parcel the two team indices correlate at r = 0.96; their top-20% shortlists overlap by 88%. Statistically they are one narrative under two names.
Pearson r 0.96 · 1,053 of 1,194 sharedFair-Share is not the teams' opposite: it correlates r = 0.83 with Team 1 and r = 0.68 with Team 2. All three lean on the same dense, populated city — they share most of their priorities.
Pearson r 0.83 / 0.68Where they part is the tail: 222 parcels — about a fifth of the Fair-Share shortlist — are sites neither team picks. They are quieter, under-served places that rank well below the teams' funding cut-off (≈ 30 / 100 on the team indices).
222 of 1,194 · team score ≈ 30| If you follow… | Fair-Share (academic) | Active city (Team 1) | Everyday-life (Team 2) |
|---|---|---|---|
| Shortlist overlap with the other two | shares 831-parcel core; adds 222 unique sites | 88% identical to Team 2; 80% inside Fair-Share | 88% identical to Team 1; 70% inside Fair-Share |
| Green goes to… | quieter, under-served, family-heavy areas missing nature | dense, central, high-traffic blocks | complete neighbourhoods with services & leisure |
| Who gains most | children & residents currently without nature | the many people already in busy areas | residents enjoying daily amenities |
| Headline win | fairness & child-health equity | visibility & footfall | everyday quality of life |
| Main risk | some chosen sites are less prominent | green-poor areas stay green-poor | same — under-served areas overlooked |
| City it builds | a more equal green city | a more concentrated green city | a more liveable green city |
fid; columns regreen_score, team1_score, team2_score).The layer is a polygon grid of 125 m × 125 m cells (cell side derived from the coordinate spacing in the data), carrying INSPIRE/EEA-style identifiers (e.g. GRD_ID “125mN6611375E595375”). Coordinates (E ≈ 595,000; N ≈ 6,611,000) are consistent with the Estonian national projected CRS L-EST97 / EPSG:3301 — to be confirmed in the layer properties. A projected, metric CRS is appropriate because the analysis is grid-based and the accessibility metric is distance-defined.
Each POI-…_W15 field is a continuous “within-15-minute-walk” accessibility value (sample values ≈ 0.0006–7.5), i.e. a distance-decayed accessibility surface rather than a raw count of reachable places; the model treats these as interval-scaled access intensities. Dwelling units (DUs_Sum) provide the demand measure. Two data-quality features shaped the implementation: several numeric fields are stored as text with “-” / blank placeholders, and field names contain Estonian diacritics (ä, õ, ü, š) and parentheses.
Each of the four terms is reduced to one value per cell (a single field, or the sum of its member attributes) and min-max normalised to [0, 1]. Because the inputs are non-negative and their minimum is effectively zero, this simplifies to division by the layer maximum:
x̂ = x / max(x) ( max taken across all cells in the dataset )Observed maxima in the working dataset were: dwellings 774, children 83.9, biodiversity-access 22.4, walkable-context 114.9. Because normalisation is relative to the dataset maximum, scores are comparable only within this extent and will shift if the study area or its maximum changes.
Weights sum to 1.0 and encode the evidence in Part 1: the biodiversity deficit carries the largest weight (0.35) as the core of a biodiversity-led, equity-driven intervention; household demand (0.25) reflects the dominance of population in green-space siting studies; children reach (0.20) operationalises the child-health objective; walkable context (0.20) keeps interventions in active neighbourhoods. The two “people” terms together (0.45) make reaching residents and children the largest combined emphasis. Weights are expert-set and should be treated as a hypothesis to be calibrated.
regreen_score = 100 × ( 0.25·D̂ + 0.20·K̂ + 0.35·(1 − Ĝ) + 0.20·Û )where D̂, K̂, Ĝ, Û are the normalised demand, children, biodiversity-access and context terms. The biodiversity term is inverted as (1 − Ĝ) so that low current access scores high. The model is linear-additive and fully compensatory (a strong term can offset a weak one).
DUs_Sum, the dwellings within each 125 m cell._W15 fields — destinations reachable within a 15-minute walk. Demand (DUs) is the cell's own dwelling count, not a walking catchment.For each of the 5,966 parcels the method outputs a 0–100 priority score, a 1–5 tier, and a per-term contribution breakdown. The actionable outputs are: (1) a prioritised map; (2) the top quintile — about 1,063 parcels — under the Fair-Share method; and (3) across all three narratives, a consensus shortlist of parcels every method ranks top-priority, recommended as the first wave of action.
For mapping, scores are reclassified into five quantile tiers (each ≈ 20% of parcels; tier 1 = top priority), in preference to fixed cut-offs because unpopulated cells carry a non-trivial deficit baseline (~35 points). Because the score is additive it is also decomposed into the points contributed by each term (c_demand 0–25, c_children 0–20, c_greendef 0–35, c_context 0–20, summing exactly to the score), and each parcel is tagged with its dominant driver — the recommended tool for auditing unexpected values.
An initial single Field Calculator expression returned NULL for every feature. Three compounding causes were identified: the source layer is read-only (delimited-text / WKT); numeric fields stored as text made arithmetic evaluate to NULL; and field references with Estonian diacritics did not resolve at evaluation time inside a nested-aggregate expression. The model was re-implemented as PyQGIS scripts that match field names via Unicode NFC normalisation, cast values tolerantly (treating “-”, blanks and text as 0), compute the layer maxima in Python, and write a new GeoPackage. The pipeline is reproducible and self-checking — each script prints its maxima, feature counts, score range and any unmatched field names.
| Output field | Type | Meaning |
|---|---|---|
| regreen_score | Double | Composite priority index, 0–100 |
| regreen_class | Integer | Quantile tier 1–5 (1 = top priority) |
| regreen_tier | Text | Tier label for the legend |
| c_demand / c_children / c_greendef / c_context | Double | Points contributed by each term (sum = score) |
| top_driver / top_share | Text / Double | Largest-contributing term and its share of the score |
The following are offered for critical review; several bear directly on interpretation.
| Issue | Effect on results | Recommended mitigation |
|---|---|---|
| Green access is a proxy | No direct vegetation / NDVI or park-area input; outdoor-recreation POI access stands in for “biodiversity”. A biologically green but POI-poor cell reads as deficit, and vice-versa. | Cross-check the deficit term against NDVI / land-cover; consider adding a vegetation index. |
| Deficit saturates in empty land | Where no outdoor-recreation POI is reachable, the full 0.35 term is awarded regardless of population, inflating remote / unpopulated cells. | Gate the score on dwellings > 0, or multiply the deficit term by a demand factor so it only counts where people are present. |
| Extent-relative normalisation | Min-max scaling is relative to the dataset maximum; scores compare only within this extent and are sensitive to outliers. | Use robust scaling (percentile capping) or fixed reference maxima for cross-study comparability. |
| Fully compensatory model | Linear addition lets a strong term mask a weak one (e.g. high reach hiding that a cell is already green). | Test a multiplicative / geometric-mean variant; always report term values beside the score. |
| Expert weighting | The 0.25 / 0.20 / 0.35 / 0.20 weights are subjective. | Sensitivity analysis (± weight perturbation) and/or AHP or stakeholder calibration. |
| Accessibility metric semantics | The exact definition of the W15 values (gravity / decay vs count) is unconfirmed, affecting interpretation. | Confirm the metric and its units with the data provider. |
| MAUP / grid effects | Results depend on the 125 m cell size and grid origin; sub-cell variation is hidden. | Test sensitivity to resolution; validate against parcel-level data where available. |
| Playground = demand | Playgrounds are treated as a children-presence signal, not green supply — a defensible but consequential choice. | Confirm with stakeholders; optionally test as partial green supply. |