NYCHA Climate Risk Index

Technical Methodology — encoding, normalization, parameters, data sources

1. Risk Formula

Risk is calculated independently for each of 4 climate domains, then averaged:

Per domain: adjusted = min_max(hazard × vulnerability) × (1 − effectiveness × mitigation)
Risk Index: = min_max( mean of 4 domain risks )

All values normalized to [0, 1]. Higher = more risk. A development with no hazard exposure in a domain scores 0 for that domain, pulling the overall average down.

2. Encoding Methods

MethodFormulaOutputUse Case
min_max(v − min) / (max − min)[0, 1]Higher raw value = higher score (solar radiation, building count below DFE)
min_max_inverse1 − (v − min) / (max − min)[0, 1]Lower raw value = higher score (groundwater depth: shallower = more hazard)
booleanTrue → 1.0, False → 0.0{0, 1}Binary presence/absence (stormwater exposure, cool roof)
categoricalLookup tableCustomOrdinal categories (Tidal onset: 2050s=1.0, 2080s=0.67, 2100=0.33, no onset=0; Equipment: Yes=1.0, Partial=0.5, none=0)

3. Normalization

StageMethodNotes
AttributesEncoded per method aboveAll output in [0, 1]
Sub-IndicesWeighted average of attributes, then min-max normalizedMissing attributes filled with 0.5 (neutral)
IndicesPer-site weighted average of sub-indices, then min-max normalizedConditional sub-indices (vulnerability/mitigation) that are N/A for a site are excluded from both the numerator and weight denominator for that site. Non-conditional missing sub-indices filled with 0.5 (neutral).
Domain RiskH × V product, min-max normalized, then mitigation appliedPer-domain before final averaging
Risk IndexMean of 4 domain risks, then min-max normalizedFinal output

Edge Cases

ConditionBehavior
All values identicalAll set to 0.5
Single valid valueSet to 0.5
No hazard exposure in a domainDomain risk = 0 (pulls overall risk down)
Hazard not present → vulnerability / mitigationNaN (displayed as “Not applicable”); excluded from composite index aggregation for that site
Groundwater depth ≥ 16 ftCapped before normalization → hazard scores 0; groundwater V and M set to N/A
Flood-exposed site with no buildings below DFEVulnerability floor of 0.1 applied (coastal, stormwater, groundwater)

4. Parameters

ParameterValuePurpose
Coastal DFE Buffer (2050)40 inches (3.33 ft)Buffer added to FEMA BFE; combined freeboard and sea level rise
Coastal DFE Buffer (2080)52 inches (4.33 ft)Buffer added to FEMA BFE; combined freeboard and sea level rise
Stormwater DFE Buffer24 inches (2.0 ft)Buffer added to stormwater BFE
Groundwater Depth Cap16 ftUSGS HERA threshold; depths ≥ 16 ft = no hazard
Mitigation Effectiveness (default)0.6Max risk reduction for heat domain
Mitigation Effectiveness (flood)0.9Enhanced for coastal, stormwater, and groundwater (engineered flood-protection measures)
Basement Depth Below Grade10 ftUniversal assumption for basement floor elevation (z_grade − 10 ft)
Vulnerability Floor0.1Minimum vulnerability for hazard-exposed properties with zero building vulnerability (applies to coastal, stormwater, and groundwater domains)
High-Rise Threshold60 ft (≈ 6 stories)Campus type classification
API Cache TTL30 daysAuto-refresh for FEMA, coastal surge, building data

5. Design Flood Elevation (DFE)

DFE = max(LiDAR DEM elevation within flood zone intersection) + freeboard buffer
DomainFlood SourceBufferOutput
Coastal (current)FEMA NFHL vector polygons40 inInhabited floors & basements below Coastal DFE (counts)
Coastal (2050)ArcGIS Coastal Surge 205040 inCoastal_DFE_2050
Coastal (2080)ArcGIS Coastal Surge 208052 inCoastal_DFE_2080
Stormwater (2050)Stormwater_2050_Moderate.tif24 inInhabited floors & basements below Stormwater DFE (counts)
Stormwater (2080)Stormwater_2080_Extreme.tif24 inStormwater_DFE_2080

