Hall-Like Transversal Stress and Sandpile Criticality on Real Production Networks

1. Énoncé du problème

Les crises économiques se propagent rarement par un canal unique ; les chocs subissent souvent une transformation transdimensionnelle. Par exemple, un choc énergétique peut se manifester sous forme d'inflation, de contraintes budgétaires ou de contraction industrielle, selon les nœuds spécifiques affectés. Les modèles de réseau traditionnels de contagion sont principalement longitudinaux, supposant que les chocs se propagent linéairement le long d'arêtes pondérées à partir d'un seul nœud. Ces modèles échouent à capturer les scénarios où un champ externe applique une contrainte simultanément sur de nombreux flux exposés, convertissant cette charge en une pression systémique d'une nature différente (contrainte transversale).

L'article comble le manque de compréhension concernant la manière dont les champs externes (par exemple, chocs géopolitiques, monétaires ou énergétiques) interagissent avec la fragilité structurelle (faible redondance, faible capacité) dans les réseaux de production réels pour déclencher des avalanches systémiques, distinctes de la criticité auto-organisée (SOC) observée dans les treillis synthétiques.

2. Méthodologie : Le modèle Tas de sable-Hall

L'auteur propose un modèle hybride combinant un mécanisme de contrainte transversale de type Hall avec une dynamique de seuil de tas de sable, calibré sur la Base de données d'entrées-sorties mondiale (WIOD) pour la période 2000–2014 (2 283 nœuds pays-secteur).

A. Le substrat réseau

Le modèle utilise la matrice de transactions intermédiaires observée du WIOD ($Z_t$) pour construire trois opérateurs de propagation alternatifs afin de garantir que les résultats ne soient pas des artefacts mécaniques de normalisations stochastiques par ligne :

1. Matrice de part de ligne ($A^{share}$) : Les lignes somment à 1 (rayon spectral $\approx 1$).
2. Matrice ajustée aux fuites ($A^{leak}$) : Intègre un facteur de fuite $\ell_{i,t}$ (part des flux sortants vers la demande finale), résultant en un rayon spectral $\rho \approx 0,33–0,34$. C'est l'opérateur principal utilisé.
3. Matrice normalisée par la ligne maximale ($A^{max}$) : Préserve l'hétérogénéité inter-nœuds sans contraintes stochastiques.

B. Contrainte transversale de type Hall

S'appuyant sur une analogie avec l'effet Hall physique, le modèle définit une contrainte transversale ($H_{i,t}$) générée lorsqu'un champ externe ($B_t$) agit sur des nœuds ayant une forte intensité de flux ($I_{i,t}$) mais de faibles tampons structurels.
$$H_{i,t} = \frac{B_t I_{i,t}}{D_{i,t} C_{i,t} + \epsilon}$$
Où :

• $B_t$ : Intensité du champ externe.
• $I_{i,t}$ : Intensité de flux relative du nœud $i$.
• $D_{i,t}$ : Redondance (substituabilité de la production).
• $C_{i,t}$ : Capacité d'absorption.
• $\epsilon$ : Constante de régularisation.

Cette contrainte agit comme une conversion multiplicative du champ externe en contrainte interne au nœud, abaissant le seuil d'activation effectif.

C. Dynamique du tas de sable

Le système évolue via une règle d'accumulation de contrainte en temps discret :
$$s_{i,t+1} = (1-\delta)s_{i,t} + \alpha x_{i,t} + \beta \sum_j A^{leak}_{ji,t} s_{j,t} + \gamma H_{i,t}$$

• Renversement (Toppling) : Le nœud $i$ se renverse si $s_{i,t} \geq \theta_i$.
• Avalanche : Lorsqu'un nœud se renverse, l'excès de contrainte est redistribué aux voisins via $A^{leak}$, tandis que le reste est dissipé.
• Ajustement du seuil : La Proposition 1 montre que la contrainte transversale abaisse efficacement le seuil d'activation : $\theta^H_i = \theta_i - \gamma H_{i,t}$.

D. Conception expérimentale

L'étude emploie des expériences de Monte Carlo contrôlées sur le substrat WIOD de 2014. Elle balaye une grille 2D de :

1. Intensité du champ externe ($\bar{B}$) : Variant de 0,25 à 2,0.
2. Contrainte de redondance ($\sigma_D$) : Une réduction contrôlée de la redondance effective ($D_{i,t} \to D_{i,t}/\sigma_D$).

3. Contributions clés

1. Analogie structurelle : Formalise un mécanisme « de type Hall » où les champs externes génèrent une contrainte transversale proportionnelle à l'exposition aux flux et inversement proportionnelle aux tampons structurels (redondance/capacité).
2. Calibration sur réseau réel : Contrairement aux modèles SOC standards utilisant des treillis synthétiques, ce modèle est ancré dans la topologie réelle de l'économie mondiale (WIOD), évitant les hypothèses idéalisées.
3. Identification de régimes : Identifie quatre régimes ordonnés distincts basés sur l'interaction conjointe de l'intensité du champ et de la contrainte de redondance, dépassant les classifications binaires « stable/instable ».
4. Criticité de taille finie : Fournit une analyse rigoureuse de taille finie montrant que, bien que les distributions de queues s'épaississent sous contrainte, elles n'exhibent pas de criticité universelle de loi de puissance (SOC) dans les réseaux réels de cette taille.

