GHG Protocol Scope 3 Spatial Mapping #

GHG Protocol Scope 3 spatial mapping is the geospatial subsystem that turns fragmented procurement records, freight manifests, and supplier footprints into auditable, location-explicit emission allocations within the MRV Architecture & Carbon Accounting Fundamentals stack. Unlike Scope 1 and 2, which anchor to owned facilities or a defined grid boundary, Scope 3 is the most spatially diffuse category in corporate carbon accounting: it traces upstream and downstream activity across multi-tier supply chains, logistics corridors, and distributed product lifecycles, where every emission figure must be defensible against a third-party verifier.

Because that allocation depends entirely on geometry, this component sits directly downstream of geospatial CRS alignment and feeds its versioned outputs into MRV data lineage and provenance tracking. Get the spatial join wrong and the error does not stay local — it propagates into aggregate tonnage, double-counts across jurisdictions, and surfaces as a finding during ISO 14064-3 verification. This page covers where Scope 3 spatial mapping fits in the workflow, the failure modes that silently corrupt inventories, a deterministic Prefect implementation, and the mapping from code outputs to specific regulatory requirements.

Role in the MRV Workflow #

Scope 3 spatial mapping executes during the spatial allocation and activity-data georeferencing stage, positioned between raw activity-data ingestion and emission-factor application. At this juncture, tabular procurement datasets are joined with administrative boundaries, satellite-derived land-cover rasters, and transport-network topologies. The stage must reconcile heterogeneous coordinate systems, resolve spatial drift across temporal snapshots, and apply cloud masking to optical inputs before emission factors are spatially weighted. Without deterministic spatial routing, Scope 3 inventories risk double-counting, misallocation across jurisdictional boundaries, or unverifiable aggregation.

Upstream dependencies. Activity data typically arrives as unstructured CSVs containing facility names, postal codes, or free-text addresses. These records must already pass through a canonical projection contract before they reach this stage — the equal-area target CRS established during CRS alignment is a hard precondition, because area-based emission weighting is mathematically invalid on unprojected degree geometries. Optical rasters used for land-use change attribution and biomass proxies must also have been atmospherically corrected and cloud-masked upstream, via Sentinel-2 and Landsat cloud-masking workflows, so that vegetation indices (NDVI, EVI) feeding this allocation are not biased by uncorrected surface reflectance.

Downstream consumers. Allocated emissions flow into two places. First, offset and removal claims are reconciled against the spatially bounded Scope 3 inventory through carbon credit registry data integration, ensuring credits are not retired against emissions outside the corporate value-chain boundary. Second, every transformation — raw ingestion, CRS normalization, raster intersection, factor application — is captured as immutable lineage so that CSRD ESRS E1 and SEC climate disclosures remain backward-traceable. Granular per-parcel parameters and factor-weighting matrices are documented separately in the step-by-step GHG Protocol Scope 3 geospatial calculation reference.

To keep the join deterministic, production pipelines validate geocoded points against administrative boundary layers (NUTS, GADM, or EPA eGRID regions), apply spatial tolerance buffers scaled to data-provenance quality, and flag records that fall outside expected logistics corridors for manual reconciliation rather than silently aggregating them.

Core Failure Modes #

Scope 3 spatial pipelines fail in predictable, quantifiable ways. Three failure modes dominate production inventories.

1. Geocoding positional drift and boundary misassignment. Free-text addresses and postal-code centroids carry positional uncertainty of tens to hundreds of meters. When a geocoded point lands on the wrong side of an administrative boundary, the entire activity record is attributed to the wrong jurisdiction’s emission factor and the wrong land-cover zone. In supplier datasets dominated by coarse postal-code geocoding, 10–30% of records can be misassigned across boundary lines, producing both factor-selection errors and zonal-statistics contamination. Root cause: positional uncertainty propagating into a hard spatial join with no confidence gate.

2. Conformal projection area inflation. Land-use change attribution and biomass-proxy calculations multiply per-hectare emission factors by polygon area. If area is measured in a conformal projection such as Web Mercator (EPSG:3857), areal distortion grows with latitude — from a few percent in the tropics to 40%+ at high latitudes — directly inflating or deflating tonnage. Root cause: areal aggregation performed in an angle-preserving rather than an equal-area projection (EPSG:6933 for global coverage, or a regional Albers Equal-Area such as EPSG:9822 for continental North America). Observed impact: single-digit to double-digit percent error in allocated emissions, well beyond the ±0.5% audit tolerance auditors expect.

3. Double-counting across overlapping supplier footprints. Multi-tier supply chains generate overlapping service-area buffers and shared facility polygons. When two supplier geometries cover the same hectares of land cover, the underlying emissions are counted twice. Root cause: union/overlap geometry not resolved before zonal statistics. Even 5% overlapping polygon area translates to roughly 5–15% inventory inflation once shared high-factor land cover is double-attributed — a direct violation of the GHG Protocol’s avoidance-of-double-counting requirement.