Calculation Steps

  1. Identify flood zone overlap: Intersect the NYCHA development polygon with the flood zone geometry (FEMA NFHL vectors for coastal; stormwater raster cells > 0 for stormwater).
  2. Sample ground elevation: Extract LiDAR DEM raster values within the intersection area. Only valid (non-nodata) cells are considered.
  3. Determine Base Flood Elevation (BFE): Take the maximum ground elevation within the flooded portion of the development. This represents the highest point where floodwater is modeled to reach.
  4. Add freeboard buffer: Add the domain-specific freeboard (40″ coastal, 24″ stormwater) to the BFE. This produces the Design Flood Elevation (DFE).
  5. Count vulnerable buildings (inhabited floor): Using the Building Elevation Subgrade dataset (z_floor = first inhabited floor elevation, NAVD88), count buildings on the development where z_floor < DFE.
  6. Count vulnerable buildings (basement): Estimate basement elevation as z_grade − 10 ft (universal assumption; configurable via BASEMENT_DEPTH_BELOW_GRADE_FT). Count buildings where basement_elevation < DFE. This captures vulnerability to building systems and infrastructure even when inhabited floors are above the flood line.

If a development has no flood zone overlap, DFE is not calculated and the building count defaults to 0. If the DEM has no valid elevation data within the intersection, DFE is also not calculated.

Raster Alignment

When the stormwater raster and DEM have different resolutions, the stormwater mask is resampled (nearest-neighbor) to match the DEM grid before combining. For coastal DFE using FEMA vectors, the development-flood zone intersection geometry is used directly as a DEM mask, avoiding raster-raster alignment.

6. Conditional Relationships

Vulnerability and mitigation sub-indices are only calculated when the corresponding hazard is present; otherwise set to NaN (“Not applicable”). When the composite hazard/vulnerability/mitigation indices are aggregated, N/A conditional sub-indices are excluded from both the numerator and the weight denominator for that site — so sites are not penalized by forced zeros for domains that do not apply to them.

ConditionTrigger (ANY)Applies To
COASTAL_TIDAL_HAZARDPFIRM=True, Surge 2050=True, Surge 2080=True, or Tidal ≠ emptyCoastal mitigation
STORMWATER_HAZARDStormwater 2020, 2050, or 2080 = TrueStormwater mitigation
GROUNDWATER_HAZARDDepth_to_Groundwater_ft < 16 ft (the no-hazard cap)Groundwater vulnerability and mitigation
HEAT_HAZARDSolarRadGroundPct > 0, bldg_height_mean > 0, or Total_Roof_Area > 0Heat vulnerability and mitigation

7. Data Sources

Auto-Refreshing (API, 30-day cache)

DataAPIFields Used
FEMA PFIRM Flood Zoneshazards.fema.gov NFHL MapServer/28FLD_ZONE, ZONE_SUBTY, SFHA_TF, STATIC_BFE
Coastal Surge 2050/2080 (NPCC projections)NYC MOCEJ ArcGIS FeatureServer, serving New York City Panel on Climate Change (NPCC) sea-level-rise and coastal flood projections (reference link shows the 2050s 100-year floodplain; 2080s and tidal-onset layers are separate NPCC datasets published by NYC Open Data)Geometry (boolean intersection)
Building Elevation SubgradeNYC Open Data SODA (bsin-59hv)the_geom, z_floor
NYC Building FootprintsNYC Open Data SODA (5zhs-2jue)the_geom, height_roof; Shape_Area computed

Manual Refresh (local rasters)

FileSourceUpdate Trigger
Stormwater_2020_Light.tifNYC Department of Environmental Protection (DEP) — NYC Stormwater Flood MapsNew DEP release or NPCC projection update
Stormwater_2050_Moderate.tifNYC Department of Environmental Protection (DEP) — NYC Stormwater Flood MapsNew DEP release or NPCC projection update
Stormwater_2080_Extreme.tifNYC Department of Environmental Protection (DEP) — NYC Stormwater Flood MapsNew DEP release or NPCC projection update
NYC_DEM_LiDAR.tifUSGS 3DEP / NYC DoITTNew LiDAR acquisition
nyc_groundwater_depth_interpolated.tifUSGS / NYC DEPNew interpolation
NYC_SolarRadiation_Merged.tifNYC / CUNY Solar MapNew analysis

NYCHA-Maintained

FileContents
NYCHA Development Data Book.csvAdministrative data, Sandy/IDA flags, completion date, building type
2025 Resident Data Book.csvDemographics (Minors %, Seniors %)
Flood Protection Components.xlsxFlood hardening inventory, critical equipment status
335 Developments Hazard Exposure desktop.xlsxTidal flooding timeline, Cloudburst designations
Developments with Replaced Roofs.xlsxCool roof installations
DFTA NYCHA Cooling Centers.xlsxDFTA senior cooling center locations (union-flagged with DYCD into OnsiteCoolingCenter)
DYCD NYCHA Cooling Centers.xlsx (optional)DYCD community/youth cooling center locations; OR-combined with DFTA into the OnsiteCoolingCenter attribute. Loads gracefully if absent.

8. Spatial Operations

