Feed 0% source
Engineering AI-generated

Catastrophic disruption cascades driven by the nonlinearity of systemic risk

Generated by a local model (nvidia/Gemma-4-26B-A4B-NVFP4) from a scientific paper, claim-checked against the full text. Provenance is open by design.

When one company fails, others can usually find a new supplier. However, if a specific group of companies fails at the same time, they might leave no alternatives. This can cause a massive chain reaction that shuts down much of the economy.

Modern supply chains are highly specialized, interconnected webs of buyer-supplier dependencies. While the failure of a single firm typically causes localized ripples, recent global events like the COVID-19 pandemic have demonstrated that these disruptions can escalate into system-wide collapses. Researchers have long tried to quantify this "systemic risk"—the danger of a local failure cascading through a network. Modeling these shifts has been difficult due to a lack of granular, firm-level data.

A new study reveals that our current understanding of risk may be incomplete. The authors report that the systemic risk caused by the simultaneous failure of a few firms can be vastly larger than the sum of their individual impacts. In the national supply chain of Ecuador, they found instances where combined failures produced risk amplifications of up to a factor of 257.

The limit of sublinear risk models

Traditional methods for quantifying systemic risk often assume a "sublinear" relationship. In finance, the authors highlight the DebtRank model. This model tracks how a bank's distress devalues the assets of its creditors. As shown in, DebtRank is inherently sublinear. This means the total loss from two banks failing simultaneously is always less than or equal to the sum of their individual losses. This happens because the model includes built-in constraints. For example, a bank cannot lose more equity than it actually possesses.

Supply chains, however, operate under different logic. The authors argue that existing models fail to account for a critical compensatory mechanism: supplier substitutability (the ability of a buyer to switch to a different vendor). In a healthy network, if a supplier fails, a customer can often pivot to an alternative provider. This acts like a shock absorber in a vehicle. It dampens the impact of a bump. But the authors demonstrate that when multiple suppliers fail at once, these "shock absorbers" vanish. This turns a manageable hiccup into a catastrophic cascade.

Breaking the substitution mechanism

To capture this nonlinearity, the authors utilize the Economic Systemic Risk Index (ESRI). This metric estimates the fraction of total production lost due to a firm's failure. Their model moves through the network in two directions. First, it calculates downstream propagation (shocks moving from suppliers to customers). This occurs via input shortages. Second, it calculates upstream propagation (shocks moving from customers to suppliers). This occurs via lost demand.

The mathematical core of their approach is the "effective market share" ($\sigma_i(t)$). This is not a static number. It is a dynamic value that changes as the network is shocked. If a firm's competitors in the same industry also fail, that firm's effective market share rises. This makes it much harder for customers to replace them. The authors illustrate this in .

Figure 1
Figure 1: Nonlinear amplification of systemic risk in a supply chain network. The left column illustrates a schematic example of nonlinear amplification. In panel A, a single supplier firm is shocked, marked by the red cross. As a result, this firm completely stops production and is shown in dark purple. The shock then propagates through the network and reduces production at downstream firms. Yellow indicates no production loss for the second supplier, blue indicates a 50% production reduction for the focal firm, and green indicates a 25% production reduction for the three final nodes. After the initial shock has cascaded through the network, total damage amounts to 40% of the total network production. Panel B shows the effect of the failure of the second supplier firm alone, which leads to a similar reduction in production among downstream firms. In panel C, both supplier firms fail simultaneously. The central firm must now stop production completely. Because the complete loss of inputs cannot be compensated for by the three final firms, they also cease production. Total damage therefore amounts to 100% of network production, more than twice the damage caused by the failure of either supplier individually. The right column shows an example from the crustacean supply chain network. In panel D, a single firm fails, marked by the red cross. As a result, some firms reduce production, shown by the green dots. In panel E, a different firm fails, leading to a much smaller cascade. In panel F, both firms fail simultaneously, triggering a large cascade in which many firms must stop production completely, shown by the dark purple dots.

A single supplier failure can be partially compensated for. However, a dual failure leaves the central firm with no alternatives. This forces a total halt in production.

To find these rare, high-risk combinations, the researchers used a "random chemistry" inspired extraction procedure. They do not check every possible combination. Instead, they first shock large, random sets of firms. They then systematically prune them. They do this until they find the smallest subset that still triggers a massive, nonlinear cascade.