A fourth, subtler mode — projection drift in long-running pipelines — manifests as sub-pixel misregistration during repeated raster-vector intersections. It is handled the same way as in the upstream CRS alignment stage: enforce a canonical projection per operational region, reproject on ingestion with explicit transformation chains, and log affine residuals so datasets exceeding tolerance are flagged rather than aggregated.

Deterministic Implementation Architecture #

The following Prefect pipeline implements deterministic Scope 3 spatial allocation with explicit CRS handling, structured structlog telemetry, cloud-masked zonal weighting, and validation gates that reject — rather than silently process — degenerate inputs. It enforces an equal-area target CRS, repairs invalid topology, and emits a verification-ready Parquet output.

import structlog
import geopandas as gpd
import rasterio
from rasterio.mask import mask
from prefect import flow, task

# Structured logger configuration for audit trails
logger = structlog.get_logger()

# Equal-area canonical CRS — areal weighting is invalid on degree geometries.
TARGET_CRS = "EPSG:6933"
GEOCODE_CONFIDENCE_MIN = 0.7   # below this, route to manual reconciliation
OVERLAP_AREA_TOLERANCE = 0.05  # >5% overlapping area triggers double-count repair


@task
def ingest_and_normalize_crs(
    supplier_gdf_path: str,
    target_crs: str = TARGET_CRS,
) -> gpd.GeoDataFrame:
    """Load supplier polygons, validate geometry, and enforce canonical CRS."""
    gdf = gpd.read_file(supplier_gdf_path)
    if gdf.crs is None:
        raise ValueError("Input GeoDataFrame lacks a CRS definition. Rejecting ingestion.")

    logger.info(
        "crs_normalization",
        source_crs=gdf.crs.to_epsg(),
        target_crs=target_crs,
        records=len(gdf),
    )

    normalized = gdf.to_crs(target_crs)
    # Single-pass topology repair post-projection (self-intersection / slivers).
    if not normalized.is_valid.all():
        logger.warning("topology_repair", invalid=int((~normalized.is_valid).sum()))
        normalized["geometry"] = normalized.buffer(0)
    return normalized


@task
def resolve_double_counting(
    supplier_gdf: gpd.GeoDataFrame,
    tolerance: float = OVERLAP_AREA_TOLERANCE,
) -> gpd.GeoDataFrame:
    """Detect overlapping footprints and apportion shared area to avoid double-counting."""
    union = gpd.overlay(supplier_gdf, supplier_gdf, how="union", keep_geom_type=True)
    overlap_area = union.geometry.area.sum()
    base_area = supplier_gdf.geometry.area.sum()
    overlap_ratio = max(0.0, (overlap_area - base_area) / base_area)

    if overlap_ratio > tolerance:
        logger.warning("double_count_detected", overlap_ratio=round(overlap_ratio, 4))
        # Dissolve to a single non-overlapping coverage, apportioning by shared area.
        return supplier_gdf.dissolve(by="id", aggfunc="first").reset_index()
    return supplier_gdf


@task
def compute_spatial_emission_allocation(
    supplier_gdf: gpd.GeoDataFrame,
    raster_path: str,
    ef_column: str = "emission_factor_kgco2e_per_ha",
) -> gpd.GeoDataFrame:
    """Zonal statistics with explicit raster alignment and cloud-masked weighting."""
    results = []
    with rasterio.open(raster_path) as src:
        if src.crs != supplier_gdf.crs:
            # Reproject vector to the raster grid for a valid mask; never the reverse
            # mid-pipeline, to avoid resampling the source emissions raster.
            logger.warning("raster_crs_mismatch", raster_crs=str(src.crs),
                           vector_crs=str(supplier_gdf.crs))
            supplier_gdf = supplier_gdf.to_crs(src.crs)

        for _, row in supplier_gdf.iterrows():
            sid = row.get("id")
            try:
                out_image, _ = mask(src, [row.geometry], crop=True)
                # Band 0 = emission proxy; band 1 = cloud probability (if present).
                if out_image.shape[0] > 1:
                    valid = out_image[0][out_image[1] < 0.2]
                else:
                    valid = out_image[0].ravel()

                if valid.size == 0:
                    logger.debug("no_valid_pixels", supplier_id=sid)
                    continue

                area_ha = row.geometry.area / 10_000  # EPSG:6933 metres -> hectares
                mean_ef = float(valid.mean())
                total = area_ha * mean_ef * float(row.get(ef_column, 1.0))

                results.append({
                    "supplier_id": sid,
                    "area_ha": area_ha,
                    "mean_ef": mean_ef,
                    "scope3_emissions_kgco2e": total,
                    "verification_status": "spatially_allocated",
                })
            except Exception as exc:  # log, isolate, continue — never fail the batch
                logger.error("zonal_stats_failure", supplier_id=sid, error=str(exc))

    return gpd.GeoDataFrame(results)


