pyAvalanches: A Python Package for Analyzing Spatiotemporal Propagation in Neuronal Avalanches
Abstract
The analysis of neuronal avalanches offers insights into brain dynamics utilizing the framework of criticality, but the reproducibility and comparability of studies are limited by the use of fragmented, lab-specific scripts. To address this issue, we introduce pyAvalanches, an open-source Python package providing a standardized, end-to-end pipeline for avalanche analysis from electrophysiological recordings (e.g., electroencephalography-EEG). Starting from the detection of neuronal avalanches the package provides their core statistical characterization, including size and duration distributions. Beyond this, the main aim of pyAvalanches is to characterize the spatiotemporal organization of activity propagation during avalanches. To this end, the core innovation of pyAvalanches is the compuation of Avalanche Transition Matrices (ATMs) to map spatiotemporal propagation patterns. Building on this, the package derives network-based metrics from the ATMs, bridging the study of the topology and organization of the underlying dynamical interactions with network neuroscience adopting the framework of neuronal avalanches. The entire workflow is encapsulated in a modular and scikit-learn compatible architecture. We demonstrate the utility of pyAvalanches through an illustrative group-level analysis on a public resting-state EEG dataset, comparing propagation patterns across different clinical populations. By providing a user-friendly, tested, and extensible tool, pyAvalanches facilitates reproducible research, enables the development of novel avalanche-based biomarkers, and makes complex avalanche analysis accessible to a broader scientific community. The package is fully documented and distributed via the Python Package Index (PyPI).
Related
- Arousal-informed neuronal avalanches operate near criticality
- Neuronal avalanches in awake macaque V1 and V4
- Neuronal Identity as an Organizational Basis for Analyzing Neural Population Dynamics
- Hierarchical connectivity sets the upper critical dimension of neural avalanches
- Eco-evolutionary dynamics create extinction avalanches in mutualistic-antagonistic networks
BibTeX
@misc{marzulli2026pyavalanches,
title = {{pyAvalanches: A Python Package for Analyzing Spatiotemporal Propagation in Neuronal Avalanches}},
author = {M. Marzulli and M. Angiolelli and C. Mannino and M. Demuru and P. Sorrentino and M. -C. Corsi},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2609.11530},
eprint = {2609.11530},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.11530},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - pyAvalanches: A Python Package for Analyzing Spatiotemporal Propagation in Neuronal Avalanches
AU - M. Marzulli
AU - M. Angiolelli
AU - C. Mannino
AU - M. Demuru
AU - P. Sorrentino
AU - M. -C. Corsi
PY - 2026
DA - 2026/10/29
PB - arXiv
DO - 10.48550/arxiv.2609.11530
UR - https://arxiv.org/abs/2609.11530
AB - The analysis of neuronal avalanches offers insights into brain dynamics utilizing the framework of criticality, but the reproducibility and comparability of studies are limited by the use of fragmented, lab-specific scripts. To address this issue, we introduce pyAvalanches, an open-source Python package providing a standardized, end-to-end pipeline for avalanche analysis from electrophysiological recordings (e.g., electroencephalography-EEG). Starting from the detection of neuronal avalanches the package provides their core statistical characterization, including size and duration distributions. Beyond this, the main aim of pyAvalanches is to characterize the spatiotemporal organization of activity propagation during avalanches. To this end, the core innovation of pyAvalanches is the compuation of Avalanche Transition Matrices (ATMs) to map spatiotemporal propagation patterns. Building on this, the package derives network-based metrics from the ATMs, bridging the study of the topology and organization of the underlying dynamical interactions with network neuroscience adopting the framework of neuronal avalanches. The entire workflow is encapsulated in a modular and scikit-learn compatible architecture. We demonstrate the utility of pyAvalanches through an illustrative group-level analysis on a public resting-state EEG dataset, comparing propagation patterns across different clinical populations. By providing a user-friendly, tested, and extensible tool, pyAvalanches facilitates reproducible research, enables the development of novel avalanche-based biomarkers, and makes complex avalanche analysis accessible to a broader scientific community. The package is fully documented and distributed via the Python Package Index (PyPI).
ER -