4. Résultats clés

A. Propriétés spectrales

• L'opérateur de part de ligne a un rayon spectral proche de 1 (artefact mécanique).
• L'opérateur ajusté aux fuites a un rayon spectral de ~0,33–0,34, confirmant que le système est intrinsèquement sous-critique et dissipatif. La dynamique d'avalanche est pilotée par l'interaction de la charge et de la topologie, et non par un opérateur de rayon unité.

B. Concentration de l'exposition

• Le potentiel de contrainte transversale ($H_{rel}$) est fortement concentré dans la queue supérieure.
• En 2014, le nœud principal (Collecte et distribution d'eau en Chine) représentait ~32 % de l'exposition relative totale.
• Les nœuds à forte exposition sont souvent des secteurs de services ou d'utilités publiques avec de forts flux entrants mais des flux sortants intermédiaires nuls enregistrés (faible redondance dans la structure des données).

C. Quatre régimes ordonnés

Les expériences contrôlées révèlent une transition de phase le long d'une frontière à pente positive dans le plan $(\bar{B}, \sigma_D)$ :

1. Absorption stable : Champ faible, redondance élevée. Taille moyenne de l'avalanche $\approx 0,08$.
2. Fragilité latente : Contrainte modérée. Taille moyenne $\approx 0,5$.
3. Transition critique : Augmentation conjointe du champ et de la contrainte. Taille moyenne $\approx 2,0$.
4. Régime d'avalanche : Champ élevé, redondance faible. Taille moyenne $\approx 5,8$.

• Constat : Ni une forte intensité de champ ni une faible redondance seules ne suffisent à déclencher une avalanche ; elles doivent augmenter conjointement.

D. Comportement des queues et criticité

• Épaississement des queues : À mesure que le système se déplace vers le régime d'avalanche, la distribution des tailles d'avalanche s'épaissit.
• Exposants de queue ($\alpha$) : Les indices de queue estimés sont $\alpha \approx 5,9–6,2$ dans les régimes actifs.
• Signification : Ces exposants sont nettement plus raides que la plage SOC canonique ($\alpha \in [1, 2]$). L'article conclut que l'économie mondiale, telle que modélisée, exhibe une dynamique d'avalanche de taille finie plutôt qu'une criticité auto-organisée universelle. Le système est sous-critique et borné par sa taille finie ($N=2\,283$).

5. Signification et implications

Signification théorique

• Au-delà de la contagion linéaire : Le modèle démontre que les systèmes économiques ne transmettent pas seulement les chocs ; ils les convertissent à travers les dimensions.
• Réfutation de la SOC universelle : Il remet en question l'application de modèles SOC universels de loi de puissance aux réseaux économiques mondiaux, suggérant que les réseaux dissipatifs réels opèrent dans un régime sous-critique de taille finie où les événements extrêmes sont rares mais dépendants du régime.

Implications politiques

• Ciblage de la conversion de contrainte : Puisque les champs externes ($B_t$) sont souvent exogènes (par exemple, prix mondiaux du pétrole), la politique doit se concentrer sur le dénominateur de l'équation de contrainte : augmenter la Redondance ($D$) et la Capacité ($C$).
• Interventions granulaires : Comme l'exposition est fortement concentrée, améliorer les tampons pour les 15 premiers nœuds (qui détiennent ~90 % du potentiel d'exposition) produit des bénéfices systémiques disproportionnés par rapport aux interventions sur les nœuds médians.
• Politique anticipative : Les politiques de résilience doivent être anticipatives. Un réseau stable dans des conditions bénignes peut se situer près de la « frontière d'avalanche » si les champs externes augmentent ; attendre que le choc se matérialise avant de construire de la redondance est insuffisant.

Limites

• Ancienneté des données : Basé sur le WIOD 2014 ; des données plus récentes pourraient modifier les caractéristiques topologiques.
• Variables proxy : La redondance et la capacité sont dérivées uniquement des flux d'entrées-sorties, manquant de données de bilan ou d'inventaire.
• Champs contrôlés : Les champs externes sont simulés, et non des chocs historiques empiriquement mesurés.

Conclusion

L'article fournit un cadre non linéaire et traitable pour comprendre comment les chocs externes activent la fragilité structurelle dans les réseaux de production réels. En combinant un mécanisme de conversion de contrainte de type Hall avec une dynamique de tas de sable, il démontre que le risque systémique découle de la configuration conjointe de la charge externe et des tampons structurels internes. Les résultats suggèrent que, bien que l'économie ne soit pas dans un état de criticité universelle, elle possède une fragilité de taille finie et dépendante du régime qui peut être gérée en ciblant les coefficients de conversion de contrainte des nœuds les plus exposés.