Rare but catastrophic amplifications

The study finds that while massive amplification is rare, it is mathematically extreme. Analyzing the Ecuadorian national supply chain, the authors report that only 0.14% of randomly sampled firm pairs exhibit an amplification factor ($\alpha$) greater than 4. This means most pairs behave predictably, but a tiny fraction causes disproportionate damage. Yet, the tail of this distribution is heavy. As shown in, the frequency of these high-amplification events follows a power law.

Figure 3
Figure 3: Number of pairs with amplification factor larger than α . We show the number of firm pairs whose amplification factor, α , exceeds a given threshold. The blue curve represents all firm pairs in the crustacean subnetwork, where amplification factors decay rapidly and reach a maximum value of α = 15. The green curve represents all firm pairs in the softdrink network. Here, amplification factors decay more slowly and reach values of up to α = 59. The light-grey curve represents two million randomly sampled firm pairs from the full supply chain network. The bulk of this distribution behaves similarly to that of the softdrink subnetwork, but its tail is substantially heavier: ten pairs exhibit amplification factors, α > 100, with a maximum value of α = 257. The dashed line shows a power law fit, following [32], with an exponent of 1.79 for amplification factors larger than 4. The black curve additionally includes the firm pairs identified by our extraction procedure, which slightly increases the weight of the tail. Note that large amplification factors can occur even when the total impact of a firm pair is small.

This implies that extreme "black swan" events occur more often than standard statistics would suggest.

The most striking result involves a specific pair of firms in the Ecuadorian network. Their individual impacts were negligible. However, their joint failure resulted in an ESRI of 0.94. This means 94% of the total network production was lost. Compared to the sum of their individual risks (0.0036), this represents an amplification factor of 257. The authors categorize these events into three structural "modes" in .

Figure 4
Figure 4: Schematic network configurations leading to strong nonlinear amplification of systemic risk. We show three representative network modes observed among firm pairs with large systemic risk. Color denotes the industry classification of a node, and its size denotes the initial market share. Red links denote buyer-supplier relationships that can completely disable the buyer if the supplier fails and cannot be replaced, i.e., high impact links with Λ ij = 1. Black links are relationships with limited impact on the buying node, i.e., Λ ij < 1. The mode in panel A consists of two nodes in the same industry, one with a large initial market share, but only a limited impact on its customers, and a small firm that serves as a high impact supplier to a "plateau-firm". In mode B, the two nodes share an industry, with a large combined market share and high impact links to many customers. Panel C shows a mode with a pair of firms in different industries. Here, one of the firms affects a large part of the other firms' industry, increasing its effective market share, σ i ( t ), see Eq. 2, thus rendering it irreplaceable.

One mode involves a small supplier whose failure becomes catastrophic only when a large player in its same industry also fails. This strips the network of its ability to substitute the small player.

Hidden vulnerabilities and model gaps

The authors note that these findings represent an upper bound on potential damage. The model uses a simplified heuristic for supplier replacement. It does not simulate active "rewiring" (the process of a firm actively seeking and signing new contracts). Therefore, the actual speed of a network's recovery might differ. The model also ignores price adjustments. In reality, wealthy firms might outbid others for scarce supplies to keep their lines running.

There are also data-driven limitations. The reconstruction relies on 2015 VAT data. This provides a snapshot of transactions. It lacks real-time information on inventory levels or production lead times. Without knowing how much stock a firm holds, it is hard to predict exactly when a supply delay becomes a production halt. Additionally, the extraction method is stochastic (randomized). While it is effective at finding destructive sets, it cannot mathematically guarantee that every possible combination has been identified.

Identifying the plateau

The research suggests that systemic risk monitoring must evolve. Instead of simply ranking firms by individual importance, managers should look for "plateau-firms." These are groups of firms whose combined failure disables the substitution mechanism for an entire industry.

The verdict depends on the application. For economic planning, this work warns that "safe" looking industries can harbor hidden, nonlinear vulnerabilities. For practitioners building resilient logistics networks, the tool is ready for prototyping. Users should supplement it with real-world data on contract flexibility and inventory buffers. Code for the extraction procedure is reportedly available; see the paper for the canonical link.

Figures from the paper