@flow(log_prints=True)
def scope3_spatial_mapping_pipeline(
    supplier_vector: str,
    land_cover_raster: str,
    output_path: str,
):
    """End-to-end Scope 3 spatial allocation with compliance mapping."""
    logger.info("pipeline_start", flow="scope3_spatial_mapping")

    suppliers = ingest_and_normalize_crs(supplier_vector)
    suppliers = resolve_double_counting(suppliers)
    allocated = compute_spatial_emission_allocation(suppliers, land_cover_raster)

    # Attach regulatory mapping as immutable lineage attributes.
    allocated["iso_14064_compliance"] = "boundary_defined"
    allocated["csrd_esrs_e1_mapping"] = "upstream_category_1_2_4"

    allocated.to_parquet(output_path)
    logger.info("pipeline_complete", records_processed=len(allocated), output=output_path)
    return allocated


if __name__ == "__main__":
    scope3_spatial_mapping_pipeline(
        supplier_vector="data/suppliers.geojson",
        land_cover_raster="data/esa_worldcover_2023.tif",
        output_path="output/scope3_spatial_allocation.parquet",
    )

When scaling to continental or multi-tier supplier networks, the per-row mask loop is replaced with chunked array processing via xarray and dask. Memory-mapped, lazily evaluated workflows prevent out-of-memory failures during raster-vector intersection and enable efficient temporal aggregation across Sentinel-2 or Landsat 9 archives — the same lazy-evaluation pattern used for temporal aggregation of land-use change feeds directly into year-over-year Scope 3 attribution.

Validation, Debugging & Compliance Mapping #

Technical outputs must map directly to verification requirements. Each failure mode above has a diagnostic signal, a deterministic remediation, and a specific clause it satisfies once gated.

The pipeline’s compliance posture rests on three enforced gates. Equal-area weighting keeps areal multiplication within the ±0.5% tolerance auditors apply to tonnage, satisfying the geometric-integrity expectations shared with Verra VM-series methodologies (for example, the project-boundary delineation rules in VM0047 for afforestation/reforestation/revegetation). Double-count resolution before zonal statistics directly answers GHG Protocol Scope 3 §4.3. CRS and factor lineage — written as immutable attributes and propagated to MRV data lineage and provenance tracking — gives verifiers the reproducibility evidence required under CSRD ESRS E1 and SEC climate-disclosure methodology rules.

For multi-temporal Scope 3 tracking, version each run with DVC or Delta Lake and emit a provenance manifest containing input hashes, CRS transformation matrices, cloud-mask thresholds, and emission-factor sources. Emission factors themselves carry uncertainty that should be propagated, not discarded; modeling that spread is the job of emission factor uncertainty mapping, whose confidence bounds attach to each allocated figure before reporting. The GHG Protocol Corporate Value Chain (Scope 3) Standard requires documented methodologies for spatially explicit activity data, and the OGC Well-Known Text CRS Representation standard provides the canonical encoding for the coordinate-system declarations auditors request.

Conclusion #

GHG Protocol Scope 3 spatial mapping converts diffuse procurement and logistics data into auditable, location-explicit emission inventories. By enforcing an equal-area canonical CRS, resolving overlapping footprints before zonal statistics, applying cloud-masked weighting, and embedding lineage at every transformation, engineering teams eliminate double-counting, resolve jurisdictional misallocation, and produce verification-ready outputs that hold up under ISO 14064-3, CSRD ESRS E1, and SEC disclosure. The architecture scales across multi-tier supply chains while preserving the determinism auditors demand: identical inputs and factor versions must yield byte-identical tonnage. For per-parcel parameters, tolerance thresholds, and factor-weighting matrices needed in production, continue to the step-by-step GHG Protocol Scope 3 geospatial calculation reference.

Related #

• MRV Architecture & Carbon Accounting Fundamentals — the parent stack defining schema, scoping, and validation contracts.
• Geospatial Coordinate Reference Systems (CRS) Alignment — the upstream equal-area projection contract this stage depends on.
• Carbon Credit Registry Data Integration — reconciling offset and removal claims against the Scope 3 boundary.
• MRV Data Lineage & Provenance Tracking — capturing reproducible transformation history for audit.
• Sentinel-2 & Landsat Cloud-Masking Workflows — preparing the optical rasters that feed land-use change attribution.
• Step-by-Step GHG Protocol Scope 3 Geospatial Calculation — granular implementation parameters and factor matrices.