Mapping Environmental Equity: A Bivariate Risk Assessment of Adelaide

Executive Summary

As urban centers face rising temperatures, standard standalone heat mapping often falls short of identifying human risk. This analysis evaluates the intersection of Extreme Thermal Stress and High Social Vulnerability across the Adelaide metropolitan area. By combining individual data vectors into a unified bivariate spatial model, we can isolate exactly where infrastructure intervention is most urgently needed to protect vulnerable communities.

---

1. The Baseline Metrics: Heat vs. Vulnerability

Before analyzing the intersection of environmental and social risks, we can look at the raw data distributions. Below is an interactive map exploring the overall Social Vulnerability Index (SVI) across Adelaide’s statistical areas.

# Render the interactive SVI baseline layer
pal_svi <- colorNumeric(palette = "YlOrRd", domain = adelaide_geo$SVI)

leaflet(adelaide_geo) %>%
  addProviderTiles(providers$CartoDB.PositronNoLabels) %>% 
  addMapPane("top_labels", zIndex = 500) %>% 
  setView(lng = 138.60, lat = -34.93, zoom = 11) %>%
  addPolygons(
    fillColor = ~pal_svi(SVI),
    stroke = TRUE, color = "white", weight = 0.4, opacity = 0.5,
    fillOpacity = 0.65, smoothFactor = 0.2,
    highlightOptions = highlightOptions(weight = 1.5, color = "#666", fillOpacity = 0.85, bringToFront = FALSE),
    label = ~paste0("SVI Score: ", round(SVI, 2))
  ) %>%
  addProviderTiles(providers$CartoDB.PositronOnlyLabels, options = pathOptions(pane = "top_labels")) %>% 
  addLegend(pal = pal_svi, values = ~SVI, title = "Vulnerability Index", position = "bottomright")

---

2. Bivariate Intersection Analysis

To uncover the true risk profile, we blend both metrics into a 3x3 classification grid. This removes the “Airport Illusion”—where massive industrial spaces look highly critical on a standard heat map but contain zero residential risk.

On this map, look directly for the Dark Purple and Deep Red zones. These represent the critical areas where maximum heat directly targets communities with the least economic or systemic capacity to cope.

leaflet(adelaide_geo) %>%
  addProviderTiles(providers$CartoDB.PositronNoLabels) %>% 
  addMapPane("top_labels", zIndex = 500) %>% 
  setView(lng = 138.60, lat = -34.93, zoom = 11) %>%
  addPolygons(
    fillColor = ~biv_pal(bi_class),
    stroke = TRUE, color = "white", weight = 0.4, opacity = 0.5,
    fillOpacity = 0.8, smoothFactor = 0.2,
    highlightOptions = highlightOptions(weight = 1.5, color = "#333", fillOpacity = 0.95, bringToFront = FALSE),
    label = ~paste0("Heat Tier: ", heat_group, " | SVI Tier: ", svi_group)
  ) %>%
  addProviderTiles(providers$CartoDB.PositronOnlyLabels, options = pathOptions(pane = "top_labels")) %>%
  addControl(
    html = '
    <div style="background: white; padding: 10px; border-radius: 5px; line-height: 15px; opacity: 0.9;">
      <strong style="font-size: 12px; display:block; margin-bottom:5px;">Risk Matrix (Heat x SVI)</strong>
      <div style="display: grid; grid-template-columns: repeat(3, 20px); gap: 2px; width: 64px; float:left; margin-right:8px;">
        <div style="background:#64acbe; width:20px; height:20px;"></div>
        <div style="background:#627f8c; width:20px; height:20px;"></div>
        <div style="background:#574249; width:20px; height:20px;"></div>
        <div style="background:#b0d5df; width:20px; height:20px;"></div>
        <div style="background:#ad9ea5; width:20px; height:20px;"></div>
        <div style="background:#985356; width:20px; height:20px;"></div>
        <div style="background:#e8e8e8; width:20px; height:20px;"></div>
        <div style="background:#e4acac; width:20px; height:20px;"></div>
        <div style="background:#c85a5a; width:20px; height:20px;"></div>
      </div>
      <div style="font-size: 10px; height: 64px; display: flex; flex-direction: column; justify-content: space-between;">
        <div>High Risk (Both)</div>
        <div>Baseline (Low)</div>
      </div>
    </div>',
    position = "bottomright"
  )

---

3. Suburb Level Aggregation

By executing a spatial join against official Australian Bureau of Statistics (ABS) suburb boundaries, we can aggregate the data to pinpoint the exact neighborhoods topping the high-risk leaderboard.

While visual maps can be biased toward large geographic areas, this statistical breakdown tracks the raw frequency of high-risk census blocks.

datatable(suburb_summary, 
          options = list(pageLength = 10, dom = 'tip'), 
          rownames = FALSE,
          class = 'cell-border stripe')

Conclusion & Key Findings

Through this analytical pipeline, several critical insights were uncovered:

• The Fringes Bear the Brunt: While central Adelaide and major industrial corridors show high physical temperatures, the most severe socioeconomic intersection is heavily concentrated in the outer residential fringes.
• The Northern Plains Corridor: Suburbs such as Penfield, Virginia, and Angle Vale form a massive, continuous cluster of vulnerability, likely exacerbated by newer developments with minimal established urban tree canopy.
• The Southern Outliers: Kangarilla represents a severe anomaly with 178 high-risk blocks, identifying it as a massive point of concern where high environmental exposure intersects vulnerable residential demographics.