Figure 2
Figure 2: Scatter plots of the ESRI for simultaneous firm failures versus the sum of individual firms' ESRIs. Panels A and B show the ESRI of simultaneous failure, ESRI( ψ i + ψ j ), against the sum of the individual ESRIs, ESRI( ψ i )+ESRI( ψ j ), for all possible pairs in the crustacean (A) and softdrink supply chain network (B) respectively. The dashed diagonal indicates equality between the ESRI of the simultaneous failure and the sum of the individual ESRIs. In the crustacean supply chain network, most pairs are clustered around the diagonal, indicating that the damage caused by simultaneous failure is comparable to the sum of the damages caused by the individual failures. Only a few pairs exhibit clearly sublinear or superlinear behavior. The red dot highlights the example shown in Fig. 1. In the softdrink supply chain network, panel B, many pairs exhibit sublinear behavior, i.e., they lie below the diagonal. Conversely, a small number of pairs lie above the diagonal: these pairs have low individual ESRIs but cause large systemic damage when they fail simultaneously. In particular, red dots mark six firm pairs with an amplification factor, α , larger than three, whose simultaneous failure affects almost the entire supply chain network. Panel C shows the corresponding results for the complete national supply chain network of Ecuador. Black dots represent a random sample of 2,000,000 firm pairs. Colored symbols indicate small firm sets identified by our extraction procedure. Our method identifies pairs (blue dots), triples (green triangles), and quadruples (purple diamonds) whose simultaneous failure can affect the failure of almost the entire network, ESRI( ∑ i ψ i ) ≈ 1, even though any individual firm has an ESRI below 10%. The inset shows a blowup of the region near the origin. There, we identify quintuples (red stars) with an amplification factor, α , close to 3. This corresponds to our detection threshold θ 1 = 3, indicated by the gray dashed line. For our extraction procedure, we use 100,000 sets of N = 500 and N = 100 initial firms, and choose as thresholds θ 1 = 3 and θ 2 = 0 . 9; see Materials and Methods.
Figure 5
Figure 5: Sensitivity analysis of the sampling method We show scatter plots of the ESRI of the simultaneous failure of identified firm sets against the sum of the ESRIs of the constituent firms. In each plot, we vary one of the parameters of the sampling scheme used to find the sets of firms, and unless otherwise noted, we use 100,000 sets of 500 firms each, with thresholds for identification chosen at θ 1 = 3 and θ 2 = 0 . 9. In panel ,A we use initial set sizes of 100 (green) and 500 firms (black). The larger initial sets are more efficient at identifying the destructive sets. In panel ,B we vary the first threshold from 3 to 6. This threshold controls how much larger the impact of a set of 500 firms has to be than the sum of its constituent parts to be considered for further examination. The larger threshold misses dots with a lower amplification factor, but still succeeds at identifying the strongly amplified sets. In panel ,C we examine the impact of the second threshold θ 2 with values of 0.9 (black), 0.75 (green), and 0.5 (blue). Almost all sets are successfully identified regardless of the choice of θ 2 . Exceptions occur for very small risk values that are missed by the more stringent threshold and a single pair of firms with a large amplification factor.
Figure 6
Figure 6: ESRI of simultaneous failure against the sum of individual ESRIs for the Softdrink subnetwork We show the value of the pair failure ESRI against the sum of individual ESRI values. Black dots are the results of an exhaustive search of all possible firm pairs. Blue dots represent node pairs identified with our extraction procedure. Note how, for the large ESRI pairs, our extraction method manages to identify all of them. We use 100,000 sets of N = 5 firms with θ 1 = 3 and θ 2 = 0 . 9.
Novelty
0.0/10
Overall
0.0/10
#systemic risk#supply chain#nonlinear dynamics#cascading failures#network science
How this was made
Generation

Model: nvidia/Gemma-4-26B-A4B-NVFP4
Persona: academic_accessible
Template: engineering_deepdive
Refinement: 0
Pipeline: forge-1.1

Verification

Evaluator: nvidia/Gemma-4-26B-A4B-NVFP4
Score: 95% (passed)
Claims verified: 15 / 15

Translation

Model: nvidia/Gemma-4-26B-A4B-NVFP4

Hardware & cost

NVIDIA GB10 · 128 GB unified · NVFP4 · 100% local · $0 cloud
Tokens: 122,475
Wall-time: 231.3s
Tokens/s: 529.5