M.S. Thesis · Data Science · UMass Dartmouth

Global Trade as a Small-World Network

A geospatial small-world analysis of 35 years of international trade, 1988-2022

I modeled global supply chains as geospatially embedded small-world networks across 188 reporting economies and nine four-year intervals - showing that world trade is efficient and tightly clustered, robust to random shocks, yet fragile to the targeted loss of hub economies.

M.S. Data Science, August 2026 · Advisor: Prof. Firas Khatib · Committee: Dr. Ashokkumar Patel, Dr. Debarun Das

35 yrs
of UN Comtrade trade modeled (1988-2022)
188
reporting economies
613K
bilateral export observations
6 / 9
hypotheses confirmed (2 partial, 1 pending)
~80%
fragmentation under 15-20% targeted hub loss

Overview

Do the structural laws that govern social and biological networks also govern the world economy? This thesis answers empirically, building nine snapshots of the global trade network - one per four-year interval from 1988 to 2022 - and testing whether each exhibits the small-world signature: short average path lengths alongside high local clustering.

The work integrates network science with geospatial analysis (the study's SCNE framework), treating each economy as a node embedded at its real-world coordinates and each bilateral trade relationship as a weighted edge. Across all nine intervals the network is unambiguously small-world; its hub structure grows more concentrated over time; and, most consequentially for supply-chain risk, it proves highly resilient to random disruption but acutely vulnerable to the targeted removal of its most connected economies.

M.S. thesis, University of Massachusetts Dartmouth (2026). Read via DOI 10.62791/20625 →

Data pipeline

Turning raw UN Comtrade dumps into nine analysis-ready networks was a substantial data-engineering task in its own right.

Raw files → analysis-ready networks
35 annual UN Comtrade files→Batch ingest + consolidate→~25M product-level records→Geo-join + harmonize + clean (~2.5 GB)→Aggregate to 613K country-pair-year rows→9 interval networks (203-238 nodes)

The raw source was 35 annual UN Comtrade HS product-level export files (1988-2022). A reproducible pandas pipeline batch-ingests and consolidates them into roughly 25 million product-level records, joins every country to its coordinates, harmonizes ISO3 codes (e.g., merging Belgium-Luxembourg into BEL), and strips self-loops and regional aggregates (EU, World, "areas nes"). Cleaning yields a geo-referenced working set (~2.5 GB), aggregated to 613,252 bilateral export observations at the country-pair-year level and assembled into nine interval networks in NetworkX.

How to read the network
188 reporting economies (+ partner-only nodes) = up to 238 nodes→Trade flows = weighted edges→Hub economies = high degree→Small-world: high clustering + short paths

Reporting economies and their trade partners are nodes; bilateral trade relationships above a USD 1M threshold are weighted, undirected edges. A handful of hub economies carry a disproportionate share of connections, giving the network short average path lengths (any two countries sit a few steps apart) alongside high local clustering - the small-world signature, and the same hub structure that makes it efficient yet vulnerable to targeted disruption.

Methods

Construction
  • Nine four-year interval snapshots (1988-1991 … 2020-2022) plus one full-period aggregate, each an undirected, weighted trade network built in NetworkX (203 nodes / 3,694 edges in 1988 growing to ~238 nodes / 16,347 edges in the aggregate)
  • Economies embedded at real geographic coordinates (GeoPandas) for spatial analysis; a USD 1M minimum trade-value threshold filters negligible links
Small-world & scale-free testing
  • Small-world coefficient σ compared against degree-preserving random baselines for every interval
  • Average clustering coefficient (C), average shortest-path length, diameter, and radius tracked across all nine intervals
  • Power-law fitting of the degree distribution by maximum-likelihood estimation (Clauset et al.)
Hub, community & resilience analysis
  • Betweenness-centrality and Herfindahl-Hirschman Index (HHI) concentration tracked over time to measure hub concentration (H5)
  • Louvain community detection for trade-bloc structure; a distance-weighted gravity model linking geography to clustering (H6)
  • Percolation resilience: random failures vs. targeted hub removals, measuring giant-component fragmentation, plus a cascade simulation and a Laplacian-spectrum (algebraic-connectivity) analysis

Key findings

Nine formal hypotheses tested: six confirmed, two partially confirmed (H2 scale-free, H9 resilience-metric convergence), one pending (H8 SIRS risk transmission).

A persistent - and shrinking - small world

Every interval is unambiguously small-world (σ > 1.2), with consistently high clustering (C in the 0.82-0.86 range). As trade deepened, the world got "smaller": average path length fell from 1.75 to 1.38, with diameter 3 and radius 2 throughout.

Structure (H1, H3, H4)
  • Small-world signature holds across all nine intervals (σ > 1.2)
  • Clustering coefficient stays high (0.82-0.86); average path length 1.75 → 1.38
  • Densification confirmed: 203 nodes / 3,694 edges (1988) → ~238 nodes / 16,347 edges (aggregate)
  • Hub betweenness centrality grew more concentrated (~+40%, 1988-2020) - routing control consolidated in a few dominant hubs (H5, confirmed)

Scale-free, but only partially

The degree distribution follows a power law, but the exponent evolved from α = 2.38 (1988) to α = 2.60 (2020) - drifting outside the tight [2.0-2.4] band the hypothesis predicted, so H2 is recorded as partially confirmed rather than a clean fit. This is a real result, reported as partial rather than rounded up.

Geography still shapes the network

Despite globalization, geography remains a strong organizing force: geographic proximity and bilateral clustering correlate at r = 0.72 (H6), and Louvain community detection recovers three geographically coherent blocs - Asia-Pacific, European, and North American. Cross-regional corridors carry global connectivity.

Resilient to chance, fragile to intent

The network absorbs 60-75% random node removal before fragmenting - but removing just 15-20% of its hub economies collapses the largest connected component from 98% to 22% (H7). A severity-0.8 cascade shock to China produces ~80% immediate trade loss with a 40-step recovery under rerouting. That asymmetry is precisely the supply-chain resilience question: it identifies which economies, if disrupted, would break global trade.

Figures

Selected results from the thesis: the shrinking, densifying world, growing hub concentration, and the rise of new hubs.

The network, then and now

The small-world signature over time

The rise of China

Why it matters

Supply-chain risk

Which nodes break the network

The hub-fragility result is a supply-chain resilience question in disguise: it pinpoints the economies whose disruption would fragment global trade - the same analysis risk and logistics teams run on their own networks.

Data engineering

Scale, handled honestly

A reproducible pipeline distilling ~25M raw records into 613K bilateral observations and nine analysis-ready networks - the kind of end-to-end data work that underpins any real analytics system.

Method transfer

Structure over surface

The same small-world tools that describe trade apply to communication, logistics, and multi-agent systems - a through-line into my work on efficient agent coordination.