SourceMethodAggregation
Building Elevation PointsSpatial join (points within polygon)min/max/mean z_floor; count below DFE
NYC Building FootprintsCentroid-in-polygon joinmean height, sum area, count
Rasters (groundwater, stormwater, solar)Polygon mask extractionMean (groundwater depth), boolean any>0 (stormwater exposure), sum of ground-area pixels > 0 (solar radiation, continuous input to min_max)
FEMA NFHL / Coastal Surge vectorsDissolved union, prepared geometry intersectionBoolean intersects + DFE via DEM sampling
Tabular (Excel, CSV)Exact name match on DEVELOPMENTFirst match or aggregated per domain rules

All spatial overlays use EPSG:2263 (NY State Plane). Vectors reprojected to match raster CRS before extraction. Invalid geometries repaired via make_valid before dissolve.

9. Ranking & Classification

All indices use descending rank (Rank 1 = highest/worst). For flood-related sub-indices, only exposed properties are ranked; zero-exposure properties show “N/A”.

third = ceil(total / 3)
Rank ≤ third → HIGH  |  Rank ≤ 2 × third → MEDIUM  |  Rank > 2 × third → LOW

10. Mitigation Action Matching

Mitigation action recommendations are tailored to each development based on its hazard exposure level. The system determines which recommendations to show using the tertile ranking of each hazard sub-index.

Decision Logic

Hazard DomainHIGH (top third)MEDIUM (middle third)LOW (bottom third)
Coastal/TidalFull mitigation packageFull mitigation packageNo section shown
StormwaterComprehensive package (H&H analysis + critical systems + site strategies + entrances + buffering)Reduced package (H&H analysis + buffering only)No section shown
GroundwaterInvestigation stepsInvestigation stepsNo section shown
Extreme HeatShown for all exposed developments (any level except N/A)

Coastal Flood Mitigation Actions

BASELINE — Protect Critical Systems (boilers, electrical switchgear, generators, elevators, fire alarm, BMS, telecom):

Passive strategies preferred. Deployable strategies may be used where passive strategies are not feasible.

Protect entrances and exits: protect or relocate ground-floor spaces below DFE; install passive barriers at entrances.

BETTERMENT — Buffering strategies to reduce flood impacts: elevated curbs or floodwalls, regraded site to direct water away from buildings, vegetated berms along waterfront edges.

Stormwater Flood Mitigation Actions — HIGH Risk

BASELINE — All coastal baseline actions plus:

BETTERMENT — Buffering per NYCHA Stormwater Design Manual: elevated curbs, adapted topography, vegetated berms, subsurface tanks and nature-based solutions.

Stormwater Flood Mitigation Actions — MEDIUM Risk

BASELINE — Build H&H Model to refine DFE per Stormwater Design Manual.

BETTERMENT — Same buffering strategies as HIGH risk.

High Groundwater Mitigation Actions

Mitigation credit is given to properties with the following measures (weights for the Groundwater Mitigation sub-index):

Recommended further investigation steps for high-groundwater sites:

  1. Review existing geotechnical surveys and soil borings
  2. Consult geotechnical engineer on groundwater impacts
  3. Install groundwater monitoring wells
  4. Evaluate active subsurface drainage strategies

Extreme Heat Mitigation Actions

BASELINE:

BETTERMENT:

11. Key Assumptions

  1. Relative ranking: Scores rank NYCHA properties against each other, not against external safety standards.
  2. Equal domain weighting: Each hazard domain contributes 0.25 to the final Risk Index.
  3. Hazard-present gating: Mitigation only assessed where hazard exists; otherwise “Not applicable.”
  4. Residual risk: Full mitigation reduces risk by 90% (flood) or 60% (heat), never to zero.
  5. Near-term weighting: Current conditions (2020) weighted more heavily than distant projections (2080).
  6. Freeboard buffers: 40 inches coastal (current, 2050 scenarios), 52 inches coastal (2080 scenario), and 24 inches stormwater — conservative margins above the computed base flood elevation (FEMA BFE for coastal, derived stormwater BFE for stormwater).
  7. Groundwater vulnerability proxy: Estimated from building age (1.0) since direct below-grade condition data is unavailable. Applied only to sites where groundwater is a concern (depth < 16 ft).
  8. Basement depth assumption: Every building is assumed to have a basement 10 ft below grade (z_grade). This universal assumption may be replaced with site-specific data when available from NYCHA engineering teams or verified subgrade observations.
  9. Vulnerability floor: Flood-exposed properties with zero building vulnerability (no inhabited floors or basements below DFE) receive a minimum vulnerability score of 0.1, reflecting site-level risks such as access disruption, grounds flooding, and infrastructure damage.